mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
Merge branch 'master' of ssh://opencfd:8007/home/dm4/OpenFOAM/repositories/OpenFOAM-dev
This commit is contained in:
Henry
@ -1,6 +1,4 @@
|
||||
{
|
||||
rho = thermo.rho();
|
||||
|
||||
// Thermodynamic density needs to be updated by psi*d(p) after the
|
||||
// pressure solution - done in 2 parts. Part 1:
|
||||
thermo.rho() -= psi*p;
|
||||
@ -54,6 +52,7 @@
|
||||
rho = thermo.rho();
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
rho.relax();
|
||||
|
||||
Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
|
||||
}
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -108,9 +108,27 @@ Foam::scalar Foam::layerAdditionRemoval::readOldThickness
|
||||
|
||||
void Foam::layerAdditionRemoval::clearAddressing() const
|
||||
{
|
||||
// Layer removal data
|
||||
deleteDemandDrivenData(pointsPairingPtr_);
|
||||
deleteDemandDrivenData(facesPairingPtr_);
|
||||
if (pointsPairingPtr_)
|
||||
{
|
||||
if (debug)
|
||||
{
|
||||
Pout<< "layerAdditionRemoval::clearAddressing()" << nl
|
||||
<< " clearing pointsPairingPtr_" << endl;
|
||||
}
|
||||
|
||||
deleteDemandDrivenData(pointsPairingPtr_);
|
||||
}
|
||||
|
||||
if (facesPairingPtr_)
|
||||
{
|
||||
if (debug)
|
||||
{
|
||||
Pout<< "layerAdditionRemoval::clearAddressing()" << nl
|
||||
<< " clearing facesPairingPtr_" << endl;
|
||||
}
|
||||
|
||||
deleteDemandDrivenData(facesPairingPtr_);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2012-2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -64,6 +64,7 @@ SourceFiles
|
||||
|
||||
#include "surfaceInterpolationScheme.H"
|
||||
#include "blendedSchemeBase.H"
|
||||
#include "surfaceInterpolate.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
@ -181,11 +182,6 @@ public:
|
||||
}
|
||||
|
||||
|
||||
//- Destructor
|
||||
virtual ~CoBlended()
|
||||
{}
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
//- Return the face-based blending factor
|
||||
@ -195,39 +191,48 @@ public:
|
||||
) const
|
||||
{
|
||||
const fvMesh& mesh = this->mesh();
|
||||
tmp<surfaceScalarField> tUflux = faceFlux_;
|
||||
|
||||
tmp<surfaceScalarField> tbf
|
||||
if (faceFlux_.dimensions() == dimDensity*dimVelocity*dimArea)
|
||||
{
|
||||
// Currently assume that the density field
|
||||
// corresponding to the mass-flux is named "rho"
|
||||
const volScalarField& rho =
|
||||
mesh.objectRegistry::template lookupObject<volScalarField>
|
||||
("rho");
|
||||
|
||||
tUflux = faceFlux_/fvc::interpolate(rho);
|
||||
}
|
||||
else if (faceFlux_.dimensions() != dimVelocity*dimArea)
|
||||
{
|
||||
FatalErrorIn
|
||||
(
|
||||
"CoBlended::blendingFactor()"
|
||||
) << "dimensions of faceFlux are not correct"
|
||||
<< exit(FatalError);
|
||||
}
|
||||
|
||||
return tmp<surfaceScalarField>
|
||||
(
|
||||
new surfaceScalarField
|
||||
(
|
||||
IOobject
|
||||
vf.name() + "BlendingFactor",
|
||||
scalar(1)
|
||||
- max
|
||||
(
|
||||
vf.name() + "BlendingFactor",
|
||||
mesh.time().timeName(),
|
||||
mesh
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("blendingFactor", dimless, 0.0)
|
||||
min
|
||||
(
|
||||
(
|
||||
mesh.time().deltaT()*mesh.deltaCoeffs()
|
||||
*mag(tUflux)/mesh.magSf()
|
||||
- Co1_
|
||||
)/(Co2_ - Co1_),
|
||||
scalar(1)
|
||||
),
|
||||
scalar(0)
|
||||
)
|
||||
)
|
||||
);
|
||||
|
||||
tbf() =
|
||||
scalar(1)
|
||||
- max
|
||||
(
|
||||
min
|
||||
(
|
||||
(
|
||||
mesh.time().deltaT()*mesh.deltaCoeffs()
|
||||
*mag(faceFlux_)/mesh.magSf()
|
||||
- Co1_
|
||||
)/(Co2_ - Co1_),
|
||||
scalar(1)
|
||||
),
|
||||
scalar(0)
|
||||
);
|
||||
|
||||
return tbf;
|
||||
}
|
||||
|
||||
|
||||
@ -240,9 +245,6 @@ public:
|
||||
{
|
||||
surfaceScalarField bf(blendingFactor(vf));
|
||||
|
||||
Info<< "weights " << max(bf).value() << " " << min(bf).value()
|
||||
<< endl;
|
||||
|
||||
return
|
||||
bf*tScheme1_().weights(vf)
|
||||
+ (scalar(1.0) - bf)*tScheme2_().weights(vf);
|
||||
|
||||
@ -323,7 +323,7 @@ void Foam::meshRefinement::getBafflePatches
|
||||
}
|
||||
|
||||
|
||||
Foam::Map<Foam::labelPair> Foam::meshRefinement::getZoneBafflePatches
|
||||
Foam::Map<Foam::labelPair> Foam::meshRefinement::getZoneBafflePatches
|
||||
(
|
||||
const bool allowBoundary,
|
||||
const labelList& globalToMasterPatch,
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -80,7 +80,7 @@ void Foam::sixDoFRigidBodyMotionConstraints::axis::constrainRotation
|
||||
pointConstraint& pc
|
||||
) const
|
||||
{
|
||||
pc.combine(pointConstraint(Tuple2<label, vector>(2, vector(0,1,0))));
|
||||
pc.combine(pointConstraint(Tuple2<label, vector>(2, axis_)));
|
||||
}
|
||||
|
||||
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -60,14 +60,8 @@ void Foam::sixDoFRigidBodyMotion::write(Ostream& os) const
|
||||
os.writeKeyword("sixDoFRigidBodyMotionRestraint")
|
||||
<< restraintType << token::END_STATEMENT << nl;
|
||||
|
||||
os.writeKeyword(word(restraintType + "Coeffs")) << nl;
|
||||
|
||||
os << indent << token::BEGIN_BLOCK << nl << incrIndent;
|
||||
|
||||
restraints_[rI].write(os);
|
||||
|
||||
os << decrIndent << indent << token::END_BLOCK << nl;
|
||||
|
||||
os << decrIndent << indent << token::END_BLOCK << endl;
|
||||
}
|
||||
|
||||
@ -91,14 +85,8 @@ void Foam::sixDoFRigidBodyMotion::write(Ostream& os) const
|
||||
|
||||
constraints_[rI].sixDoFRigidBodyMotionConstraint::write(os);
|
||||
|
||||
os.writeKeyword(word(constraintType + "Coeffs")) << nl;
|
||||
|
||||
os << indent << token::BEGIN_BLOCK << nl << incrIndent;
|
||||
|
||||
constraints_[rI].write(os);
|
||||
|
||||
os << decrIndent << indent << token::END_BLOCK << nl;
|
||||
|
||||
os << decrIndent << indent << token::END_BLOCK << endl;
|
||||
}
|
||||
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2013-2014 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -105,7 +105,7 @@ Foam::sixDoFRigidBodyMotionSolver::sixDoFRigidBodyMotionSolver
|
||||
{
|
||||
if (rhoName_ == "rhoInf")
|
||||
{
|
||||
rhoInf_ = readScalar(dict.lookup("rhoInf"));
|
||||
rhoInf_ = readScalar(coeffDict().lookup("rhoInf"));
|
||||
}
|
||||
|
||||
// Calculate scaling factor everywhere
|
||||
@ -202,8 +202,16 @@ void Foam::sixDoFRigidBodyMotionSolver::solve()
|
||||
|
||||
f.calcForcesMoment();
|
||||
|
||||
uniformDimensionedVectorField g =
|
||||
db().lookupObject<uniformDimensionedVectorField>("g");
|
||||
dimensionedVector g("g", dimAcceleration, vector::zero);
|
||||
|
||||
if (db().foundObject<uniformDimensionedVectorField>("g"))
|
||||
{
|
||||
g = db().lookupObject<uniformDimensionedVectorField>("g");
|
||||
}
|
||||
else if (coeffDict().found("g"))
|
||||
{
|
||||
coeffDict().lookup("g") >> g;
|
||||
}
|
||||
|
||||
// scalar ramp = min(max((this->db().time().value() - 5)/10, 0), 1);
|
||||
scalar ramp = 1.0;
|
||||
|
||||
@ -28,7 +28,7 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
default none;
|
||||
div(phi,U) Gauss filteredLinear2V 0.2 0.05;
|
||||
div(phi,U) Gauss LUST grad(U);
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,K) Gauss limitedLinear 1;
|
||||
div(phi,Yi_h) Gauss multivariateSelection
|
||||
|
||||
@ -30,7 +30,7 @@ boundaryField
|
||||
|
||||
outlet
|
||||
{
|
||||
type inletOutlet;
|
||||
type pressureInletOutletVelocity;
|
||||
inletValue uniform (0 0 0);
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
|
||||
@ -28,14 +28,12 @@ boundaryField
|
||||
outlet
|
||||
{
|
||||
type waveTransmissive;
|
||||
field p;
|
||||
phi phi;
|
||||
rho rho;
|
||||
psi thermo:psi;
|
||||
gamma 1.3;
|
||||
fieldInf 1e5;
|
||||
lInf 0.3;
|
||||
value uniform 1e5;
|
||||
value $internalField;
|
||||
}
|
||||
|
||||
upperWall
|
||||
|
||||
@ -15,34 +15,15 @@ FoamFile
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
gasThermoTypes
|
||||
{
|
||||
constant
|
||||
{
|
||||
type hePsiThermo;
|
||||
mixture pureMixture;
|
||||
transport const;
|
||||
thermo hConst;
|
||||
equationOfState perfectGas;
|
||||
specie specie;
|
||||
energy sensibleEnthalpy;
|
||||
}
|
||||
|
||||
janaf
|
||||
{
|
||||
type hePsiThermo;
|
||||
mixture pureMixture;
|
||||
transport sutherland;
|
||||
thermo janaf;
|
||||
equationOfState perfectGas;
|
||||
specie specie;
|
||||
energy sensibleEnthalpy;
|
||||
}
|
||||
}
|
||||
|
||||
thermoType
|
||||
{
|
||||
${:gasThermoTypes.constant};
|
||||
type hePsiThermo;
|
||||
mixture pureMixture;
|
||||
transport const;
|
||||
thermo eConst;
|
||||
equationOfState perfectGas;
|
||||
specie specie;
|
||||
energy sensibleInternalEnergy;
|
||||
}
|
||||
|
||||
mixture
|
||||
@ -54,7 +35,7 @@ mixture
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Cp 1000;
|
||||
Cv 712;
|
||||
Hf 0;
|
||||
}
|
||||
transport
|
||||
|
||||
@ -17,7 +17,7 @@ FoamFile
|
||||
|
||||
application rhoPimpleFoam;
|
||||
|
||||
startFrom latestTime;
|
||||
startFrom startTime;
|
||||
|
||||
startTime 0;
|
||||
|
||||
|
||||
@ -28,9 +28,10 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
default none;
|
||||
div(phi,U) Gauss filteredLinear2V 0.2 0;
|
||||
div(phi,h) Gauss filteredLinear2 0.2 0;
|
||||
div(phi,U) Gauss LUST grad(U);
|
||||
div(phi,e) Gauss LUST grad(e);
|
||||
div(phi,K) Gauss linear;
|
||||
div(phiv,p) Gauss linear;
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,B) Gauss limitedLinear 1;
|
||||
div(phi,muTilda) Gauss limitedLinear 1;
|
||||
|
||||
@ -31,15 +31,15 @@ solvers
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
"(U|h|k|nuTilda)"
|
||||
"(U|e|k|nuTilda)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-6;
|
||||
relTol 0.01;
|
||||
}
|
||||
|
||||
"(U|h|k|nuTilda)Final"
|
||||
"(U|e|k|nuTilda)Final"
|
||||
{
|
||||
$U;
|
||||
relTol 0;
|
||||
@ -63,7 +63,7 @@ relaxationFactors
|
||||
}
|
||||
equations
|
||||
{
|
||||
"(U|h|k|epsilon|omega).*" 1;
|
||||
".*" 1;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -23,7 +23,7 @@ solvers
|
||||
tolerance 1e-08;
|
||||
relTol 0.05;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 20;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -23,7 +23,7 @@ solvers
|
||||
tolerance 1e-08;
|
||||
relTol 0.05;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 20;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -23,7 +23,7 @@ solvers
|
||||
tolerance 1e-08;
|
||||
relTol 0.05;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 20;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -15,6 +15,9 @@ runApplication blockMesh
|
||||
runApplication surfaceFeatureExtract
|
||||
|
||||
runApplication snappyHexMesh -overwrite
|
||||
|
||||
runApplication renumberMesh -overwrite
|
||||
|
||||
# force removal of fields generated by snappy
|
||||
\rm -rf 0
|
||||
cp -rf 0.org 0
|
||||
|
||||
@ -5,24 +5,17 @@
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
|
||||
// Calculate Q (from U)
|
||||
Q
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object motionProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
// Where to load it from
|
||||
functionObjectLibs ("libutilityFunctionObjects.so");
|
||||
|
||||
dynamicFvMesh dynamicMotionSolverFvMesh;
|
||||
type Q;
|
||||
|
||||
motionSolverLibs ("libfvMotionSolvers.so" "libsixDoFRigidBodyMotion.so");
|
||||
|
||||
solver displacementLaplacian;
|
||||
|
||||
displacementLaplacianCoeffs
|
||||
{
|
||||
diffusivity inverseDistance (movingBlock);
|
||||
// Output every
|
||||
outputControl outputTime;
|
||||
}
|
||||
|
||||
|
||||
@ -58,4 +58,12 @@ adjustTimeStep yes;
|
||||
|
||||
maxCo 2;
|
||||
|
||||
functions
|
||||
{
|
||||
#include "readFields"
|
||||
#include "Q"
|
||||
#include "surfaces"
|
||||
#include "forces"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
||||
@ -5,34 +5,28 @@
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
|
||||
forces
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object topoSetDict;
|
||||
type forces;
|
||||
|
||||
functionObjectLibs ( "libforces.so" );
|
||||
|
||||
outputControl timeStep;
|
||||
timeInterval 1;
|
||||
|
||||
log yes;
|
||||
|
||||
patches ( "propeller.*" );
|
||||
pName p;
|
||||
UName U;
|
||||
rhoName rhoInf; // Indicates incompressible
|
||||
log true;
|
||||
rhoInf 1; // Redundant for incompressible
|
||||
|
||||
CofR (0 0 0); // Rotation around centre line of propeller
|
||||
pitchAxis (0 1 0);
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
actions
|
||||
(
|
||||
{
|
||||
name c0;
|
||||
type cellSet;
|
||||
action new;
|
||||
source boxToCell;
|
||||
sourceInfo
|
||||
{
|
||||
box (0.35 0.35 0.44) (0.65 0.65 0.56);
|
||||
}
|
||||
}
|
||||
|
||||
{
|
||||
name c0;
|
||||
type cellSet;
|
||||
action invert;
|
||||
}
|
||||
);
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -5,45 +5,16 @@
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object fvSchemes;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
ddtSchemes
|
||||
// Make sure all fields for functionObjects are loaded. Prevents any
|
||||
// problems running with execFlowFunctionObjects.
|
||||
readFields
|
||||
{
|
||||
}
|
||||
// Where to load it from (if not already in solver)
|
||||
functionObjectLibs ("libfieldFunctionObjects.so");
|
||||
|
||||
gradSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
}
|
||||
|
||||
divSchemes
|
||||
{
|
||||
}
|
||||
|
||||
laplacianSchemes
|
||||
{
|
||||
default Gauss linear corrected;
|
||||
}
|
||||
|
||||
interpolationSchemes
|
||||
{
|
||||
default linear;
|
||||
}
|
||||
|
||||
snGradSchemes
|
||||
{
|
||||
}
|
||||
|
||||
fluxRequired
|
||||
{
|
||||
type readFields;
|
||||
fields (p U k);
|
||||
}
|
||||
|
||||
|
||||
@ -0,0 +1,52 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
surfaces
|
||||
{
|
||||
type surfaces;
|
||||
functionObjectLibs ("libsampling.so");
|
||||
outputControl outputTime;
|
||||
|
||||
surfaceFormat vtk;
|
||||
fields ( p U Q);
|
||||
|
||||
interpolationScheme cellPoint;
|
||||
|
||||
surfaces
|
||||
(
|
||||
zNormal
|
||||
{
|
||||
type cuttingPlane;
|
||||
planeType pointAndNormal;
|
||||
pointAndNormalDict
|
||||
{
|
||||
basePoint (0 0 0);
|
||||
normalVector (0 0 1);
|
||||
}
|
||||
interpolate true;
|
||||
}
|
||||
|
||||
isoQ
|
||||
{
|
||||
type isoSurface;
|
||||
isoField Q;
|
||||
isoValue 1000;
|
||||
interpolate true;
|
||||
}
|
||||
|
||||
propeller
|
||||
{
|
||||
type patch;
|
||||
patches ( "propeller.*" );
|
||||
interpolate true;
|
||||
}
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -23,79 +23,7 @@ boundaryField
|
||||
{
|
||||
wing
|
||||
{
|
||||
type sixDoFRigidBodyDisplacement;
|
||||
mass 22.9;
|
||||
centreOfMass (0.4974612746 -0.01671895744 0.125);
|
||||
momentOfInertia (1.958864357 3.920839234 2.057121362);
|
||||
orientation
|
||||
(
|
||||
0.9953705935 0.09611129781 0
|
||||
-0.09611129781 0.9953705935 0
|
||||
0 0 1
|
||||
);
|
||||
velocity (0 0 0);
|
||||
acceleration (0 0 0);
|
||||
angularMomentum (0 0 -2.0);
|
||||
torque (0 0 0);
|
||||
rhoName rhoInf;
|
||||
rhoInf 1;
|
||||
g (0 -9.81 0);
|
||||
report on;
|
||||
constraints
|
||||
{
|
||||
maxIterations 500;
|
||||
|
||||
yLine
|
||||
{
|
||||
sixDoFRigidBodyMotionConstraint line;
|
||||
lineCoeffs
|
||||
{
|
||||
refPoint (0.25 0.007 0.125);
|
||||
direction (0 1 0);
|
||||
}
|
||||
tolerance 1e-9;
|
||||
relaxationFactor 0.7;
|
||||
}
|
||||
|
||||
zAxis
|
||||
{
|
||||
sixDoFRigidBodyMotionConstraint axis;
|
||||
axisCoeffs
|
||||
{
|
||||
axis (0 0 1);
|
||||
}
|
||||
tolerance 1e-06;
|
||||
relaxationFactor 0.7;
|
||||
}
|
||||
}
|
||||
restraints
|
||||
{
|
||||
verticalSpring
|
||||
{
|
||||
sixDoFRigidBodyMotionRestraint linearSpring;
|
||||
|
||||
linearSpringCoeffs
|
||||
{
|
||||
anchor (0.25 0.007 0.125);
|
||||
refAttachmentPt (0.25 0.007 0.125);
|
||||
stiffness 4000;
|
||||
damping 2;
|
||||
restLength 0;
|
||||
}
|
||||
}
|
||||
axialSpring
|
||||
{
|
||||
sixDoFRigidBodyMotionRestraint linearAxialAngularSpring;
|
||||
|
||||
linearAxialAngularSpringCoeffs
|
||||
{
|
||||
axis (0 0 1);
|
||||
stiffness 700;
|
||||
damping 0.5;
|
||||
referenceOrientation $orientation;
|
||||
}
|
||||
}
|
||||
}
|
||||
type calculated;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
|
||||
|
||||
@ -10,19 +10,75 @@ FoamFile
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object motionProperties;
|
||||
object dynamicMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dynamicFvMesh dynamicMotionSolverFvMesh;
|
||||
|
||||
motionSolverLibs ("libfvMotionSolvers.so" "libsixDoFRigidBodyMotion.so");
|
||||
motionSolverLibs ("libsixDoFRigidBodyMotion.so");
|
||||
|
||||
solver displacementLaplacian;
|
||||
solver sixDoFRigidBodyMotion;
|
||||
|
||||
displacementLaplacianCoeffs
|
||||
sixDoFRigidBodyMotionCoeffs
|
||||
{
|
||||
diffusivity inverseDistance (wing);
|
||||
patches (wing);
|
||||
innerDistance 0.3;
|
||||
outerDistance 1;
|
||||
|
||||
mass 22.9;
|
||||
centreOfMass (0.4974612746 -0.01671895744 0.125);
|
||||
momentOfInertia (1.958864357 3.920839234 2.057121362);
|
||||
orientation
|
||||
(
|
||||
0.9953705935 0.09611129781 0
|
||||
-0.09611129781 0.9953705935 0
|
||||
0 0 1
|
||||
);
|
||||
angularMomentum (0 0 -2);
|
||||
g (0 -9.81 0);
|
||||
rhoName rhoInf;
|
||||
rhoInf 1;
|
||||
report on;
|
||||
|
||||
constraints
|
||||
{
|
||||
yLine
|
||||
{
|
||||
sixDoFRigidBodyMotionConstraint line;
|
||||
direction (0 1 0);
|
||||
}
|
||||
|
||||
zAxis
|
||||
{
|
||||
sixDoFRigidBodyMotionConstraint axis;
|
||||
axis (0 0 1);
|
||||
}
|
||||
}
|
||||
|
||||
restraints
|
||||
{
|
||||
verticalSpring
|
||||
{
|
||||
sixDoFRigidBodyMotionRestraint linearSpring;
|
||||
|
||||
anchor (0.25 0.007 0.125);
|
||||
refAttachmentPt (0.25 0.007 0.125);
|
||||
stiffness 4000;
|
||||
damping 2;
|
||||
restLength 0;
|
||||
}
|
||||
|
||||
axialSpring
|
||||
{
|
||||
sixDoFRigidBodyMotionRestraint linearAxialAngularSpring;
|
||||
|
||||
axis (0 0 1);
|
||||
stiffness 700;
|
||||
damping 0.5;
|
||||
referenceOrientation $orientation;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
|
||||
@ -21,40 +21,40 @@ FoamFile
|
||||
{
|
||||
type patch;
|
||||
nFaces 72;
|
||||
startFace 22911;
|
||||
startFace 25179;
|
||||
}
|
||||
inlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 40;
|
||||
startFace 22983;
|
||||
startFace 25251;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 62;
|
||||
startFace 23023;
|
||||
startFace 25291;
|
||||
}
|
||||
front
|
||||
{
|
||||
type empty;
|
||||
inGroups 1(empty);
|
||||
nFaces 11431;
|
||||
startFace 23085;
|
||||
nFaces 12565;
|
||||
startFace 25353;
|
||||
}
|
||||
back
|
||||
{
|
||||
type empty;
|
||||
inGroups 1(empty);
|
||||
nFaces 11431;
|
||||
startFace 34516;
|
||||
nFaces 12565;
|
||||
startFace 37918;
|
||||
}
|
||||
wing
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(wall);
|
||||
nFaces 378;
|
||||
startFace 45947;
|
||||
startFace 50483;
|
||||
}
|
||||
)
|
||||
|
||||
|
||||
@ -17,7 +17,7 @@ FoamFile
|
||||
|
||||
application pimpleDyMFoam;
|
||||
|
||||
startFrom latestTime;
|
||||
startFrom startTime;
|
||||
|
||||
startTime 0;
|
||||
|
||||
@ -49,28 +49,4 @@ adjustTimeStep yes;
|
||||
|
||||
maxCo 0.9;
|
||||
|
||||
libs
|
||||
(
|
||||
"libOpenFOAM.so"
|
||||
"libforces.so"
|
||||
);
|
||||
|
||||
functions
|
||||
{
|
||||
forces
|
||||
{
|
||||
type forces;
|
||||
functionObjectLibs ("libforces.so");
|
||||
outputControl timeStep;
|
||||
outputInterval 10;
|
||||
patches (wing);
|
||||
pName p;
|
||||
UName U;
|
||||
rhoName rhoInf;
|
||||
log true;
|
||||
rhoInf 1;
|
||||
CofR (0.25 0.007 0.125);
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
||||
@ -21,40 +21,40 @@ FoamFile
|
||||
{
|
||||
type patch;
|
||||
nFaces 72;
|
||||
startFace 22911;
|
||||
startFace 25179;
|
||||
}
|
||||
inlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 40;
|
||||
startFace 22983;
|
||||
startFace 25251;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 62;
|
||||
startFace 23023;
|
||||
startFace 25291;
|
||||
}
|
||||
front
|
||||
{
|
||||
type empty;
|
||||
inGroups 1(empty);
|
||||
nFaces 11431;
|
||||
startFace 23085;
|
||||
nFaces 12565;
|
||||
startFace 25353;
|
||||
}
|
||||
back
|
||||
{
|
||||
type empty;
|
||||
inGroups 1(empty);
|
||||
nFaces 11431;
|
||||
startFace 34516;
|
||||
nFaces 12565;
|
||||
startFace 37918;
|
||||
}
|
||||
wing
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(wall);
|
||||
nFaces 378;
|
||||
startFace 45947;
|
||||
startFace 50483;
|
||||
}
|
||||
)
|
||||
|
||||
|
||||
@ -28,7 +28,7 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
default none;
|
||||
div(phi,U) Gauss filteredLinear2V 0.2 0;
|
||||
div(phi,U) Gauss LUST grad(U);
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,B) Gauss limitedLinear 1;
|
||||
div(phi,nuTilda) Gauss limitedLinear 1;
|
||||
|
||||
@ -23,7 +23,7 @@ solvers
|
||||
tolerance 1e-08;
|
||||
relTol 0.05;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 20;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -23,7 +23,7 @@ solvers
|
||||
tolerance 1e-08;
|
||||
relTol 0.05;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 20;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -23,7 +23,7 @@ solvers
|
||||
tolerance 1e-08;
|
||||
relTol 0.05;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 20;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
@ -32,7 +32,7 @@ solvers
|
||||
"(U|k|epsilon)"
|
||||
{
|
||||
solver smoothSolver;
|
||||
smoother GaussSeidel;
|
||||
smoother symGaussSeidel;
|
||||
nSweeps 2;
|
||||
tolerance 1e-07;
|
||||
relTol 0.1;
|
||||
|
||||
@ -43,16 +43,16 @@ solvers
|
||||
|
||||
"(U.air|k|omega)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
"(U.air|k|omega)Final"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
@ -29,12 +29,12 @@ divSchemes
|
||||
{
|
||||
default none;
|
||||
|
||||
div(phi,U) Gauss limitedLinearV 1;
|
||||
div(phi,U) Gauss limitedLinear 1;
|
||||
div(phi,Yi_h) Gauss limitedLinear 1;
|
||||
div(phi,K) Gauss limitedLinear 1;
|
||||
div(phid,p) Gauss limitedLinear 1;
|
||||
div(phi,epsilon) Gauss limitedLinear 1;
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
|
||||
div((muEff*dev2(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
@ -24,10 +24,16 @@ solvers
|
||||
|
||||
p
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
tolerance 1e-6;
|
||||
relTol 0.1;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
|
||||
pFinal
|
||||
@ -38,8 +44,8 @@ solvers
|
||||
|
||||
"(U|h|k|epsilon)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-6;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
@ -23,7 +23,7 @@ startTime 0;
|
||||
|
||||
stopAt endTime;
|
||||
|
||||
endTime 10;
|
||||
endTime 5;
|
||||
|
||||
deltaT 2e-4;
|
||||
|
||||
|
||||
@ -43,26 +43,30 @@ solvers
|
||||
|
||||
"(U.air|k|omega)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
"(U.air|k|omega)Final"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
kinematicCloud:alpha
|
||||
{
|
||||
solver PCG; //PBiCG;
|
||||
preconditioner DIC; //DILU;
|
||||
solver GAMG;
|
||||
tolerance 1e-06;
|
||||
relTol 0;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -19,10 +19,16 @@ solvers
|
||||
{
|
||||
p
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
tolerance 1e-06;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.01;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
|
||||
pFinal
|
||||
@ -33,24 +39,24 @@ solvers
|
||||
|
||||
U.air
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
U.airFinal
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
kinematicCloud:alpha
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-06;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
@ -43,8 +43,8 @@ solvers
|
||||
|
||||
"(U|k|omega).*"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
@ -41,7 +41,7 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
@ -56,8 +56,8 @@ solvers
|
||||
|
||||
"(U|h|R|k|epsilon)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
@ -71,16 +71,12 @@ solvers
|
||||
|
||||
"(Yi|CO2|O2|N2|CH4|H2|H2O|CO)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
tolerance 1e-06;
|
||||
relTol 0.1;
|
||||
$h;
|
||||
}
|
||||
|
||||
"(Yi|CO2|O2|N2|CH4|H2|H2O|CO)Final"
|
||||
{
|
||||
$Yi;
|
||||
preconditioner DILU;
|
||||
tolerance 1e-06;
|
||||
relTol 0;
|
||||
}
|
||||
@ -107,4 +103,5 @@ relaxationFactors
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
}
|
||||
@ -54,10 +54,16 @@ solvers
|
||||
|
||||
p_rgh
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
|
||||
p_rghFinal
|
||||
@ -73,6 +79,7 @@ solvers
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
PIMPLE
|
||||
{
|
||||
nCorrectors 2;
|
||||
@ -80,6 +87,7 @@ PIMPLE
|
||||
momentumPredictor yes;
|
||||
}
|
||||
|
||||
|
||||
relaxationFactors
|
||||
{
|
||||
fields
|
||||
|
||||
@ -19,15 +19,15 @@ solvers
|
||||
{
|
||||
"(Uf|hf|deltaf\*rhof)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-10;
|
||||
relTol 0;
|
||||
}
|
||||
deltaf
|
||||
{
|
||||
solver PBiCG; // PCG;
|
||||
preconditioner DILU; // DIC;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-10;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h|k|epsilon)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
}
|
||||
@ -49,7 +49,7 @@ solvers
|
||||
|
||||
p_rgh
|
||||
{
|
||||
solver GAMG;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
smoother DICGaussSeidel;
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
}
|
||||
@ -54,10 +54,16 @@ solvers
|
||||
|
||||
p_rgh
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
|
||||
p_rghFinal
|
||||
|
||||
@ -19,15 +19,15 @@ solvers
|
||||
{
|
||||
"(Uf|hf|deltaf\*rhof)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-10;
|
||||
relTol 0;
|
||||
}
|
||||
deltaf
|
||||
{
|
||||
solver PBiCG; // PCG;
|
||||
preconditioner DILU; // DIC;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-10;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
}
|
||||
@ -54,10 +54,16 @@ solvers
|
||||
|
||||
p_rgh
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
|
||||
p_rghFinal
|
||||
|
||||
@ -19,15 +19,15 @@ solvers
|
||||
{
|
||||
"(Uf|hf|deltaf\*rhof)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-10;
|
||||
relTol 0;
|
||||
}
|
||||
deltaf
|
||||
{
|
||||
solver PBiCG; // PCG;
|
||||
preconditioner DILU; // DIC;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-10;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h|k|epsilon)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-06;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h|k|epsilon)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
@ -34,8 +34,8 @@ solvers
|
||||
|
||||
"(U|h|k|omega)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-06;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
@ -34,18 +34,24 @@ solvers
|
||||
|
||||
"(U|h)"
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
p
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
tolerance 1e-06;
|
||||
solver GAMG;
|
||||
tolerance 0;
|
||||
relTol 0.1;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
|
||||
pFinal
|
||||
|
||||
@ -6,6 +6,8 @@ cd ${0%/*} || exit 1 # run from this directory
|
||||
|
||||
cleanCase
|
||||
|
||||
rm -rf 0 > /dev/null 2>&1
|
||||
|
||||
rm -rf constant/triSurface/blob.stl.gz > /dev/null 2>&1
|
||||
|
||||
rm -rf constant/cellAlignments > /dev/null 2>&1
|
||||
|
||||
@ -30,7 +30,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression 1;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phid1,p_rgh) Gauss upwind;
|
||||
div(phid2,p_rgh) Gauss upwind;
|
||||
div(rhoPhi,T) Gauss linearUpwind unlimited;
|
||||
|
||||
@ -23,6 +23,9 @@ runApplication surfaceBooleanFeatures intersection \
|
||||
runApplication foamyHexMesh
|
||||
|
||||
runApplication collapseEdges -latestTime -collapseFaceSet indirectPatchFaces
|
||||
mv log.collapseEdges log.collapseFaceSet
|
||||
|
||||
runApplication collapseEdges -latestTime -collapseFaces
|
||||
|
||||
|
||||
# ----------------------------------------------------------------- end-of-file
|
||||
|
||||
@ -1,66 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class polyBoundaryMesh;
|
||||
location "constant/polyMesh";
|
||||
object boundary;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
7
|
||||
(
|
||||
letters_text
|
||||
{
|
||||
type patch;
|
||||
nFaces 1064;
|
||||
startFace 132063;
|
||||
}
|
||||
box_bottom
|
||||
{
|
||||
type patch;
|
||||
nFaces 324;
|
||||
startFace 133127;
|
||||
}
|
||||
box_top
|
||||
{
|
||||
type patch;
|
||||
nFaces 324;
|
||||
startFace 133451;
|
||||
}
|
||||
box_inlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 207;
|
||||
startFace 133775;
|
||||
}
|
||||
box_outlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 207;
|
||||
startFace 133982;
|
||||
}
|
||||
front
|
||||
{
|
||||
type empty;
|
||||
inGroups 1(empty);
|
||||
nFaces 64792;
|
||||
startFace 134189;
|
||||
}
|
||||
back
|
||||
{
|
||||
type empty;
|
||||
inGroups 1(empty);
|
||||
nFaces 64792;
|
||||
startFace 198981;
|
||||
}
|
||||
)
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -1,80 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class pointVectorField;
|
||||
location "0.01";
|
||||
object pointDisplacement;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 1 0 0 0 0 0];
|
||||
|
||||
internalField uniform (0 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
stationaryWalls
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
movingBlock
|
||||
{
|
||||
type uncoupledSixDoFRigidBodyDisplacement;
|
||||
centreOfMass (0.5 0.5 0.5);
|
||||
momentOfInertia (0.1052 0.1052 0.1778);
|
||||
mass 9.6;
|
||||
velocity (0 0 0);
|
||||
acceleration (0 0 0);
|
||||
angularMomentum (0 0 0);
|
||||
torque (0 0 0);
|
||||
gravity (0 0 0);
|
||||
rhoInf 1;
|
||||
report on;
|
||||
restraints
|
||||
{
|
||||
topSpring
|
||||
{
|
||||
sixDoFRigidBodyMotionRestraint linearSpring;
|
||||
|
||||
linearSpringCoeffs
|
||||
{
|
||||
anchor (0.5 0.5 1);
|
||||
refAttachmentPt $centreOfMass;
|
||||
stiffness 5000;
|
||||
damping 50;
|
||||
restLength 0.4;
|
||||
}
|
||||
}
|
||||
}
|
||||
constraints
|
||||
{
|
||||
fixedOrientation
|
||||
{
|
||||
sixDoFRigidBodyMotionConstraint orientation;
|
||||
}
|
||||
|
||||
zLine
|
||||
{
|
||||
sixDoFRigidBodyMotionConstraint line;
|
||||
|
||||
lineCoeffs
|
||||
{
|
||||
direction (0 0 1);
|
||||
}
|
||||
}
|
||||
}
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -1,11 +0,0 @@
|
||||
#!/bin/sh
|
||||
cd ${0%/*} || exit 1 # run from this directory
|
||||
|
||||
# Source tutorial clean functions
|
||||
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
|
||||
|
||||
rm -rf 0 t_vs_cm t_vs_lv shm.eps > /dev/null 2>&1
|
||||
|
||||
cleanCase
|
||||
|
||||
# ----------------------------------------------------------------- end-of-file
|
||||
@ -1,19 +0,0 @@
|
||||
#!/bin/sh
|
||||
|
||||
cd ${0%/*} || exit 1 # run from this directory
|
||||
|
||||
# Source tutorial run functions
|
||||
. $WM_PROJECT_DIR/bin/tools/RunFunctions
|
||||
|
||||
# Set application name
|
||||
application=`getApplication`
|
||||
|
||||
runApplication blockMesh
|
||||
runApplication topoSet
|
||||
runApplication subsetMesh -overwrite c0 -patch movingBlock
|
||||
cp -r 0.org 0 > /dev/null 2>&1
|
||||
runApplication $application
|
||||
./extractData log.$application
|
||||
gnuplot shm.gnuplot
|
||||
|
||||
# ----------------------------------------------------------------- end-of-file
|
||||
@ -1,66 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object blockMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
convertToMeters 1;
|
||||
|
||||
vertices
|
||||
(
|
||||
(0 0 0)
|
||||
(1 0 0)
|
||||
(1 1 0)
|
||||
(0 1 0)
|
||||
(0 0 1)
|
||||
(1 0 1)
|
||||
(1 1 1)
|
||||
(0 1 1)
|
||||
);
|
||||
|
||||
blocks
|
||||
(
|
||||
hex (0 1 2 3 4 5 6 7) (9 9 14) simpleGrading (1 1 1)
|
||||
);
|
||||
|
||||
edges
|
||||
(
|
||||
);
|
||||
|
||||
boundary
|
||||
(
|
||||
stationaryWalls
|
||||
{
|
||||
type patch;
|
||||
faces
|
||||
(
|
||||
(0 3 2 1)
|
||||
(2 6 5 1)
|
||||
(1 5 4 0)
|
||||
(3 7 6 2)
|
||||
(0 4 7 3)
|
||||
(4 5 6 7)
|
||||
);
|
||||
}
|
||||
movingBlock
|
||||
{
|
||||
type patch;
|
||||
faces ();
|
||||
}
|
||||
);
|
||||
|
||||
mergePatchPairs
|
||||
(
|
||||
);
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -1,42 +0,0 @@
|
||||
#!/bin/sh
|
||||
#------------------------------------------------------------------------------
|
||||
# ========= |
|
||||
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
# \\ / O peration |
|
||||
# \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
|
||||
# \\/ M anipulation |
|
||||
#-------------------------------------------------------------------------------
|
||||
# License
|
||||
# This file is part of OpenFOAM.
|
||||
#
|
||||
# OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
# under the terms of the GNU General Public License as published by
|
||||
# the Free Software Foundation, either version 3 of the License, or
|
||||
# (at your option) any later version.
|
||||
#
|
||||
# OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
# for more details.
|
||||
#
|
||||
# You should have received a copy of the GNU General Public License
|
||||
# along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
#
|
||||
# Script
|
||||
# extractData
|
||||
#
|
||||
# Description
|
||||
# Extracts motion data from a simple harmonic motion dynamicMesh case
|
||||
#
|
||||
#------------------------------------------------------------------------------
|
||||
|
||||
grep "Centre of mass" $1 | cut -d ":" -f 2 | cut -d " " -f 4 | tr -d ")" > cM
|
||||
grep "Linear velocity" $1 | cut -d ":" -f 2 | cut -d " " -f 4 | tr -d ")" > lV
|
||||
grep -e "^Time = " $1 | cut -d " " -f 3 > times
|
||||
|
||||
paste times cM > t_vs_cm
|
||||
paste times lV > t_vs_lv
|
||||
|
||||
rm cM lV times
|
||||
|
||||
#------------------------------------------------------------------------------
|
||||
@ -1,76 +0,0 @@
|
||||
#------------------------------------------------------------------------------
|
||||
# ========= |
|
||||
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
# \\ / O peration |
|
||||
# \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
|
||||
# \\/ M anipulation |
|
||||
#-------------------------------------------------------------------------------
|
||||
# License
|
||||
# This file is part of OpenFOAM.
|
||||
#
|
||||
# OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
# under the terms of the GNU General Public License as published by
|
||||
# the Free Software Foundation, either version 3 of the License, or
|
||||
# (at your option) any later version.
|
||||
#
|
||||
# OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
# for more details.
|
||||
#
|
||||
# You should have received a copy of the GNU General Public License
|
||||
# along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
#
|
||||
# Script
|
||||
# shm.gnuplot
|
||||
#
|
||||
# Description
|
||||
# Creates an .eps graph of OpenFOAM results vs analytical solution
|
||||
# for a simple harmonic motion dynamicMesh case
|
||||
#
|
||||
#------------------------------------------------------------------------------
|
||||
|
||||
reset
|
||||
|
||||
set samples 2000
|
||||
|
||||
k = 5000.0
|
||||
m = 9.6
|
||||
c = 50.0
|
||||
a = -0.1
|
||||
|
||||
omega = sqrt(k/m)
|
||||
zeta = c/(2.0*m*omega)
|
||||
|
||||
phi = atan((sqrt(1.0 - zeta**2))/zeta)
|
||||
A = a/sin(phi)
|
||||
|
||||
pos(A, t, omega, phi, zeta) = A*exp(-zeta*omega*t)*sin(sqrt(1-zeta**2)*omega*t + phi)
|
||||
vel(A, t, omega, phi, zeta) = \
|
||||
A*exp(-zeta*omega*t)*\
|
||||
( \
|
||||
sqrt(1-zeta**2)*omega*cos(sqrt(1-zeta**2)*omega*t + phi) \
|
||||
- zeta*omega*sin(sqrt(1-zeta**2)*omega*t + phi) \
|
||||
)
|
||||
|
||||
set xlabel "Time/[s]"
|
||||
set ylabel "Position"
|
||||
|
||||
set ytics nomirror
|
||||
set y2tics
|
||||
|
||||
set yrange [-0.1:0.1]
|
||||
set y2range [-2:2]
|
||||
|
||||
set xzeroaxis
|
||||
|
||||
set terminal postscript eps color enhanced solid
|
||||
set output "shm.eps"
|
||||
|
||||
plot \
|
||||
"t_vs_cm" u 1:($2 - 0.6) w l t "Simulation, centre of mass relative to start", \
|
||||
pos(A, x, omega, phi, zeta) w l t "Analytical solution, centre of mass", \
|
||||
"t_vs_lv" u 1:2 w l axes x1y2 t "Simulation, vertical velocity", \
|
||||
vel(A, x, omega, phi, zeta) w l axes x1y2 t "Analytical solution, vertical velocity"
|
||||
|
||||
#------------------------------------------------------------------------------
|
||||
@ -1,61 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object controlDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
application moveDynamicMesh;
|
||||
|
||||
startFrom startTime;
|
||||
|
||||
startTime 0;
|
||||
|
||||
stopAt endTime;
|
||||
|
||||
endTime 4;
|
||||
|
||||
deltaT 0.002;
|
||||
|
||||
writeControl adjustableRunTime;
|
||||
|
||||
writeInterval 0.25;
|
||||
|
||||
purgeWrite 0;
|
||||
|
||||
writeFormat ascii;
|
||||
|
||||
writePrecision 12;
|
||||
|
||||
writeCompression uncompressed;
|
||||
|
||||
timeFormat general;
|
||||
|
||||
timePrecision 6;
|
||||
|
||||
runTimeModifiable yes;
|
||||
|
||||
adjustTimeStep yes;
|
||||
|
||||
maxCo 0.2;
|
||||
|
||||
maxDeltaT 0.025;
|
||||
|
||||
libs
|
||||
(
|
||||
"libOpenFOAM.so"
|
||||
"libincompressibleRASModels.so"
|
||||
"libforces.so"
|
||||
);
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -1,34 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object fvSolution;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
solvers
|
||||
{
|
||||
cellDisplacement
|
||||
{
|
||||
solver GAMG;
|
||||
tolerance 1e-06;
|
||||
relTol 0;
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,omega) Gauss upwind;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linear;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss limitedLinearV 1;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -28,7 +28,7 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression 1;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
|
||||
div(rhoPhi,U) Gauss vanLeerV;
|
||||
div(phi,thermo:rho.water) Gauss linear;
|
||||
|
||||
@ -28,7 +28,7 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression 1;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,thermo:rho.water) Gauss upwind;
|
||||
|
||||
@ -41,7 +41,7 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -28,7 +28,7 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression 1;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,thermo:rho.water) Gauss upwind;
|
||||
|
||||
@ -36,7 +36,7 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
"div\(phi,.*rho.*\)" Gauss upwind;
|
||||
div(rhoPhi,T) Gauss upwind;
|
||||
div(rhoPhi,K) Gauss upwind;
|
||||
|
||||
@ -29,10 +29,10 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration off;
|
||||
cacheAgglomeration on;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 2;
|
||||
mergeLevels 1;
|
||||
}
|
||||
tolerance 1e-05;
|
||||
relTol 0;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss vanLeerV;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,epsilon) Gauss upwind;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
|
||||
@ -30,7 +30,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression 1;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phid1,p_rgh) Gauss upwind;
|
||||
div(phid2,p_rgh) Gauss upwind;
|
||||
div(rhoPhi,T) Gauss linearUpwind unlimited;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -27,9 +27,9 @@ gradSchemes
|
||||
|
||||
divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss limitedLinearV 1;
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linear;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,B) Gauss limitedLinear 1;
|
||||
div(B) Gauss linear;
|
||||
|
||||
@ -36,7 +36,7 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linear;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,epsilon) Gauss upwind;
|
||||
div(phi,R) Gauss upwind;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linear;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,epsilon) Gauss upwind;
|
||||
div(phi,R) Gauss upwind;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,omega) Gauss upwind;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
|
||||
@ -36,7 +36,7 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -31,7 +31,7 @@ divSchemes
|
||||
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
|
||||
"div\(phi,(k|omega)\)" Gauss upwind;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
|
||||
@ -44,7 +44,7 @@ solvers
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss linear;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,epsilon) Gauss upwind;
|
||||
div(phi,R) Gauss upwind;
|
||||
|
||||
@ -29,7 +29,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss limitedLinearV 1;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -31,7 +31,7 @@ divSchemes
|
||||
default none;
|
||||
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
UD Gauss upwind;
|
||||
div(rhoPhi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,k) Gauss upwind;
|
||||
|
||||
@ -31,7 +31,7 @@ divSchemes
|
||||
div(phi,omega) Gauss linearUpwind grad(omega);
|
||||
div(phi,k) Gauss linearUpwind grad(k);
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -81,7 +81,7 @@ solvers
|
||||
nPostSweeps 2;
|
||||
nFinestSweeps 2;
|
||||
|
||||
cacheAgglomeration false;
|
||||
cacheAgglomeration true;
|
||||
nCellsInCoarsestLevel 10;
|
||||
agglomerator faceAreaPair;
|
||||
mergeLevels 1;
|
||||
|
||||
@ -31,7 +31,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
@ -31,7 +31,7 @@ divSchemes
|
||||
{
|
||||
div(rhoPhi,U) Gauss upwind;
|
||||
div(phi,alpha) Gauss vanLeer;
|
||||
div(phirb,alpha) Gauss interfaceCompression;
|
||||
div(phirb,alpha) Gauss linear;
|
||||
div((muEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
|
||||
Reference in New Issue
Block a user