mus_particle_config_module Module

mus_particle_config_module contains routines for loading the data for LBM-DEM simulations of particles in a flow from the lua script. All particle data is read from the "particles" table in the lua script. This table should look as follows: particles = { nParticles = 1, -- Total number of particles in the simulation -- particle kind, can be: -- * 'MEM' for fully-resolved particles using the Momentum-Exchange Method -- * 'DPS' for four-way coupled unresolved particles based on the -- Volume-Averaged Navier-Stokes equations -- * 'DPS_twoway' for two-way coupled unresolved particles, neglecting the -- effect of volume fraction -- * 'DPS_oneway' for one-way coupled unresolved particles. That is, particles -- do not affect the flow kind = 'DPS',

-- Describe how particles should interact with boundaries. They can either bounce -- off them ('wall') or have periodic boundaries ('periodic'). Currently only -- interactions with prismatic boundaries are supported. If not specified, -- particles will be removed from the domain once they hit a boundary of the -- fluid domain (also works for arbitrary geometries). boundaries = { domain_bnd = { 0.0, L, 0.0, W, 0.0, H }, bnd_kind = {'wall', 'wall', 'wall', 'wall', 'wall', 'wall'} },

-- Number of DEM subcycles per LBM time step to use nDEMsubcycles = 50,

-- How often to log the particle data (1 = log every LBM time step) particleLogInterval = 1,

-- Size of the buffers in the particle communication routines. Should be -- sufficient to hold all particles anticipated to be on a domain boundary -- at any time during the simulation. Safest option is to make this equal -- to or greater than nParticles. particleBufferSize = 100,

-- Collision time (in physical units) of particle-particle collisions particle_collision_time = 0.1*dt,

-- Collision tolerance, i.e. what should the gap between particles be -- before we call it a collision particle_collision_tol = 0.0,

-- Size of the gap before capping the calculation of lubrication forces critical_gap_Flub = 0.01,

-- Particle positions. Should contain nParticles tables of the form -- {x,y,z,rx,ry,rz} with x, y, z the translational positions and -- rx, ry, rz the rotation angles about the corresponding axes. position = { { x_p_phy, y_p_phy, z_p_phy, 0.0, 0.0, 0.0} },

-- Particle initial velocities. Should contain nParticles tables of the form -- {ux, uy, uz, urx, ury, urz} velocity = { {0.0, 0.0, 0.0, 0.0, 0.0, 0.0} },

-- External forces on particles. Should contain nParticles tables of the form -- {Fx, Fy, Fz, Frx, Fry, Frz} force = { {0.0, 0.0, -F_gravity + F_buoyancy, 0.0, 0.0, 0.0} },

-- Particle radii, should have nParticles entries radius = { 0.5*Dia_p_phy },

-- Particle masses, should have nParticles entries mass = { m_p_phy } }