00001 const fvMesh& mesh = fluidRegions[i]; 00002 00003 basicRhoThermo& thermo = thermoFluid[i]; 00004 volScalarField& rho = rhoFluid[i]; 00005 volScalarField& K = KFluid[i]; 00006 volVectorField& U = UFluid[i]; 00007 surfaceScalarField& phi = phiFluid[i]; 00008 00009 compressible::turbulenceModel& turb = turbulence[i]; 00010 volScalarField& DpDt = DpDtFluid[i]; 00011 00012 volScalarField& p = thermo.p(); 00013 const volScalarField& psi = thermo.psi(); 00014 volScalarField& h = thermo.h(); 00015 00016 volScalarField& p_rgh = p_rghFluid[i]; 00017 const volScalarField& gh = ghFluid[i]; 00018 const surfaceScalarField& ghf = ghfFluid[i]; 00019 00020 const dimensionedScalar initialMass 00021 ( 00022 "initialMass", 00023 dimMass, 00024 initialMassFluid[i] 00025 );