$$$$ FIMP NOTICE CHAT 11/09/12 21:16:12 7124 DATE 11/09/12 Operateur FIMP Voir aussi : NAVI -------------- I) SYNTAXE (EQEX) : Cf operateur EQEX _______________ 'ZONE' $paroi 'OPER' FIMP COEF 'INCO' 'TN' Description : _____________ FIMP operator discretizes a flux density or a specific source. IN 2D 2D shell elements (SEG2 ou SEG3) -> Flux (in K/ms) volume elements (TRI3 TRI7 etc) -> specific source (in K/m2s) IN 3D 3D line elements (SEG2 ou SEG3) -> non sense !! 3D shell elements (TRI3 TRI7 etc) -> Flux (in K/m2s) volumic elements (CUB8 CU27 etc) -> specific source (in K/m3s) Commentaires : ______________ $paroi MMODEL object 'NAVIER_STOKES' type associated to the surface (or the volume) where the flux (or source) is prescribed. COEF flux density (or specific source) CHPOINT SCAL CENTRE or CHPOINT SCAL SOMMET or FLOTTANT or MOT (conventionally the inlet flux is positive) TN Temperature field CHPOINT SCAL SOMMET Remarks : ___________ This operator computes a source term on a surface or inside a volume regarding the nature of the geometrical support, for the unknown TN. This unknown must be a CHPOINT SCAL SOMMET. However an extension is possible if we want to add a volumetric source to the equation: Div U = 0 In this case, the operator is applied to the pressure (INCO PRES) and the support of the pressure is specified in the option 'INCOD'. Only volumetric geometrical supports are allowed. ex : 'OPTI' 'INCOD' KPRES 'ZONE' $MT 'OPER' 'FIMP' COEF 'INCO' 'PRES' where KPRES = 'CENTRE' or 'CENTREP1' or 'MSOMMET' When COEF is of type MOT, we look for data in INCO table at the index COEF. OPTION : (EQEX) ________ Finite element formulation EF or EFM1 (available) Dual unknown OPTI INCOD SOMMET (default option) II Numerical discretisation of the Euler Equations ___________________________________________________ IIa : perfect mono-component polytropic gas ________________________________________________ Finite-Volume "cell-centered" discretisation of the Euler equations of gas dynamics for perfect mono-component polytropic gas Unknowns: -------- density, momentum, total energy per unit volume (conserved variables) One can compute : IIa.1 The contribution of the gravity force to the residuum IIa.2 The contribution of the gravity force to the Jacobian IIa.1 Residuum ________________ RCHRES = 'FIMP' 'VF' 'GRAVMONO' 'RESI' LISTINCO CHPRN CHPGN CHPGRAV ; LISTINCO : LISTMOTS object Names of the components of the resultant vector (RCHRES) They are named in the following order: name of the density, name of the momentum, name of the total energy per unit volume CHPRN : CHPOINT object that contains the total mass density (the name of its component is 'SCAL'). CHPGN : CHPOINT object that contains the momentum (two components in 2D, 'UX ','UY ', same geometrical support (SGP) as CHPRN). CHPGRAV : CHPOINT object that contains the gravity (two components in 2D, 'UX ','UY ', same (SGP) as CHPRN). RCHRES : CHPOINT object (components = LISTINCO, same SPG as CHPRN) IIa.2 Jacobian ______________ RJAC = 'FIMP' 'VF' 'GRAVMONO' 'JACOCONS' LISTINCO CHPRN CHPGN CHPGRAV ; LISTINCO : LISTMOTS object They are named in the following order: name of the density, name of the momentum, name of the total energy per unit volume CHPRN : CHPOINT object that contains the total mass density (the name of its component is 'SCAL'). CHPGN : CHPOINT object that contains the momentum (two components in 2D, 'UX ','UY ', same geometrical support (SGP) as CHPRN). CHPGRAV : CHPOINT object that contains the gravity (two components in 2D, 'UX ','UY ', same (SGP) as CHPRN). RJAC : MATRIK object (same SPG as CHPRN) (primal variables = dual variables = LISTINCO) It contains the Jacobian of the residuum with respect to the conservative varriables. III Discretisation of the Navier-Stokes equations with k-epsilon turbulent closure model __________________________________________________________ IIIa : gaz multi-component with specific heats depending on temperature ___________________________________________________ Discretisation of the Navier-Stokes multi-component equations with k-epsilon turbulent closure model Unknowns: density, momentum, total energy per unit volume, densities of the species which are in (TABGAS.'ESPEULE'), turbulent kinetic energy per unit volume, rate of dissipation of the turbulent energy per unit volume (conservative variables) One can compute: IIIa.1 the contribution to the residual due to the: gravitational force, the source terms of the species conservation eqs., the source terms of the equations for the turbulent kinetic energy and for the rate of dissipation of the turbulent energy The Residual _____________ RCHRES = 'FIMP' 'VF' 'GRAVKEPS' 'RESI' TABGAS LISTINCO CHPRN CHPGN CHPGRAV CHPRKA CHPREPS GRADR GRADV ; TABGAS : table with the gas properties; here we need: a) turbulent constant \sigma_t (see gravity term in k-\eps equations) name: TABGAS . 'SIGT' b) turbulent constant C_{\eps 1} (see \eps-equation) name: TABGAS . 'CEPS1' c) turbulent constant C_{\eps 2} (see \eps-equation) name: TABGAS . 'CEPS2' d) rate of productions for the species (source terms); type 'CHPOINT' (SPG = 'CENTRE', the number of components is equal to the number of species in TABGAS . 'ESPEULE') e) TABGAS . 'ESPEULE' - the names of the species which are explicitely in the equations f) TABGAS . 'MUTURB' - turbulent viscosity (SPG = 'CENTRE'; number of compon. = 1) g) the formatio enthalpies of the species at the reference temperature T=0K LISTINCO : object of the type LISTMOTS. Names of the components of the resultant residual vector (RCHRES). They are given in the following order: the name of the density, of the momentum, of the total energy per unit volume, the species which are in (TABGAS.'ESPEULE'), the kinetic turbulent energy per unit volume, the rate of dissipation of the turbulent energy per unit volume CHPRN : CHPOINT which contains the density values (one component, 'SCAL'). CHPGN : CHPOINT which contains the values of the momentum (two components in 2D, 'UX ', 'UY ', the same Geometrical Support (GS) as for CHPRN). CHPGRAV : CHPOINT which contains the gravity vectors (two components in 2D, 'UX ', 'UY ', the same Geometrical Support as for CHPRN). CHPRKA : CHPOINT which contains kinetic turbulent energy per unit volume (one component, the same Geometrical Support as for CHPRN). CHPREPS : CHPOINT which contains the rate of dissipation of the turbulent energy per unit volume (one component, the same Geometrical Support as for CHPRN). GRADR : gradient of the density (GS = 'CENTRE', components: 'P1DX','P1DY',('P1DZ')) GRADV : gradient of the velocity vector (GS = 'CENTRE', components: 'P1DX','P1DY', ('P1DZ') 'P2DX','P2DY', ('P2DZ') ('P3DX','P3DY', ('P3DZ'))) RCHRES : object of the type CHPOINT (components = LISTINCO, the same Geometrical Support as for CHPRN).
© Cast3M 2003 - All rights reserved.
Disclaimer