SDInjection Class Reference

Super-droplets initialization structure. More...

#include <ERF_SDInitialization.H>

Inheritance diagram for SDInjection:

Collaboration diagram for SDInjection:

Public Member Functions
virtual	~SDInjection ()=default
void	updateDt (const amrex::Real a_dt)
	Update time-dependent quantities for particle injection. More...
void	readInputs (const std::string &, const std::string &, const amrex::Geometry &, const MatVec &, const MatVec &) override
void	readInputs (const std::string &a_prefix, const amrex::Geometry &a_geom, const MatVec &a_species_mat, const MatVec &a_aerosol_mat, const amrex::Real a_dt)
	Read super-droplet injection parameters from input file. More...
void	printParameters (const MatVec &, const MatVec &) const override
int	numSDPerCell (const amrex::Real a_dv) const override
	Compute number of super-droplets to inject per grid cell. More...
amrex::Real	numParticlesPerCell (const amrex::Real a_dv) const override
Public Member Functions inherited from SDInitProperties
virtual void	setDefaults (const amrex::Geometry &a_geom, const MatVec &a_species_mat, const MatVec &a_aerosol_mat)
	Set default values for initialization parameters. More...
void	getDistribution (amrex::Vector< amrex::Real > &a_mass, const int a_np, const amrex::Real a_density, const std::string &a_init_type, const amrex::Real a_mass_min, const amrex::Real a_mass_max, const amrex::Real a_mass_mean, const amrex::Real a_radius_min, const amrex::Real a_radius_max, const amrex::Real a_radius_mean, const amrex::Real a_radius_gstd, std::mt19937 &a_rng) const
	Get a distribution with constant multiplicity. More...
void	getDistribution (amrex::Vector< amrex::Real > &a_mass, amrex::Vector< amrex::Real > &a_mult, const amrex::Real a_dV, const int a_np, const amrex::Real a_density, const std::string &a_init_type, const amrex::Real a_mass_min, const amrex::Real a_mass_max, const amrex::Real a_mass_mean, const amrex::Real a_radius_min, const amrex::Real a_radius_max, const amrex::Real a_radius_mean, const amrex::Real a_radius_gstd, std::mt19937 &a_rng) const
	Get a distribution with sampled multiplicity. More...
void	getAerosolDistribution (amrex::Vector< amrex::Real > &a_mass, const int a_idx, const int a_np, const amrex::Real a_density, std::mt19937 &a_rng) const
	Compute the aerosol mass distribution. More...
void	getAerosolDistribution (amrex::Vector< amrex::Real > &a_mass, amrex::Vector< amrex::Real > &a_mult, amrex::Real a_dV, int a_idx, int a_np, amrex::Real a_density, std::mt19937 &a_rng) const
	Compute the aerosol mass distribution with sampled multiplicity. More...
void	getSpeciesDistribution (amrex::Vector< amrex::Real > &a_mass, const int a_idx, const int a_np, const amrex::Real a_density, std::mt19937 &a_rng) const
	Compute the species mass distribution with constant multiplicity. More...
void	getSpeciesDistribution (amrex::Vector< amrex::Real > &a_mass, amrex::Vector< amrex::Real > &a_mult, amrex::Real a_dV, int a_idx, int a_np, amrex::Real a_density, std::mt19937 &a_rng) const
	Compute the species mass distribution with sampled multiplicity. More...
bool	sampledMultiplicity () const
	Determine whether multiplicity is sampled or constant. More...
amrex::Real	volume () const
	Calculate the volume of the particle domain. More...

Public Attributes
amrex::Real	m_inj_rate = 0
amrex::Real	m_sd_inj_rate = 0
amrex::Real	m_numdens_sd = -1
amrex::Vector< amrex::Real >	m_domain_vel = {0.0,0.0,0.0}
amrex::Real	m_tstart = 0.0
amrex::Real	m_tstop = 1.0e17
Public Attributes inherited from SDInitProperties
int	m_ppc = 1
SDInitShape	m_type = SDInitShape::uniform
amrex::Real	m_numdens = -1
amrex::Real	m_max_multiplicity = 1000000
std::vector< amrex::Real >	m_mass_species_min
std::vector< amrex::Real >	m_mass_species_max
std::vector< amrex::Real >	m_mass_species_mean
std::vector< amrex::Real >	m_radius_species_min
std::vector< amrex::Real >	m_radius_species_max
std::vector< amrex::Real >	m_radius_species_mean
std::vector< amrex::Real >	m_radius_species_geom_std
std::vector< std::string >	m_species_init_type
std::vector< amrex::Real >	m_mass_aerosol_min
std::vector< amrex::Real >	m_mass_aerosol_max
std::vector< amrex::Real >	m_mass_aerosol_mean
std::vector< amrex::Real >	m_radius_aerosol_min
std::vector< amrex::Real >	m_radius_aerosol_max
std::vector< amrex::Real >	m_radius_aerosol_mean
std::vector< amrex::Real >	m_radius_aerosol_geom_std
std::vector< std::string >	m_aerosol_init_type
int	m_num_species = 0
int	m_num_aerosols = 0
amrex::RealBox	m_particle_domain
amrex::Vector< amrex::Real >	m_init_particle_p1
amrex::Vector< amrex::Real >	m_init_particle_p2
SDMultiplicityType	m_mult_type

Additional Inherited Members
Public Types inherited from SDInitProperties
using	MatVec = std::vector< std::unique_ptr< MaterialProperties > >

Detailed Description

Super-droplets initialization structure.

Constructor & Destructor Documentation

~SDInjection()

virtual SDInjection::~SDInjection ( )

virtualdefault

Injection rate (number of physical particles per meter^3 second)

Member Function Documentation

numParticlesPerCell()

amrex::Real SDInjection::numParticlesPerCell ( const amrex::Real  a_dv ) const

inlineoverridevirtual

Compute number of physical particles per grid cell

Implements SDInitProperties.

        {
            amrex::Real num_par_per_cell = 0.0;
            if (this->m_numdens >= 0) {
                num_par_per_cell = std::ceil(this->m_numdens*a_dv);
            } else {
                num_par_per_cell = 1;
            }
            return num_par_per_cell;
        }

numSDPerCell()

int SDInjection::numSDPerCell ( const amrex::Real  a_dv ) const

inlineoverridevirtual

Compute number of super-droplets to inject per grid cell.

Determines the number of super-droplets to inject based on either:

1. The super-droplet number density and cell volume, or
2. A constant number per cell (m_ppc) if number density is not set

Parameters

[in]

a_dv

Volume of the grid cell

Returns

Number of super-droplets to inject in this cell

Implements SDInitProperties.

        {
            int num_sd_per_cell = 0;
            if (m_numdens_sd >= 0) {
                num_sd_per_cell = static_cast<int>(std::ceil(m_numdens_sd*a_dv));
            } else {
                num_sd_per_cell = this->m_ppc;
            }
            return num_sd_per_cell;
        }

printParameters()

void SDInjection::printParameters ( const MatVec &  a_species_mat, const MatVec &  a_aerosol_mat  ) const

overridevirtual

print super-droplets injection parameters to screen

Reimplemented from SDInitProperties.

{
    using namespace amrex;
    Print() << "    Injection rate (# m^{-3} s^{-1}): " << m_inj_rate << "\n"
            << "    Injection domain velocity [m/s]: "
            << m_domain_vel[0] << ","
            << m_domain_vel[1] << ","
            << m_domain_vel[2] << "\n"
            << "    Time (start, stop) [s]: " << m_tstart << ", " << m_tstop << "\n"
            << "    SD injection rate: " << m_sd_inj_rate << "\n";
    SDInitProperties::printParameters(a_species_mat, a_aerosol_mat);
}

Here is the call graph for this function:

readInputs() [1/2]

void SDInjection::readInputs ( const std::string &  , const std::string &  , const amrex::Geometry &  , const MatVec &  , const MatVec &    )

inlineoverridevirtual

read super-droplets injection parameters

Reimplemented from SDInitProperties.

        {
            amrex::Abort("SDInjection::readInputs(): Do not use this interface");
        }

readInputs() [2/2]

void SDInjection::readInputs	(	const std::string &	a_prefix,
		const amrex::Geometry &	a_geom,
		const MatVec &	a_species_mat,
		const MatVec &	a_aerosol_mat,
		const amrex::Real	a_dt
	)

Read super-droplet injection parameters from input file.

This function reads injection-specific parameters like injection rate, timing parameters, and domain velocity.

Parameters

[in]	a_prefix	Prefix for parameter parser
[in]	a_geom	Simulation geometry information
[in]	a_species_mat	Vector of species material properties
[in]	a_aerosol_mat	Vector of aerosol material properties
[in]	a_dt	Current timestep size

{
    BL_PROFILE("SDInjection::readInputs");
 
    SDInitProperties::readInputs( a_prefix, "", a_geom, a_species_mat, a_aerosol_mat);
 
    amrex::ignore_unused(a_geom);
    using namespace amrex;
 
    amrex::ParmParse pp(a_prefix);
    pp.query("rate", m_inj_rate);
    pp.query("sd_rate", m_sd_inj_rate);
    pp.query("t_start", m_tstart);
    pp.query("t_stop", m_tstop);
    pp.queryarr("domain_velocity", m_domain_vel);
 
    this->m_numdens = m_inj_rate * a_dt;
    m_numdens_sd = (m_sd_inj_rate > 0 ? std::max(m_sd_inj_rate*a_dt, 1.0) : -1);
}

Here is the call graph for this function:

updateDt()

void SDInjection::updateDt ( const amrex::Real  a_dt )

inline

Update time-dependent quantities for particle injection.

This function updates:

1. Number density based on injection rate and timestep
2. Super-droplet number density based on SD injection rate
3. Position of the injection domain based on domain velocity

Parameters

[in]

a_dt

Timestep size in seconds

        {
            this->m_numdens = m_inj_rate * a_dt;
            m_numdens_sd = (m_sd_inj_rate>0 ? std::max(m_sd_inj_rate*a_dt, 1.0) : -1);
 
            if (this->m_type == SDInitShape::uniform) {
                amrex::Vector<amrex::Real> lo = {0.0, 0.0, 0.0};
                amrex::Vector<amrex::Real> hi = {0.0, 0.0, 0.0};
                for (int dir = 0; dir < AMREX_SPACEDIM; dir++) {
                    lo[dir] = this->m_particle_domain.lo(dir) + m_domain_vel[dir] * a_dt;
                    hi[dir] = this->m_particle_domain.hi(dir) + m_domain_vel[dir] * a_dt;
                }
                this->m_particle_domain.setLo(lo);
                this->m_particle_domain.setHi(hi);
            } else if (m_type == SDInitShape::bubble) {
                amrex::Vector<amrex::Real> center = {0.0, 0.0, 0.0};
                for (int dir = 0; dir < AMREX_SPACEDIM; dir++) {
                    center[dir] = this->m_particle_domain.lo(dir) + m_domain_vel[dir] * a_dt;
                }
                this->m_particle_domain.setLo(center);
            }
        }

Member Data Documentation

m_domain_vel

amrex::Vector<amrex::Real> SDInjection::m_domain_vel = {0.0,0.0,0.0}

Start time for injection

Referenced by updateDt().

m_inj_rate

amrex::Real SDInjection::m_inj_rate = 0

Injection rate (number of SDs per meter^3 second)

m_numdens_sd

amrex::Real SDInjection::m_numdens_sd = -1

Injection domain velocity

Referenced by numSDPerCell(), and updateDt().

m_sd_inj_rate

amrex::Real SDInjection::m_sd_inj_rate = 0

Initial number density of super-droplets

Referenced by updateDt().

m_tstart

amrex::Real SDInjection::m_tstart = 0.0

Stop time for injection

m_tstop

amrex::Real SDInjection::m_tstop = 1.0e17

---

The documentation for this class was generated from the following files:

• Source/Particles/ERF_SDInitialization.H
• Source/Particles/ERF_SDInitialization.cpp