programmer's documentation
Data Types | Modules | Functions/Subroutines | Variables
lagran.f90 File Reference

Module for Lagrangian model. More...

Data Types

interface  cs_f_lagr_pointers
 
interface  lagr_part_copy
 Copy attributes from one particle to another. More...
 
interface  lagr_current_to_previous
 Copy current attributes to previous attributes. More...
 
interface  cs_lagr_resize_particle_set
 Resize particle set buffers if needed. More...
 
interface  lagr_restart_read_particle_data
 Read particle data from checkpoint. More...
 
interface  lagr_restart_write_particle_data
 Write particle data to checkpoint. More...
 

Modules

module  lagran
 

Functions/Subroutines

subroutine init_zone_class_param (ii, jj)
 
subroutine lagr_init_zone_class_param (i_cz_params, r_cz_params)
 Initialize Lagrangian module parameters for a given zone and class. More...
 
subroutine lagr_define_zone_class_param (class_id, zone_id, i_cz_params, r_cz_params)
 Define Lagrangian module parameters for a given zone and class. More...
 
subroutine lagr_update_pointers ()
 Update Lagrangian pointers. More...
 
integer function lagr_resize_particle_set (n_min_particles)
 Ensure particle data buffers are large enough. More...
 
subroutine lagr_status (model, restart, frozen)
 Return Lagrangian model status. More...
 
subroutine init_lagr_arrays
 
subroutine finalize_lagr_arrays
 

Variables

integer, save iilagr
 
integer, save isuila
 activation (=1) or not (=0) of a Lagrangian calculation restart. The calculation restart file read when this option is activated only contains the data related to the particles (see also isuist) the global calculation must also be a restart calculation More...
 
integer, save isttio
 
integer(c_int), pointer, save nbpart
 number of current local particles Always useful, but initialized and updated without intervention of the user More...
 
integer, save nbpnew
 number of new entering particles More...
 
integer(c_int), pointer, save nbperr
 number of in error deleted particles More...
 
integer, save nbptot
 total number of injected particles, since the beginning, including calculation restarts More...
 
integer(c_int), pointer, save nbpout
 contains normally exited particles and in detection error exited particles More...
 
integer, save nbpert
 
integer, save ndepot
 number of deposed particles for ever, but keeped for post-processing More...
 
integer(c_int), pointer, save nbpdep
 number of deposed particles More...
 
integer, save nbpres
 number of re-entrained particles More...
 
real(c_double), pointer, save dnbpar
 real value for nbpart More...
 
double precision, save dnbpnw
 real value for nbpnew More...
 
real(c_double), pointer, save dnbper
 real value for nbperr More...
 
real(c_double), pointer, save dnbpou
 real value for nbpout More...
 
real(c_double), pointer, save dnbdep
 real value for nbpdep More...
 
double precision, save dnbres
 real value for nbres More...
 
integer, save npclon
 number of cloned new particles More...
 
integer, save npkill
 number of killed by error particles More...
 
integer, save npcsup
 number of cloned particles More...
 
double precision, save dnpclo
 real value for npclon More...
 
double precision, save dnpkil
 real value for npkill More...
 
double precision, save dnpcsu
 real value for npcsup More...
 
integer, dimension(:,:), pointer, save ipepa
 
double precision, dimension(:,:), pointer, save pepa
 
double precision, dimension(:,:), pointer, save eptp
 
double precision, dimension(:,:), pointer, save eptpa
 
double precision, dimension(:,:), pointer, save ptsvar
 
integer, save iphyla
 activates (>0) or deactivates (=0) the physical models associated to the particles: More...
 
integer, save itpvar
 activation (=1) or not (=0) of an evolution equation on the particle temperature (in degrees Celsius). Useful if iphyla=1 and if there is a thermal scalar associated with the continuous phase More...
 
integer, save idpvar
 activation (=1) or not (=0) of an evolution equation on the particle diameter. Useful if iphyla = 1 More...
 
integer, save impvar
 activation (=1) or not (=0) of an evolution equation on the particle mass Useful if si iphyla = 1 More...
 
double precision, save tpart
 initialization temperature (in degree Celsius) for the particles already present in the calculation domain when an evolution equation on the particle temperature is activated during a calculation (iphyla = 1 and itpvar = 1). Useful if isuila = 1 and itpvar = 0 in the previous calculation More...
 
double precision, save cppart
 initialization value for the specific heat ( $ J.kg^{-1}.K^{-1} $) of the particles already present in the calculation domain when an evolution equation on the particle temperature is activated during a calculation (iphyla = 1 and itpvar = 1). Useful if isuila = 1 and itpvar = 0 in the previous calculation More...
 
integer, save idepst
 
integer, save idlvo
 
integer ngeol
 geometric parameters stored More...
 
integer, save jimark
 Additional pointer in the IPEPA array (contains the particule state) More...
 
integer, save jdfac
 Additional pointer in the IPEPA array (contains the particule state) More...
 
integer, save jryplu
 Additional pointer in the IPEPA array (contains the particule state) More...
 
integer, save jrinpf
 Additional pointer in the IPEPA array (contains the particule state) More...
 
integer, save jord1
 Additional pointer in the IPEPA array (contains the particule state) More...
 
integer, save ireent
 
integer, save irough
 
integer, save iflow
 
integer, save jroll
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
integer, save jnbasg
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
integer, save jnbasp
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
integer, save jdepo
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
integer, save jfadh
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
integer, save jmfadh
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
integer, save jndisp
 Additional pointer in IPEPA and PEPA arrays (contains particule state) More...
 
double precision, save espasg
 Parameter of the particle resuspension model. More...
 
double precision, save denasp
 Parameter of the particle resuspension model. More...
 
double precision, save modyeq
 Parameter of the particle resuspension model. More...
 
double precision, save rayasp
 Parameter of the particle resuspension model. More...
 
double precision, save rayasg
 Parameter of the particle resuspension model. More...
 
integer, save ipreci
 
integer, save nbrclas
 number of particle classes More...
 
double precision, save dprec
 diameter of particles formed by precipitation More...
 
double precision, save rho_preci
 density of particles formed by precipitation More...
 
integer, save iclogst
 
integer, save jnbpoi
 Additional pointer in ITEPA and TEPA arrays (contains particule state) More...
 
integer, save jrtdep
 Additional pointer in ITEPA and TEPA arrays (contains particule state) More...
 
integer, save jrhcon
 Additional pointer in ITEPA and TEPA arrays (contains particule state) More...
 
double precision, save jamlim
 Parameter of the particle clogging model. More...
 
double precision, save mporos
 
double precision, save csthpp
 
integer, save iplas
 absolute iteration number (including the restarts) in the Lagrangian module ( i.e. Lagrangian time step number) More...
 
integer, save iplar
 relative iteration number (including the restarts) in the Lagrangian module More...
 
double precision, save dtp
 duration of a Lagrangian iteration More...
 
double precision, save ttclag
 physical time of the Lagrangian simulation More...
 
integer, save ierr
 error indicator More...
 
integer, dimension(:), allocatable nbprec
 
double precision, dimension(:), allocatable solub
 
double precision, dimension(:,:), allocatable mp_diss
 
integer, save jxp
 pointer to particle X coordinate for pointer eptp More...
 
integer, save jyp
 pointer to particle Y coordinate for pointer eptp More...
 
integer, save jzp
 pointer to particle Z coordinate for pointer eptp More...
 
integer, save jup
 pointer to particle X velocity component for pointer eptp More...
 
integer, save jvp
 pointer to particle Y velocity component for pointer eptp More...
 
integer, save jwp
 pointer to particle Z velocity component for pointer eptp More...
 
integer, save juf
 pointer to locally undisturbed X fluid velocity component for pointer eptp More...
 
integer, save jvf
 pointer to locally undisturbed Y fluid velocity component for pointer eptp More...
 
integer, save jwf
 pointer to locally undisturbed Z fluid velocity component for pointer eptp More...
 
integer, save jmp
 pointer to particle mass for pointer eptp More...
 
integer, save jdp
 pointer to particle diameter for pointer eptp More...
 
integer, save jdp2
 pointer to particle extent for array eptp More...
 
integer, save jtp
 pointer to particle and locally undisturbed fluid flow temperature (Celsius) for pointer eptp More...
 
integer, save jtf
 pointer to particle and locally undisturbed fluid flow temperature (Celsius) for pointer eptp More...
 
integer, save jcp
 pointer to particle specific heat for pointer eptp More...
 
integer, dimension(nlayer), save jhp
 pointer to coal particle temperature ( $K$) for pointer eptp More...
 
integer, save jmwat
 pointer to water mass (for coal) for pointer eptp More...
 
integer, dimension(nlayer), save jmch
 pointer to mass of reactive coal of the coal particle for pointer eptp More...
 
integer, dimension(nlayer), save jmck
 pointer to mass of coke of the coal particle for pointer eptp More...
 
integer, dimension(nusvar), save jvls
 pointer to additional user variable for pointer eptp More...
 
integer, save jtaux
 pointer to work array for the second order in time for pointer pepa More...
 
integer, dimension(3), save jbx1
 pointer to turbulence characteristics for second order in time for array pepa More...
 
integer, dimension(3), save jtsup
 pointer to velocity prediction for second order in time for array pepa More...
 
integer, dimension(3), save jtsuf
 pointer to seen velocity prediction for second order in time for array pepa More...
 
integer, save jrval
 pointer to random number associated with a particle for array pepa More...
 
integer, save jrtsp
 pointer to particle residence time for array pepa More...
 
integer, save jrpoi
 pointer to particle statistic weight for array pepa More...
 
integer, save jreps
 pointer to particle emissivity for array pepa More...
 
integer, save jrd0p
 pointer to coal particle initial diameter for array pepa More...
 
integer, save jrdck
 pointer to coal particle shrinking core diameter for array pepa More...
 
integer, dimension(nlayer), save jrhock
 pointer to coal density for array pepa More...
 
integer, save jisor
 pointer to number of the current cell containing the particle for ipepa array; this number is re-actualised during the trajectography step More...
 
integer, save jisora
 pointer to number of the previous cell containing the particle for ipepa array More...
 
integer, save jirka
 pointer to number of the previous rank containing the particle for ipepa array More...
 
integer, save jinch
 pointer to number of the coal particle for ipepa array More...
 
integer, save jclst
 pointer to class of the particle for ipepa array More...
 
integer, save nvls
 pointer to number of additional variables related to the particles for ipepa array. The additional variables can be accessed in the arrays eptp and eptpa by means of the pointer jvls: eptp(nbpt,jvls(ii)) and eptpa(nbpt,jvls(ii)) (nbpt is the index-number of the treated particle, and ii an integer between 1 and nvls) More...
 
integer, save nfrlag
 number of boundary zones More...
 
integer, save injcon
 activates (=1) or not (=0) the continuous injection of particles. this option allows to inject particles continuously during the duration of the Lagrangian time step dtp rather than only once at the beginning of the Lagrangian iteration. It helps avoiding the fractioning of the particles cloud close to the injection areas. More...
 
integer, dimension(nflagm), save ilflag
 list of number of boundary zones More...
 
integer, dimension(nflagm), save iusncl
 for all the nfrlag boundary zones previously identified, the number of classes nbclas (a class is a set of particles sharing the same physical properties and the same characteristics concerning the injection in the calculation domain) of entering particles is given: iusncl(izone) = nbclas. By default, the number of particle classes is zero. The maximum number of classes is nflagm (parameter stored in lagpar, whose default value is 20). More...
 
integer, dimension(nflagm), save iusclb
 for all the nfrlag boundary zones previously identified, a particle boundary condition type is given. The categories of particle boundary condition types are marked out by the key words ientrl, isortl, irebol, idepo1, idepo2, iencrl. More...
 
integer, dimension(nflagm), save iusmoy
 mean over a zone (if mean per zones is activated) More...
 
integer, dimension(:,:,:), allocatable iuslag
 Some pieces of information must be given for each particle class associated with an injection zone. The first part consists in integers contained in the array iuslag. There are at the most ndlaim integers. This data must be provided for each class iclas and each particle injection zone izone. They are marked out by means of "pointers": More...
 
double precision, dimension(nflagm), save deblag
 massic flow rate for a boudary zone More...
 
integer, save ijnbp
 number of particles per class and per boudary zone More...
 
integer, save ijfre
 injection frequency (if < 0 : particle are introduced only at first iteration More...
 
integer, save ijuvw
 velocity condition type: More...
 
integer, save ijprpd
 
integer, save ijprtp
 
integer, save ijprdp
 type of user profiles in uslag2: More...
 
integer, save irawcl
 type of coal initial composition (if iphyla=2) More...
 
integer, save inuchl
 coal number of the particle (if iphyla=2) More...
 
integer, save iclst
 number of the statistics group More...
 
double precision, dimension(:,:,:), allocatable ruslag
 some pieces of information must be given for each particle class associated with an injection zone. The second and last part consists in real numbers contained in the array ruslag. There are at the most ndlagm such real numbers. This data must be provided for each class iclas and each particle injection zone izone. They are marked out by means of "pointers": More...
 
integer, save iuno
 particle velocity magnitude More...
 
integer, save iupt
 particle u component by class and zone More...
 
integer, save ivpt
 particle v component by class and zone More...
 
integer, save iwpt
 particle w component by class and zone More...
 
integer, save itpt
 particle temperature More...
 
integer, save idpt
 particle diameter More...
 
integer, save ivdpt
 particle diameter variance More...
 
integer, save iropt
 density More...
 
integer, save icpt
 particle specific heat More...
 
integer, save ipoit
 particle weight More...
 
integer, save idebt
 flow rate More...
 
integer, save iepsi
 particle emissivity More...
 
integer, dimension(nlayer), save ihpt
 particle temperature More...
 
integer, save ifrmwt
 water mass fraction in coal particles More...
 
integer, dimension(nlayer), save ifrmch
 active coal mass fraction in coal particles More...
 
integer, dimension(nlayer), save ifrmck
 coke mass fraction in coal particles More...
 
integer, save irdck
 diameter of shrinking core More...
 
integer, save ird0p
 initial particle diameter (for coal particles) More...
 
integer, dimension(nlayer), save irhock0
 coke density after pyrolysis (for coal particles) More...
 
integer, save ilvx
 mean dispersed phase velocity X component More...
 
integer, save ilvy
 mean dispersed phase velocity Y component More...
 
integer, save ilvz
 mean dispersed phase velocity Z component More...
 
integer, save ilpd
 sum of the statistical weights More...
 
integer, save ilfv
 dispersed phase volumetric concentration More...
 
integer, save ilts
 recidence time More...
 
integer, save iltp
 phase temperature ( $\text{\degresC}$) More...
 
integer, save ildp
 dispersed phase mean diameter More...
 
integer, save ilmp
 dispersed phase mean mass More...
 
integer, dimension(nlayer), save ilhp
 temperature of the coal particle cloud ( $K$) More...
 
integer, save ilmwat
 water mass More...
 
integer, dimension(nlayer), save ilmch
 mass of reactive coal of the coal particle cloud More...
 
integer, dimension(nlayer), save ilmck
 mass of coke of the coal particle cloud More...
 
integer, save ildck
 shriking core diameter More...
 
integer, dimension(nussta), save ilvu
 supplementary user volumetric statistics More...
 
integer, save iactfv
 
integer, save iactvx
 
integer, save iactvy
 
integer, save iactvz
 
integer, save iactts
 
integer, save istala
 activation (=1) or not (=0) of the calculation of the volume statistics related to the dispersed phase. if istala = 1, the calculation of the statistics is activated starting from the absolute iteration (including the restarts) idstnt. by default, the statistics are not steady (reset to zero at every Lagrangian iteration). But if isttio=1, since the flow is steady, the statistics will be averaged over the different time steps. the statistics represent the significant results on the particle cloud More...
 
integer, save isuist
 
integer, save nvlsts
 number of additional user volume statistic the additional statistics (or their cumulated value in the steady case) can be accessed in the array statis by means of the pointer ilvu: statis(iel,ilvu(ii)) (iel is the cell index-number and ii an integer between 1 and nvlsts). Useful if istala = 1 More...
 
integer, save idstnt
 absolute Lagrangian iteration number (includings the restarts) after which the calculation of the volume statistics is activated. Useful if istala = 1 More...
 
integer, save nstist
 
integer, save npst
 number of iterations during which steady volume statistics have been cumulated. Useful if istala=1, isttio=1 and if nstist is inferior or equal to the current Lagrangian iteration. npst is initialized and updated automatically by the code, its value is not to be modified by the user More...
 
integer, save npstt
 number of iterations during which volume statistics have been calculated (the potential iterations during which unsteady statistics have been calculated are counted in npstt). Useful if istala=1. npstt is initialized and updated automatically by the code, its value is not to be modified by the user More...
 
double precision, save tstat
 if the volume statistics are calculated in a steady way, tstat represents the physical time during which the statistics have been cumulated. if the volume statistics are calculated in a unsteady way, then tstat=dtp (it is the Lagrangian time step, because the statistics are reset to zero at every iteration). Useful if istala=1. tstat is initialized and updated automatically by the code, its value is not to be modified by the user More...
 
double precision, save seuil
 every cell of the calculation domain contains a certain quantity of particles, representing a certain statistical weight (sum of the statistical weights of all the particles present in the cell). seuil is the limit statistical weight value, below which the contribution of the cell in term of statistical weight is not taken into account in the volume statistics (for the complete turbulent dispersion model, in the Poisson's equation used to correct the mean velocities or in the listing and post-processing outputs). Useful if istala = 1 More...
 
character(len=32), dimension(nvplmx), save nomlag
 name of the volumetric statistics, displayed in the listing and the post-processing files. The default value is given above, with "XXXX" representing a four digit number (for instance 0001, 0011 ...). Useful if istala = 1. Warning: this name is also used to reference information in the restart file (isuist =1). If the name of a variable is changed between two calculations, it will not be possible to read its value from the restart file More...
 
character(len=32), dimension(nvplmx), save nomlav
 
integer, dimension(nvplmx), save ihslag
 historic statistics options More...
 
integer, save nbclst
 statistic per zone and per class More...
 
integer, save ltsdyn
 activation (=1) or not (=0) of the two-way coupling on the dynamics of the continuous phase. Useful if iilagr = 2 and iccvfg = 0 More...
 
integer, save ltsmas
 activation (=1) or not (=0) of the two-way coupling on the mass. Useful if iilagr = 2, iphyla = 1 and impvar = 1 More...
 
integer, save ltsthe
 if iphyla = 1 and itpvar = 1, ltsthe activates (=1) or not (=0) the two-way coupling on temperature. if iphyla = 2, ltsthe activates (=1) or not (=0) the two-way coupling on the eulerian variables related to pulverised coal combustion. Useful if iilagr = 2 More...
 
integer, save itsvx
 explicit source term for the continuous phase X velocity More...
 
integer, save itsvy
 explicit source term for the continuous phase Y velocity More...
 
integer, save itsvz
 explicit source term for the continuous phase Z velocity More...
 
integer, save itsli
 implicit source term for the continuous phase velocity and for the turbulent energy if the $k-\varepsilon$ model is used More...
 
integer, save itske
 explicit source term for the turbulent dissipation and the turbulent energy if the $k-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itsr11
 source term for the Reynolds stress and the turbulent dissipation if the $R_{ij}-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itsr12
 source term for the Reynolds stress and the turbulent dissipation if the $R_{ij}-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itsr13
 source term for the Reynolds stress and the turbulent dissipation if the $R_{ij}-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itsr22
 source term for the Reynolds stress and the turbulent dissipation if the $R_{ij}-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itsr23
 source term for the Reynolds stress and the turbulent dissipation if the $R_{ij}-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itsr33
 source term for the Reynolds stress and the turbulent dissipation if the $R_{ij}-\varepsilon$ turbulence model is used for the continuous phase More...
 
integer, save itste
 explicit thermal source term for the thermal scalar of the continuous phase More...
 
integer, save itsti
 implicit thermal source term for the thermal scalar of the continuous phase More...
 
integer, save itsmas
 mass source term More...
 
integer, dimension(ncharm2), save itsmv1
 source term for the light volatile matters More...
 
integer, dimension(ncharm2), save itsmv2
 source term for the heavy volatile matters More...
 
integer, save itsco
 source term for the carbon released during heterogeneous combustion More...
 
integer, save itsfp4
 Variance of the air scalar. More...
 
integer, save nstits
 
integer, save npts
 number of time steps for source terms accumulations More...
 
integer, save ntxerr
 nomber of cells, whose vulumetric rate DODO (concentration ?)is greather than 0.8 More...
 
double precision, save vmax
 maximum volumetric concentration reached More...
 
double precision, save tmamax
 maximum massic concentration reached More...
 
integer, save iroule
 activates (=1) or not (=0) of the particle cloning/fusion technique (option also called "Russian roulette"). when iroule = 1, the importance function must be specified via the array croule in the user subroutine uslaru More...
 
integer, save iencra
 activates (=1) or not (=0) the option of coal particle fouling. It then is necessary to specify the domain boundaries on which fouling may take place. Useful if iphyla = 2 More...
 
integer(c_int), pointer, save npencr
 encrustation data More...
 
double precision, dimension(ncharm2), save enc1
 encrustation data More...
 
double precision, dimension(ncharm2), save enc2
 encrustation data More...
 
double precision, dimension(ncharm2), save tprenc
 limit temperature (in degree Celsius) below which the coal particles do not cause any fouling (if the fouling model is activated). Useful if iphyla = 2 and iencra = 1 More...
 
double precision, dimension(ncharm2), save visref
 ash critical viscosity in $ kg.m^{-1}.s^{-1} $, in the fouling model cf J.D. Watt et T. Fereday (J.Inst.Fuel, Vol.42-p99). Useful if iphyla = 2 and iencra = 1 More...
 
real(c_double), pointer, save dnpenc
 encrustation data More...
 
double precision, save cstham
 Hamaker constant for the particle/fluid/substrate system. More...
 
double precision, save epseau
 Dielectric constant of the fluid. More...
 
double precision, save phi_p
 Electrokinetic potential of the first solid - particle. More...
 
double precision, save phi_s
 Electrokinetic potential of the second solid - surface. More...
 
double precision, save valen
 Valency of ions in the solution (used for EDL forces) More...
 
double precision, save fion
 Ionic force. More...
 
double precision cstfar
 Faraday constant (C/mol) More...
 
double precision epsvid
 Vacuum permittivity (F/m) More...
 
double precision kboltz
 Boltzmann constant (J/K) More...
 
double precision dcutof
 Cut-off distance for adhesion forces (assumed to be the Born distance) (m) More...
 
double precision lambwl
 Characteristic retardation wavelength (m) for Hamaker constant. More...
 
integer, save lamvbr
 brownnian motion activation More...
 
integer, save nor
 number of lagrangian under step (1 or 2) More...
 
integer, save nordre
 order of integration for the stochastic differential equations More...
 
integer, save modcpl
 activates (>0) or not (=0) the complete turbulent dispersion model. When modcpl is strictly positive, its value is interpreted as the absolute Lagrangian time step number (including restarts) after which the complete model is applied. Since the complete model uses volume statistics, modcpl must either be 0 or be larger than idstnt. Useful if istala = 1 More...
 
integer, save idirla
 direction (1=x, 2=y, 3=z) of the complete model. it corresponds to the main directions of the flow. Useful if modcpl > 0 More...
 
integer, save idistu
 activation (=1) or not (=0) of the particle turbulent dispersion. The turbulent dispersion is compatible only with the RANS turbulent models ( $k-\varepsilon$, $R_{ij}-\varepsilon$, v2f or $k-\omega$). (iturb=20, 21, 30, 31, 50 or 60). More...
 
integer, save idiffl
 idiffl=1 suppresses the crossing trajectory effect, making turbulent dispersion for the particles identical to the turbulent diffusion of fluid particles. Useful if idistu=1 More...
 
integer, save ilapoi
 activation (=1) or not (=0) of the solution of a Poisson's equation for the correction of the particle instantaneous velocities (in order to obtain a null divergence). this option is not validated and reserved to the development team. Do not change the default value More...
 
integer, save iadded_mass
 activation (=1) or not (=0) of the added-mass term.

\[ \DP{u_p} = - \dfrac{1}{\rho_p} \grad P + \dfrac{u_s-u_p}{\tau_p} + g +1/2 C_A \dfrac{\rho_f}{\rho_p} \left( \dfrac{Du}{Dt}-\DP{u_p} \right) \]

and

\[ \rho_f \dfrac{Du}{Dt} \simeq - \grad P + \rho_f g \]

with $ C_A = 1$. Then

\[ \DP{u_p} = - \dfrac{1}{\rho_p} \dfrac{1+C_A/2} {1+C_A/2\dfrac{\rho_f}{\rho_p}} \grad P + \dfrac{u_s-u_p}{\widetilde{\tau}_p} + g \]

with

\[ \widetilde{\tau_p} = (1 + C_A /2 \dfrac{\rho_f}{\rho_p}) \tau_p \]

More...
 
double precision, save added_mass_const
 Added-mass constant ( $ C_A = 1$) More...
 
integer, save nusbor
 number additional user data to record for the calculation of additional boundary statistics in parbor. Useful if iensi3=1 More...
 
integer, save nstbor
 
integer, save npstf
 number of iterations during which steady boundary statistics have been cumulated. Useful if iensi3=1, isttio=1 and nstbor inferior or equal to the current Lagrangian iteration. npstf is initialized and updated automatically by the code, its value is not to be modified by the user More...
 
integer, save npstft
 number of iterations during which boundary statistics have been calculated (the potential iterations during which unsteady statistics have been calculated are counted in npstft). Useful if iensi3=1. npstft is initialized and updated automatically by the code, its value is not to be modified by the user More...
 
integer, save inbrbd
 activation (=1) or not (=0) of the recording of the number of particle/boundary interactions, and of the calculation of the associated boundary statistics. inbrbd = 1 is a compulsory condition to use the particulate average imoybr = 2. Useful if iensi3=1 More...
 
integer, save iflmbd
 activation (=1) or not (=0) of the recording of the particulate mass flow related to the particle/boundary interactions, and of the calculation of the associated boundary statistics. inbrbd = 1 is a compulsory condition to use iflmbd=1. Useful if iensi3=1 and inbrbd=1 More...
 
integer, save iangbd
 activation (=1) or not (=0) of the recording of the angle between a particle trajectory and a boundary face involved in a particle/boundary interaction, and of the calculation of the associated boundary statistics. Useful if iensi3=1 More...
 
integer, save ivitbd
 activation (=1) or not (=0) of the recording of the velocity of a particle involved in a particle/boundary interaction, and of the calculation of the associated boundary statistics. Useful if iensi3=1 More...
 
integer, save iclgst
 activation (=1) or not (=0) of the recording of clogging parameters involved in a particle/boundary interaction, and of the calculation of the associated boundary statistics. Useful if iensi3=1 More...
 
integer, save iencnbbd
 
integer, save iencmabd
 
integer, save iencdibd
 
integer, save iencckbd
 
integer, save inbr
 number of particle/boundary interactions More...
 
integer, save iflm
 particle mass flow at the boundary faces More...
 
integer, save iang
 mean interaction angle with the boundary faces More...
 
integer, save ivit
 mean interaction velocity with the boundary faces More...
 
integer, save ires
 number of resuspended particles More...
 
integer, save iflres
 mass flow of resuspended particles at the boundary faces More...
 
integer, save iencnb
 
integer, save iencma
 
integer, save iencdi
 
integer, save iencck
 
integer, dimension(nusbrd), save iusb
 supplementary user boundary statistics More...
 
integer, dimension(nusbrd+10), save imoybr
 the recordings in parbor at every particle/boundary interaction are cumulated values (possibly reset to zero at every iteration in the unsteady case). They must therefore be divided by a quantity to get boundary statistics. The user can choose between two average types: More...
 
integer, save inclg
 
integer, save inclgt
 
integer, save iclogt
 
integer, save iclogh
 
integer, save iscovc
 
integer, save ihdepm
 
integer, save ihdepv
 
integer, save ihdiam
 
integer, save ihsum
 
double precision, save tstatp
 if the recording of the boundary statistics is steady, tstatp contains the cumulated physical duration of the recording of the boundary statistics. if the recording of the boundary statisticss is unsteady, then tstat=dtp (it is the Lagrangian time step, because the statistics are reset to zero at every time step). Useful if iensi3=1 More...
 
double precision, save seuilf
 every boundary face of the mesh undergoes a certain number of interactions with particles, expressed in term of statistical weight (sum of the statistical weights of all the particles which have interacted with the boundary face). seuilf is the limit statistical weight value, below which the contribution of the face is not taken into account in the statistics at the boundaries for post-processing. Useful if iensi3=1 More...
 
character(len=50), dimension(nvplmx), save nombrd
 name of the boundary statistics, displayed in the listing and the post-processing files. Useful if iensi3=1. Warning: this name is also used to reference information in the restart file (isuist =1). If the name of a variable is changed between two calculations, it will not be possible to read its value from the restart file More...
 
integer, save iifrla
 ifrlag pointer in ia to identify boundary zones and faces More...
 
integer, save iensi3
 
integer, save ivisv1
 associates (=1) or not (=0) the variable "velocity of the locally undisturbed fluid flow field" with the output of particles or trajectories. More...
 
integer, save ivisv2
 associates (=1) or not (=0) the variable "particle velocity" with the output of particles or trajectories. More...
 
integer, save ivistp
 associates (=1) or not (=0) the variable "residence time" with the output of particles or trajectories. More...
 
integer, save ivisdm
 associates (=1) or not (=0) the variable "particle diameter" with the output of particles or trajectories. More...
 
integer, save iviste
 associates (=1) or not (=0) the variable "particle temperature" with the output of particles or trajectories. More...
 
integer, save ivismp
 associates (=1) or not (=0) the variable "particle mass" with the output of particles or trajectories. More...
 
integer, save ivisdk
 associates (=1) or not (=0) the variable "shrinking core diameter of the coal particles" with the output of particles or trajectories. useful only if iphyla = 2 More...
 
integer, save ivisch
 associates (=1) or not (=0) the variable "mass of reactive coal of the coal particles" with the output of particles or trajectories. useful only if iphyla = 2 More...
 
integer, save ivisck
 associates (=1) or not (=0) the variable "mass of coal of the coal particles" with the output of particles or trajectories. useful only if iphyla = 2 More...
 
integer, save iviswat
 

Detailed Description

Module for Lagrangian model.