From 7b9f550505b9739cd145b213eea9f7070b48e472 Mon Sep 17 00:00:00 2001 From: pkuLmq <1900011062@pku.edu.cn> Date: Mon, 11 Jul 2022 14:21:48 +0800 Subject: [PATCH 1/6] add diffusion velocity correction for dfHighSpeedFoam, correct alpha in rhoEEqn.H and clean the code --- .../solvers/dfHighSpeedFoam/Make/options | 4 ++ .../solvers/dfHighSpeedFoam/createFields.H | 43 +++++++++++++++++- .../solvers/dfHighSpeedFoam/dfHighSpeedFoam.C | 18 ++++++-- .../solvers/dfHighSpeedFoam/rhoEEqn.H | 31 ++++++------- .../solvers/dfHighSpeedFoam/rhoYEqn.H | 45 ++++++++++++------- .../det_1D_H2_01mm/system/blockMeshDict | 2 - .../det_1D_H2_01mm/system/fvSchemes | 26 +++++------ 7 files changed, 114 insertions(+), 55 deletions(-) diff --git a/applications/solvers/dfHighSpeedFoam/Make/options b/applications/solvers/dfHighSpeedFoam/Make/options index f5f42f89..e8d6d66c 100644 --- a/applications/solvers/dfHighSpeedFoam/Make/options +++ b/applications/solvers/dfHighSpeedFoam/Make/options @@ -1,4 +1,7 @@ +-include $(GENERAL_RULES)/mplibType + EXE_INC = -std=c++14 \ + $(PFLAGS) $(PINC) \ -I$(FOAM_APP)/solvers/compressible/rhoCentralFoam/BCs/lnInclude \ -I$(LIB_SRC)/finiteVolume/cfdTools \ -I$(LIB_SRC)/finiteVolume/lnInclude \ @@ -8,6 +11,7 @@ EXE_INC = -std=c++14 \ -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \ -I$(LIB_SRC)/sampling/lnInclude \ -I$(LIB_SRC)/dynamicFvMesh/lnInclude \ + -I$(LIB_SRC)/Pstream/mpi \ -I$(LIB_SRC)/meshTools/lnInclude \ -I$(DF_SRC)/CanteraMixture/lnInclude \ -I$(DF_SRC)/dfChemistryModel \ diff --git a/applications/solvers/dfHighSpeedFoam/createFields.H b/applications/solvers/dfHighSpeedFoam/createFields.H index 563eef8f..5dded90b 100644 --- a/applications/solvers/dfHighSpeedFoam/createFields.H +++ b/applications/solvers/dfHighSpeedFoam/createFields.H @@ -126,4 +126,45 @@ forAll(Y, i) { fields.add(Y[i]); } -fields.add(thermo.he()); \ No newline at end of file +fields.add(thermo.he()); + +const scalar Sct = chemistry.lookupOrDefault("Sct", 1.); +volScalarField diffAlphaD +( + IOobject + ( + "diffAlphaD", + runTime.timeName(), + mesh, + IOobject::NO_READ, + IOobject::NO_WRITE + ), + mesh, + dimensionedScalar(dimEnergy/dimTime/dimVolume, 0) +); +volVectorField hDiffCorrFlux +( + IOobject + ( + "hDiffCorrFlux", + runTime.timeName(), + mesh, + IOobject::NO_READ, + IOobject::NO_WRITE + ), + mesh, + dimensionedVector(dimensionSet(1,0,-3,0,0,0,0), Zero) +); +volVectorField sumYDiffError +( + IOobject + ( + "sumYDiffError", + runTime.timeName(), + mesh, + IOobject::NO_READ, + IOobject::NO_WRITE + ), + mesh, + dimensionedVector("sumYDiffError", dimDynamicViscosity/dimLength, Zero) +); diff --git a/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C b/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C index 8654cb95..042d1a6a 100644 --- a/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C +++ b/applications/solvers/dfHighSpeedFoam/dfHighSpeedFoam.C @@ -42,8 +42,7 @@ Description #include "directionInterpolate.H" #include "localEulerDdtScheme.H" #include "fvcSmooth.H" -//add -//#include "fvOptions.H" +#include "PstreamGlobals.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -57,7 +56,11 @@ int main(int argc, char *argv[]) #include "createDynamicFvMesh.H" #include "createFields.H" #include "createTimeControls.H" - //#include "createFvOptions.H" + + double time_monitor_flow; + double time_monitor_chem; + double time_monitor_Y; + clock_t start, end; turbulence->validate(); @@ -200,11 +203,14 @@ int main(int argc, char *argv[]) volTensorField tauMC("tauMC", muEff*dev2(Foam::T(fvc::grad(U)))); // --- Solve density - //solve(fvm::ddt(rho) + fvc::div(phi)); #include "rhoEqn.H" + start = std::clock(); // --- Solve momentum #include "rhoUEqn.H" + end = std::clock(); + time_monitor_flow += double(end - start) / double(CLOCKS_PER_SEC); + // --- Solve species #include "rhoYEqn.H" @@ -216,6 +222,10 @@ int main(int argc, char *argv[]) runTime.write(); + Info<< "MonitorTime_chem = " << time_monitor_chem << " s" << nl << endl; + Info<< "MonitorTime_Y = " << time_monitor_Y << " s" << nl << endl; + Info<< "MonitorTime_flow = " << time_monitor_flow << " s" << nl << endl; + Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << nl << endl; diff --git a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H index 23eaf01b..c6f8c3e8 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H @@ -8,13 +8,22 @@ surfaceScalarField sigmaDotU & (a_pos*U_pos + a_neg*U_neg) ); -solve +fvScalarMatrix EEqn ( - fvm::ddt(rhoE) - + fvc::div(phiEp) - - fvc::div(sigmaDotU) + fvm::ddt(rhoE) + + fvc::div(phiEp) + - fvc::div(sigmaDotU) + // alpha in deepflame is considered to calculate h by default (kappa/Cp), so multiply gamma to correct alpha + - fvc::laplacian(turbulence->alphaEff()*thermo.gamma(), e) + + diffAlphaD + == + fvc::div(hDiffCorrFlux) ); +EEqn.relax(); + +EEqn.solve(); + e = rhoE/rho - 0.5*magSqr(U); e.correctBoundaryConditions(); @@ -23,20 +32,6 @@ chemistry.correctThermo(); // before this, we must update ha = e + p/rho rhoE.boundaryFieldRef() == rho.boundaryField()*(e.boundaryField() + 0.5*magSqr(U.boundaryField())); -if (!inviscid) -{ - fvScalarMatrix eEqn - ( - fvm::ddt(rho, e) - fvc::ddt(rho, e) - - fvc::laplacian(turbulence->alphaEff(), e) - ); - eEqn.solve("e"); - - ha = e + p/rho; - chemistry.correctThermo(); - rhoE = rho*(e + 0.5*magSqr(U)); -} - Info<< "min/max(T) = " << min(T).value() << ", " << max(T).value() << endl; diff --git a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H index 84221586..8360e195 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H @@ -1,3 +1,7 @@ +hDiffCorrFlux = Zero; +diffAlphaD = Zero; +sumYDiffError = Zero; + tmp> mvConvection ( fv::convectionScheme::New @@ -9,39 +13,47 @@ tmp> mvConvection ) ); +forAll(Y, i) { + sumYDiffError += chemistry.rhoD(i)*fvc::grad(Y[i]); +} +const surfaceScalarField phiUc = linearInterpolate(sumYDiffError) & mesh.Sf(); + +{ + start = std::clock(); chemistry.solve(mesh.time().deltaTValue()); + label flag_mpi_init; + MPI_Initialized(&flag_mpi_init); + if(flag_mpi_init) MPI_Barrier(PstreamGlobals::MPI_COMM_FOAM); + end = std::clock(); + time_monitor_chem += double(end - start) / double(CLOCKS_PER_SEC); volScalarField Yt(0*Y[0]); + start = std::clock(); forAll(Y, i) { + volScalarField& Yi = Y[i]; + hDiffCorrFlux += chemistry.hai(i)*(chemistry.rhoD(i)*fvc::grad(Yi) - Yi*sumYDiffError); + diffAlphaD += fvc::laplacian(thermo.alpha()*chemistry.hai(i), Yi); + if (i != inertIndex) { - volScalarField& Yi = Y[i]; - + tmp DEff = chemistry.rhoD(i) + turbulence->mut()/Sct; fvScalarMatrix YiEqn ( - fvm::ddt(rho, Yi) - + mvConvection->fvmDiv(phi, Yi) + fvm::ddt(rho, Yi) + + mvConvection->fvmDiv(phi, Yi) + - fvm::laplacian(DEff(), Yi) + + mvConvection->fvmDiv(phiUc, Yi) == - chemistry.RR(i) - //+ fvOptions(rho, Yi) + chemistry.RR(i) ); - if (!inviscid) - { - YiEqn -= fvm::laplacian(turbulence->muEff(), Yi); - } - YiEqn.relax(); - //fvOptions.constrain(YiEqn); - YiEqn.solve("Yi"); - //fvOptions.correct(Yi); - Yi.max(0.0); Yt += Yi; } @@ -49,4 +61,7 @@ tmp> mvConvection Y[inertIndex] = scalar(1) - Yt; Y[inertIndex].max(0.0); + + end = std::clock(); + time_monitor_Y += double(end - start) / double(CLOCKS_PER_SEC); } diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict b/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict index 8199aa9c..1a73e764 100755 --- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict +++ b/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/blockMeshDict @@ -41,7 +41,6 @@ boundary ( Left { - // type patch; type symmetryPlane; faces ( @@ -50,7 +49,6 @@ boundary } Right { - // type patch; type symmetryPlane; faces ( diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSchemes b/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSchemes index f4f06e9a..3f654d84 100755 --- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSchemes +++ b/examples/dfHighSpeedFoam/det_1D_H2_01mm/system/fvSchemes @@ -29,33 +29,29 @@ gradSchemes divSchemes { - default none; - div(tauMC) Gauss linear; - div(phi,O2) Gauss vanLeer01; - div(phi,H2) Gauss vanLeer01; - div(phi,H2O) Gauss vanLeer01; - div(phi,AR) Gauss vanLeer01; - div(phi,Yi) Gauss vanLeer01; - div(phi,Yi_h) Gauss vanLeer; + default none; + div(tauMC) Gauss linear; + div(hDiffCorrFlux) Gauss cubic; + div(phi,Yi_h) Gauss vanLeer; } laplacianSchemes { - default Gauss linear uncorrected; + default Gauss linear uncorrected; } interpolationSchemes { - default linear; - reconstruct(rho) Minmod; - reconstruct(U) MinmodV; - reconstruct(T) Minmod; - reconstruct(Yi) Minmod; + default linear; + reconstruct(rho) Minmod; + reconstruct(U) MinmodV; + reconstruct(T) Minmod; + reconstruct(Yi) Minmod; } snGradSchemes { - default uncorrected; + default uncorrected; } From 313182773dc44b51ea0b40066e5f1a8935bc183d Mon Sep 17 00:00:00 2001 From: pkuLmq <1900011062@pku.edu.cn> Date: Tue, 19 Jul 2022 11:02:41 +0800 Subject: [PATCH 2/6] add diffusion velocity correction and timing function for dfHighSpeedFoam, correct alpha in rhoEEqn.H and clean the code --- .../solvers/dfHighSpeedFoam/rhoEEqn.H | 24 +++++++++++++------ .../solvers/dfHighSpeedFoam/rhoYEqn.H | 7 +++++- .../constant/CanteraTorchProperties | 5 ++++ 3 files changed, 28 insertions(+), 8 deletions(-) diff --git a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H index c6f8c3e8..b80e1d17 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H @@ -8,22 +8,16 @@ surfaceScalarField sigmaDotU & (a_pos*U_pos + a_neg*U_neg) ); -fvScalarMatrix EEqn +solve ( fvm::ddt(rhoE) + fvc::div(phiEp) - fvc::div(sigmaDotU) - // alpha in deepflame is considered to calculate h by default (kappa/Cp), so multiply gamma to correct alpha - - fvc::laplacian(turbulence->alphaEff()*thermo.gamma(), e) + diffAlphaD == fvc::div(hDiffCorrFlux) ); -EEqn.relax(); - -EEqn.solve(); - e = rhoE/rho - 0.5*magSqr(U); e.correctBoundaryConditions(); @@ -32,6 +26,22 @@ chemistry.correctThermo(); // before this, we must update ha = e + p/rho rhoE.boundaryFieldRef() == rho.boundaryField()*(e.boundaryField() + 0.5*magSqr(U.boundaryField())); +if (!inviscid) +{ + fvScalarMatrix eEqn + ( + fvm::ddt(rho, e) - fvc::ddt(rho, e) + // alpha in deepflame is considered to calculate h by default (kappa/Cp), so multiply gamma to correct alpha + - fvc::laplacian(turbulence->alphaEff()*thermo.gamma(), e) + ); + + eEqn.solve("e"); + + ha = e + p/rho; + chemistry.correctThermo(); + rhoE = rho*(e + 0.5*magSqr(U)); +} + Info<< "min/max(T) = " << min(T).value() << ", " << max(T).value() << endl; diff --git a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H index 8360e195..5493bd39 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H @@ -39,7 +39,7 @@ const surfaceScalarField phiUc = linearInterpolate(sumYDiffError) & mesh.Sf(); if (i != inertIndex) { - tmp DEff = chemistry.rhoD(i) + turbulence->mut()/Sct; + tmp DEff = chemistry.rhoD(i); fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) @@ -50,6 +50,11 @@ const surfaceScalarField phiUc = linearInterpolate(sumYDiffError) & mesh.Sf(); chemistry.RR(i) ); + if (!inviscid) + { + YiEqn -= fvm::laplacian(turbulence->muEff(), Yi); + } + YiEqn.relax(); YiEqn.solve("Yi"); diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties b/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties index d4b25919..78517df2 100755 --- a/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties +++ b/examples/dfHighSpeedFoam/det_1D_H2_01mm/constant/CanteraTorchProperties @@ -19,6 +19,11 @@ FoamFile chemistry on; CanteraMechanismFile "H2_Li.xml"; transportModel "UnityLewis"; +loadbalancing +{ + active false; + //log true; +} odeCoeffs {} inertSpecie "N2"; From e0bb5e6cf2c3beac5210d98b48f5bd7e52341713 Mon Sep 17 00:00:00 2001 From: pkuLmq <1900011062@pku.edu.cn> Date: Tue, 19 Jul 2022 11:55:12 +0800 Subject: [PATCH 3/6] add diffusion velocity correction and timing function for dfHighSpeedFoam, correct alpha in rhoEEqn.H and clean the code --- applications/solvers/dfHighSpeedFoam/rhoEEqn.H | 12 ++++++------ applications/solvers/dfHighSpeedFoam/rhoYEqn.H | 6 ++---- examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz | Bin 0 -> 374 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz | Bin 0 -> 371 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz | Bin 0 -> 373 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/T.orig.gz | Bin 0 -> 372 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/U.orig.gz | Bin 0 -> 377 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz | Bin 0 -> 377 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz | Bin 0 -> 391 bytes 9 files changed, 8 insertions(+), 10 deletions(-) create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/T.orig.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.orig.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz diff --git a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H index b80e1d17..c1ba18d6 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoEEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoEEqn.H @@ -13,9 +13,6 @@ solve fvm::ddt(rhoE) + fvc::div(phiEp) - fvc::div(sigmaDotU) - + diffAlphaD - == - fvc::div(hDiffCorrFlux) ); e = rhoE/rho - 0.5*magSqr(U); @@ -30,11 +27,14 @@ if (!inviscid) { fvScalarMatrix eEqn ( - fvm::ddt(rho, e) - fvc::ddt(rho, e) + fvm::ddt(rho, e) - fvc::ddt(rho, e) // alpha in deepflame is considered to calculate h by default (kappa/Cp), so multiply gamma to correct alpha - - fvc::laplacian(turbulence->alphaEff()*thermo.gamma(), e) + - fvc::laplacian(turbulence->alphaEff()*thermo.gamma(), e) + + diffAlphaD + == + fvc::div(hDiffCorrFlux) ); - + eEqn.solve("e"); ha = e + p/rho; diff --git a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H index 5493bd39..f0a50aec 100644 --- a/applications/solvers/dfHighSpeedFoam/rhoYEqn.H +++ b/applications/solvers/dfHighSpeedFoam/rhoYEqn.H @@ -39,20 +39,18 @@ const surfaceScalarField phiUc = linearInterpolate(sumYDiffError) & mesh.Sf(); if (i != inertIndex) { - tmp DEff = chemistry.rhoD(i); + tmp DEff = chemistry.rhoD(i) + turbulence->mut()/Sct; fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi) - - fvm::laplacian(DEff(), Yi) - + mvConvection->fvmDiv(phiUc, Yi) == chemistry.RR(i) ); if (!inviscid) { - YiEqn -= fvm::laplacian(turbulence->muEff(), Yi); + YiEqn -= fvm::laplacian(DEff(), Yi) - mvConvection->fvmDiv(phiUc, Yi); } YiEqn.relax(); diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz new file mode 100644 index 0000000000000000000000000000000000000000..a09d4b08708063b58a3429719b83315db118ffda GIT binary patch literal 374 zcmV-+0g3(}iwFp6lA&S%14uFetx~aW!Y~k>b6;^{3<9S}RV$(nh$`w(L{(5bG#yMh zz?EZ1jzNtI|6a#QThs*xo)_ot-Mi|Hlan>#~H zP>@{yb#YpdN~WX95a%26UpUiSlM#-Sdz@=67y8p_v?5Ldr&y5TR%q2~ywai+mWGHM zYm0Kon07iV&GWuC>w4Hv*Rwc6wLG?mFJww5<%VxO(@{MUC&-m3wzA!{hxCjz&PK=} zaRw713BLT|U1P-W9`vKw7gyTk+%{f$nF-O9Gs(-ceZo%5S;i$dNgZv|mpc0ms1tRX zqF#JZ)@nTrVmJU}7+&l!2Cx)4srnVIxwi<9&U*|XlqI9MOxpiot3>@qjQY{w_WCMD zfQ8lQq?o2m`B0M@&$&Hw-a literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz new file mode 100644 index 0000000000000000000000000000000000000000..cb797eb47bdc8a180bcf0aead79473163c350a20 GIT binary patch literal 371 zcmV-(0gV11iwFp6lA&S%15Ppktx~aW!Y~k>b6;^{3<4w63W=x#qRLcCRZu%L9ZWdD zm19SagBlh7y^fQ%s0$1{FV5Y&ch65@VebsW9!|Sm?42b#Y66INKfwPuxNed(cZL|F zAi4hQ;yfo6&nL4nE_dX=aG^_+5>DfLTxu;>`qOE&B2EJ5n3LgFXw_=G(y9=a#)vy> z^J2`Hb~+o)v%WUlde~3bvp7MuKDCE0WC|zchDV<1sGf*3;5i*FK z!B|LwFTZ%-81cJ@{Xyi58*MUf8?U@bg=osDb6;^{1c4DlD-m@-RGEsX3TlVa!ATBq z<=Bzqphks%uj8a`R8@k3=f$~u_wM;AtUo-1{t)Mb0S?a+9d!eUci+IjIJ)j;)7%+i zgo0%H+r?#0Dq7AKQ(W)Kf8k1(CLx?h_qf(t#`?o)v?5Ldmza~`R%q2~yii;SOH;(X zwRtgROgo*e=GjP_T|Mll>sg$kO1t*(nM~oN-0-buI;tn)0=W`-Dcem)*j$jt#ppAL zoWV#)f=|Er&=~Q%C!=xbi(73nZX2(>NQ7w0iR4AmKH;F{D&dlwsE)SjOP%}#)EoDj zqK@AwX|1;e5!D55{26NTu9AT2)!$K*MsSoQBOsoaWphC?IS2B3KU{-U3 zBY=%DRXw6qS)Y|+a%pD(t&%-Y!!Qtr=lEBgIzg(`DGO2~hN>tNRS}dj z5C$i?wnmN}IZj1H``>Y#02NC+^!4K0J@?^Pp?+{eegL!K5Q3ASMH2#W_apePz0D>` zcc%a`a1{Ap78WJScrl$%V7W#21<$J18Nw{Show>?RiCy;-^8|n1(c|n;Y#*BUNNn> z!3n_Dn6jEss%)RNV#P@5ONzQQN zZ26SLcEOknBwuduu3N;-z8Qs)3$B$enCW(9RmOQ2&IGHf{s}uJR+t&3$RiL6v@3t}A+p$bnD_{>2WZgyg(H**>D*?Ja~c1Y-!h^Bxhxr9my3i2MJbmb`r< zTsh*}5~-?X&UE7l>y7t-YePJH@k~=9R83K!(aob^5+mXri1)-d+Y<;DrD+_}F6O8s XxvT>Zgg9R=rPh~kD3$)Ba037UX7aSp literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz new file mode 100644 index 0000000000000000000000000000000000000000..e400169bc66ed0f4986bb9f2cc34944b709c3408 GIT binary patch literal 377 zcmV-<0fzn`iwFp6lA&S%16gEcW?^+~bO5bVu};G<5S`<%ICV_bI$}X;#84GwqAH>l z2oVM+xwJ-(9XU=#MEmd9PJoJ~9s0aDzIX55J%xsoGYBWJ7>{9cmT0Ld0l51C{Kvso z)3m=+fCM;-;;#y;8fCJYFK4jXq5p!lYITmVNbX^yl*rVl)fhyq1gxM&%?wv^(0IYD z;Ra^_J7el*MyaworDD}o>8)$_-L4nrAd6Fb@Qk{#QnvZlGcDB>o>|D`tAHtjRg=pL4;QW>{gb#5!k!>BL3b*M-V| z0~$r6{!wKgm=9X71(vK8CMHKhDTG%$loC?#3Z>g+u=p0i$$5_n;nJX%OeDkK&`Rz$ zg~;tqGS$k0>CRVeI>7^$h6MH!n65@>x~f8>yGOw!#w0kB;E8W{XDC?Jrt=i%$+J%2 XtDbly#7A;1xJbSM!FMQwasvPWa2B=R literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz new file mode 100644 index 0000000000000000000000000000000000000000..7576a206b4f6ac218ea4b555e941a3c9b24d070c GIT binary patch literal 391 zcmV;20eJo&iwFp6lA&S%18^>Ha%pD(tx{2I!Y~ki&#yRqR5aBw29`aDqfZmZl&xw+Fy)>})hm z{hk6Oz)>uIqA)E{CezVm2=g_%6U@{`=LqBE4(3XUOnq35AYvt83MFc0xROER8MBHT z90IJ3DXSr+%IYi?D|$+=8nR4l22E7>l<;k2fZkEmrJD!W_+)#$~pJfxnNZlM%XGb%b8$0X{_~eq4IBo zT2aeis_YH(pjF0XPbh`<#TKQ6EO>#k86x!d(uc171`zr68WX~$K`og`!b8|d-V6!- zsDCxMj;#(D1sQ^}hP-QTjC~S>WNIT9OxKPU9J|M*aT|YZVCoW~s*3`Rt{()G=x5)% l|07?m|Hr_hH1%JMBe literal 0 HcmV?d00001 From 624ca1516567406f36d5b1572ab67523c738dd0f Mon Sep 17 00:00:00 2001 From: pkuLmq <1900011062@pku.edu.cn> Date: Tue, 19 Jul 2022 11:56:00 +0800 Subject: [PATCH 4/6] add diffusion velocity correction and timing function for dfHighSpeedFoam, correct alpha in rhoEEqn.H and clean the code --- examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz | Bin 374 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz | Bin 371 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz | Bin 373 -> 0 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/T.orig.gz | Bin 372 -> 0 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/U.orig.gz | Bin 377 -> 0 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz | Bin 377 -> 0 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz | Bin 391 -> 0 bytes 7 files changed, 0 insertions(+), 0 deletions(-) delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/T.orig.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.orig.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz deleted file mode 100644 index a09d4b08708063b58a3429719b83315db118ffda..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 374 zcmV-+0g3(}iwFp6lA&S%14uFetx~aW!Y~k>b6;^{3<9S}RV$(nh$`w(L{(5bG#yMh zz?EZ1jzNtI|6a#QThs*xo)_ot-Mi|Hlan>#~H zP>@{yb#YpdN~WX95a%26UpUiSlM#-Sdz@=67y8p_v?5Ldr&y5TR%q2~ywai+mWGHM zYm0Kon07iV&GWuC>w4Hv*Rwc6wLG?mFJww5<%VxO(@{MUC&-m3wzA!{hxCjz&PK=} zaRw713BLT|U1P-W9`vKw7gyTk+%{f$nF-O9Gs(-ceZo%5S;i$dNgZv|mpc0ms1tRX zqF#JZ)@nTrVmJU}7+&l!2Cx)4srnVIxwi<9&U*|XlqI9MOxpiot3>@qjQY{w_WCMD zfQ8lQq?o2m`B0M@&$&Hw-a diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz deleted file mode 100644 index cb797eb47bdc8a180bcf0aead79473163c350a20..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 371 zcmV-(0gV11iwFp6lA&S%15Ppktx~aW!Y~k>b6;^{3<4w63W=x#qRLcCRZu%L9ZWdD zm19SagBlh7y^fQ%s0$1{FV5Y&ch65@VebsW9!|Sm?42b#Y66INKfwPuxNed(cZL|F zAi4hQ;yfo6&nL4nE_dX=aG^_+5>DfLTxu;>`qOE&B2EJ5n3LgFXw_=G(y9=a#)vy> z^J2`Hb~+o)v%WUlde~3bvp7MuKDCE0WC|zchDV<1sGf*3;5i*FK z!B|LwFTZ%-81cJ@{Xyi58*MUf8?U@bg=osDb6;^{1c4DlD-m@-RGEsX3TlVa!ATBq z<=Bzqphks%uj8a`R8@k3=f$~u_wM;AtUo-1{t)Mb0S?a+9d!eUci+IjIJ)j;)7%+i zgo0%H+r?#0Dq7AKQ(W)Kf8k1(CLx?h_qf(t#`?o)v?5Ldmza~`R%q2~yii;SOH;(X zwRtgROgo*e=GjP_T|Mll>sg$kO1t*(nM~oN-0-buI;tn)0=W`-Dcem)*j$jt#ppAL zoWV#)f=|Er&=~Q%C!=xbi(73nZX2(>NQ7w0iR4AmKH;F{D&dlwsE)SjOP%}#)EoDj zqK@AwX|1;e5!D55{26NTu9AT2)!$K*MsSoQBOsoaWphC?IS2B3KU{-U3 zBY=%DRXw6qS)Y|+a%pD(t&%-Y!!Qtr=lEBgIzg(`DGO2~hN>tNRS}dj z5C$i?wnmN}IZj1H``>Y#02NC+^!4K0J@?^Pp?+{eegL!K5Q3ASMH2#W_apePz0D>` zcc%a`a1{Ap78WJScrl$%V7W#21<$J18Nw{Show>?RiCy;-^8|n1(c|n;Y#*BUNNn> z!3n_Dn6jEss%)RNV#P@5ONzQQN zZ26SLcEOknBwuduu3N;-z8Qs)3$B$enCW(9RmOQ2&IGHf{s}uJR+t&3$RiL6v@3t}A+p$bnD_{>2WZgyg(H**>D*?Ja~c1Y-!h^Bxhxr9my3i2MJbmb`r< zTsh*}5~-?X&UE7l>y7t-YePJH@k~=9R83K!(aob^5+mXri1)-d+Y<;DrD+_}F6O8s XxvT>Zgg9R=rPh~kD3$)Ba037UX7aSp diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz deleted file mode 100644 index e400169bc66ed0f4986bb9f2cc34944b709c3408..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 377 zcmV-<0fzn`iwFp6lA&S%16gEcW?^+~bO5bVu};G<5S`<%ICV_bI$}X;#84GwqAH>l z2oVM+xwJ-(9XU=#MEmd9PJoJ~9s0aDzIX55J%xsoGYBWJ7>{9cmT0Ld0l51C{Kvso z)3m=+fCM;-;;#y;8fCJYFK4jXq5p!lYITmVNbX^yl*rVl)fhyq1gxM&%?wv^(0IYD z;Ra^_J7el*MyaworDD}o>8)$_-L4nrAd6Fb@Qk{#QnvZlGcDB>o>|D`tAHtjRg=pL4;QW>{gb#5!k!>BL3b*M-V| z0~$r6{!wKgm=9X71(vK8CMHKhDTG%$loC?#3Z>g+u=p0i$$5_n;nJX%OeDkK&`Rz$ zg~;tqGS$k0>CRVeI>7^$h6MH!n65@>x~f8>yGOw!#w0kB;E8W{XDC?Jrt=i%$+J%2 XtDbly#7A;1xJbSM!FMQwasvPWa2B=R diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.orig.gz deleted file mode 100644 index 7576a206b4f6ac218ea4b555e941a3c9b24d070c..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 391 zcmV;20eJo&iwFp6lA&S%18^>Ha%pD(tx{2I!Y~ki&#yRqR5aBw29`aDqfZmZl&xw+Fy)>})hm z{hk6Oz)>uIqA)E{CezVm2=g_%6U@{`=LqBE4(3XUOnq35AYvt83MFc0xROER8MBHT z90IJ3DXSr+%IYi?D|$+=8nR4l22E7>l<;k2fZkEmrJD!W_+)#$~pJfxnNZlM%XGb%b8$0X{_~eq4IBo zT2aeis_YH(pjF0XPbh`<#TKQ6EO>#k86x!d(uc171`zr68WX~$K`og`!b8|d-V6!- zsDCxMj;#(D1sQ^}hP-QTjC~S>WNIT9OxKPU9J|M*aT|YZVCoW~s*3`Rt{()G=x5)% l|07?m|Hr_hH1%JMBe From 058166d7243ce7f4f35fe9f71e044d60ecdc07c1 Mon Sep 17 00:00:00 2001 From: pkuLmq <1900011062@pku.edu.cn> Date: Tue, 19 Jul 2022 12:17:32 +0800 Subject: [PATCH 5/6] correct det_1D_H2_01mm --- examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz | Bin 0 -> 374 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz | Bin 0 -> 371 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz | Bin 0 -> 373 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/T.gz | Bin 0 -> 456 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz | Bin 0 -> 368 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz | Bin 0 -> 377 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz | Bin 0 -> 504 bytes .../det_1D_H2_01mm/0_orig/{T.orig.gz => T.gz} | Bin .../det_1D_H2_01mm/0_orig/{U.orig.gz => U.gz} | Bin .../det_1D_H2_01mm/0_orig/{p.orig.gz => p.gz} | Bin .../det_1D_H2_01mm/system/controlDict | 2 +- 11 files changed, 1 insertion(+), 1 deletion(-) create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/T.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz create mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz rename examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/{T.orig.gz => T.gz} (100%) rename examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/{U.orig.gz => U.gz} (100%) rename examples/dfHighSpeedFoam/det_1D_H2_01mm/0_orig/{p.orig.gz => p.gz} (100%) diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz new file mode 100644 index 0000000000000000000000000000000000000000..a09d4b08708063b58a3429719b83315db118ffda GIT binary patch literal 374 zcmV-+0g3(}iwFp6lA&S%14uFetx~aW!Y~k>b6;^{3<9S}RV$(nh$`w(L{(5bG#yMh zz?EZ1jzNtI|6a#QThs*xo)_ot-Mi|Hlan>#~H zP>@{yb#YpdN~WX95a%26UpUiSlM#-Sdz@=67y8p_v?5Ldr&y5TR%q2~ywai+mWGHM zYm0Kon07iV&GWuC>w4Hv*Rwc6wLG?mFJww5<%VxO(@{MUC&-m3wzA!{hxCjz&PK=} zaRw713BLT|U1P-W9`vKw7gyTk+%{f$nF-O9Gs(-ceZo%5S;i$dNgZv|mpc0ms1tRX zqF#JZ)@nTrVmJU}7+&l!2Cx)4srnVIxwi<9&U*|XlqI9MOxpiot3>@qjQY{w_WCMD zfQ8lQq?o2m`B0M@&$&Hw-a literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz new file mode 100644 index 0000000000000000000000000000000000000000..cb797eb47bdc8a180bcf0aead79473163c350a20 GIT binary patch literal 371 zcmV-(0gV11iwFp6lA&S%15Ppktx~aW!Y~k>b6;^{3<4w63W=x#qRLcCRZu%L9ZWdD zm19SagBlh7y^fQ%s0$1{FV5Y&ch65@VebsW9!|Sm?42b#Y66INKfwPuxNed(cZL|F zAi4hQ;yfo6&nL4nE_dX=aG^_+5>DfLTxu;>`qOE&B2EJ5n3LgFXw_=G(y9=a#)vy> z^J2`Hb~+o)v%WUlde~3bvp7MuKDCE0WC|zchDV<1sGf*3;5i*FK z!B|LwFTZ%-81cJ@{Xyi58*MUf8?U@bg=osDb6;^{1c4DlD-m@-RGEsX3TlVa!ATBq z<=Bzqphks%uj8a`R8@k3=f$~u_wM;AtUo-1{t)Mb0S?a+9d!eUci+IjIJ)j;)7%+i zgo0%H+r?#0Dq7AKQ(W)Kf8k1(CLx?h_qf(t#`?o)v?5Ldmza~`R%q2~yii;SOH;(X zwRtgROgo*e=GjP_T|Mll>sg$kO1t*(nM~oN-0-buI;tn)0=W`-Dcem)*j$jt#ppAL zoWV#)f=|Er&=~Q%C!=xbi(73nZX2(>NQ7w0iR4AmKH;F{D&dlwsE)SjOP%}#)EoDj zqK@AwX|AJVE zI={ZVs0e3wU3Ff6M(gpWn>U}jJd+mQ|KwEu>pLeZzePT~C+_`!^;tue%-|C~TEVXq z_DUDkn4dnzW-NHKCBSsD%E#P!TdzOuT%W(|fyCXTr)KWB`uOQd^I-GbRc1eGWp48I zGCvL4via&ZiL_n8-Rg6WO+Kv=7`gkv88NrNeZu#@ad#HapSR@ljUC@gS32(xyZhwY zx3iTCXQ$mbv{rvzQS1w=*L~aG>(~7I-q3ex@xF`t^|NYDn9g6d*Goup?~TbSUTbqQ zW43oD=3ZY}a(;Hia{Faok>WmYOD6KxGe5t3?PlKeqQu!B7Is?yRCyIUYuT%*R&VS6 z6@L|j1S@v(gUHjxgH$tC%ENwlzPRDmec|6x! wpZeT;cGu$Ho_%)n;IcJ`-Z>P=6w7Z3*SdG^d(hI|bKmoev&_HyEs%i$08&}(Z~y=R literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz new file mode 100644 index 0000000000000000000000000000000000000000..2cb6c7e49efd9748c48e9a7f290445660dd352c6 GIT binary patch literal 368 zcmV-$0gwJ4iwFP!000001Fcd)Pr@)1zV}zWaMF>{;l%{vfnau$U=;D-JapRwYPxpm zI^rVyce{2d@d7ctCQbYLecyZUOKA|Eh#-V%zYpO_poJy`;O-FcH+J?lN&J}t#K2MH zzppSaQO5JhYy^u9x*I&JT4xB;_#PHYiBx@BjZVZ$z#K}{%y1<;jaN)7Zg2#!F{Z3W zlq##UQmhy#y>7{l-Ev_9GCyt)mZ&Q$Wy!Y=X`v1XGhmXJwdgACNOFc7XGcIuY#WTZ zK=S1ZZ+#=K_RSy~yW~pgf*Ic{t1`}gIuop_x)HWYEHN`m$E~&gU8wBGpkCDTNtJ$B zoF~1E$(~RO!NnG(gyg(H*$&b5_7=eq!XZR{y^jgu(x8@1#N9urC2wbhD~GtYh*Z@w zXS#8S&YK5Z8F%p4HcV6E@kgtgqClgYN5Lfef%k5I;_J=7$yk)8IYD-q18Gl)^Yv2E O`uYtL88pVQ0{{TY>ajHd literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz new file mode 100644 index 0000000000000000000000000000000000000000..e400169bc66ed0f4986bb9f2cc34944b709c3408 GIT binary patch literal 377 zcmV-<0fzn`iwFp6lA&S%16gEcW?^+~bO5bVu};G<5S`<%ICV_bI$}X;#84GwqAH>l z2oVM+xwJ-(9XU=#MEmd9PJoJ~9s0aDzIX55J%xsoGYBWJ7>{9cmT0Ld0l51C{Kvso z)3m=+fCM;-;;#y;8fCJYFK4jXq5p!lYITmVNbX^yl*rVl)fhyq1gxM&%?wv^(0IYD z;Ra^_J7el*MyaworDD}o>8)$_-L4nrAd6Fb@Qk{#QnvZlGcDB>o>|D`tAHtjRg=pL4;QW>{gb#5!k!>BL3b*M-V| z0~$r6{!wKgm=9X71(vK8CMHKhDTG%$loC?#3Z>g+u=p0i$$5_n;nJX%OeDkK&`Rz$ zg~;tqGS$k0>CRVeI>7^$h6MH!n65@>x~f8>yGOw!#w0kB;E8W{XDC?Jrt=i%$+J%2 XtDbly#7A;1xJbSM!FMQwasvPWa2B=R literal 0 HcmV?d00001 diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz new file mode 100644 index 0000000000000000000000000000000000000000..0eafccb139677bdac2ecdb7a177611e779fbf333 GIT binary patch literal 504 zcmb2|=3oE==C_v|vu_)S9Qe4~?&#v>x~fwr3%X6t4X~=oIvH`gETXJHX6eNZs>#-= z3(uCX*WDd0*)!$hVZQ&*C(PH~`&~uy^}!DgJZ(jFH9CizJAGd{b?&L>`Dy%a_RXV< zc3)T87%_Qm-T6~{4D{yBkM(<6m3(M&%+K7JQkMPpr?!5-dFg*n)feAcVn)*sP4Ud# z`?|Wu^z^A~X+Nj2?D`LRnqr%fsRowfVtndSTAE(Dfb;r;cLfAY7zT{CZKuTA}sH2s>_+MK+Y z>W;en`))=#b;+H(UjO#de4Q1%OQria|G~=brQ2+m$J|&|7v8$LM!>Xlk;mSk<%=)A ztPuLUy>HRQ*&kUT;Y$(u?^w2JaMsK1c~htCzA-Cj>ANW=g^AZncf6Fly`sEzqI%%l zPj{|7eHT2 Date: Tue, 19 Jul 2022 12:18:34 +0800 Subject: [PATCH 6/6] correct det_1D_H2_01mm --- examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz | Bin 374 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz | Bin 371 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz | Bin 373 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/T.gz | Bin 456 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz | Bin 368 -> 0 bytes .../dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz | Bin 377 -> 0 bytes examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz | Bin 504 -> 0 bytes 7 files changed, 0 insertions(+), 0 deletions(-) delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/O2.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/T.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz delete mode 100644 examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/H2.gz deleted file mode 100644 index a09d4b08708063b58a3429719b83315db118ffda..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 374 zcmV-+0g3(}iwFp6lA&S%14uFetx~aW!Y~k>b6;^{3<9S}RV$(nh$`w(L{(5bG#yMh zz?EZ1jzNtI|6a#QThs*xo)_ot-Mi|Hlan>#~H zP>@{yb#YpdN~WX95a%26UpUiSlM#-Sdz@=67y8p_v?5Ldr&y5TR%q2~ywai+mWGHM zYm0Kon07iV&GWuC>w4Hv*Rwc6wLG?mFJww5<%VxO(@{MUC&-m3wzA!{hxCjz&PK=} zaRw713BLT|U1P-W9`vKw7gyTk+%{f$nF-O9Gs(-ceZo%5S;i$dNgZv|mpc0ms1tRX zqF#JZ)@nTrVmJU}7+&l!2Cx)4srnVIxwi<9&U*|XlqI9MOxpiot3>@qjQY{w_WCMD zfQ8lQq?o2m`B0M@&$&Hw-a diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/N2.gz deleted file mode 100644 index cb797eb47bdc8a180bcf0aead79473163c350a20..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 371 zcmV-(0gV11iwFp6lA&S%15Ppktx~aW!Y~k>b6;^{3<4w63W=x#qRLcCRZu%L9ZWdD zm19SagBlh7y^fQ%s0$1{FV5Y&ch65@VebsW9!|Sm?42b#Y66INKfwPuxNed(cZL|F zAi4hQ;yfo6&nL4nE_dX=aG^_+5>DfLTxu;>`qOE&B2EJ5n3LgFXw_=G(y9=a#)vy> z^J2`Hb~+o)v%WUlde~3bvp7MuKDCE0WC|zchDV<1sGf*3;5i*FK z!B|LwFTZ%-81cJ@{Xyi58*MUf8?U@bg=osDb6;^{1c4DlD-m@-RGEsX3TlVa!ATBq z<=Bzqphks%uj8a`R8@k3=f$~u_wM;AtUo-1{t)Mb0S?a+9d!eUci+IjIJ)j;)7%+i zgo0%H+r?#0Dq7AKQ(W)Kf8k1(CLx?h_qf(t#`?o)v?5Ldmza~`R%q2~yii;SOH;(X zwRtgROgo*e=GjP_T|Mll>sg$kO1t*(nM~oN-0-buI;tn)0=W`-Dcem)*j$jt#ppAL zoWV#)f=|Er&=~Q%C!=xbi(73nZX2(>NQ7w0iR4AmKH;F{D&dlwsE)SjOP%}#)EoDj zqK@AwX|AJVE zI={ZVs0e3wU3Ff6M(gpWn>U}jJd+mQ|KwEu>pLeZzePT~C+_`!^;tue%-|C~TEVXq z_DUDkn4dnzW-NHKCBSsD%E#P!TdzOuT%W(|fyCXTr)KWB`uOQd^I-GbRc1eGWp48I zGCvL4via&ZiL_n8-Rg6WO+Kv=7`gkv88NrNeZu#@ad#HapSR@ljUC@gS32(xyZhwY zx3iTCXQ$mbv{rvzQS1w=*L~aG>(~7I-q3ex@xF`t^|NYDn9g6d*Goup?~TbSUTbqQ zW43oD=3ZY}a(;Hia{Faok>WmYOD6KxGe5t3?PlKeqQu!B7Is?yRCyIUYuT%*R&VS6 z6@L|j1S@v(gUHjxgH$tC%ENwlzPRDmec|6x! wpZeT;cGu$Ho_%)n;IcJ`-Z>P=6w7Z3*SdG^d(hI|bKmoev&_HyEs%i$08&}(Z~y=R diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/U.gz deleted file mode 100644 index 2cb6c7e49efd9748c48e9a7f290445660dd352c6..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 368 zcmV-$0gwJ4iwFP!000001Fcd)Pr@)1zV}zWaMF>{;l%{vfnau$U=;D-JapRwYPxpm zI^rVyce{2d@d7ctCQbYLecyZUOKA|Eh#-V%zYpO_poJy`;O-FcH+J?lN&J}t#K2MH zzppSaQO5JhYy^u9x*I&JT4xB;_#PHYiBx@BjZVZ$z#K}{%y1<;jaN)7Zg2#!F{Z3W zlq##UQmhy#y>7{l-Ev_9GCyt)mZ&Q$Wy!Y=X`v1XGhmXJwdgACNOFc7XGcIuY#WTZ zK=S1ZZ+#=K_RSy~yW~pgf*Ic{t1`}gIuop_x)HWYEHN`m$E~&gU8wBGpkCDTNtJ$B zoF~1E$(~RO!NnG(gyg(H*$&b5_7=eq!XZR{y^jgu(x8@1#N9urC2wbhD~GtYh*Z@w zXS#8S&YK5Z8F%p4HcV6E@kgtgqClgYN5Lfef%k5I;_J=7$yk)8IYD-q18Gl)^Yv2E O`uYtL88pVQ0{{TY>ajHd diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/Ydefault.gz deleted file mode 100644 index e400169bc66ed0f4986bb9f2cc34944b709c3408..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 377 zcmV-<0fzn`iwFp6lA&S%16gEcW?^+~bO5bVu};G<5S`<%ICV_bI$}X;#84GwqAH>l z2oVM+xwJ-(9XU=#MEmd9PJoJ~9s0aDzIX55J%xsoGYBWJ7>{9cmT0Ld0l51C{Kvso z)3m=+fCM;-;;#y;8fCJYFK4jXq5p!lYITmVNbX^yl*rVl)fhyq1gxM&%?wv^(0IYD z;Ra^_J7el*MyaworDD}o>8)$_-L4nrAd6Fb@Qk{#QnvZlGcDB>o>|D`tAHtjRg=pL4;QW>{gb#5!k!>BL3b*M-V| z0~$r6{!wKgm=9X71(vK8CMHKhDTG%$loC?#3Z>g+u=p0i$$5_n;nJX%OeDkK&`Rz$ zg~;tqGS$k0>CRVeI>7^$h6MH!n65@>x~f8>yGOw!#w0kB;E8W{XDC?Jrt=i%$+J%2 XtDbly#7A;1xJbSM!FMQwasvPWa2B=R diff --git a/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz b/examples/dfHighSpeedFoam/det_1D_H2_01mm/0/p.gz deleted file mode 100644 index 0eafccb139677bdac2ecdb7a177611e779fbf333..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 504 zcmb2|=3oE==C_v|vu_)S9Qe4~?&#v>x~fwr3%X6t4X~=oIvH`gETXJHX6eNZs>#-= z3(uCX*WDd0*)!$hVZQ&*C(PH~`&~uy^}!DgJZ(jFH9CizJAGd{b?&L>`Dy%a_RXV< zc3)T87%_Qm-T6~{4D{yBkM(<6m3(M&%+K7JQkMPpr?!5-dFg*n)feAcVn)*sP4Ud# z`?|Wu^z^A~X+Nj2?D`LRnqr%fsRowfVtndSTAE(Dfb;r;cLfAY7zT{CZKuTA}sH2s>_+MK+Y z>W;en`))=#b;+H(UjO#de4Q1%OQria|G~=brQ2+m$J|&|7v8$LM!>Xlk;mSk<%=)A ztPuLUy>HRQ*&kUT;Y$(u?^w2JaMsK1c~htCzA-Cj>ANW=g^AZncf6Fly`sEzqI%%l zPj{|7eHT2