zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

removing old averaging - replaced by averaging function object

Browse Source
This commit is contained in:
andy
2008-06-26 18:01:47 +01:00
committed by graham
parent a2ba29bcb4
commit 5de6dec9a2
7 changed files with 123 additions and 172 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
applications/solvers/multiphase/twoPhaseEulerFoam/averaging/calculateAverages.H
View File

@ -1,10 +0,0 @@
nAveragingSteps++;
scalar nm1Coeff = scalar(nAveragingSteps - 1)/nAveragingSteps;
scalar nCoeff = 1.0/nAveragingSteps;
UaMean = nm1Coeff*UaMean + nCoeff*Ua;
UbMean = nm1Coeff*UbMean + nCoeff*Ub;
alphaMean = nm1Coeff*alphaMean + nCoeff*alpha;
pMean = nm1Coeff*pMean + nCoeff*p;

138
applications/solvers/multiphase/twoPhaseEulerFoam/averaging/createAverages.H
View File

@ -1,138 +0,0 @@
label nAveragingSteps = 1;
IFstream nAveragingStepsFile
(
runTime.path()/runTime.timeName()/"nAveragingSteps.raw"
);
volVectorField* UaMeanPtr;
volVectorField* UbMeanPtr;
volScalarField* alphaMeanPtr;
volScalarField* pMeanPtr;
if (nAveragingStepsFile.good())
{
nAveragingStepsFile >> nAveragingSteps;
Info<< "Reading field UaMean\n" << endl;
UaMeanPtr = new volVectorField
(
IOobject
(
"UaMean",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Reading field UbMean\n" << endl;
UbMeanPtr = new volVectorField
(
IOobject
(
"UbMean",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Reading field alphaMean\n" << endl;
alphaMeanPtr = new volScalarField
(
IOobject
(
"alphaMean",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Reading field pMean\n" << endl;
pMeanPtr = new volScalarField
(
IOobject
(
"pMean",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
}
else
{
Info<< "Creating field UaMean\n" << endl;
UaMeanPtr = new volVectorField
(
IOobject
(
"UaMean",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua
);
Info<< "Creating field UbMean\n" << endl;
UbMeanPtr = new volVectorField
(
IOobject
(
"UbMean",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ub
);
Info<< "Creating field alphaMean\n" << endl;
alphaMeanPtr = new volScalarField
(
IOobject
(
"alphaMean",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
alpha
);
Info<< "Creating field pMean\n" << endl;
pMeanPtr = new volScalarField
(
IOobject
(
"pMean",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
p
);
}
volVectorField& UaMean = *UaMeanPtr;
volVectorField& UbMean = *UbMeanPtr;
volScalarField& alphaMean = *alphaMeanPtr;
volScalarField& pMean = *pMeanPtr;

21
applications/solvers/multiphase/twoPhaseEulerFoam/averaging/writeNaveragingSteps.H
View File

@ -1,21 +0,0 @@
if (runTime.outputTime())
{
OFstream nAveragingStepsFile
(
runTime.path()/runTime.timeName()/"nAveragingSteps.raw"
);
if (nAveragingStepsFile.good())
{
nAveragingStepsFile << nAveragingSteps << endl;
}
else
{
FatalErrorIn(args.executable())
<< "Cannot open file "
<< runTime.path()/runTime.timeName()/"nAveragingSteps.raw"
<< abort(FatalError);
}
mkDir(runTime.path()/runTime.timeName()/"data");
}

3
applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
View File

@ -54,7 +54,6 @@ int main(int argc, char *argv[])
# include "createMesh.H"
# include "readEnvironmentalProperties.H"
# include "createFields.H"
# include "createAverages.H"
# include "readPPProperties.H"
# include "initContinuityErrs.H"
# include "readTimeControls.H"
@ -100,9 +99,7 @@ int main(int argc, char *argv[])
kineticTheory.solve();
nuEffa += kineticTheory.mua()/rhoa;
}
# include "calculateAverages.H"
# include "write.H"
# include "writeNaveragingSteps.H"
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"

41
tutorials/twoPhaseEulerFoam/bed/system/controlDict
View File

@ -50,4 +50,45 @@ maxCo 0.5;
maxDeltaT 1;
functions
(
fieldAverage1
{
// Type of functionObject
type fieldAverage;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldAverage.so");
// Fields to be averaged - runTime modifiable
fields
(
Ua
{
mean on;
prime2Mean off;
base time;
}
Ub
{
mean on;
prime2Mean off;
base time;
}
alpha
{
mean on;
prime2Mean off;
base time;
}
p
{
mean on;
prime2Mean off;
base time;
}
);
}
);
// ************************************************************************* //

41
tutorials/twoPhaseEulerFoam/bed2/system/controlDict
View File

@ -50,4 +50,45 @@ maxCo 0.9;
maxDeltaT 1e-5;
functions
(
fieldAverage1
{
// Type of functionObject
type fieldAverage;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldAverage.so");
// Fields to be averaged - runTime modifiable
fields
(
Ua
{
mean on;
prime2Mean off;
base time;
}
Ub
{
mean on;
prime2Mean off;
base time;
}
alpha
{
mean on;
prime2Mean off;
base time;
}
p
{
mean on;
prime2Mean off;
base time;
}
);
}
);
// ************************************************************************* //

41
tutorials/twoPhaseEulerFoam/bubbleColumn/system/controlDict
View File

@ -50,4 +50,45 @@ maxCo 0.5;
maxDeltaT 1;
functions
(
fieldAverage1
{
// Type of functionObject
type fieldAverage;
// Where to load it from (if not already in solver)
functionObjectLibs ("libfieldAverage.so");
// Fields to be averaged - runTime modifiable
fields
(
Ua
{
mean on;
prime2Mean off;
base time;
}
Ub
{
mean on;
prime2Mean off;
base time;
}
alpha
{
mean on;
prime2Mean off;
base time;
}
p
{
mean on;
prime2Mean off;
base time;
}
);
}
);
// ************************************************************************* //