diff --git a/applications/solvers/dfLowMachFoam/EEqn.H b/applications/solvers/dfLowMachFoam/EEqn.H index f8fc95351..0dae9da88 100644 --- a/applications/solvers/dfLowMachFoam/EEqn.H +++ b/applications/solvers/dfLowMachFoam/EEqn.H @@ -1,5 +1,13 @@ { volScalarField& he = thermo.he(); + tmp thcSource(new fvScalarMatrix(he, dimEnergy/dimTime)); + fvScalarMatrix& hcSource = thcSource.ref(); + forAll(Y, i) + { + hcSource.source() -= parcels.rhoTrans(i)*chemistry->mixture().Hc(i)/runTime.deltaT(); + } + //hSource.source() -= parcels.hsTrans()/runTime.deltaT();// equivalent to parcels.Sh(he) + start1 = std::clock(); fvScalarMatrix EEqn ( @@ -15,10 +23,16 @@ fvm::laplacian(turbulence->alpha(), he) - diffAlphaD + fvc::div(hDiffCorrFlux) + // + rho*(U&g) + + parcels.Sh(he) + + hcSource ) : ( fvm::laplacian(turbulence->alphaEff(), he) + // + rho*(U&g) + + parcels.Sh(he) + + hcSource ) ) ); diff --git a/applications/solvers/dfLowMachFoam/Make/options b/applications/solvers/dfLowMachFoam/Make/options index a0dc9ba1c..49835a766 100644 --- a/applications/solvers/dfLowMachFoam/Make/options +++ b/applications/solvers/dfLowMachFoam/Make/options @@ -6,18 +6,33 @@ EXE_INC = -std=c++14 \ -Wno-unused-variable \ -Wno-unused-but-set-variable \ -Wno-old-style-cast \ + -I. \ $(PFLAGS) $(PINC) \ $(if $(LIBTORCH_ROOT),-DUSE_LIBTORCH,) \ $(if $(PYTHON_INC_DIR),-DUSE_PYTORCH,) \ + -I$(FOAM_APP)/solvers/lagrangian/reactingParcelFoam \ + -I$(FOAM_APP)/solvers/compressible/rhoPimpleFoam \ -I$(LIB_SRC)/transportModels/compressible/lnInclude \ -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \ -I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \ -I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \ + -I$(FOAM_APP)/solvers/lagrangian/reactingParcelFoam \ -I$(LIB_SRC)/finiteVolume/cfdTools \ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude \ -I$(LIB_SRC)/sampling/lnInclude \ + -I$(LIB_SRC)/regionModels/regionModel/lnInclude \ + -I$(LIB_SRC)/regionModels/surfaceFilmModels/lnInclude \ -I$(LIB_SRC)/dynamicFvMesh/lnInclude \ + -I$(LIB_SRC)/lagrangian/basic/lnInclude \ + -I$(DF_SRC)/lagrangian/intermediate/lnInclude \ + -I$(LIB_SRC)/lagrangian/intermediate/lnInclude \ + -I$(DF_SRC)/lagrangian/spray/lnInclude \ + -I$(LIB_SRC)/lagrangian/spray/lnInclude \ + -I$(LIB_SRC)/lagrangian/distributionModels/lnInclude \ + -I$(DF_SRC)/thermophysicalModels/thermophysicalProperties/lnInclude \ + -I$(LIB_SRC)/thermophysicalModels/thermophysicalProperties/lnInclude \ + -I$(DF_SRC)/thermophysicalModels/SLGThermo/lnInclude \ -I$(LIB_SRC)/Pstream/mpi \ -I$(DF_SRC)/dfCanteraMixture/lnInclude \ -I$(DF_SRC)/dfChemistryModel/lnInclude \ @@ -35,12 +50,21 @@ EXE_INC = -std=c++14 \ EXE_LIBS = \ -lcompressibleTransportModels \ -lturbulenceModels \ + -llagrangian \ + -lregionModels \ + -ldfSurfaceFilmModels \ -lfiniteVolume \ + -ltopoChangerFvMesh \ -lmeshTools \ -lsampling \ -L$(DF_LIBBIN) \ -ldfFluidThermophysicalModels \ -ldfCompressibleTurbulenceModels \ + -ldfThermophysicalProperties \ + -ldfSLGThermo \ + -ldfLagrangianIntermediate \ + -ldfLagrangianTurbulence \ + -ldfLagrangianSpray \ -ldfCanteraMixture \ -ldfChemistryModel \ -ldfCombustionModels \ diff --git a/applications/solvers/dfLowMachFoam/UEqn.H b/applications/solvers/dfLowMachFoam/UEqn.H index 23a9bc7ba..77ad8087e 100644 --- a/applications/solvers/dfLowMachFoam/UEqn.H +++ b/applications/solvers/dfLowMachFoam/UEqn.H @@ -4,6 +4,9 @@ ( fvm::ddt(rho, U) + fvm::div(phi, U) + turbulence->divDevRhoReff(U) + == + // rho()*g + parcels.SU(U) ); fvVectorMatrix& UEqn = tUEqn.ref(); end1 = std::clock(); diff --git a/applications/solvers/dfLowMachFoam/YEqn.H b/applications/solvers/dfLowMachFoam/YEqn.H index a80452998..6203bbef9 100644 --- a/applications/solvers/dfLowMachFoam/YEqn.H +++ b/applications/solvers/dfLowMachFoam/YEqn.H @@ -107,8 +107,8 @@ == ( splitting - ? fvm::laplacian(DEff(), Yi) - : (fvm::laplacian(DEff(), Yi) + combustion->R(Yi)) + ? (fvm::laplacian(DEff(), Yi) + parcels.SYi(i, Yi) ) + : (fvm::laplacian(DEff(), Yi) + parcels.SYi(i, Yi) + combustion->R(Yi)) ) ); diff --git a/applications/solvers/dfLowMachFoam/createClouds.H b/applications/solvers/dfLowMachFoam/createClouds.H new file mode 100644 index 000000000..ee0985ff7 --- /dev/null +++ b/applications/solvers/dfLowMachFoam/createClouds.H @@ -0,0 +1,9 @@ +Info<< "\nConstructing reacting cloud" << endl; +basicSprayCloud parcels +( + "sprayCloud", + rho, + U, + g, + slgThermo +); diff --git a/applications/solvers/dfLowMachFoam/createFields.H b/applications/solvers/dfLowMachFoam/createFields.H index 42214ea22..444882e72 100644 --- a/applications/solvers/dfLowMachFoam/createFields.H +++ b/applications/solvers/dfLowMachFoam/createFields.H @@ -1,4 +1,5 @@ #include "createRDeltaT.H" +#include "readGravitationalAcceleration.H" Info<< "Reading thermophysical properties\n" << endl; @@ -7,6 +8,8 @@ fluidThermo* pThermo = new heRhoThermo(mesh, word::nu fluidThermo& thermo = *pThermo; // thermo.validate(args.executable(), "ha"); +SLGThermo slgThermo(mesh, thermo); + const volScalarField& psi = thermo.psi(); volScalarField& p = thermo.p(); volScalarField& T = thermo.T(); @@ -115,6 +118,9 @@ volScalarField K("K", 0.5*magSqr(U)); multivariateSurfaceInterpolationScheme::fieldTable fields; +#include "createMRF.H" +#include "createClouds.H" + if(combModelName!="ESF" && combModelName!="flareFGM" && combModelName!="DeePFGM" && combModelName!="FSD") { forAll(Y, i) diff --git a/applications/solvers/dfLowMachFoam/dfLowMachFoam.C b/applications/solvers/dfLowMachFoam/dfLowMachFoam.C index 590836ce0..b0bbd3288 100644 --- a/applications/solvers/dfLowMachFoam/dfLowMachFoam.C +++ b/applications/solvers/dfLowMachFoam/dfLowMachFoam.C @@ -64,6 +64,9 @@ Description #include "basicThermo.H" #include "CombustionModel.H" +#include "basicSprayCloud.H" +#include "SLGThermo.H" + //#define GPUSolver_ // #define TIME // #define DEBUG_ @@ -153,6 +156,7 @@ int main(int argc, char *argv[]) double time_monitor_Y = 0; double time_monitor_E = 0; double time_monitor_p = 0; + double time_monitor_parcels=0; double time_monitor_chemistry_correctThermo = 0; double time_monitor_turbulence_correct = 0; double time_monitor_chem = 0; // combustion correct @@ -261,6 +265,19 @@ int main(int argc, char *argv[]) runTime++; Info<< "Time = " << runTime.timeName() << nl << endl; + + if (parcels.solution().active()) + { + start = std::clock(); + + // Store the particle positions + parcels.storeGlobalPositions(); + + parcels.evolve(); + + end = std::clock(); + time_monitor_parcels += double(end - start) / double(CLOCKS_PER_SEC); + } // store old time fields #ifdef GPUSolver_ @@ -522,6 +539,10 @@ int main(int argc, char *argv[]) Info<< "Y Equations = " << time_monitor_Y - time_monitor_chem << " s" << endl; Info<< "E Equations = " << time_monitor_E << " s" << endl; Info<< "p Equations = " << time_monitor_p << " s" << endl; + if (parcels.solution().active()) + { + Info<< "calculate parcels = " << time_monitor_parcels << " s" << endl; + } Info<< "chemistry correctThermo = " << time_monitor_chemistry_correctThermo << " s" << endl; Info<< "turbulence correct = " << time_monitor_turbulence_correct << " s" << endl; Info<< "combustion correct(in Y) = " << time_monitor_chem << " s" << endl; diff --git a/applications/solvers/dfLowMachFoam/rhoEqn.H b/applications/solvers/dfLowMachFoam/rhoEqn.H index 7f7d7c15b..65f119319 100644 --- a/applications/solvers/dfLowMachFoam/rhoEqn.H +++ b/applications/solvers/dfLowMachFoam/rhoEqn.H @@ -34,6 +34,8 @@ Description ( fvm::ddt(rho) + fvc::div(phi) + == + parcels.Srho(rho) ); end1 = std::clock(); time_monitor_rhoEqn += double(end1 - start1) / double(CLOCKS_PER_SEC); diff --git a/examples/dfLowMachFoam/fgm/twoD_SandiaD_DeePFGM/constant/sprayCloudProperties b/examples/dfLowMachFoam/fgm/twoD_SandiaD_DeePFGM/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/fgm/twoD_SandiaD_DeePFGM/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/fgm/twoD_SandiaD_flareFGM/constant/g b/examples/dfLowMachFoam/fgm/twoD_SandiaD_flareFGM/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/fgm/twoD_SandiaD_flareFGM/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/fgm/twoD_SandiaD_flareFGM/constant/sprayCloudProperties b/examples/dfLowMachFoam/fgm/twoD_SandiaD_flareFGM/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/fgm/twoD_SandiaD_flareFGM/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/N2.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/N2.gz new file mode 100644 index 000000000..8b0f6005f Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/N2.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/O2.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/O2.gz new file mode 100644 index 000000000..9d47ae21a Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/O2.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/T.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/T.gz new file mode 100644 index 000000000..533967017 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/T.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/U.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/U.gz new file mode 100644 index 000000000..12cfcbeff Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/U.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/Ydefault.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/Ydefault.gz new file mode 100644 index 000000000..8000cee23 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/Ydefault.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/alphat.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/alphat.gz new file mode 100644 index 000000000..10e23cc04 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/alphat.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/epsilon.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/epsilon.gz new file mode 100644 index 000000000..d1307e06a Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/epsilon.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/k.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/k.gz new file mode 100644 index 000000000..8126dc7c3 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/k.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/nut.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/nut.gz new file mode 100644 index 000000000..b3022cbd4 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/nut.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/p.gz b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/p.gz new file mode 100644 index 000000000..33aac7d27 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/0/p.gz differ diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/Allclean b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/Allclean new file mode 100755 index 000000000..2cd9e2246 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/Allclean @@ -0,0 +1,9 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # Run from this directory + +echo "Cleaning log.*" +rm log.* +echo "Cleaning processor*" +rm -r processor* +echo "Cleaning polyMesh/" +rm -r constant/polyMesh diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/Allrun b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/Allrun new file mode 100755 index 000000000..0a0119f71 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/Allrun @@ -0,0 +1,11 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # Run from this directory + +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +application=dfSprayFoam + +runApplication blockMesh +runApplication decomposePar +runApplication mpirun -np 4 --allow-run-as-root $application -parallel diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/C7_oneStep.yaml b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/C7_oneStep.yaml new file mode 120000 index 000000000..c1698cb6c --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/C7_oneStep.yaml @@ -0,0 +1 @@ +../../../mechanisms/C7H16/C7_oneStep.yaml \ No newline at end of file diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/CanteraTorchProperties b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/CanteraTorchProperties new file mode 100644 index 000000000..9ffa283cf --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/CanteraTorchProperties @@ -0,0 +1,51 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object CanteraTorchProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +chemistry on; + +CanteraMechanismFile "C7_oneStep.yaml"; + +transportModel "Mix"; + +odeCoeffs +{ + "relTol" 1e-6; + "absTol" 1e-10; +} + +inertSpecie "N2"; + +splittingStrategy off; + +TorchSettings +{ + torch off; + GPU off; + log off; + torchModel ""; + coresPerNode 4; +} + +loadbalancing +{ + active true; + log false; + algorithm allAverage;//headTail; +} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/combustionProperties b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/combustionProperties new file mode 100644 index 000000000..7418500d0 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/combustionProperties @@ -0,0 +1,20 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object combustionProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +combustionModel laminar; + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/g b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/sprayCloudProperties new file mode 100644 index 000000000..505d69eb0 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/sprayCloudProperties @@ -0,0 +1,234 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active true; + coupled true; + transient yes; + cellValueSourceCorrection on; + maxCo 0.3; + + sourceTerms + { + schemes + { + rho explicit 1; + U explicit 1; + Yi explicit 1; + h explicit 1; + radiation explicit 1; + } + } + + interpolationSchemes + { + rho cell; + U cellPoint; + thermo:mu cell; + T cell; + Cp cell; + kappa cell; + p cell; + } + + integrationSchemes + { + U Euler; + T analytical; + } +} + + +constantProperties +{ + T0 320; + + // place holders for rho0 and Cp0 + // - reset from liquid properties using T0 + rho0 1000; + Cp0 4187; + + constantVolume false; +} + + +subModels +{ + particleForces + { + sphereDrag; + } + + injectionModels + { + model1 + { + type coneInjection; + SOI 0; + massTotal 6.0e-6; + parcelBasisType mass; + injectionMethod disc; + flowType flowRateAndDischarge; + dInner 0; + dOuter 1.9e-4; + duration 1.25e-3; + position (0 0.0995 0); + direction (0 -1 0); + parcelsPerSecond 20000000; + flowRateProfile table + ( + (0 0.1272) + (4.16667e-05 6.1634) + (8.33333e-05 9.4778) + (0.000125 9.5806) + (0.000166667 9.4184) + (0.000208333 9.0926) + (0.00025 8.7011) + (0.000291667 8.2239) + (0.000333333 8.0401) + (0.000375 8.8450) + (0.000416667 8.9174) + (0.000458333 8.8688) + (0.0005 8.8882) + (0.000541667 8.6923) + (0.000583333 8.0014) + (0.000625 7.2582) + (0.000666667 7.2757) + (0.000708333 6.9680) + (0.00075 6.7608) + (0.000791667 6.6502) + (0.000833333 6.7695) + (0.000875 5.5774) + (0.000916667 4.8649) + (0.000958333 5.0805) + (0.001 4.9547) + (0.00104167 4.5613) + (0.00108333 4.4536) + (0.001125 5.2651) + (0.00116667 5.2560) + (0.00120833 5.1737) + (0.00125 3.9213) + (0.001251 0.0000) + (1000 0.0000) + ); + + Cd constant 0.9; + + thetaInner constant 0.0; + thetaOuter constant 10.0; + + sizeDistribution + { + type RosinRammler; + + RosinRammlerDistribution + { + minValue 1e-06; + maxValue 0.00015; + d 0.00015; + n 3; + } + } + } + } + + dispersionModel none; + + patchInteractionModel standardWallInteraction; + + heatTransferModel RanzMarshall; + + compositionModel singlePhaseMixture; + + phaseChangeModel liquidEvaporationBoil; + + surfaceFilmModel none; + + atomizationModel none; + + breakupModel ReitzDiwakar; // ReitzKHRT; + + stochasticCollisionModel none; + + radiation off; + + standardWallInteractionCoeffs + { + type rebound; + } + + RanzMarshallCoeffs + { + BirdCorrection true; + } + + singlePhaseMixtureCoeffs + { + phases + ( + liquid + { + C7H16 1; + } + ); + } + + liquidEvaporationBoilCoeffs + { + enthalpyTransfer enthalpyDifference; + + activeLiquids ( C7H16 ); + } + + ReitzDiwakarCoeffs + { + solveOscillationEq yes; + Cbag 6; + Cb 0.785; + Cstrip 0.5; + Cs 10; + } + +/* + ReitzKHRTCoeffs + { + solveOscillationEq yes; + B0 0.61; + B1 40; + Ctau 1; + CRT 0.1; + msLimit 0.2; + WeberLimit 6; + } +*/ + TABCoeffs + { + y0 0; + yDot0 0; + Cmu 10; + Comega 8; + WeCrit 12; + } +} + + +cloudFunctions +{} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/thermophysicalProperties b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/thermophysicalProperties new file mode 100644 index 000000000..06be802bb --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/thermophysicalProperties @@ -0,0 +1,24 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +liquids +{ + C7H16; +} + +solids +{} +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/turbulenceProperties b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/turbulenceProperties new file mode 100644 index 000000000..762189e6a --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/constant/turbulenceProperties @@ -0,0 +1,30 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object turbulenceProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType RAS; + +RAS +{ + RASModel kEpsilon; + + turbulence on; + + printCoeffs on; +} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/blockMeshDict b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/blockMeshDict new file mode 100644 index 000000000..b7bfbe91a --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/blockMeshDict @@ -0,0 +1,57 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 0.001; + +vertices +( + (-10 0 -10) + (-10 0 10) + (10 0 10) + (10 0 -10) + (-10 100 -10) + (-10 100 10) + (10 100 10) + (10 100 -10) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (41 41 100) simpleGrading (1 1 1) +); + +edges +( +); + +patches +( + wall walls + ( + (2 6 5 1) + (0 4 7 3) + (0 1 5 4) + (4 5 6 7) + (7 6 2 3) + (3 2 1 0) + ) +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/controlDict b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/controlDict new file mode 100644 index 000000000..d646345d5 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/controlDict @@ -0,0 +1,53 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application sprayFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 1e-4; + +deltaT 2.5e-06; + +writeControl adjustableRunTime; + +writeInterval 1e-03; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression on; + +timeFormat general; + +timePrecision 6; + +adjustTimeStep yes; + +maxCo 0.1; + +runTimeModifiable yes; + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/decomposeParDict b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/decomposeParDict new file mode 100644 index 000000000..aa67a0627 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/decomposeParDict @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object decomposeParDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 4; + +method scotch; + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/fvSchemes b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/fvSchemes new file mode 100644 index 000000000..dd991a49b --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/fvSchemes @@ -0,0 +1,59 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + div(phi,U) Gauss upwind; + div(phid,p) Gauss upwind; + div(phi,K) Gauss linear; + div(phi,k) Gauss upwind; + div(phi,epsilon) Gauss upwind; + div(U) Gauss linear; + div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; + div(phi,Yi_h) Gauss upwind; + div(hDiffCorrFlux) Gauss cubic; +} + +laplacianSchemes +{ + default Gauss linear orthogonal; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default orthogonal; +} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/fvSolution b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/fvSolution new file mode 100644 index 000000000..5790575b0 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_aachenBomb/system/fvSolution @@ -0,0 +1,100 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + rho + { + solver PCG; + preconditioner DIC; + tolerance 1e-05; + relTol 0.1; + } + + + rhoFinal + { + $rho; + tolerance 1e-05; + relTol 0; + } + + "(U|k|epsilon)" + { + solver smoothSolver; + smoother symGaussSeidel; + tolerance 1e-06; + relTol 0.1; + } + + p + { + solver GAMG; + tolerance 0; + relTol 0.1; + smoother GaussSeidel; + } + + pFinal + { + $p; + tolerance 1e-06; + relTol 0; + } + + "(U|k|epsilon)Final" + { + $U; + tolerance 1e-06; + relTol 0; + } + + "(ha|Yi|O2|N2|H2O)" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-6; + relTol 0.1; + } + + "(ha|Yi|O2|N2|H2O)Final" + { + $Yi; + relTol 0; + } + + + +} + +PIMPLE +{ + transonic no; + nCorrectors 2; + nNonOrthogonalCorrectors 0; + momentumPredictor yes; +} + +relaxationFactors +{ + equations + { + ".*" 1; + } +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_reactingTGV/H2/cvodeIntegrator/constant/g b/examples/dfLowMachFoam/notorch/threeD_reactingTGV/H2/cvodeIntegrator/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_reactingTGV/H2/cvodeIntegrator/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/threeD_reactingTGV/H2/cvodeIntegrator/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/threeD_reactingTGV/H2/cvodeIntegrator/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/threeD_reactingTGV/H2/cvodeIntegrator/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_HIT_flame/CH4/constant/g b/examples/dfLowMachFoam/notorch/twoD_HIT_flame/CH4/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_HIT_flame/CH4/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_HIT_flame/CH4/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_HIT_flame/CH4/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_HIT_flame/CH4/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_SandiaD_EDC_Trans/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_SandiaD_EDC_Trans/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_SandiaD_EDC_Trans/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_reactingTGV/CH4/constant/g b/examples/dfLowMachFoam/notorch/twoD_reactingTGV/CH4/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_reactingTGV/CH4/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_reactingTGV/CH4/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_reactingTGV/CH4/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_reactingTGV/CH4/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_supercirticalTGV/constant/g b/examples/dfLowMachFoam/notorch/twoD_supercirticalTGV/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_supercirticalTGV/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_supercirticalTGV/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_supercirticalTGV/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_supercirticalTGV/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/C2H5OH.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/C2H5OH.gz new file mode 100644 index 000000000..5fda47883 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/C2H5OH.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/CO2.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/CO2.gz new file mode 100644 index 000000000..2c6d7a110 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/CO2.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/H2O.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/H2O.gz new file mode 100644 index 000000000..081e153db Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/H2O.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/N2.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/N2.gz new file mode 100644 index 000000000..b0f607672 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/N2.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/O2.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/O2.gz new file mode 100644 index 000000000..b87d87b9c Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/O2.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/T.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/T.gz new file mode 100644 index 000000000..e63ef0567 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/T.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/U.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/U.gz new file mode 100644 index 000000000..b949ecc75 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/U.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/Ydefault.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/Ydefault.gz new file mode 100644 index 000000000..f168cce92 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/Ydefault.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/alphat.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/alphat.gz new file mode 100644 index 000000000..8a36ba1be Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/alphat.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/epsilon.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/epsilon.gz new file mode 100644 index 000000000..17beb3229 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/epsilon.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/k.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/k.gz new file mode 100644 index 000000000..358ce9180 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/k.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/nut.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/nut.gz new file mode 100644 index 000000000..e5aeb2a80 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/nut.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/p.gz b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/p.gz new file mode 100644 index 000000000..e9d947dd8 Binary files /dev/null and b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/0/p.gz differ diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/Allclean b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/Allclean new file mode 100755 index 000000000..2cd9e2246 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/Allclean @@ -0,0 +1,9 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # Run from this directory + +echo "Cleaning log.*" +rm log.* +echo "Cleaning processor*" +rm -r processor* +echo "Cleaning polyMesh/" +rm -r constant/polyMesh diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/Allrun b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/Allrun new file mode 100755 index 000000000..206b5a7f5 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/Allrun @@ -0,0 +1,11 @@ +#!/bin/sh +cd ${0%/*} || exit 1 # Run from this directory + +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +application=dfSprayFoam + +runApplication blockMesh +runApplication decomposePar +runApplication mpirun -np 4 --allow-run-as-root $application -parallel diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/CanteraTorchProperties b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/CanteraTorchProperties new file mode 100644 index 000000000..5ef6c6c54 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/CanteraTorchProperties @@ -0,0 +1,51 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object CanteraTorchProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +chemistry on; + +CanteraMechanismFile "ethanol66.yaml"; + +transportModel "Mix"; + +odeCoeffs +{ + "relTol" 1e-6; + "absTol" 1e-10; +} + +inertSpecie "N2"; + +splittingStrategy off; + +TorchSettings +{ + torch off; + GPU off; + log off; + torchModel ""; + coresPerNode 4; +} + +loadbalancing +{ + active true; + log false; + algorithm allAverage;//headTail; +} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/combustionProperties b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/combustionProperties new file mode 100644 index 000000000..7a5f8b79b --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/combustionProperties @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object combustionProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +combustionModel laminar;//EDC + +active on; + +EDCCoeffs +{ + version v2005; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/g b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/g new file mode 100644 index 000000000..a23db6c71 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 0 -9.81); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/radiationProperties b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/radiationProperties new file mode 100644 index 000000000..f151acfff --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/radiationProperties @@ -0,0 +1,20 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object radiationProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +radiationModel none; + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/sprayCloudProperties new file mode 100644 index 000000000..ecd2b6878 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/sprayCloudProperties @@ -0,0 +1,193 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active true; + coupled true; + transient yes; + cellValueSourceCorrection on; + maxCo 0.1; + + sourceTerms + { + schemes + { + rho explicit 1; + U explicit 1; + Yi explicit 1; + h explicit 1; + radiation explicit 1; + } + } + + interpolationSchemes + { + rho cell; + U cellPoint; + thermo:mu cell; + T cell; + Cp cell; + kappa cell; + p cell; + } + + integrationSchemes + { + U Euler; + T analytical; + } +} + + +constantProperties +{ + T0 267; + + // place holders for rho0 and Cp0 + // - reset from liquid properties using T0 + rho0 1000; + Cp0 4187; + + constantVolume false; +} + + +subModels +{ + particleForces + { + sphereDrag; + gravity; + } + + injectionModels + { + model1 + { + type patchInjection; + patchName inletFuel; + duration 1; + parcelsPerSecond 1100000; //55000 + massTotal 0.24e-3; + parcelBasisType mass; + //massFlowRate 0.24e-03;//kg/s + flowRateProfile constant 1; + //nParticle 1; + SOI 0.0; + U0 (0 0 24); + //Cd constant 0.9; + sizeDistribution + { + type RosinRammler; + + RosinRammlerDistribution + { + minValue 5e-06; + maxValue 0.0001; + d 0.00002; + n 3; + } + } + } + } + + dispersionModel stochasticDispersionRAS; + + patchInteractionModel standardWallInteraction; + + heatTransferModel RanzMarshall; + + compositionModel singlePhaseMixture; + + phaseChangeModel liquidEvaporation; //liquidEvaporationBoil; + + surfaceFilmModel none; + + atomizationModel none; + + breakupModel ReitzDiwakar; // ReitzKHRT; + + stochasticCollisionModel none; + + radiation off; + + standardWallInteractionCoeffs + { + type rebound; + } + + RanzMarshallCoeffs + { + BirdCorrection true; + } + + singlePhaseMixtureCoeffs + { + phases + ( + liquid + { + C2H5OH 1; + } + ); + } + + liquidEvaporationCoeffs + { + enthalpyTransfer enthalpyDifference; + + activeLiquids ( C2H5OH ); + } + + ReitzDiwakarCoeffs + { + solveOscillationEq yes; + Cbag 6; + Cb 0.785; + Cstrip 0.5; + Cs 10; + } + +/* + ReitzKHRTCoeffs + { + solveOscillationEq yes; + B0 0.61; + B1 40; + Ctau 1; + CRT 0.1; + msLimit 0.2; + WeberLimit 6; + } +*/ + TABCoeffs + { + y0 0; + yDot0 0; + Cmu 10; + Comega 8; + WeCrit 12; + } +} + + +cloudFunctions +{} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/thermophysicalProperties b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/thermophysicalProperties new file mode 100644 index 000000000..c2ad7aa4e --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/thermophysicalProperties @@ -0,0 +1,40 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "constant"; + object thermophysicalProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +thermoType +{ + type hePsiThermo; + mixture reactingMixture; + transport sutherland; + thermo janaf; + energy sensibleEnthalpy; + equationOfState perfectGas; + specie specie; +} + +inertSpecie N2; + +liquids +{ + C2H5OH; +} + +solids +{} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/turbulenceProperties b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/turbulenceProperties new file mode 100644 index 000000000..278cffffc --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/constant/turbulenceProperties @@ -0,0 +1,37 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object turbulenceProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +simulationType RAS; + +RAS +{ + RASModel kEpsilon; + + turbulence on; + + printCoeffs on; + + kEpsilonCoeffs + { + C1 1.6; + C2 1.92; + } +} + + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/ethanol66.yaml b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/ethanol66.yaml new file mode 120000 index 000000000..1e2c35bb2 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/ethanol66.yaml @@ -0,0 +1 @@ +../../../mechanisms/C2H5OH/ethanol66.yaml \ No newline at end of file diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/blockMeshDict b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/blockMeshDict new file mode 100644 index 000000000..c9f105810 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/blockMeshDict @@ -0,0 +1,176 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 0.001; + +vertices +( + (0 0 -50) // 0 + + (5.25 -0.18333403983 -50) // 1 + (5.25 0.18333403983 -50) // 2 + + (5.75 -0.20079442457 -50) // 3 + (5.75 0.20079442457 -50) // 4 + + (12.5 -0.43650961864 -50) // 5 + (12.5 0.43650961864 -50) // 6 + + (0 0 0) // 7 + + (5.25 -0.18333403983 0) // 8 + (5.25 0.18333403983 0) // 9 + + (5.75 -0.20079442457 0) // 10 + (5.75 0.20079442457 0) // 11 + + (12.5 -0.43650961864 0) // 12 + (12.5 0.43650961864 0) // 13 + + + (104 -3.63176002714 0) // 14 + (104 3.63176002714 0) // 15 + + (0 0 500) // 16 + + (5.25 -0.18333403983 500) // 17 + (5.25 0.18333403983 500) // 18 + + (5.75 -0.20079442457 500) // 19 + (5.75 0.20079442457 500) // 20 + + (12.5 -0.43650961864 500) // 21 + (12.5 0.43650961864 500) // 22 + + (104 -3.63176002714 500) // 23 + (104 3.63176002714 500) // 24 +); + +blocks +( + hex ( 0 1 2 0 7 8 9 7) (11 1 20) simpleGrading (1 1 1) + hex ( 3 5 6 4 10 12 13 11) (13 1 20) simpleGrading (1 1 1) + + hex ( 7 8 9 7 16 17 18 16) (11 1 200) simpleGrading (1 1 2) + hex ( 8 10 11 9 17 19 20 18) (1 1 200) simpleGrading (1 1 2) + hex (10 12 13 11 19 21 22 20) (13 1 200) simpleGrading (1 1 2) + hex (12 14 15 13 21 23 24 22) (100 1 200) simpleGrading (3 1 2) +); + +boundary +( + inletFuel + { + type patch; + faces + ( + (1 0 0 2) + ); + } + + wallOutside + { + type wall; + faces + ( + (14 15 24 23) + ); + } + + wallTube + { + type wall; + faces + ( + (1 2 9 8) + (10 11 9 8) + (4 3 10 11) + (5 6 13 12) + ); + } + + inletPilot + { + type patch; + faces + ( + (5 3 4 6) + ); + } + + inletAir + { + type patch; + faces + ( + (14 12 13 15) + ); + } + + outlet + { + type patch; + faces + ( + (16 17 18 16) + (17 19 20 18) + (19 21 22 20) + (21 22 24 23) + ); + } + + axis + { + type empty; + faces + ( + (0 7 7 0) + (7 16 16 7) + ); + } + + frontAndBack_pos + { + type wedge; + faces + ( + (2 0 7 9) + (6 4 11 13) + + (9 7 16 18) + (11 9 18 20) + (13 11 20 22) + (15 13 22 24) + ); + } + + frontAndBack_neg + { + type wedge; + faces + ( + (0 1 8 7) + (3 5 12 10) + + (7 8 17 16) + (8 10 19 17) + (10 12 21 19) + (12 14 23 21) + ); + } +); + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/controlDict b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/controlDict new file mode 100644 index 000000000..d3be94a1d --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/controlDict @@ -0,0 +1,62 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application sprayFoam; + +startFrom latestTime; + +startTime 0.0; + +stopAt endTime; + +endTime 0.15; + +deltaT 2e-07; + +writeControl adjustableRunTime; + +writeInterval 0.002; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 12; + +writeCompression on; + +timeFormat general; + +timePrecision 6; + +adjustTimeStep yes; + +maxCo 0.3; + +runTimeModifiable true; + +OptimisationSwitches +{ + fileHandler collated; + //- collated: thread buffer size for queued file writes. + // If set to 0 or not sufficient for the file size, threading is not used. + // Default: 1e9 + maxThreadFileBufferSize 0;//1e9; +} + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/decomposeParDict b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/decomposeParDict new file mode 100644 index 000000000..9a79d19bd --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/decomposeParDict @@ -0,0 +1,35 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 4.x | +| \\ / A nd | Web: www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object decomposeParDict; +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +numberOfSubdomains 4; + +method simple; + +simpleCoeffs +{ + n (1 1 4 ); + delta 0.001; +} + + + +coeffs +{ + n (1 1 40); +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/fvSchemes b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/fvSchemes new file mode 100644 index 000000000..55bfcf337 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/fvSchemes @@ -0,0 +1,70 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; +} + +divSchemes +{ + default none; + + div(phi,U) Gauss limitedLinearV 1; + div(phi,Yi) Gauss limitedLinear01 1; + div(phi,ha) Gauss limitedLinear 1; + div(phi,K) Gauss limitedLinear 1; + div(phid,p) Gauss limitedLinear 1; + div(phi,epsilon) Gauss limitedLinear 1; + div(phi,Yi_h) Gauss limitedLinear01 1; + div(phi,k) Gauss limitedLinear 1; + div(hDiffCorrFlux) Gauss linear; + div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; +} + +laplacianSchemes +{ + default Gauss linear orthogonal; +} + +interpolationSchemes +{ + default linear; +} + +snGradSchemes +{ + default orthogonal; +} + +fluxRequired +{ + default no; + p ; +} + +wallDist +{ + method meshWave; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/fvSolution b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/fvSolution new file mode 100644 index 000000000..4c8ac09e9 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_sydneySprayBurner/system/fvSolution @@ -0,0 +1,99 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + location "system"; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + "rho.*" + { + solver diagonal; + } + +/* p + { + solver PCG; + preconditioner DIC; + tolerance 1e-6; + relTol 0.01; + } + + pFinal + { + $p; + relTol 0; + } +*/ + p + { + solver GAMG; + tolerance 1e-6; + relTol 0.01; + smoother GaussSeidel; + } + + pFinal + { + solver GAMG; + tolerance 1e-6; + relTol 0; + smoother GaussSeidel; + } + + + "(U|ha|k|epsilon)" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-6; + relTol 0.01; + } + + "(U|ha|k|epsilon)Final" + { + $U; + relTol 0; + } + + "Yi.*" + { + solver PBiCGStab; + preconditioner DILU; + tolerance 1e-6; + relTol 0.01; + } +} + +PIMPLE +{ + transonic no; + momentumPredictor yes; + nOuterCorrectors 1; + nCorrectors 2; + nNonOrthogonalCorrectors 1; + +// maxDeltaT 1e-5; +// maxCo 0.25; +} + +relaxationFactors +{ + equations + { + ".*" 1; + } +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/CH4/constant/g b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/CH4/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/CH4/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/CH4/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/CH4/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_temporally_evolving_jet/CH4/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_tripleFlame/H2/constant/g b/examples/dfLowMachFoam/notorch/twoD_tripleFlame/H2/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_tripleFlame/H2/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/examples/dfLowMachFoam/notorch/twoD_tripleFlame/H2/constant/sprayCloudProperties b/examples/dfLowMachFoam/notorch/twoD_tripleFlame/H2/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/examples/dfLowMachFoam/notorch/twoD_tripleFlame/H2/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/test/Tu500K-Phi1/constant/g b/test/Tu500K-Phi1/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/test/Tu500K-Phi1/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/test/Tu500K-Phi1/constant/sprayCloudProperties b/test/Tu500K-Phi1/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/test/Tu500K-Phi1/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* // diff --git a/test/corrtest.cpp b/test/corrtest.cpp index f97ea9f83..2ae8a26a5 100644 --- a/test/corrtest.cpp +++ b/test/corrtest.cpp @@ -38,7 +38,6 @@ float T10 = readSandia(10,"2DSandia/data_T.xy"); float T11 = readSandia(11,"2DSandia/data_T.xy"); - TEST(corrtest,dfHighSpeedFoam){ EXPECT_NEAR(v,1979.33,19.79); // within 1% of the theroetical value } diff --git a/test/dfLowMachFoam/twoD_reactingTGV/H2/cvodeSolver/constant/g b/test/dfLowMachFoam/twoD_reactingTGV/H2/cvodeSolver/constant/g new file mode 100644 index 000000000..28cfc2369 --- /dev/null +++ b/test/dfLowMachFoam/twoD_reactingTGV/H2/cvodeSolver/constant/g @@ -0,0 +1,22 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class uniformDimensionedVectorField; + location "constant"; + object g; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -2 0 0 0 0]; +value (0 -9.81 0); + + +// ************************************************************************* // diff --git a/test/dfLowMachFoam/twoD_reactingTGV/H2/cvodeSolver/constant/sprayCloudProperties b/test/dfLowMachFoam/twoD_reactingTGV/H2/cvodeSolver/constant/sprayCloudProperties new file mode 100644 index 000000000..2010a8858 --- /dev/null +++ b/test/dfLowMachFoam/twoD_reactingTGV/H2/cvodeSolver/constant/sprayCloudProperties @@ -0,0 +1,26 @@ +/*--------------------------------*- C++ -*----------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | Website: https://openfoam.org + \\ / A nd | Version: 7 + \\/ M anipulation | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format binary; + class dictionary; + location "constant"; + object SprayCloudProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solution +{ + active false; + coupled false; + transient no; + cellValueSourceCorrection off; +} + +// ************************************************************************* //