Initialization of liquid droplets

[le-snuapl]

Dear all,

So far I have no experience in using the liquid droplets in PeleLM.
I want to initialize N droplets of diameter D at random locations in the domain. I dont need any continuous injection.
From looking at the SprayTest I think I understand how to set the thermodynamic data in the input file.
But how do I define the location size and number of droplets at the beginning of the simulation?

Thank you very much for your time already.

[baperry2]

SprayTest is a good example to follow. Should be a pretty simple modification to the SprayParticleContainer::uniformSprayInit function defined in PelePhysics/Source/Spray/SprayInjection.cpp

[le-snuapl]

Thank you for the fast reply! I wanted to check how the uniformSprayInit works first in a simple 2D case without any modifications.
I prepared a 2D case with just 1 particle but it keeps exiting with the "Bad Particle!" notification.
I tried modifying the amount of particles as well as the particle parameter data but it doesnt help.
Probably its something really obvious but I dont really see it. I have attached the settings and the SprayParticlesInitInsert which I copied from the SprayTest.

input.2d:
---------------------- DOMAIN DEFINITION ------------------------
geometry.is_periodic = 0 0 # For each dir, 0: non-perio, 1: periodic
geometry.coord_sys = 0 # 0 => cart, 1 => RZ
geometry.prob_lo = -0.1 -0.1 0.0 # x_lo y_lo (z_lo)
geometry.prob_hi = 0.1 0.1 0.0 # x_hi y_hi (z_hi)

#---------------------- BC FLAGS ---------------------------------

Interior, Inflow, Outflow, Symmetry,

SlipWallAdiab, NoSlipWallAdiab, SlipWallIsotherm, NoSlipWallIsotherm

peleLM.lo_bc = Outflow Outflow # bc in x_lo y_lo (z_lo)
peleLM.hi_bc = Outflow Outflow # bc in x_hi y_hi (z_hi)

#---------------------- AMR CONTROL ------------------------------
amr.n_cell = 128 128 # Level 0 number of cells in each direction
amr.max_level = 0 # maximum level number allowed

#---------------------- Problem ----------------------------------

parameters for the spray

prob.num_particles = (1, 1, 1)
prob.part_temp = 230.0
prob.part_dia = 20.E-6
prob.part_vel = 0.0 0.0 0.0

InitSprayParticles function:
SprayParticleContainer::InitSprayParticles(const bool init_parts)
{
if (!init_parts) {
return;
}
const int level = 0;
amrex::ParmParse pp("prob");
int numRedist = -1;
amrex::IntVect partNum(AMREX_D_DECL(100, 100, 100));
// Find the number of redistributions during particle initialization
pp.query("init_redist", numRedist);
pp.query("num_particles", partNum);
amrex::RealVect partVel;
std::array<amrex::Real, AMREX_SPACEDIM> pvel = {0.0};
pp.queryamrex::Real("part_vel", pvel);
for (int dir = 0; dir < AMREX_SPACEDIM; ++dir) {
partVel[dir] = pvel[dir];
}
amrex::Real partDia, partTemp;
pp.get("part_dia", partDia);
pp.get("part_temp", partTemp);
std::array<amrex::Real, SPRAY_FUEL_NUM> partY = {0.0};
partY[0] = 1.;
uniformSprayInit(
partNum, partVel, partDia, partTemp, partY.begin(), level, numRedist);
}

[Soobin120596]

I think you can refer to the file initspraydata_3d.dat in RegTests\Spray-EB to initialize the droplets. The droplet initialization file is specified in spray-eb.inp using particles.init_file = "initspraydata_3d.dat".