Boundary conditions

| programmer's documentation

- Main Page
- Related Pages
- Modules
- Data Structures
- Files

All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Variables
Boundary conditions
Module for Lagrangian model | General module for specific physics
Collaboration diagram for Boundary conditions:

## Variables
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 nclagm (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. These pieces of information 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. These pieces of information 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, dimension(nozppm), save | iqimp
| imposed flow zone indicator in a way which is similar to the process described in the framework of the EBU module, the user chooses for every inlet face to impose the mass flow or not (iqimp(izone)=1 or 0). If the mass flow is imposed, the user must set the air mass flow value qimpat(izone), its direction in rcodcl(ifac,iu), rcodcl(ifac,iv) and rcodcl(ifac,iw) and the incoming air temperature timpat(izone) in Kelvin. If the velocity is imposed, he has to set rcodcl(ifac,iu), rcodcl(ifac,iv) and rcodcl(ifac,iw). More...
integer, dimension(nozppm), save | icalke
| condition type turbulence indicator More...
double precision, dimension(nozppm), save | xintur
| turbulent intensity (k=1.5(uref*xintur)**2) More...
double precision, dimension(nozppm), save | dh
| hydraulic diameter More...
integer, save | nozapm
| index of maximum reached boundary zone More...
integer, save | nzfppp
| number of boundary zones on current process More...
integer, dimension(nbzppm), save | ilzppp
| list of bondary zones index More...

## Detailed Description

## Variable Documentation

double precision, dimension(nflagm), save deblag

massic flow rate for a boudary zone

double precision, dimension(nozppm), save dh

hydraulic diameter

integer, dimension(nozppm), save icalke

condition type turbulence indicator

- 0 : given by the user
- 1 : automatic, from hydraulic diameter and input velocity performed.
- 2 : automatic, from turbulent intensity and input velocity performed.

integer, save iclst

number of the statistics group

integer, save icpt

particle specific heat

integer, save idebt

flow rate

integer, save idpt

particle diameter

integer, save iepsi

particle emissivity

integer, dimension(nlayer), save ifrmch

active coal mass fraction in coal particles

integer, dimension(nlayer), save ifrmck

coke mass fraction in coal particles

integer, save ifrmwt

water mass fraction in coal particles

integer, dimension(nlayer), save ihpt

particle temperature

integer, save ijfre

injection frequency (if < 0 : particle are introduced only at first iteration

integer, save ijnbp

number of particles per class and per boudary zone

integer, save ijprdp

type of user profiles in uslag2:

- 1: flat profile of diameter given in uslag2
- 2: user profile to be given

integer, save ijprpd

- 1 uniform distribution,
- 2 presence rate profile given by user.

integer, save ijprtp

- 1 constant temperature profile given in uslag2
- 2 temperature profile given by the user

integer, save ijuvw

velocity condition type:

- -1 imposed fluid velocity
- 0 imposed fluid velocity along the normal of the bondary face, with iuno norm.
- 1 imposed velocity: iupt ivpt iwpt must be given.
- 2 velocity profile given by user.

integer, dimension(nflagm), save ilflag

list of number of boundary zones

integer, dimension(nbzppm), save ilzppp

list of bondary zones index

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.

integer, save inuchl

coal number of the particle (if iphyla=2)

integer, save ipoit

particle weight

integer, dimension(nozppm), save iqimp

imposed flow zone indicator in a way which is similar to the process described in the framework of the EBU module, the user chooses for every inlet face to impose the mass flow or not (iqimp(izone)=1 or 0). If the mass flow is imposed, the user must set the air mass flow value qimpat(izone), its direction in rcodcl(ifac,iu), rcodcl(ifac,iv) and rcodcl(ifac,iw) and the incoming air temperature timpat(izone) in Kelvin. If the velocity is imposed, he has to set rcodcl(ifac,iu), rcodcl(ifac,iv) and rcodcl(ifac,iw).

integer, save irawcl

type of coal initial composition (if iphyla=2)

- 1: coal initial composition is given by DP_FCP
- 0: user profile to be given in uslag2

integer, save ird0p

initial particle diameter (for coal particles)

integer, save irdck

diameter of shrinking core

integer, dimension(nlayer), save irhock0

coke density after pyrolysis (for coal particles)

integer, save iropt

density

integer, save itpt

particle temperature

integer, save iuno

particle velocity magnitude

integer, save iupt

particle u component by class and zone

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.

- if iusclb(izone) = ientrl, izone is a particle injection zone. For each particle class associated with this zone, information must be provided (see below). If a particle trajectory may cross an injection zone, then this particle leaves the calculation domain. - if iusclb(izone) = isortl, the particles interacting with the zone izone leave definitely the calculation domain. - if iusclb(izone) = irebol, the particles undergo an elastic rebound on the boundary zone izone. - if iusclb(izone) = idepo1, the particles settle definitely on the boundary zone izone. These particles leave the calculation domain and are definitely erased from the calculation - if iusclb(izone) = idepo2, the particles settle definitevely on the boundary zone izone and they are kept in the calculation
- if iusclb(izone) = idepo2, the particles settle definitevely on the boundary zone izone and they are kept in the calculation domain: the particles do not disappear after touching the boundary zone. However, using idepo2 type zones necessitates more memory than using idepo1 type zones.
- if iusclb(izone) = iencrl, the particles which are coal particles (if iphyla = 2) can become fouled up on the zone izone. The slagging is a idepo1 type deposit of the coal particle if a certain criterion is respected. Otherwise, the coal particle rebounds (irebol type behaviour). This boundary condition type is available if iencra = 1. A limit temperature tprenc, a critical viscosity visref and the coal composition in mineral matters must be given in the subroutine uslag1.

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. These pieces of information must be provided for each class iclas and each particle injection zone izone. They are marked out by means of "pointers":

- iuslag(iclas,izone,ijnbp): number of particles to inject in the calculation domain per class and per zone.
- iuslag(iclas,izone,ijfre): injection period (expressed in number of time steps). If the period is null, then there is injection only at the first absolute Lagrangian time step (including the restart calculations).
- iuslag(iclas,izone,ijuvw): type of velocity condition:
- if iuslag(iclas,izone,ijuvw) = 1, the particle velocity vector is imposed, and its components must be given in the array ruslag (see below).
- if iuslag(iclas,izone,ijuvw) = 0, the particle velocity is imposed perpendicular to the injection boundary face and with the norm ruslag(iclas,izone,iuno).
- if iuslag(iclas,izone,ijuvw) = -1, the particle injection velocity is equal to the fluid velocity at the center of the cell neighbouring the injection boundary face.
- iuslag(iclas,izone,inuchl): when the particles are coal particles (iphyla = 2), this part of the array contains the coal index-number, between 1 and ncharb (defined by the user in the thermochemical file dp_FCP, with ncharb <= ncharm = 3).

integer, dimension(nflagm), save iusmoy

mean over a zone (if mean per zones is activated)

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 nclagm (parameter stored in lagpar, whose default value is 20).

integer, save ivdpt

particle diameter variance

integer, save ivpt

particle v component by class and zone

integer, save iwpt

particle w component by class and zone

integer, save nfrlag

number of boundary zones

integer, save nozapm

index of maximum reached boundary zone

integer, save nzfppp

number of boundary zones on current process

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. These pieces of information must be provided for each class iclas and each particle injection zone izone. They are marked out by means of "pointers":

double precision, dimension(nozppm), save xintur

turbulent intensity (k=1.5(uref*xintur)**2)