!------------------------------------------------------------------------------- ! Code_Saturne version 2.0.0-rc2 ! -------------------------- ! This file is part of the Code_Saturne Kernel, element of the ! Code_Saturne CFD tool. ! Copyright (C) 1998-2009 EDF S.A., France ! contact: saturne-support@edf.fr ! The Code_Saturne Kernel is free software; you can redistribute it ! and/or modify it under the terms of the GNU General Public License ! as published by the Free Software Foundation; either version 2 of ! the License, or (at your option) any later version. ! The Code_Saturne Kernel is distributed in the hope that it will be ! useful, but WITHOUT ANY WARRANTY; without even the implied warranty ! of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the ! GNU General Public License for more details. ! You should have received a copy of the GNU General Public License ! along with the Code_Saturne Kernel; if not, write to the ! Free Software Foundation, Inc., ! 51 Franklin St, Fifth Floor, ! Boston, MA 02110-1301 USA !------------------------------------------------------------------------------- subroutine uscpcl & !================ ( idbia0 , idbra0 , & ndim , ncelet , ncel , nfac , nfabor , nfml , nprfml , & nnod , lndfac , lndfbr , ncelbr , & nvar , nscal , nphas , & nideve , nrdeve , nituse , nrtuse , & ifacel , ifabor , ifmfbr , ifmcel , iprfml , maxelt , lstelt , & ipnfac , nodfac , ipnfbr , nodfbr , & icodcl , itrifb , itypfb , izfppp , & idevel , ituser , ia , & xyzcen , surfac , surfbo , cdgfac , cdgfbo , xyznod , volume , & dt , rtp , rtpa , propce , propfa , propfb , & coefa , coefb , rcodcl , & w1 , w2 , w3 , w4 , w5 , w6 , coefu , & rdevel , rtuser , ra ) !=============================================================================== ! PURPOSE : ! -------- ! USER'S ROUTINE FOR EXTENDED PHYSIC ! PULVERISED COAL COMBUSTION ! ALLOCATION OF BOUNDARY CONDITIONS (ICODCL,RCODCL) ! FOR VARIABLES UNKNOWN DURING USCLIM ! ! Introduction ! ============ ! Here we define boundary conditions on a per-face basis. ! Boundary faces may be identified using the 'getfbr' subroutine. ! getfbr(string, nelts, eltlst) : ! - string is a user-supplied character string containing ! selection criteria; ! - nelts is set by the subroutine. It is an integer value ! corresponding to the number of boundary faces verifying the ! selection criteria; ! - lstelt is set by the subroutine. It is an integer array of ! size nelts containing the list of boundary faces verifying ! the selection criteria. ! string may contain: ! - references to colors (ex.: 1, 8, 26, ... ! - references to groups (ex.: inlet, group1, ...) ! - geometric criteria (ex. x < 0.1, y >= 0.25, ...) ! These criteria may be combined using logical operators ! ('and', 'or') and parentheses. ! Example: '1 and (group2 or group3) and y < 1' will select boundary ! faces of color 1, belonging to groups 'group2' or 'group3' and ! with face center coordinate y less than 1. ! Boundary condition types ! ======================== ! Boundary conditions may be assigned in two ways. ! For "standard" boundary conditions: ! ----------------------------------- ! (inlet, free outlet, wall, symmetry), we define a code ! in the 'itypfb' array (of dimensions number of boundary faces, ! number of phases). This code will then be used by a non-user ! subroutine to assign the following conditions (scalars in ! particular will receive the conditions of the phase to which ! they are assigned). Thus: ! Code | Boundary type ! -------------------------- ! ientre | Inlet ! isolib | Free outlet ! isymet | Symmetry ! iparoi | Wall (smooth) ! iparug | Rough wall ! Integers ientre, isolib, isymet, iparoi, iparug ! are defined elsewhere (param.h). Their value is greater than ! or equal to 1 and less than or equal to ntypmx ! (value fixed in paramx.h) ! In addition, some values must be defined: ! - Inlet (more precisely, inlet/outlet with prescribed flow, as ! the flow may be prescribed as an outflow): ! -> Dirichlet conditions on variables ! other than pressure are mandatory if the flow is incoming, ! optional if the flow is outgoing (the code assigns 0 flux ! if no Dirichlet is specified); thus, ! at face 'ifac', for the variable 'ivar': rcodcl(ifac, ivar, 1) ! - Smooth wall: (= impermeable solid, with smooth friction) ! -> Velocity value for sliding wall if applicable ! at face ifac, rcodcl(ifac, iu, 1) ! rcodcl(ifac, iv, 1) ! rcodcl(ifac, iw, 1) ! -> Specific code and prescribed temperature value ! at wall, if applicable: ! at face ifac, icodcl(ifac, ivar) = 5 ! rcodcl(ifac, ivar, 1) = prescribed temperature ! -> Specific code and prescribed flux value ! at wall, if applicable: ! at face ifac, icodcl(ifac, ivar) = 3 ! rcodcl(ifac, ivar, 3) = prescribed flux ! = ! Note that the default condition for scalars ! (other than k and epsilon) is homogeneous Neumann. ! - Rough wall: (= impermeable solid, with rough friction) ! -> Velocity value for sliding wall if applicable ! at face ifac, rcodcl(ifac, iu, 1) ! rcodcl(ifac, iv, 1) ! rcodcl(ifac, iw, 1) ! -> Value of the dynamic roughness height to specify in ! rcodcl(ifac, iu, 3) (value for iv et iw not used) ! -> Specific code and prescribed temperature value ! at rough wall, if applicable: ! at face ifac, icodcl(ifac, ivar) = 6 ! rcodcl(ifac, ivar, 1) = prescribed temperature ! rcodcl(ifac, ivar, 3) = dynamic roughness height ! -> Specific code and prescribed flux value ! at rough wall, if applicable: ! at face ifac, icodcl(ifac, ivar) = 3 ! rcodcl(ifac, ivar, 3) = prescribed flux ! = ! Note that the default condition for scalars ! (other than k and epsilon) is homogeneous Neumann. ! - Symmetry (= impermeable frictionless wall): ! -> Nothing to specify ! - Free outlet (more precisely free inlet/outlet with prescribed pressure) ! -> Nothing to prescribe for pressure and velocity ! For scalars and turbulent values, a Dirichlet value may optionally ! be specified. The behavior is as follows: ! * pressure is always handled as a Dirichlet condition ! * if the mass flow is inflowing: ! we retain the velocity at infinity ! Dirichlet condition for scalars and turbulent values ! (or zero flux if the user has not specified a ! Dirichlet value) ! if the mass flow is outflowing: ! we prescribe zero flux on the velocity, the scalars, ! and turbulent values ! Note that the pressure will be reset to P0 ! on the first free outlet face found ! For "non-standard" conditions: ! ------------------------------ ! Other than (inlet, free outlet, wall, symmetry), we define ! - on one hand, for each face: ! -> an admissible 'itypfb' value ! (i.e. greater than or equal to 1 and less than or equal to ! ntypmx; see its value in paramx.h). ! The values predefined in paramx.h: ! 'ientre', 'isolib', 'isymet', 'iparoi', 'iparug' are in ! this range, and it is preferable not to assign one of these ! integers to 'itypfb' randomly or in an inconsiderate manner. ! To avoid this, we may use 'iindef' if we wish to avoid ! checking values in paramx.h. 'iindef' is an admissible ! value to which no predefined boundary condition is attached. ! Note that the 'itypfb' array is reinitialized at each time ! step to the non-admissible value of 0. If we forget to ! modify 'typfb' for a given face, the code will stop. ! - and on the other hand, for each face and each variable: ! -> a code icodcl(ifac, ivar) ! -> three real values rcodcl(ifac, ivar, 1) ! rcodcl(ifac, ivar, 2) ! rcodcl(ifac, ivar, 3) ! The value of 'icodcl' is taken from the following: ! 1: Dirichlet (usable for any variable) ! 3: Neumann (usable for any variable) ! 4: Symmetry (usable only for the velocity and ! components of the Rij tensor) ! 5: Smooth wall (usable for any variable except for pressure) ! 6: Rough wall (usable for any variable except for pressure) ! 9: Free outlet (usable only for velocity) ! The values of the 3 'rcodcl' components are ! rcodcl(ifac, ivar, 1): ! Dirichlet for the variable if icodcl(ifac, ivar) = 1 ! wall value (sliding velocity, temp) if icodcl(ifac, ivar) = 5 ! The dimension of rcodcl(ifac, ivar, 1) is that of the ! resolved variable: ex U (velocity in m/s), ! T (temperature in degrees) ! H (enthalpy in J/kg) ! F (passive scalar in -) ! rcodcl(ifac, ivar, 2): ! "exterior" exchange coefficient (between the prescribed value ! and the value at the domain boundary) ! rinfin = infinite by default ! For velocities U, in kg/(m2 s): ! rcodcl(ifac, ivar, 2) = (viscl+visct) / d ! For the pressure P, in s/m: ! rcodcl(ifac, ivar, 2) = dt / d ! For temperatures T, in Watt/(m2 degres): ! rcodcl(ifac, ivar, 2) = Cp*(viscls+visct/sigmas) / d ! For enthalpies H, in kg /(m2 s): ! rcodcl(ifac, ivar, 2) = (viscls+visct/sigmas) / d ! For other scalars F in: ! rcodcl(ifac, ivar, 2) = (viscls+visct/sigmas) / d ! (d has the dimension of a distance in m) ! ! rcodcl(ifac, ivar, 3) if icodcl(ifac, ivar) <> 6: ! Flux density (< 0 if gain, n outwards-facing normal) ! if icodcl(ifac, ivar)= 3 ! For velocities U, in kg/(m s2) = J: ! rcodcl(ifac, ivar, 3) = -(viscl+visct) * (grad U).n ! For pressure P, en kg/(m2 s): ! rcodcl(ifac, ivar, 3) = -dt * (grad P).n ! For temperatures T, in Watt/m2: ! rcodcl(ifac, ivar, 3) = -Cp*(viscls+visct/sigmas) * (grad T).n ! For enthalpies H, in Watt/m2: ! rcodcl(ifac, ivar, 3) = -(viscls+visct/sigmas) * (grad H).n ! For other scalars F in : ! rcodcl(ifac, ivar, 3) = -(viscls+visct/sigmas) * (grad F).n ! rcodcl(ifac, ivar, 3) if icodcl(ifac, ivar) = 6: ! Roughness for the rough wall law ! For velocities U, dynamic roughness ! rcodcl(ifac, ivar, 3) = rugd ! For other scalars, thermal roughness ! rcodcl(ifac, ivar, 3) = rugt ! Note that if the user assigns a value to itypfb equal to ! ientre, isolib, isymet, iparoi, or iparug ! and does not modify icodcl (zero value by default), ! itypfb will define the boundary condition type. ! To the contrary, if the user prescribes ! icodcl(ifac, ivar) (nonzero), ! the values assigned to rcodcl will be used for the considered ! face and variable (if rcodcl values are not set, the default ! values will be used for the face and variable, so: ! rcodcl(ifac, ivar, 1) = 0.d0 ! rcodcl(ifac, ivar, 2) = rinfin ! rcodcl(ifac, ivar, 3) = 0.d0) ! Especially, we may have for example: ! -> set itypfb(ifac, iphas) = iparoi ! which prescribes default wall conditions for all variables at ! face ifac, ! -> and define IN ADDITION for variable ivar on this face ! specific conditions by specifying ! icodcl(ifac, ivar) and the 3 rcodcl values. ! The user may also assign to itypfb a value not equal to ! ientre, isolib, isymet, iparoi, iparug, iindef ! but greater than or equal to 1 and less than or equal to ! ntypmx (see values in param.h) to distinguish ! groups or colors in other subroutines which are specific ! to the case and in which itypfb is accessible. ! In this case though it will be necessary to ! prescribe boundary conditions by assigning values to ! icodcl and to the 3 rcodcl fields (as the value of itypfb ! will not be predefined in the code). ! Consistency rules ! ================= ! A few consistency rules between 'icodcl' codes for ! variables with non-standard boundary conditions: ! Codes for velocity components must be identical ! Codes for Rij components must be identical ! If code (velocity or Rij) = 4 ! we must have code (velocity and Rij) = 4 ! If code (velocity or turbulence) = 5 ! we must have code (velocity and turbulence) = 5 ! If code (velocity or turbulence) = 6 ! we must have code (velocity and turbulence) = 6 ! If scalar code (except pressure or fluctuations) = 5 ! we must have velocity code = 5 ! If scalar code (except pressure or fluctuations) = 6 ! we must have velocity code = 6 ! Remarks ! ======= ! Caution: to prescribe a flux (nonzero) to Rij, ! the viscosity to take into account is viscl ! even if visct exists (visct=rho cmu k2/epsilon) ! We have the ordering array for boundary faces from the ! previous time step (except for the fist time step, ! where 'itrifb' has not been set yet). ! The array of boundary face types 'itypfb' has been ! reset before entering the subroutine. ! Note how to access some variables: ! Cell values ! Let iel = ifabor(ifac) ! * Density phase iphas, cell iel: ! propce(iel, ipproc(irom(iphas))) ! * Dynamic molecular viscosity phase iphas, cell iel: ! propce(iel, ipproc(iviscl(iphas))) ! * Turbulent viscosity dynamique phase iphas, cell iel: ! propce(iel, ipproc(ivisct(iphas))) ! * Specific heat phase iphas, cell iel: ! propce(iel, ipproc(icp(iphasl)) ! * Diffusivity: lambda scalaire iscal, cell iel: ! propce(iel, ipproc(ivisls(iscal))) ! Boundary face values ! * Density phase iphas, boundary face ifac : ! propfb(ifac, ipprob(irom(iphas))) ! * Mass flow relative to variable ivar, boundary face ifac: ! (i.e. the mass flow used for convecting ivar) ! propfb(ifac, pprob(ifluma(ivar ))) ! * For other values at boundary face ifac: ! take as an approximation the value in the adjacent cell iel ! i.e. as above with iel = ifabor(ifac). !------------------------------------------------------------------------------- ! Arguments !__________________.____._____.________________________________________________. ! name !type!mode ! role ! !__________________!____!_____!________________________________________________! ! idbia0 ! i ! <-- ! number of first free position in ia ! ! idbra0 ! i ! <-- ! number of first free position in ra ! ! ndim ! i ! <-- ! spatial dimension ! ! ncelet ! i ! <-- ! number of extended (real + ghost) cells ! ! ncel ! i ! <-- ! number of cells ! ! nfac ! i ! <-- ! number of interior faces ! ! nfabor ! i ! <-- ! number of boundary faces ! ! nfml ! i ! <-- ! number of families (group classes) ! ! nprfml ! i ! <-- ! number of properties per family (group class) ! ! nnod ! i ! <-- ! number of vertices ! ! lndfac ! i ! <-- ! size of nodfac indexed array ! ! lndfbr ! i ! <-- ! size of nodfbr indexed array ! ! ncelbr ! i ! <-- ! number of cells with faces on boundary ! ! nvar ! i ! <-- ! total number of variables ! ! nscal ! i ! <-- ! total number of scalars ! ! nphas ! i ! <-- ! number of phases ! ! nideve, nrdeve ! i ! <-- ! sizes of idevel and rdevel arrays ! ! nituse, nrtuse ! i ! <-- ! sizes of ituser and rtuser arrays ! ! ifacel(2, nfac) ! ia ! <-- ! interior faces -> cells connectivity ! ! ifabor(nfabor) ! ia ! <-- ! boundary faces -> cells connectivity ! ! ifmfbr(nfabor) ! ia ! <-- ! boundary face family numbers ! ! ifmcel(ncelet) ! ia ! <-- ! cell family numbers ! ! iprfml ! ia ! <-- ! property numbers per family ! ! (nfml, nprfml) ! ! ! ! ! maxelt ! e ! <-- ! max number of cells and faces (int/boundary) ! ! lstelt(maxelt) ! ia ! --- ! work array ! ! ipnfac(nfac+1) ! ia ! <-- ! interior faces -> vertices index (optional) ! ! nodfac(lndfac) ! ia ! <-- ! interior faces -> vertices list (optional) ! ! ipnfbr(nfabor+1) ! ia ! <-- ! boundary faces -> vertices index (optional) ! ! nodfac(lndfbr) ! ia ! <-- ! boundary faces -> vertices list (optional) ! ! icodcl ! ia ! --> ! boundary condition code ! ! (nfabor, nvar) ! ! ! = 1 -> Dirichlet ! ! ! ! ! = 2 -> flux density ! ! ! ! ! = 4 -> sliding wall and u.n=0 (velocity) ! ! ! ! ! = 5 -> friction and u.n=0 (velocity) ! ! ! ! ! = 6 -> roughness and u.n=0 (velocity) ! ! ! ! ! = 9 -> free inlet/outlet (velocity) ! ! ! ! ! inflowing possibly blocked ! ! itrifb(nfabor ! ia ! <-- ! indirection for boundary faces ordering) ! ! (nfabor, nphas) ! ! ! ! ! itypfb ! ia ! --> ! boundary face types ! ! (nfabor, nphas) ! ! ! ! ! idevel(nideve) ! ia ! <-- ! integer work array for temporary developpement ! ! ituser(nituse ! ia ! <-- ! user-reserved integer work array ! ! ia(*) ! ia ! --- ! main integer work array ! ! xyzcen ! ra ! <-- ! cell centers ! ! (ndim, ncelet) ! ! ! ! ! surfac ! ra ! <-- ! interior faces surface vectors ! ! (ndim, nfac) ! ! ! ! ! surfbo ! ra ! <-- ! boundary faces surface vectors ! ! (ndim, nfavor) ! ! ! ! ! cdgfac ! ra ! <-- ! interior faces centers of gravity ! ! (ndim, nfac) ! ! ! ! ! cdgfbo ! ra ! <-- ! boundary faces centers of gravity ! ! (ndim, nfabor) ! ! ! ! ! xyznod ! ra ! <-- ! vertex coordinates (optional) ! ! (ndim, nnod) ! ! ! ! ! volume(ncelet) ! ra ! <-- ! cell volumes ! ! dt(ncelet) ! ra ! <-- ! time step (per cell) ! ! rtp, rtpa ! ra ! <-- ! calculated variables at cell centers ! ! (ncelet, *) ! ! ! (at current and preceding time steps) ! ! propce(ncelet, *)! ra ! <-- ! physical properties at cell centers ! ! propfa(nfac, *) ! ra ! <-- ! physical properties at interior face centers ! ! propfb(nfabor, *)! ra ! <-- ! physical properties at boundary face centers ! ! coefa, coefb ! ra ! <-- ! boundary conditions ! ! (nfabor, *) ! ! ! ! ! rcodcl ! ra ! --> ! boundary condition values ! ! ! ! ! rcodcl(1) = Dirichlet value ! ! ! ! ! rcodcl(2) = exterior exchange coefficient ! ! ! ! ! (infinite if no exchange) ! ! ! ! ! rcodcl(3) = flux density value ! ! ! ! ! (negative for gain) in w/m2 or ! ! ! ! ! roughness height (m) if icodcl=6 ! ! ! ! ! for velocities ( vistl+visct)*gradu ! ! ! ! ! for pressure dt*gradp ! ! ! ! ! for scalars cp*(viscls+visct/sigmas)*gradt ! ! w1,2,3,4,5,6 ! ra ! --- ! work arrays ! ! (ncelet) ! ! ! (computation of pressure gradient) ! ! coefu ! ra ! --- ! tab de trav ! ! (nfabor, 3) ! ! ! (computation of pressure gradient) ! ! rdevel(nrdeve) ! ra ! <-> ! tab reel complementaire developemt ! ! rdevel(nideve) ! ra ! <-- ! real work array for temporary developpement ! ! rtuser(nituse ! ra ! <-- ! user-reserved real work array ! ! ra(*) ! ra ! --- ! main real work array ! !__________________!____!_____!________________________________________________! ! Type: i (integer), r (real), s (string), a (array), l (logical), ! and composite types (ex: ra real array) ! mode: <-- input, --> output, <-> modifies data, --- work array !=============================================================================== implicit none !=============================================================================== ! COMMON OF DATA !=============================================================================== include "ihmpre.h" include "paramx.h" include "pointe.h" include "numvar.h" include "optcal.h" include "cstphy.h" include "cstnum.h" include "entsor.h" include "parall.h" include "period.h" include "ppppar.h" include "ppthch.h" include "coincl.h" include "cpincl.h" include "ppincl.h" !=============================================================================== ! Arguments integer idbia0 , idbra0 integer ndim , ncelet , ncel , nfac , nfabor integer nfml , nprfml integer nnod , lndfac , lndfbr , ncelbr integer nvar , nscal , nphas integer nideve , nrdeve , nituse , nrtuse integer ifacel(2,nfac) , ifabor(nfabor) integer ifmfbr(nfabor) , ifmcel(ncelet) integer iprfml(nfml,nprfml) integer maxelt, lstelt(maxelt) integer ipnfac(nfac+1), nodfac(lndfac) integer ipnfbr(nfabor+1), nodfbr(lndfbr) integer icodcl(nfabor,nvar) integer itrifb(nfabor,nphas), itypfb(nfabor,nphas) integer izfppp(nfabor) integer idevel(nideve), ituser(nituse), ia(*) double precision xyzcen(ndim,ncelet) double precision surfac(ndim,nfac), surfbo(ndim,nfabor) double precision cdgfac(ndim,nfac), cdgfbo(ndim,nfabor) double precision xyznod(ndim,nnod), volume(ncelet) double precision dt(ncelet), rtp(ncelet,*), rtpa(ncelet,*) double precision propce(ncelet,*) double precision propfa(nfac,*), propfb(nfabor,*) double precision coefa(nfabor,*), coefb(nfabor,*) double precision rcodcl(nfabor,nvar,3) double precision w1(ncelet),w2(ncelet),w3(ncelet) double precision w4(ncelet),w5(ncelet),w6(ncelet) double precision coefu(nfabor,ndim) double precision rdevel(nrdeve), rtuser(nrtuse), ra(*) ! LOCAL VARIABLES integer idebia, idebra integer ifac, iphas, ii integer izone integer icha, iclapc integer ilelt, nlelt double precision uref2, d2s3 double precision xkent, xeent !=============================================================================== !=============================================================================== ! 1. INITIALISATIONS !=============================================================================== idebia = idbia0 idebra = idbra0 d2s3 = 2.d0/3.d0 !=============================================================================== ! 2. ALLOCATION OF BOUNDARY CONDITIONS VECTOR ! LOOP ON BOUNDARY FACE ! FAMILY AND PROPERTIES ARE DETERMINED ! BOUNDARY CONDITION ALLOCATED ! ! BOUNDARY CONDITIONS ON BONDARY FACE HAVE TO BE ALLOCATED HERE ! ! USER'S WORK TO DO ! !=============================================================================== iphas = 1 ! ---- BOUNDARY FACE corresponding to AIR INLET ! e.g. : secondary or tertiary air CALL GETFBR('AIR1.1A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 1 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.576d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) !rcodcl(ifac,iu(iphas),1) = 3.83d0 !rcodcl(ifac,iv(iphas),1) = 3.214d0 !rcodcl(ifac,iw(iphas),1) = 0.d0 do ifac = 1,nfac rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) rcodcl(ifac, iw(iphas), 1) = -surfbo(3, ifac) enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.4d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR1.2A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 2 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.576d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -3.214d0 rcodcl(ifac,iv(iphas),1) = 3.83d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.4d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR1.3A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 3 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.576d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -3.83d0 rcodcl(ifac,iv(iphas),1) = -3.214d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.4d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR1.4A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 4 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.576d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 3.214d0 rcodcl(ifac,iv(iphas),1) = -3.83d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.4d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.1A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 5 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 3.83d0 rcodcl(ifac,iv(iphas),1) = 3.214d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.1A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 6 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 3.83d0 rcodcl(ifac,iv(iphas),1) = 3.214d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.1B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 7 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 30.642d0 rcodcl(ifac,iv(iphas),1) = 25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.1B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 8 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 30.642d0 rcodcl(ifac,iv(iphas),1) = 25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.1C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 9 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 30.642d0 rcodcl(ifac,iv(iphas),1) = 25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.1C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 10 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 30.642d0 rcodcl(ifac,iv(iphas),1) = 25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.1D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 11 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 30.642d0 rcodcl(ifac,iv(iphas),1) = 25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.1D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 12 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 30.642d0 rcodcl(ifac,iv(iphas),1) = 25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('BNOFA.1',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 13 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.863d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 19.151d0 rcodcl(ifac,iv(iphas),1) = 16.07d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.2d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.2A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 14 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -3.214d0 rcodcl(ifac,iv(iphas),1) = 3.83d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.2A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 15 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -3.214d0 rcodcl(ifac,iv(iphas),1) = 3.83d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.2B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 16 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -25.712d0 rcodcl(ifac,iv(iphas),1) = 30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.2B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 17 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -25.712d0 rcodcl(ifac,iv(iphas),1) = 30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.2C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 18 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -25.712d0 rcodcl(ifac,iv(iphas),1) = 30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.2C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 19 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -25.712d0 rcodcl(ifac,iv(iphas),1) = 30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.2D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 20 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -25.712d0 rcodcl(ifac,iv(iphas),1) = 30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.2D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 21 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -25.712d0 rcodcl(ifac,iv(iphas),1) = 30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('BNOFA.2',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 22 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.863d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -16.07d0 rcodcl(ifac,iv(iphas),1) = 19.252d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.2d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.3A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 23 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -3.83d0 rcodcl(ifac,iv(iphas),1) = -3.214d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.3A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 24 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -3.83d0 rcodcl(ifac,iv(iphas),1) = -3.214d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.3B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 25 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -30.642d0 rcodcl(ifac,iv(iphas),1) = -25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.3B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 26 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -30.642d0 rcodcl(ifac,iv(iphas),1) = -25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.3C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 27 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -30.642d0 rcodcl(ifac,iv(iphas),1) = -25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.3C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 28 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -30.642d0 rcodcl(ifac,iv(iphas),1) = -25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.3D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 29 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -30.642d0 rcodcl(ifac,iv(iphas),1) = -25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.3D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 30 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -30.642d0 rcodcl(ifac,iv(iphas),1) = -25.712d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('BNOFA.3',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 31 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.863d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -19.151d0 rcodcl(ifac,iv(iphas),1) = -16.07d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.2d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.4A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 32 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 3.214d0 rcodcl(ifac,iv(iphas),1) = -3.83d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.4A',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 33 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.287d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 3.214d0 rcodcl(ifac,iv(iphas),1) = -3.83d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.4B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 34 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 25.712d0 rcodcl(ifac,iv(iphas),1) = -30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.4B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 35 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 25.712d0 rcodcl(ifac,iv(iphas),1) = -30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.4C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 36 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 25.712d0 rcodcl(ifac,iv(iphas),1) = -30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.4C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 37 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 25.712d0 rcodcl(ifac,iv(iphas),1) = -30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2L.4D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 38 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 25.712d0 rcodcl(ifac,iv(iphas),1) = -30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('AIR2H.4D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 39 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 2.3d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 25.712d0 rcodcl(ifac,iv(iphas),1) = -30.642d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.3d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('BNOFA.4',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 40 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 0.863d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 16.07d0 rcodcl(ifac,iv(iphas),1) = -19.151d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.2d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) ! (ITYTUR est un indicateur qui vaut ITURB/10) if (itytur(iphas).eq.2) then rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent elseif(itytur(iphas).eq.3) then rcodcl(ifac,ir11(iphas),1) = d2s3*xkent rcodcl(ifac,ir22(iphas),1) = d2s3*xkent rcodcl(ifac,ir33(iphas),1) = d2s3*xkent rcodcl(ifac,ir12(iphas),1) = 0.d0 rcodcl(ifac,ir13(iphas),1) = 0.d0 rcodcl(ifac,ir23(iphas),1) = 0.d0 rcodcl(ifac,iep(iphas),1) = xeent elseif (iturb(iphas).eq.50) then rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent rcodcl(ifac,iphi(iphas),1) = d2s3 rcodcl(ifac,ifb(iphas),1) = 0.d0 elseif (iturb(iphas).eq.60) then rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iomg(iphas),1) = xeent/cmu/xkent endif endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAL.1',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 41 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 28.53d0 rcodcl(ifac,iv(iphas),1) = 28.036d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAH.1',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 42 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 28.53d0 rcodcl(ifac,iv(iphas),1) = 28.036d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAL.2',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 43 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -28.036d0 rcodcl(ifac,iv(iphas),1) = 28.53d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAH.2',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 44 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -28.036d0 rcodcl(ifac,iv(iphas),1) = 28.53d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAL.3',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 45 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -28.53d0 rcodcl(ifac,iv(iphas),1) = -28.036d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAH.3',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 46 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -28.53d0 rcodcl(ifac,iv(iphas),1) = -28.036d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAL.4',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 47 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 28.036d0 rcodcl(ifac,iv(iphas),1) = -28.53d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo CALL GETFBR('OFAH.4',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 48 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientat(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer mass flow rate in kg/s qimpat(izone) = 1.687d0 ! - Oxidizer's Temperature in K timpat(izone) = 550.d0 ! - The color 12 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 28.036d0 rcodcl(ifac,iv(iphas),1) = -28.53d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo ! ------ Turbulence treatment ! Boundary conditions of turbulence icalke(izone) = 1 dh(izone) = 0.25d0 xintur(izone) = 0.05d0 if(icalke(izone).eq.0) then ! Calcul de k et epsilon en entree (XKENT et XEENT) a partir ! l'intensite turbulente et de lois standards en conduite ! circulaire (leur initialisation est inutile mais plus ! propre) uref2 = rcodcl(ifac,iu(iphas),1)**2 & +rcodcl(ifac,iv(iphas),1)**2 & +rcodcl(ifac,iw(iphas),1)**2 uref2 = max(uref2,1.d-12) xkent = epzero xeent = epzero call keenin & !========== ( uref2, xintur(izone), dh(izone), cmu, xkappa, & xkent, xeent ) rcodcl(ifac,ik(iphas),1) = xkent rcodcl(ifac,iep(iphas),1) = xeent endif ! Exemple de cas ou ICALKE(IZONE) = 0 : FIN ! ------ Traitement des scalaires physiques particulieres ! Ils sont traites automatiquement ! ------ Traitement des scalaires utilisateurs ! Exemple : On traite les scalaires rattaches a la phase courante : DEBUT ! Eliminer ces lignes pour la clarte s'il n'y en a pas if ( (nscal-nscapp).gt.0 ) then do ii = 1, (nscal-nscapp) if(iphsca(ii).eq.iphas) then rcodcl(ifac,isca(ii),1) = 1.d0 endif enddo endif enddo ! ---- BOUNDARY FACE for pulverised COAL & primary air INLET CALL GETFBR('AIR1.1B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 49 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 15.321d0 rcodcl(ifac,iv(iphas),1) = 12.856d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.1C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 50 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 15.321d0 rcodcl(ifac,iv(iphas),1) = 12.856d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.1D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 51 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 15.321d0 rcodcl(ifac,iv(iphas),1) = 12.856d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.2B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 52 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -12.856d0 rcodcl(ifac,iv(iphas),1) = 15.321d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.2C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 53 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -12.856d0 rcodcl(ifac,iv(iphas),1) = 15.321d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.2D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 54 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -12.856d0 rcodcl(ifac,iv(iphas),1) = 15.321d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = -surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.3B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 55 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -15.321d0 rcodcl(ifac,iv(iphas),1) = -12.856d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.3C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 56 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -15.321d0 rcodcl(ifac,iv(iphas),1) = -12.856d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.3D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 57 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = -15.321d0 rcodcl(ifac,iv(iphas),1) = -12.856d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.4B',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 58 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 12.856d0 rcodcl(ifac,iv(iphas),1) = -15.321d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.4C',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 59 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 12.856d0 rcodcl(ifac,iv(iphas),1) = -15.321d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo CALL GETFBR('AIR1.4D',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = ientre ! zone's number (from 1 to n) izone = 60 ! - Allocation of the zone's number to the face izfppp(ifac) = izone ! - For theses inlet faces, mass flux is fixed ientcp(izone) = 1 iqimp(izone) = 1 ! - Oxidizer's number (1 to 3) inmoxy(izone) = 1 ! - Oxidizer's mass flow rate in kg/s qimpat(izone) = 3.349d0 ! - Oxidizer's Temperature in K timpat(izone) = 379.d0 ! Coal inlet, initialisation do icha = 1, ncharm qimpcp(izone,icha) = zero timpcp(izone,icha) = zero do iclapc = 1, ncpcmx distch(izone,icha,iclapc) = zero enddo enddo ! Code_Saturne deals with NCHA different coals (component of blend) ! every coal is described by NCLPCH(icha) class of particles ! (each of them described by an inlet diameter) ! ! - Treatment for the first coal icha = 1 ! - Coal mass flow rate in kg/s qimpcp(izone,icha) = 0.93d0 ! - PERCENTAGE mass fraction of each granulometric class do iclapc = 1, nclpch(icha) distch(izone,icha,iclapc) = 100.d0/dble(nclpch(icha)) enddo ! - Inlet temperature for coal & primary air timpcp(izone,icha) = 379.d0 ! - The color 11 becomes an fixed flow rate inlet ! The user gives speed vector direction ! (speed vector norm is irrelevent) rcodcl(ifac,iu(iphas),1) = 12.856d0 rcodcl(ifac,iv(iphas),1) = -15.321d0 rcodcl(ifac,iw(iphas),1) = 0.d0 !do ifac = 1,nfac !rcodcl(ifac, iu(iphas), 1) = -surfbo(1, ifac) !rcodcl(ifac, iv(iphas), 1) = surfbo(2, ifac) !rcodcl(ifac, iw(iphas), 1) = 0.d0 !enddo icalke(izone) = 1 ! ! - If ICALKE = 0 the boundary conditions of turbulence at ! the inlet are calculated as follows: dh(izone) = 0.4d0 xintur(izone) = 0.05d0 enddo ! The color 15 become a WALL CALL GETFBR('WALL',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! WALL : NUL MASS FLUX (PRESSURE FLUX is zero valued) ! FRICTION FOR SPEED (& TURBULENCE) ! NUL SCALAR FLUX ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = iparoi ! zone's number (from 1 to n) izone = 61 ! - Allocation of the zone's number to the face izfppp(ifac) = izone enddo CALL GETFBR('PG1',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! WALL : NUL MASS FLUX (PRESSURE FLUX is zero valued) ! FRICTION FOR SPEED (& TURBULENCE) ! NUL SCALAR FLUX ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = iparoi ! zone's number (from 1 to n) izone = 62 ! - Allocation of the zone's number to the face izfppp(ifac) = izone enddo CALL GETFBR('PG2',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! WALL : NUL MASS FLUX (PRESSURE FLUX is zero valued) ! FRICTION FOR SPEED (& TURBULENCE) ! NUL SCALAR FLUX ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = iparoi ! zone's number (from 1 to n) izone = 63 ! - Allocation of the zone's number to the face izfppp(ifac) = izone enddo ! The color 19 becomes an OUTLET CALL GETFBR('EXIT',NLELT,LSTELT) !========== do ilelt = 1, nlelt ifac = lstelt(ilelt) ! OUTLET : NUL FLUX for SPEED & SCALARS, FIXED PRESSURE ! Kind of boundary conditions for standard variables itypfb(ifac,iphas) = isolib ! zone's number (from 1 to n) izone = 64 ! - Allocation of the zone's number to the face izfppp(ifac) = izone enddo ! The color 14 becomes a symetry plane !---- ! END !---- return end subroutine