From 7d1e727372cbbe37d375efa35d885a657dd2aff0 Mon Sep 17 00:00:00 2001 From: andy Date: Tue, 1 Mar 2011 17:52:10 +0000 Subject: [PATCH 1/8] ENH: pass args by reference instead of copy --- .../IntegrationScheme/Analytical/Analytical.C | 4 ++-- .../IntegrationScheme/Analytical/Analytical.H | 4 ++-- .../intermediate/IntegrationScheme/Euler/Euler.C | 4 ++-- .../intermediate/IntegrationScheme/Euler/Euler.H | 4 ++-- .../IntegrationScheme/IntegrationScheme.C | 8 ++++---- .../IntegrationScheme/IntegrationScheme.H | 4 ++-- 6 files changed, 14 insertions(+), 14 deletions(-) diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C index fb237b0d04..89c72685a2 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C @@ -59,9 +59,9 @@ template typename Foam::IntegrationScheme::integrationResult Foam::Analytical::integrate ( - const Type phi, + const Type& phi, const scalar dt, - const Type alpha, + const Type& alpha, const scalar beta ) const { diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H index 15697136b2..91a188dc3a 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H @@ -81,9 +81,9 @@ public: //- Perform the integration virtual typename IntegrationScheme::integrationResult integrate ( - const Type phi, + const Type& phi, const scalar dt, - const Type alpha, + const Type& alpha, const scalar beta ) const; }; diff --git a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C index 8e7e1aea2b..1ebaa9dc29 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C @@ -59,9 +59,9 @@ template typename Foam::IntegrationScheme::integrationResult Foam::Euler::integrate ( - const Type phi, + const Type& phi, const scalar dt, - const Type alpha, + const Type& alpha, const scalar beta ) const { diff --git a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H index 1140713b13..de58579a9d 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H +++ b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H @@ -78,9 +78,9 @@ public: //- Perform the integration virtual typename IntegrationScheme::integrationResult integrate ( - const Type phi, + const Type& phi, const scalar dt, - const Type alpha, + const Type& alpha, const scalar beta ) const; }; diff --git a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C index a91e45d90e..d3cac7f4e0 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C +++ b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C @@ -60,9 +60,9 @@ template typename Foam::IntegrationScheme::integrationResult Foam::IntegrationScheme::integrate ( - const Type phi, + const Type& phi, const scalar dt, - const Type alpha, + const Type& alpha, const scalar beta ) const { @@ -71,9 +71,9 @@ Foam::IntegrationScheme::integrate "Foam::IntegrationScheme::integrationResult" "Foam::IntegrationScheme::integrate" "(" - "const Type, " + "const Type&, " "const scalar, " - "const Type, " + "const Type&, " "const scalar" ") const" ); diff --git a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H index 93726d261f..d7f0d0ef07 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H +++ b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H @@ -183,9 +183,9 @@ public: //- Perform the Integration virtual integrationResult integrate ( - const Type phi, + const Type& phi, const scalar dt, - const Type alpha, + const Type& alpha, const scalar beta ) const; }; From b2de414a997a46be989a7d190e9da25519680e9b Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 11:00:13 +0000 Subject: [PATCH 2/8] ENH: Updated steadyReactingParcelFoam -> LTSReactingParcelFoam --- .../LTSReactingParcelFoam.C} | 46 ++++++++----------- .../LTSReactingParcelFoam/Make/files | 3 ++ .../Make/options | 0 .../UEqn.H | 0 .../YEqn.H | 0 .../chemistry.H | 0 .../createClouds.H | 0 .../createExplicitSources.H | 0 .../createFields.H | 4 +- .../createPorousZones.H | 0 .../hsEqn.H | 0 .../pEqn.H | 3 -- .../readAdditionalSolutionControls.H | 0 .../readChemistryProperties.H | 0 .../readTimeControls.H | 2 +- .../rhoEqn.H | 2 +- .../setPressureWork.H | 0 .../timeScales.H | 2 +- .../steadyReactingParcelFoam/Make/files | 3 -- .../counterFlowFlame2D/0/CH4 | 0 .../counterFlowFlame2D/0/CO2 | 0 .../counterFlowFlame2D/0/H2O | 0 .../counterFlowFlame2D/0/N2 | 0 .../counterFlowFlame2D/0/O2 | 0 .../counterFlowFlame2D/0/T | 0 .../counterFlowFlame2D/0/U | 0 .../counterFlowFlame2D/0/alphat | 0 .../counterFlowFlame2D/0/epsilon | 0 .../counterFlowFlame2D/0/k | 0 .../counterFlowFlame2D/0/mut | 0 .../counterFlowFlame2D/0/p | 0 .../counterFlowFlame2D/0/rho | 0 .../counterFlowFlame2D/Allrun | 2 +- .../counterFlowFlame2D/constant/RASProperties | 0 .../constant/chemistryProperties | 0 .../constant/energySourceProperties | 0 .../counterFlowFlame2D/constant/g | 0 .../constant/massSourceProperties | 0 .../constant/momentumSourceProperties | 0 .../constant/polyMesh/blockMeshDict | 0 .../constant/polyMesh/boundary | 0 .../constant/radiationProperties | 0 .../constant/reactingCloud1Properties | 0 .../counterFlowFlame2D/constant/reactions | 0 .../constant/thermo.icoPoly8 | 0 .../constant/thermophysicalProperties | 0 .../constant/turbulenceProperties | 0 .../counterFlowFlame2D/system/controlDict | 4 ++ .../counterFlowFlame2D/system/fvSchemes | 0 .../counterFlowFlame2D/system/fvSolution | 3 +- .../verticalChannel/0.org/H2O | 0 .../verticalChannel/0.org/T | 0 .../verticalChannel/0.org/U | 0 .../verticalChannel/0.org/air | 0 .../verticalChannel/0.org/alphat | 0 .../verticalChannel/0.org/k | 0 .../verticalChannel/0.org/mut | 0 .../verticalChannel/0.org/omega | 0 .../verticalChannel/0.org/p | 0 .../verticalChannel/0/H2O | 0 .../verticalChannel/0/T | 0 .../verticalChannel/0/U | 0 .../verticalChannel/0/air | 0 .../verticalChannel/0/alphat | 0 .../verticalChannel/0/k | 0 .../verticalChannel/0/mut | 0 .../verticalChannel/0/omega | 0 .../verticalChannel/0/p | 0 .../verticalChannel/Allclean | 0 .../verticalChannel/Allrun | 2 +- .../verticalChannel/constant/RASProperties | 0 .../constant/chemistryProperties | 0 .../constant/energySourceProperties | 0 .../verticalChannel/constant/g | 0 .../constant/massSourceProperties | 0 .../constant/momentumSourceProperties | 0 .../constant/particleTrackDict | 0 .../constant/polyMesh/blockMeshDict | 0 .../constant/polyMesh/boundary | 0 .../verticalChannel/constant/porousZones | 0 .../constant/radiationProperties | 0 .../constant/reactingCloud1Properties | 0 .../verticalChannel/constant/reactions | 0 .../constant/thermo.incompressiblePoly | 0 .../constant/thermophysicalProperties | 0 .../constant/turbulenceProperties | 0 .../verticalChannel/system/controlDict | 0 .../verticalChannel/system/fvSchemes | 0 .../verticalChannel/system/fvSolution | 3 +- 89 files changed, 38 insertions(+), 41 deletions(-) rename applications/solvers/lagrangian/{steadyReactingParcelFoam/steadyReactingParcelFoam.C => LTSReactingParcelFoam/LTSReactingParcelFoam.C} (78%) create mode 100644 applications/solvers/lagrangian/LTSReactingParcelFoam/Make/files rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/Make/options (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/UEqn.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/YEqn.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/chemistry.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/createClouds.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/createExplicitSources.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/createFields.H (96%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/createPorousZones.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/hsEqn.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/pEqn.H (95%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/readAdditionalSolutionControls.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/readChemistryProperties.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/readTimeControls.H (97%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/rhoEqn.H (95%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/setPressureWork.H (100%) rename applications/solvers/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/timeScales.H (99%) delete mode 100644 applications/solvers/lagrangian/steadyReactingParcelFoam/Make/files rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/CH4 (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/CO2 (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/H2O (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/N2 (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/O2 (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/T (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/U (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/alphat (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/epsilon (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/k (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/mut (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/p (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/0/rho (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/Allrun (82%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/RASProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/chemistryProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/energySourceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/g (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/massSourceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/momentumSourceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/polyMesh/blockMeshDict (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/polyMesh/boundary (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/radiationProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/reactingCloud1Properties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/reactions (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/thermo.icoPoly8 (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/thermophysicalProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/constant/turbulenceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/system/controlDict (96%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/system/fvSchemes (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/counterFlowFlame2D/system/fvSolution (98%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/H2O (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/T (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/U (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/air (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/alphat (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/k (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/mut (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/omega (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0.org/p (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/H2O (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/T (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/U (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/air (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/alphat (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/k (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/mut (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/omega (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/0/p (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/Allclean (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/Allrun (87%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/RASProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/chemistryProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/energySourceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/g (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/massSourceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/momentumSourceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/particleTrackDict (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/polyMesh/blockMeshDict (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/polyMesh/boundary (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/porousZones (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/radiationProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/reactingCloud1Properties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/reactions (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/thermo.incompressiblePoly (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/thermophysicalProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/constant/turbulenceProperties (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/system/controlDict (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/system/fvSchemes (100%) rename tutorials/lagrangian/{steadyReactingParcelFoam => LTSReactingParcelFoam}/verticalChannel/system/fvSolution (98%) diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/steadyReactingParcelFoam.C b/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C similarity index 78% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/steadyReactingParcelFoam.C rename to applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C index a6aa1290a3..1160760bb6 100644 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/steadyReactingParcelFoam.C +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/LTSReactingParcelFoam.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2008-2010 OpenCFD Ltd. + \\ / A nd | Copyright (C) 2010-2011 OpenCFD Ltd. \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -22,17 +22,13 @@ License along with OpenFOAM. If not, see . Application - porousExplicitSourceReactingParcelFoam + LTSReactingParcelFoam Description - Transient PISO solver for compressible, laminar or turbulent flow with - reacting multiphase Lagrangian parcels for porous media, including explicit - sources for mass, momentum and energy - - The solver includes: - - reacting multiphase parcel cloud - - porous media - - mass, momentum and energy sources + Local time stepping (LTS) solver for steady, compressible, laminar or + turbulent reacting and non-reacting flow with multiphase Lagrangian + parcels and porous media, including explicit sources for mass, momentum + and energy Note: ddtPhiCorr not used here when porous zones are active - not well defined for porous calculations @@ -75,7 +71,7 @@ int main(int argc, char *argv[]) while (runTime.run()) { - #include "readSIMPLEControls.H" + #include "readPISOControls.H" #include "readChemistryProperties.H" #include "readAdditionalSolutionControls.H" #include "readTimeControls.H" @@ -84,26 +80,24 @@ int main(int argc, char *argv[]) Info<< "Time = " << runTime.timeName() << nl << endl; - p.storePrevIter(); + parcels.evolve(); - // --- Pressure-velocity corrector + #include "chemistry.H" + #include "timeScales.H" + + #include "rhoEqn.H" + #include "UEqn.H" + #include "YEqn.H" + #include "hsEqn.H" + + // --- PISO loop + for (int corr=0; corrcorrect(); } + turbulence->correct(); + if (runTime.write()) { chemistry.dQ()().write(); diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/Make/files b/applications/solvers/lagrangian/LTSReactingParcelFoam/Make/files new file mode 100644 index 0000000000..1efa8dad41 --- /dev/null +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/Make/files @@ -0,0 +1,3 @@ +LTSReactingParcelFoam.C + +EXE = $(FOAM_APPBIN)/LTSReactingParcelFoam diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/Make/options b/applications/solvers/lagrangian/LTSReactingParcelFoam/Make/options similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/Make/options rename to applications/solvers/lagrangian/LTSReactingParcelFoam/Make/options diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/UEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/UEqn.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/UEqn.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/YEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/YEqn.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/YEqn.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/YEqn.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/chemistry.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/chemistry.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/chemistry.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/chemistry.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/createClouds.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/createClouds.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/createClouds.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/createClouds.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/createExplicitSources.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/createExplicitSources.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/createExplicitSources.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/createExplicitSources.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/createFields.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/createFields.H similarity index 96% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/createFields.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/createFields.H index 68734ab9a7..5ca55b84f6 100644 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/createFields.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/createFields.H @@ -73,12 +73,12 @@ dimensionedScalar rhoMax ( - mesh.solutionDict().subDict("SIMPLE").lookup("rhoMax") + mesh.solutionDict().subDict("PISO").lookup("rhoMax") ); dimensionedScalar rhoMin ( - mesh.solutionDict().subDict("SIMPLE").lookup("rhoMin") + mesh.solutionDict().subDict("PISO").lookup("rhoMin") ); Info<< "Creating turbulence model\n" << endl; diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/createPorousZones.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/createPorousZones.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/createPorousZones.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/createPorousZones.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/hsEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/hsEqn.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/hsEqn.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/hsEqn.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H similarity index 95% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H index c9c4890b01..fccb98c8cb 100644 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/pEqn.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H @@ -48,9 +48,6 @@ } } - // Explicitly relax pressure for momentum corrector - p.relax(); - Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl; // Second part of thermodynamic density update diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/readAdditionalSolutionControls.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/readAdditionalSolutionControls.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/readAdditionalSolutionControls.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/readAdditionalSolutionControls.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/readChemistryProperties.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/readChemistryProperties.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/readChemistryProperties.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/readChemistryProperties.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/readTimeControls.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H similarity index 97% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/readTimeControls.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H index 0ab2be096e..420eb70c7b 100644 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/readTimeControls.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/readTimeControls.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2010-2010 OpenCFD Ltd. + \\ / A nd | Copyright (C) 2010-2011 OpenCFD Ltd. \\/ M anipulation | ------------------------------------------------------------------------------- License diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/rhoEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/rhoEqn.H similarity index 95% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/rhoEqn.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/rhoEqn.H index e42db72399..f8f1d15ab7 100644 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/rhoEqn.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/rhoEqn.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2008-2010 OpenCFD Ltd. + \\ / A nd | Copyright (C) 2008-2011 OpenCFD Ltd. \\/ M anipulation | ------------------------------------------------------------------------------- License diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/setPressureWork.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/setPressureWork.H similarity index 100% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/setPressureWork.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/setPressureWork.H diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/timeScales.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H similarity index 99% rename from applications/solvers/lagrangian/steadyReactingParcelFoam/timeScales.H rename to applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H index 0e26fcdac1..b44271ab8c 100644 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/timeScales.H +++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/timeScales.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2010-2010 OpenCFD Ltd. + \\ / A nd | Copyright (C) 2010-2011 OpenCFD Ltd. \\/ M anipulation | ------------------------------------------------------------------------------- License diff --git a/applications/solvers/lagrangian/steadyReactingParcelFoam/Make/files b/applications/solvers/lagrangian/steadyReactingParcelFoam/Make/files deleted file mode 100644 index a64e1ee786..0000000000 --- a/applications/solvers/lagrangian/steadyReactingParcelFoam/Make/files +++ /dev/null @@ -1,3 +0,0 @@ -steadyReactingParcelFoam.C - -EXE = $(FOAM_APPBIN)/steadyReactingParcelFoam diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/CH4 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/CH4 similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/CH4 rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/CH4 diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/CO2 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/CO2 similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/CO2 rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/CO2 diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/H2O b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/H2O similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/H2O rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/H2O diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/N2 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/N2 similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/N2 rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/N2 diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/O2 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/O2 similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/O2 rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/O2 diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/T b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/T rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/T diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/U b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/U rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/U diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/alphat b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/alphat similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/alphat rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/alphat diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/epsilon b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/epsilon rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/epsilon diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/k b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/k rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/k diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/mut b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/mut similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/mut rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/mut diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/p b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/p rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/p diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/rho b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/rho similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/0/rho rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/0/rho diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/Allrun b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun similarity index 82% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/Allrun rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun index 4263158119..ce144f896a 100755 --- a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/Allrun +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/Allrun @@ -8,4 +8,4 @@ cd ${0%/*} || exit 1 # run from this directory runApplication blockMesh # run the solver -runApplication steadyReactingParcelFoam +runApplication LTSReactingParcelFoam diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/RASProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/RASProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/RASProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/RASProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/chemistryProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/energySourceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/energySourceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/energySourceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/energySourceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/g b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/g similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/g rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/g diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/massSourceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/massSourceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/massSourceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/massSourceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/momentumSourceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/momentumSourceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/momentumSourceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/momentumSourceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/blockMeshDict b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/blockMeshDict similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/blockMeshDict rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/blockMeshDict diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/polyMesh/boundary diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/radiationProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/radiationProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/radiationProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/radiationProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/reactingCloud1Properties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/reactingCloud1Properties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/reactingCloud1Properties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/reactingCloud1Properties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/reactions b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/reactions similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/reactions rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/reactions diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermo.icoPoly8 diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/thermophysicalProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/constant/turbulenceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/controlDict b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict similarity index 96% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/controlDict rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict index 9cb7a51bef..2d1233ffbe 100644 --- a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/controlDict +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/controlDict @@ -51,4 +51,8 @@ maxCo 0.5; maxDeltaT 1; +alphaTauSpecie 1e10; + +alphaTauTemp 0.005; + // ************************************************************************* // diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/fvSchemes b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/fvSchemes rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSchemes diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution similarity index 98% rename from tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/fvSolution rename to tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution index 14096cac1e..04d104577b 100755 --- a/tutorials/lagrangian/steadyReactingParcelFoam/counterFlowFlame2D/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution @@ -53,8 +53,9 @@ solvers } } -SIMPLE +PISO { + nCorrectors 2; nNonOrthogonalCorrectors 0; momentumPredictor yes; rhoMin rhoMin [1 -3 0 0 0] 0.1; diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/H2O b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/H2O similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/H2O rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/H2O diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/T b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/T similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/T rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/T diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/U b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/U similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/U rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/U diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/air b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/air similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/air rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/air diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/alphat b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/alphat similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/alphat rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/alphat diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/k b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/k similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/k rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/k diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/mut b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/mut similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/mut rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/mut diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/omega b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/omega similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/omega rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/omega diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/p b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/p similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0.org/p rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0.org/p diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/H2O b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/H2O similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/H2O rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/H2O diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/T b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/T similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/T rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/T diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/U b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/U similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/U rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/U diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/air b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/air similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/air rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/air diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/alphat b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/alphat similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/alphat rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/alphat diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/k b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/k similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/k rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/k diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/mut b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/mut similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/mut rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/mut diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/omega b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/omega similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/omega rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/omega diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/p b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/p similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/0/p rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/0/p diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/Allclean b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/Allclean similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/Allclean rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/Allclean diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/Allrun b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/Allrun similarity index 87% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/Allrun rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/Allrun index ab18f7f797..192c30fd62 100755 --- a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/Allrun +++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/Allrun @@ -8,6 +8,6 @@ cd ${0%/*} || exit 1 # run from this directory runApplication blockMesh # run the solver -runApplication steadyReactingParcelFoam +runApplication LTSReactingParcelFoam # ----------------------------------------------------------------- end-of-file diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/RASProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/RASProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/RASProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/RASProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/chemistryProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/chemistryProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/chemistryProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/chemistryProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/energySourceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/energySourceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/energySourceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/energySourceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/g b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/g similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/g rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/g diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/massSourceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/massSourceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/massSourceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/massSourceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/momentumSourceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/momentumSourceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/momentumSourceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/momentumSourceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/particleTrackDict b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/particleTrackDict similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/particleTrackDict rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/particleTrackDict diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/polyMesh/blockMeshDict b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/polyMesh/blockMeshDict similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/polyMesh/blockMeshDict rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/polyMesh/blockMeshDict diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/polyMesh/boundary b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/polyMesh/boundary similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/polyMesh/boundary rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/polyMesh/boundary diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/porousZones b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/porousZones similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/porousZones rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/porousZones diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/radiationProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/radiationProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/radiationProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/radiationProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactingCloud1Properties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/reactions b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactions similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/reactions rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/reactions diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/thermo.incompressiblePoly b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/thermo.incompressiblePoly similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/thermo.incompressiblePoly rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/thermo.incompressiblePoly diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/thermophysicalProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/thermophysicalProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/thermophysicalProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/thermophysicalProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/turbulenceProperties b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/turbulenceProperties similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/constant/turbulenceProperties rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/constant/turbulenceProperties diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/controlDict b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/controlDict similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/controlDict rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/controlDict diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/fvSchemes b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSchemes similarity index 100% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/fvSchemes rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSchemes diff --git a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution similarity index 98% rename from tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/fvSolution rename to tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution index 25043f8651..a78ea6bfef 100644 --- a/tutorials/lagrangian/steadyReactingParcelFoam/verticalChannel/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution @@ -61,8 +61,9 @@ solvers }; } -SIMPLE +PISO { + nCorrectors 2; nNonOrthogonalCorrectors 0; // 10; momentumPredictor yes; rhoMin rhoMin [ 1 -3 0 0 0 ] 0; // 0.5; From b91eb061b2687222a05dd752a59459e08a908b5b Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 11:01:47 +0000 Subject: [PATCH 3/8] ENH: Clouds - only read parcel fields when solution is active --- .../clouds/Templates/CollidingCloud/CollidingCloud.C | 8 ++++---- .../clouds/Templates/KinematicCloud/KinematicCloud.C | 8 ++++---- .../clouds/Templates/ReactingCloud/ReactingCloud.C | 10 +++++----- .../ReactingMultiphaseCloud/ReactingMultiphaseCloud.C | 8 ++++---- .../clouds/Templates/ThermoCloud/ThermoCloud.C | 8 ++++---- 5 files changed, 21 insertions(+), 21 deletions(-) diff --git a/src/lagrangian/intermediate/clouds/Templates/CollidingCloud/CollidingCloud.C b/src/lagrangian/intermediate/clouds/Templates/CollidingCloud/CollidingCloud.C index 925468ad61..10f7f5c1e8 100644 --- a/src/lagrangian/intermediate/clouds/Templates/CollidingCloud/CollidingCloud.C +++ b/src/lagrangian/intermediate/clouds/Templates/CollidingCloud/CollidingCloud.C @@ -114,11 +114,11 @@ Foam::CollidingCloud::CollidingCloud if (this->solution().active()) { setModels(); - } - if (readFields) - { - parcelType::readFields(*this); + if (readFields) + { + parcelType::readFields(*this); + } } } diff --git a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C index 4275cf78f8..0b9b1f7931 100644 --- a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C +++ b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C @@ -360,11 +360,11 @@ Foam::KinematicCloud::KinematicCloud if (solution_.active()) { setModels(); - } - if (readFields) - { - parcelType::readFields(*this); + if (readFields) + { + parcelType::readFields(*this); + } } if (solution_.resetSourcesOnStartup()) diff --git a/src/lagrangian/intermediate/clouds/Templates/ReactingCloud/ReactingCloud.C b/src/lagrangian/intermediate/clouds/Templates/ReactingCloud/ReactingCloud.C index 1c465fc450..eb1f66e85a 100644 --- a/src/lagrangian/intermediate/clouds/Templates/ReactingCloud/ReactingCloud.C +++ b/src/lagrangian/intermediate/clouds/Templates/ReactingCloud/ReactingCloud.C @@ -117,6 +117,11 @@ Foam::ReactingCloud::ReactingCloud if (this->solution().active()) { setModels(); + + if (readFields) + { + parcelType::readFields(*this, this->composition()); + } } // Set storage for mass source fields and initialise to zero @@ -142,11 +147,6 @@ Foam::ReactingCloud::ReactingCloud ); } - if (readFields) - { - parcelType::readFields(*this, this->composition()); - } - if (this->solution().resetSourcesOnStartup()) { resetSourceTerms(); diff --git a/src/lagrangian/intermediate/clouds/Templates/ReactingMultiphaseCloud/ReactingMultiphaseCloud.C b/src/lagrangian/intermediate/clouds/Templates/ReactingMultiphaseCloud/ReactingMultiphaseCloud.C index 89374ecb48..0290eacb7c 100644 --- a/src/lagrangian/intermediate/clouds/Templates/ReactingMultiphaseCloud/ReactingMultiphaseCloud.C +++ b/src/lagrangian/intermediate/clouds/Templates/ReactingMultiphaseCloud/ReactingMultiphaseCloud.C @@ -94,11 +94,11 @@ Foam::ReactingMultiphaseCloud::ReactingMultiphaseCloud if (this->solution().active()) { setModels(); - } - if (readFields) - { - parcelType::readFields(*this, this->composition()); + if (readFields) + { + parcelType::readFields(*this, this->composition()); + } } if (this->solution().resetSourcesOnStartup()) diff --git a/src/lagrangian/intermediate/clouds/Templates/ThermoCloud/ThermoCloud.C b/src/lagrangian/intermediate/clouds/Templates/ThermoCloud/ThermoCloud.C index d2ce89b5c4..2dbe453580 100644 --- a/src/lagrangian/intermediate/clouds/Templates/ThermoCloud/ThermoCloud.C +++ b/src/lagrangian/intermediate/clouds/Templates/ThermoCloud/ThermoCloud.C @@ -135,11 +135,11 @@ Foam::ThermoCloud::ThermoCloud if (this->solution().active()) { setModels(); - } - if (readFields) - { - parcelType::readFields(*this); + if (readFields) + { + parcelType::readFields(*this); + } } if (this->solution().resetSourcesOnStartup()) From ad5a84941660f6dda675b713c6c7cd370de86304 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 11:42:37 +0000 Subject: [PATCH 4/8] ENH: Tutoral updates --- .../counterFlowFlame2D/system/fvSolution | 5 ----- .../LTSReactingParcelFoam/verticalChannel/system/fvSolution | 5 ----- 2 files changed, 10 deletions(-) diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution index 04d104577b..799d7c9376 100755 --- a/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/counterFlowFlame2D/system/fvSolution @@ -69,10 +69,5 @@ additional solveSpecies true; } -relaxationFactors -{ - p 0.3; -} - // ************************************************************************* // diff --git a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution index a78ea6bfef..1d4e95dae0 100644 --- a/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution +++ b/tutorials/lagrangian/LTSReactingParcelFoam/verticalChannel/system/fvSolution @@ -70,11 +70,6 @@ PISO rhoMax rhoMax [ 1 -3 0 0 0 ] 100; // 1.5; } -relaxationFactors -{ - p 0.3; -} - additional { eWork true; From a5ce23af338effdde7a42982aa62e68c110a2ba9 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 14:31:00 +0000 Subject: [PATCH 5/8] ENH: Improvements to particle integration methods --- .../IntegrationScheme/Analytical/Analytical.C | 17 +++++++++++++---- .../IntegrationScheme/Analytical/Analytical.H | 2 +- .../IntegrationScheme/Euler/Euler.C | 4 ++-- .../IntegrationScheme/Euler/Euler.H | 2 +- .../IntegrationScheme/IntegrationScheme.C | 2 +- .../IntegrationScheme/IntegrationScheme.H | 2 +- .../Templates/KinematicParcel/KinematicParcel.C | 7 +++---- .../Templates/ThermoParcel/ThermoParcel.C | 2 +- 8 files changed, 23 insertions(+), 15 deletions(-) diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C index 89c72685a2..d7995dd405 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C @@ -61,7 +61,7 @@ Foam::Analytical::integrate ( const Type& phi, const scalar dt, - const Type& alpha, + const Type& alphaBeta, const scalar beta ) const { @@ -69,9 +69,18 @@ Foam::Analytical::integrate const scalar expTerm = exp(min(50, -beta*dt)); - retValue.average() = - alpha + (phi - alpha)*(1 - expTerm)/(beta*dt + ROOTVSMALL); - retValue.value() = alpha + (phi - alpha)*expTerm; + if (beta > ROOTVSMALL) + { + const Type alpha = alphaBeta/beta; + retValue.average() = alpha + (phi - alpha)*(1 - expTerm)/(beta*dt); + retValue.value() = alpha + (phi - alpha)*expTerm; + } + else + { + retValue.average() = phi; + retValue.value() = phi; + } + return retValue; } diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H index 91a188dc3a..5eca4f17a8 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.H @@ -83,7 +83,7 @@ public: ( const Type& phi, const scalar dt, - const Type& alpha, + const Type& alphaBeta, const scalar beta ) const; }; diff --git a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C index 1ebaa9dc29..ba9cc4c072 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C @@ -61,12 +61,12 @@ Foam::Euler::integrate ( const Type& phi, const scalar dt, - const Type& alpha, + const Type& alphaBeta, const scalar beta ) const { typename IntegrationScheme::integrationResult retValue; - retValue.value() = (phi + beta*dt*alpha)/(1.0 + beta*dt); + retValue.value() = (phi + alphaBeta*dt)/(1.0 + beta*dt); retValue.average() = 0.5*(phi + retValue.value()); return retValue; diff --git a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H index de58579a9d..3737d6e902 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H +++ b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.H @@ -80,7 +80,7 @@ public: ( const Type& phi, const scalar dt, - const Type& alpha, + const Type& alphaBeta, const scalar beta ) const; }; diff --git a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C index d3cac7f4e0..a8ea549fce 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C +++ b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.C @@ -62,7 +62,7 @@ Foam::IntegrationScheme::integrate ( const Type& phi, const scalar dt, - const Type& alpha, + const Type& alphaBeta, const scalar beta ) const { diff --git a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H index d7f0d0ef07..941185dcd9 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H +++ b/src/lagrangian/intermediate/IntegrationScheme/IntegrationScheme/IntegrationScheme.H @@ -185,7 +185,7 @@ public: ( const Type& phi, const scalar dt, - const Type& alpha, + const Type& alphaBeta, const scalar beta ) const; }; diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C index 58390134b3..9a83e69908 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C @@ -192,21 +192,20 @@ const Foam::vector Foam::KinematicParcel::calcVelocity const parcelType& p = static_cast(*this); const forceSuSp Fcp = forces.calcCoupled(p, dt, mass, Re, mu); const forceSuSp Fncp = forces.calcNonCoupled(p, dt, mass, Re, mu); - forceSuSp Feff = Fcp + Fncp; - Feff.Sp() += ROOTVSMALL; + const forceSuSp Feff = Fcp + Fncp; // New particle velocity //~~~~~~~~~~~~~~~~~~~~~~ // Update velocity - treat as 3-D - const vector ap = Uc_ + (Feff.Su() + Su)/Feff.Sp(); + const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/mass; const scalar bp = Feff.Sp()/mass; Spu = Feff.Sp(); IntegrationScheme::integrationResult Ures = - td.cloud().UIntegrator().integrate(U, dt, ap, bp); + td.cloud().UIntegrator().integrate(U, dt, abp, bp); vector Unew = Ures.value(); diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C index 622ef85fcf..064f945ad1 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C @@ -333,7 +333,7 @@ Foam::scalar Foam::ThermoParcel::calcHeatTransfer // Integrate to find the new parcel temperature IntegrationScheme::integrationResult Tres = - td.cloud().TIntegrator().integrate(T, dt, ap, bp); + td.cloud().TIntegrator().integrate(T, dt, ap*bp, bp); scalar Tnew = max(Tres.value(), td.cloud().constProps().TMin()); From ce29cba31a1768495afed354411d23f1c99263f6 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 14:31:27 +0000 Subject: [PATCH 6/8] ENH: use solution dt instead of time dt --- .../clouds/Templates/KinematicCloud/KinematicCloud.C | 2 +- .../InjectionModel/InjectionModel/InjectionModel.C | 2 +- .../PatchInteractionModel/PatchInteractionModel.C | 6 ++++-- 3 files changed, 6 insertions(+), 4 deletions(-) diff --git a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C index 0b9b1f7931..3e49d29cd5 100644 --- a/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C +++ b/src/lagrangian/intermediate/clouds/Templates/KinematicCloud/KinematicCloud.C @@ -504,7 +504,7 @@ void Foam::KinematicCloud::checkParcelProperties parcel.rho() = constProps_.rho0(); } - const scalar carrierDt = this->db().time().deltaTValue(); + const scalar carrierDt = solution_.deltaT(); parcel.stepFraction() = (carrierDt - lagrangianDt)/carrierDt; parcel.typeId() = constProps_.parcelTypeId(); } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C index 94b5a3d91b..5bf1c140c3 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C @@ -520,7 +520,7 @@ void Foam::InjectionModel::inject(TrackData& td) } const scalar time = this->owner().db().time().value(); - const scalar carrierDt = this->owner().db().time().deltaTValue(); + const scalar carrierDt = this->owner().solution().deltaT(); const polyMesh& mesh = this->owner().mesh(); // Prepare for next time step diff --git a/src/lagrangian/intermediate/submodels/Kinematic/PatchInteractionModel/PatchInteractionModel/PatchInteractionModel.C b/src/lagrangian/intermediate/submodels/Kinematic/PatchInteractionModel/PatchInteractionModel/PatchInteractionModel.C index e2e1c67d24..3236774c18 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/PatchInteractionModel/PatchInteractionModel/PatchInteractionModel.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/PatchInteractionModel/PatchInteractionModel/PatchInteractionModel.C @@ -290,7 +290,7 @@ void Foam::PatchInteractionModel::patchData } else { - Up = (Cf - Cf00)/mesh.time().deltaTValue(); + Up = (Cf - Cf00)/this->owner().solution().deltaT(); } if (mag(dn) > SMALL) @@ -312,7 +312,9 @@ void Foam::PatchInteractionModel::patchData // magOmega = sin(angle between unit normals) // Normalise omega vector by magOmega, then multiply by // angle/dt to give the correct angular velocity vector. - omega *= Foam::asin(magOmega)/(magOmega*mesh.time().deltaTValue()); + omega *= + Foam::asin(magOmega) + /(magOmega*this->owner().solution().deltaT()); // Project position onto face and calculate this position // relative to the face centre. From 73ef3c4dffe230939077825ccd4827e698154ab4 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 15:40:58 +0000 Subject: [PATCH 7/8] ENH: #include clean-up --- .../intermediate/IntegrationScheme/Analytical/Analytical.C | 1 - src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C | 1 - 2 files changed, 2 deletions(-) diff --git a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C index d7995dd405..be7539c6f3 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Analytical/Analytical.C @@ -24,7 +24,6 @@ License \*---------------------------------------------------------------------------*/ #include "Analytical.H" -#include "addToRunTimeSelectionTable.H" // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // diff --git a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C index ba9cc4c072..0ba07fb2af 100644 --- a/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C +++ b/src/lagrangian/intermediate/IntegrationScheme/Euler/Euler.C @@ -24,7 +24,6 @@ License \*---------------------------------------------------------------------------*/ #include "Euler.H" -#include "addToRunTimeSelectionTable.H" // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // From 432f72a88a0810d68fbf9eee74719dcaa721cd31 Mon Sep 17 00:00:00 2001 From: andy Date: Wed, 2 Mar 2011 15:41:24 +0000 Subject: [PATCH 8/8] ENH: upadtes to semi-implicit cloud source calcs --- .../parcels/Templates/KinematicParcel/KinematicParcel.C | 2 +- .../intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C index 9a83e69908..0a3e05d200 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C @@ -202,7 +202,7 @@ const Foam::vector Foam::KinematicParcel::calcVelocity const vector abp = (Feff.Sp()*Uc_ + (Feff.Su() + Su))/mass; const scalar bp = Feff.Sp()/mass; - Spu = Feff.Sp(); + Spu = Feff.Sp()*dt/td.cloud().solution().deltaT(); IntegrationScheme::integrationResult Ures = td.cloud().UIntegrator().integrate(U, dt, abp, bp); diff --git a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C index 064f945ad1..9df054c2d0 100644 --- a/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/ThermoParcel/ThermoParcel.C @@ -339,7 +339,7 @@ Foam::scalar Foam::ThermoParcel::calcHeatTransfer dhsTrans += dt*htc*As*(0.5*(T + Tnew) - Tc_); - Cuh = bp; + Cuh = bp*dt/td.cloud().solution().deltaT(); return Tnew; }