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 00011 volScalarField& p = thermo.p(); 00012 const volScalarField& psi = thermo.psi(); 00013 volScalarField& h = thermo.h(); 00014 00015 const dimensionedScalar initialMass 00016 ( 00017 "initialMass", 00018 dimMass, 00019 initialMassFluid[i] 00020 ); 00021 00022 const label pRefCell = pRefCellFluid[i]; 00023 const scalar pRefValue = pRefValueFluid[i]; 00024 00025 volScalarField& p_rgh = p_rghFluid[i]; 00026 const volScalarField& gh = ghFluid[i]; 00027 const surfaceScalarField& ghf = ghfFluid[i];