diff --git a/src/OSspecific/POSIX/POSIX.C b/src/OSspecific/POSIX/POSIX.C
index 48e4b3a4dfe05fb0dca7ddc19053b3aaff0e075e..dba1e73f3a536110fe98f6aee6822ea9f36f475e 100644
--- a/src/OSspecific/POSIX/POSIX.C
+++ b/src/OSspecific/POSIX/POSIX.C
@@ -57,6 +57,10 @@ Description
#include <netinet/in.h>
+#ifdef USE_RANDOM
+# include <climits>
+#endif
+
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
defineTypeNameAndDebug(Foam::POSIX, 0);
@@ -68,16 +72,19 @@ pid_t Foam::pid()
return ::getpid();
}
+
pid_t Foam::ppid()
{
return ::getppid();
}
+
pid_t Foam::pgid()
{
return ::getpgrp();
}
+
bool Foam::env(const word& envName)
{
return ::getenv(envName.c_str()) != NULL;
@@ -890,7 +897,6 @@ bool Foam::mvBak(const fileName& src, const std::string& ext)
}
-
// Remove a file, returning true if successful otherwise false
bool Foam::rm(const fileName& file)
{
@@ -1221,4 +1227,34 @@ Foam::fileNameList Foam::dlLoaded()
}
+void Foam::osRandomSeed(const label seed)
+{
+#ifdef USE_RANDOM
+ srandom((unsigned int)seed);
+#else
+ srand48(seed);
+#endif
+}
+
+
+Foam::label Foam::osRandomInteger()
+{
+#ifdef USE_RANDOM
+ return random();
+#else
+ return lrand48();
+#endif
+}
+
+
+Foam::scalar Foam::osRandomDouble()
+{
+#ifdef USE_RANDOM
+ return (scalar)random();
+#else
+ return drand48();
+#endif
+}
+
+
// ************************************************************************* //
diff --git a/src/OpenFOAM/include/OSspecific.H b/src/OpenFOAM/include/OSspecific.H
index f01ad296476fa90f4d4ae67d34be44c8b122c1f5..b77499733974b480a3b6277cd484e5a5756728dc 100644
--- a/src/OpenFOAM/include/OSspecific.H
+++ b/src/OpenFOAM/include/OSspecific.H
@@ -200,6 +200,18 @@ bool dlSymFound(void* handle, const std::string& symbol);
fileNameList dlLoaded();
+// Low level random numbers. Use Random class instead.
+
+//- Seed random number generator.
+void osRandomSeed(const label seed);
+
+//- Return random integer (uniform distribution between 0 and 2^31)
+label osRandomInteger();
+
+//- Return random double precision (uniform distribution between 0 and 1)
+scalar osRandomDouble();
+
+
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
diff --git a/src/OpenFOAM/primitives/random/Random/Random.C b/src/OpenFOAM/primitives/random/Random/Random.C
index 704e5ef118119c71d9555c44a0d6526cf53fda7e..e2496ee648625f6bf45d4833ae1298f558143fd4 100644
--- a/src/OpenFOAM/primitives/random/Random/Random.C
+++ b/src/OpenFOAM/primitives/random/Random/Random.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2004-2010 OpenCFD Ltd.
+ \\ / A nd | Copyright (C) 2004-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -24,6 +24,7 @@ License
\*---------------------------------------------------------------------------*/
#include "Random.H"
+#include "OSspecific.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -37,11 +38,6 @@ namespace Foam
# error "The random number generator may not work!"
#endif
-#ifdef USE_RANDOM
-# include <climits>
-#else
-# include <cstdlib>
-#endif
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
@@ -57,22 +53,13 @@ Random::Random(const label seed)
Seed = 1;
}
-# ifdef USE_RANDOM
- srandom((unsigned int)Seed);
-# else
- srand48(Seed);
-# endif
-
+ osRandomSeed(Seed);
}
int Random::bit()
{
-# ifdef USE_RANDOM
- if (random() > INT_MAX/2)
-# else
- if (lrand48() > INT_MAX/2)
-# endif
+ if (osRandomInteger() > INT_MAX/2)
{
return 1;
}
@@ -85,11 +72,7 @@ int Random::bit()
scalar Random::scalar01()
{
-# ifdef USE_RANDOM
- return (scalar)random()/INT_MAX;
-# else
- return drand48();
-# endif
+ return osRandomDouble();
}
@@ -140,11 +123,7 @@ tensor Random::tensor01()
label Random::integer(const label lower, const label upper)
{
-# ifdef USE_RANDOM
- return lower + (random() % (upper+1-lower));
-# else
- return lower + (lrand48() % (upper+1-lower));
-# endif
+ return lower + (osRandomInteger() % (upper+1-lower));
}
diff --git a/src/OpenFOAM/primitives/random/cachedRandom/cachedRandom.C b/src/OpenFOAM/primitives/random/cachedRandom/cachedRandom.C
index b48c1f4e63269f892bde020164f005367a03f740..ac0c2705030ba3ef67e2aa9ae6decad8517a71ae 100644
--- a/src/OpenFOAM/primitives/random/cachedRandom/cachedRandom.C
+++ b/src/OpenFOAM/primitives/random/cachedRandom/cachedRandom.C
@@ -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
@@ -24,7 +24,7 @@ License
\*---------------------------------------------------------------------------*/
#include "cachedRandom.H"
-#include <cstdlib>
+#include "OSspecific.H"
#if INT_MAX != 2147483647
# error "INT_MAX != 2147483647"
@@ -37,7 +37,7 @@ Foam::scalar Foam::cachedRandom::scalar01()
{
if (sampleI_ < 0)
{
- return drand48();
+ return osRandomDouble();
}
if (sampleI_ == samples_.size() - 1)
@@ -76,7 +76,7 @@ Foam::cachedRandom::cachedRandom(const label seed, const label count)
}
// Initialise samples
- srand48(seed_);
+ osRandomSeed(seed_);
forAll(samples_, i)
{
samples_[i] = drand48();
@@ -98,7 +98,7 @@ Foam::cachedRandom::cachedRandom(const cachedRandom& cr, const bool reset)
) << "Copy constructor called, but samples not being cached. "
<< "This may lead to non-repeatable behaviour" << endl;
- srand48(seed_);
+ osRandomSeed(seed_);
}
else if (reset)
{
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairModel/PairSpringSliderDashpot/PairSpringSliderDashpot.C b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairModel/PairSpringSliderDashpot/PairSpringSliderDashpot.C
index 812296e6f5b1bb161dc732388d2d105ae56a6323..5580d3cb05a195e63d92a4ad627e6c566d6f3295 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairModel/PairSpringSliderDashpot/PairSpringSliderDashpot.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairModel/PairSpringSliderDashpot/PairSpringSliderDashpot.C
@@ -214,7 +214,8 @@ void Foam::PairSpringSliderDashpot<CloudType>::evaluatePair
rHat_AB
*(kN*pow(normalOverlapMag, b_) - etaN*(U_AB & rHat_AB));
- // Cohesion force
+ // Cohesion force, energy density multiplied by the area of
+ // particle-particle overlap
if (cohesion_)
{
fN_AB +=
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.C b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.C
index 9ae09dd98f639553aa51fc1cdf98f13c7896094e..946ec3adb89ad85c17cad24ecf842bb46f8717b3 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.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) 2008-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -76,7 +76,8 @@ void Foam::WallLocalSpringSliderDashpot<CloudType>::evaluateWall
typename CloudType::parcelType& p,
const point& site,
const WallSiteData<vector>& data,
- scalar pREff
+ scalar pREff,
+ bool cohesion
) const
{
// wall patch index
@@ -88,14 +89,18 @@ void Foam::WallLocalSpringSliderDashpot<CloudType>::evaluateWall
scalar alpha = alpha_[wPI];
scalar b = b_[wPI];
scalar mu = mu_[wPI];
+ scalar cohesionEnergyDensity = cohesionEnergyDensity_[wPI];
+ cohesion = cohesion && cohesion_[wPI];
vector r_PW = p.position() - site;
vector U_PW = p.U() - data.wallData();
- scalar normalOverlapMag = max(pREff - mag(r_PW), 0.0);
+ scalar r_PW_mag = mag(r_PW);
- vector rHat_PW = r_PW/(mag(r_PW) + VSMALL);
+ scalar normalOverlapMag = max(pREff - r_PW_mag, 0.0);
+
+ vector rHat_PW = r_PW/(r_PW_mag + VSMALL);
scalar kN = (4.0/3.0)*sqrt(pREff)*Estar;
@@ -105,6 +110,16 @@ void Foam::WallLocalSpringSliderDashpot<CloudType>::evaluateWall
rHat_PW
*(kN*pow(normalOverlapMag, b) - etaN*(U_PW & rHat_PW));
+ // Cohesion force, energy density multiplied by the area of wall/particle
+ // overlap
+ if (cohesion)
+ {
+ fN_PW +=
+ -cohesionEnergyDensity
+ *mathematical::pi*(sqr(pREff) - sqr(r_PW_mag))
+ *rHat_PW;
+ }
+
p.f() += fN_PW;
vector USlip_PW =
@@ -168,6 +183,8 @@ Foam::WallLocalSpringSliderDashpot<CloudType>::WallLocalSpringSliderDashpot
alpha_(),
b_(),
mu_(),
+ cohesionEnergyDensity_(),
+ cohesion_(),
patchMap_(),
maxEstarIndex_(-1),
collisionResolutionSteps_
@@ -212,6 +229,8 @@ Foam::WallLocalSpringSliderDashpot<CloudType>::WallLocalSpringSliderDashpot
alpha_.setSize(nWallPatches);
b_.setSize(nWallPatches);
mu_.setSize(nWallPatches);
+ cohesionEnergyDensity_.setSize(nWallPatches);
+ cohesion_.setSize(nWallPatches);
scalar maxEstar = -GREAT;
@@ -238,6 +257,13 @@ Foam::WallLocalSpringSliderDashpot<CloudType>::WallLocalSpringSliderDashpot
mu_[wPI] = readScalar(patchCoeffDict.lookup("mu"));
+ cohesionEnergyDensity_[wPI] = readScalar
+ (
+ patchCoeffDict.lookup("cohesionEnergyDensity")
+ );
+
+ cohesion_[wPI] = (mag(cohesionEnergyDensity_[wPI]) > VSMALL);
+
if (Estar_[wPI] > maxEstar)
{
maxEstarIndex_ = wPI;
@@ -325,20 +351,22 @@ void Foam::WallLocalSpringSliderDashpot<CloudType>::evaluateWall
p,
flatSitePoints[siteI],
flatSiteData[siteI],
- pREff
+ pREff,
+ true
);
}
forAll(sharpSitePoints, siteI)
{
- // Treating sharp sites like flat sites
+ // Treating sharp sites like flat sites, except suppress cohesion
evaluateWall
(
p,
sharpSitePoints[siteI],
sharpSiteData[siteI],
- pREff
+ pREff,
+ false
);
}
}
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.H b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.H
index 700ce321abdfac0f021a1da37e91228a83f52773..c4ad801b13afaa239a0a3baef207b545b899bf6a 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.H
+++ b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallLocalSpringSliderDashpot/WallLocalSpringSliderDashpot.H
@@ -65,6 +65,12 @@ class WallLocalSpringSliderDashpot
//- Coefficient of friction in for tangential sliding
scalarList mu_;
+ //- Cohesion energy density [J/m^3]
+ scalarList cohesionEnergyDensity_;
+
+ // Switch cohesion on and off
+ boolList cohesion_;
+
//- Mapping the patch index to the model data
labelList patchMap_;
@@ -115,7 +121,8 @@ class WallLocalSpringSliderDashpot
typename CloudType::parcelType& p,
const point& site,
const WallSiteData<vector>& data,
- scalar pREff
+ scalar pREff,
+ bool cohesion
) const;
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.C b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.C
index 5f14b8e274fec2dfd8dccd4b8b0702fa10b48adc..1f547f0d5d760e4d34912a64c4483dbe1dca37c2 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.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) 2008-2011 OpenCFD Ltd.
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -77,16 +77,19 @@ void Foam::WallSpringSliderDashpot<CloudType>::evaluateWall
const point& site,
const WallSiteData<vector>& data,
scalar pREff,
- scalar kN
+ scalar kN,
+ bool cohesion
) const
{
vector r_PW = p.position() - site;
vector U_PW = p.U() - data.wallData();
- scalar normalOverlapMag = max(pREff - mag(r_PW), 0.0);
+ scalar r_PW_mag = mag(r_PW);
- vector rHat_PW = r_PW/(mag(r_PW) + VSMALL);
+ scalar normalOverlapMag = max(pREff - r_PW_mag, 0.0);
+
+ vector rHat_PW = r_PW/(r_PW_mag + VSMALL);
scalar etaN = alpha_*sqrt(p.mass()*kN)*pow025(normalOverlapMag);
@@ -94,6 +97,16 @@ void Foam::WallSpringSliderDashpot<CloudType>::evaluateWall
rHat_PW
*(kN*pow(normalOverlapMag, b_) - etaN*(U_PW & rHat_PW));
+ // Cohesion force, energy density multiplied by the area of wall/particle
+ // overlap
+ if (cohesion)
+ {
+ fN_PW +=
+ -cohesionEnergyDensity_
+ *mathematical::pi*(sqr(pREff) - sqr(r_PW_mag))
+ *rHat_PW;
+ }
+
p.f() += fN_PW;
vector USlip_PW =
@@ -157,6 +170,11 @@ Foam::WallSpringSliderDashpot<CloudType>::WallSpringSliderDashpot
alpha_(readScalar(this->coeffDict().lookup("alpha"))),
b_(readScalar(this->coeffDict().lookup("b"))),
mu_(readScalar(this->coeffDict().lookup("mu"))),
+ cohesionEnergyDensity_
+ (
+ readScalar(this->coeffDict().lookup("cohesionEnergyDensity"))
+ ),
+ cohesion_(false),
collisionResolutionSteps_
(
readScalar
@@ -183,6 +201,8 @@ Foam::WallSpringSliderDashpot<CloudType>::WallSpringSliderDashpot
Estar_ = 1/((1 - sqr(pNu))/pE + (1 - sqr(nu))/E);
Gstar_ = 1/(2*((2 + pNu - sqr(pNu))/pE + (2 + nu - sqr(nu))/E));
+
+ cohesion_ = (mag(cohesionEnergyDensity_) > VSMALL);
}
@@ -266,13 +286,14 @@ void Foam::WallSpringSliderDashpot<CloudType>::evaluateWall
flatSitePoints[siteI],
flatSiteData[siteI],
pREff,
- kN
+ kN,
+ cohesion_
);
}
forAll(sharpSitePoints, siteI)
{
- // Treating sharp sites like flat sites
+ // Treating sharp sites like flat sites, except suppress cohesion
evaluateWall
(
@@ -280,7 +301,8 @@ void Foam::WallSpringSliderDashpot<CloudType>::evaluateWall
sharpSitePoints[siteI],
sharpSiteData[siteI],
pREff,
- kN
+ kN,
+ false
);
}
}
diff --git a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.H b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.H
index 10010c60c31ce5238720894e82ef9c42424ee22c..e6e93ddb24623d89c5c843b726c433a7398b5ed3 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.H
+++ b/src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/WallModel/WallSpringSliderDashpot/WallSpringSliderDashpot.H
@@ -65,6 +65,12 @@ class WallSpringSliderDashpot
//- Coefficient of friction in for tangential sliding
scalar mu_;
+ //- Cohesion energy density [J/m^3]
+ scalar cohesionEnergyDensity_;
+
+ // Switch cohesion on and off
+ bool cohesion_;
+
//- The number of steps over which to resolve the minimum
// harmonic approximation of the collision period
scalar collisionResolutionSteps_;
@@ -110,7 +116,8 @@ class WallSpringSliderDashpot
const point& site,
const WallSiteData<vector>& data,
scalar pREff,
- scalar kN
+ scalar kN,
+ bool cohesion
) const;
diff --git a/src/postProcessing/functionObjects/IO/partialWrite/partialWrite.C b/src/postProcessing/functionObjects/IO/partialWrite/partialWrite.C
index 327c7b804ad0d93fc54b8f600780cd00d5475126..1ea7ae253d131cc38176cbc717f9f9e741ab5272 100644
--- a/src/postProcessing/functionObjects/IO/partialWrite/partialWrite.C
+++ b/src/postProcessing/functionObjects/IO/partialWrite/partialWrite.C
@@ -27,6 +27,7 @@ License
#include "dictionary.H"
#include "Time.H"
#include "IOobjectList.H"
+#include "polyMesh.H"
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@@ -123,14 +124,26 @@ void Foam::partialWrite::write()
else
{
// Delete all but marked objects
+ fileName dbDir;
+ if (isA<polyMesh>(obr_))
+ {
+ dbDir = dynamic_cast<const polyMesh&>(obr_).dbDir();
+ }
+
IOobjectList objects(obr_, obr_.time().timeName());
forAllConstIter(HashPtrTable<IOobject>, objects, iter)
{
if (!objectNames_.found(iter()->name()))
{
- const fileName f = obr_.time().timePath()/iter()->name();
- //Pout<< " rm " << f << endl;
+ const fileName f =
+ obr_.time().timePath()
+ /dbDir
+ /iter()->name();
+ if (debug)
+ {
+ Pout<< " rm " << f << endl;
+ }
rm(f);
}
}
diff --git a/tutorials/combustion/fireFoam/les/oppositeBurningPanels/Allrun b/tutorials/combustion/fireFoam/les/oppositeBurningPanels/Allrun
index ae53653faa373895e5c55c08070127365a133e34..026f614b53623028cfabc159623d6e6392fa55f4 100755
--- a/tutorials/combustion/fireFoam/les/oppositeBurningPanels/Allrun
+++ b/tutorials/combustion/fireFoam/les/oppositeBurningPanels/Allrun
@@ -32,4 +32,8 @@ decomposePar -region panelRegion > log.decomposParPanelRegion.log 2>&1
runParallel `getApplication` 6
+paraFoam -touch
+paraFoam -touch -region panelRegion
+
+
# ----------------------------------------------------------------- end-of-file
diff --git a/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperEmptying/constant/kinematicCloudProperties b/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperEmptying/constant/kinematicCloudProperties
index 7422e40c4404239739ae9f0e11a19ade0b78adca..ccd5d89ef86aee30868917ca107055d81d7623b1 100644
--- a/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperEmptying/constant/kinematicCloudProperties
+++ b/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperEmptying/constant/kinematicCloudProperties
@@ -112,6 +112,7 @@ subModels
alpha 0.12;
b 1.5;
mu 0.43;
+ cohesionEnergyDensity 0;
}
frontAndBack
{
@@ -120,6 +121,7 @@ subModels
alpha 0.12;
b 1.5;
mu 0.1;
+ cohesionEnergyDensity 0;
}
};
}
diff --git a/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperInitialState/constant/kinematicCloudProperties b/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperInitialState/constant/kinematicCloudProperties
index 36448f7b3072cd013d7d88f700374341ec2868ac..3adb6ab77122dd63a356966d9d22de6d7eb0f55e 100644
--- a/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperInitialState/constant/kinematicCloudProperties
+++ b/tutorials/lagrangian/icoUncoupledKinematicParcelFoam/hopper/hopperInitialState/constant/kinematicCloudProperties
@@ -121,6 +121,7 @@ subModels
alpha 0.12;
b 1.5;
mu 0.43;
+ cohesionEnergyDensity 0;
}
frontAndBack
{
@@ -129,6 +130,7 @@ subModels
alpha 0.12;
b 1.5;
mu 0.1;
+ cohesionEnergyDensity 0;
}
};
}
diff --git a/tutorials/lagrangian/reactingParcelFilmFoam/cylinder/Allrun.pre b/tutorials/lagrangian/reactingParcelFilmFoam/cylinder/Allrun.pre
index 747433b184cd7c4e4f4813cbc0ff1f43757ea9d2..bd6fc4f90f93c116105bd5bd534fcf3a42212481 100755
--- a/tutorials/lagrangian/reactingParcelFilmFoam/cylinder/Allrun.pre
+++ b/tutorials/lagrangian/reactingParcelFilmFoam/cylinder/Allrun.pre
@@ -8,3 +8,6 @@ runApplication setSet -batch wallFilmRegion.setSet
mv log.setSet log.wallFilmRegion.setSet
runApplication extrudeToRegionMesh -overwrite
+
+paraFoam -touch
+paraFoam -touch -region wallFilmRegion
diff --git a/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/Allrun.pre b/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/Allrun.pre
index 55834022011a73d1e1e9d86dbc86df1e234db9e7..b5423180b6ab7fd33db08e26a3e36d8ff2145b5f 100755
--- a/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/Allrun.pre
+++ b/tutorials/lagrangian/reactingParcelFilmFoam/hotBoxes/Allrun.pre
@@ -32,3 +32,6 @@ cp -r system/wallFilmRegion.org system/wallFilmRegion
find ./0 -maxdepth 1 -type f -exec \
sed -i "s/wallFilm/\"(region0_to.*)\"/g" {} \;
+
+paraFoam -touch
+paraFoam -touch -region wallFilmRegion
diff --git a/tutorials/lagrangian/reactingParcelFilmFoam/rivuletPanel/Allrun.pre b/tutorials/lagrangian/reactingParcelFilmFoam/rivuletPanel/Allrun.pre
index 0e3c1209c5bc981b7e0be9cfaf2192a70b2b46dd..d9b17ad0979f5a80eabeaabbb0d3cdbc402d3c88 100755
--- a/tutorials/lagrangian/reactingParcelFilmFoam/rivuletPanel/Allrun.pre
+++ b/tutorials/lagrangian/reactingParcelFilmFoam/rivuletPanel/Allrun.pre
@@ -10,3 +10,5 @@ mv log.setSet log.wallFilmRegion.setSet
runApplication extrudeToRegionMesh -overwrite
+paraFoam -touch
+paraFoam -touch -region wallFilmRegion
diff --git a/tutorials/lagrangian/reactingParcelFilmFoam/splashPanel/Allrun.pre b/tutorials/lagrangian/reactingParcelFilmFoam/splashPanel/Allrun.pre
index 14da52453ed735ef3b94371277f95c0421540811..bbae3db32c4fa9e323206b4901aa7a27c0526603 100755
--- a/tutorials/lagrangian/reactingParcelFilmFoam/splashPanel/Allrun.pre
+++ b/tutorials/lagrangian/reactingParcelFilmFoam/splashPanel/Allrun.pre
@@ -14,3 +14,6 @@ runApplication setSet -region wallFilmRegion -batch createWallFilmRegionPatches.
mv log.setSet log.createWallFilmRegionPatches.setSet
runApplication createPatch -region wallFilmRegion -overwrite
+
+paraFoam -touch
+paraFoam -touch -region wallFilmRegion