updated

tutorials/Lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/Make/files

+rhoPisoTwinParcelFoam.C
+
+EXE = $(FOAM_USER_APPBIN)/rhoPisoTwinParcelFoam

tutorials/Lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/createClouds.H

+    Info<< "Constructing thermoCloud1" << endl;
+    basicThermoCloud thermoCloud1
+    (
+        "thermoCloud1",
+        rho,
+        U,
+        g,
+        thermo()
+    );
+
+    Info<< "Constructing kinematicCloud1" << endl;
+    basicKinematicCloud kinematicCloud1
+    (
+        "kinematicCloud1",
+        rho,
+        U,
+        thermo().mu(),
+        g
+    );
+

tutorials/Lagrangian/rhoTurbTwinParcelFoam/rhoTurbTwinParcelFoam/createFields.H → tutorials/Lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/createFields.H

@@ -43,30 +43,16 @@
     Info<< "Creating turbulence model\n" << endl;
     autoPtr<compressible::turbulenceModel> turbulence
     (
-        compressible::turbulenceModel::New(rho, U, phi, thermo())
+        compressible::turbulenceModel::New
+        (
+            rho,
+            U,
+            phi,
+            thermo()
+        )
     );
 
+
     Info<< "Creating field DpDt\n" << endl;
     volScalarField DpDt =
         fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);
-
-    Info<< "Constructing thermoCloud1" << endl;
-    basicThermoCloud thermoCloud1
-    (
-        "thermoCloud1",
-        rho,
-        U,
-        g,
-        thermo()
-    );
-
-    Info<< "Constructing kinematicCloud1" << endl;
-    basicKinematicCloud kinematicCloud1
-    (
-        "kinematicCloud1",
-        rho,
-        U,
-        thermo().mu(),
-        g
-    );
-

tutorials/Lagrangian/rhoTurbTwinParcelFoam/rhoTurbTwinParcelFoam/rhoTurbTwinParcelFoam.C → tutorials/Lagrangian/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam/rhoPisoTwinParcelFoam.C

@@ -23,7 +23,7 @@ License
     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
 
 Application
-    rhoTurbFoam
+    rhoPisoTwinParcelFoam
 
 Description
     Transient solver for compressible, turbulent flow with two thermo-clouds.
@@ -48,6 +48,7 @@ int main(int argc, char *argv[])
 #   include "createMesh.H"
 #   include "readEnvironmentalProperties.H"
 #   include "createFields.H"
+#   include "createClouds.H"
 #   include "readPISOControls.H"
 #   include "initContinuityErrs.H"
 #   include "readTimeControls.H"