From c4f89ef1a4546ceb6dc4ab28522ec17e8add67c1 Mon Sep 17 00:00:00 2001 From: pkuLmq <1900011062@pku.edu.cn> Date: Thu, 19 Sep 2024 17:19:09 +0800 Subject: [PATCH] do a little adjustment --- .../solvers/dfHighSpeedFoam/createFields.H | 29 ------------------- .../solvers/dfHighSpeedFoam/rhoEEqn.H | 4 --- applications/solvers/dfHighSpeedFoam/rhoEqn.H | 8 ----- 3 files changed, 41 deletions(-) diff --git a/applications/solvers/dfHighSpeedFoam/createFields.H b/applications/solvers/dfHighSpeedFoam/createFields.H index f5c7375a..08f1dfcd 100644 --- a/applications/solvers/dfHighSpeedFoam/createFields.H +++ b/applications/solvers/dfHighSpeedFoam/createFields.H @@ -375,32 +375,3 @@ volScalarField pMax ), thermo.p() ); - -// use to analysis -volScalarField::Internal Sm -( - IOobject - ( - "Sm", - runTime.timeName(), - mesh, - IOobject::NO_READ, - IOobject::AUTO_WRITE - ), - mesh, - dimensionedScalar(dimMass/dimTime/dimVolume, Zero) -); - -volScalarField::Internal Se -( - IOobject - ( - "Se", - runTime.timeName(), - mesh, - IOobject::NO_READ, - IOobject::AUTO_WRITE - ), - mesh, - dimensionedScalar(dimEnergy/dimTime/dimVolume, Zero) -); \ No newline at end of file diff --git a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H index c2bc5368..f9d29cf5 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H @@ -18,15 +18,11 @@ if((ddtSchemes == "RK2SSP") || (ddtSchemes == "RK3SSP")) } else { - // use to analysis - Se = parcels.hsTrans()/runTime.deltaT()/mesh.V(); - tmp thcSource(new fvScalarMatrix(rhoE, dimEnergy/dimTime)); fvScalarMatrix& hcSource = thcSource.ref(); forAll(Y, i) { hcSource.source() -= parcels.rhoTrans(i)*chemistry->mixture().Hc(i)/runTime.deltaT(); - Se.field() += parcels.rhoTrans(i)*chemistry->mixture().Hc(i)/runTime.deltaT()/mesh.V(); } solve diff --git a/applications/solvers/dfHighSpeedFoam/rhoEqn.H b/applications/solvers/dfHighSpeedFoam/rhoEqn.H index 3b2105c4..c484cb2c 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoEqn.H @@ -10,14 +10,6 @@ if ((ddtSchemes == "RK2SSP") || (ddtSchemes == "RK3SSP")) } else { - // use to analysis - volScalarField::Internal massSource(0.0*Sm); - forAll(parcels.rhoTrans(),i) - { - massSource += parcels.rhoTrans(i)/runTime.deltaT()/mesh.V(); - } - Sm = massSource; - solve ( fvm::ddt(rho)