mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-11-28 03:28:01 +00:00
Merge branch 'master' of /home/dm4/OpenFOAM/OpenFOAM-dev
This commit is contained in:
@ -5,5 +5,6 @@ set -x
|
||||
wmake
|
||||
wmake rhoPimplecFoam
|
||||
wmake rhoLTSPimpleFoam
|
||||
wmake rhoPimpleDyMFoam
|
||||
|
||||
# ----------------------------------------------------------------- end-of-file
|
||||
|
||||
@ -0,0 +1,3 @@
|
||||
rhoPimpleDyMFoam.C
|
||||
|
||||
EXE = $(FOAM_APPBIN)/rhoPimpleDyMFoam
|
||||
@ -0,0 +1,27 @@
|
||||
EXE_INC = \
|
||||
-I.. \
|
||||
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
|
||||
-I$(LIB_SRC)/turbulenceModels/compressible/turbulenceModel \
|
||||
-I$(LIB_SRC)/finiteVolume/cfdTools \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude \
|
||||
-I$(LIB_SRC)/fvOptions/lnInclude \
|
||||
-I$(LIB_SRC)/dynamicFvMesh/lnInclude \
|
||||
-I$(LIB_SRC)/dynamicMesh/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfluidThermophysicalModels \
|
||||
-lspecie \
|
||||
-lcompressibleTurbulenceModel \
|
||||
-lcompressibleRASModels \
|
||||
-lcompressibleLESModels \
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lsampling \
|
||||
-lfvOptions \
|
||||
-ldynamicFvMesh \
|
||||
-ltopoChangerFvMesh \
|
||||
-ldynamicMesh \
|
||||
-lmeshTools
|
||||
@ -0,0 +1,59 @@
|
||||
{
|
||||
if (mesh.changing())
|
||||
{
|
||||
forAll(U.boundaryField(), patchi)
|
||||
{
|
||||
if (U.boundaryField()[patchi].fixesValue())
|
||||
{
|
||||
U.boundaryField()[patchi].initEvaluate();
|
||||
}
|
||||
}
|
||||
|
||||
forAll(U.boundaryField(), patchi)
|
||||
{
|
||||
if (U.boundaryField()[patchi].fixesValue())
|
||||
{
|
||||
U.boundaryField()[patchi].evaluate();
|
||||
|
||||
phi.boundaryField()[patchi] =
|
||||
rho.boundaryField()[patchi]
|
||||
*(
|
||||
U.boundaryField()[patchi]
|
||||
& mesh.Sf().boundaryField()[patchi]
|
||||
);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
volScalarField pcorr
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"pcorr",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("pcorr", p.dimensions(), 0.0),
|
||||
pcorrTypes
|
||||
);
|
||||
|
||||
dimensionedScalar Dp("Dp", dimTime, 1.0);
|
||||
|
||||
while (pimple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix pcorrEqn
|
||||
(
|
||||
fvm::laplacian(Dp, pcorr) == fvc::div(phi) - divrhoU
|
||||
);
|
||||
|
||||
pcorrEqn.solve();
|
||||
|
||||
if (pimple.finalNonOrthogonalIter())
|
||||
{
|
||||
phi -= pcorrEqn.flux();
|
||||
}
|
||||
}
|
||||
}
|
||||
@ -0,0 +1,13 @@
|
||||
wordList pcorrTypes
|
||||
(
|
||||
p.boundaryField().size(),
|
||||
zeroGradientFvPatchScalarField::typeName
|
||||
);
|
||||
|
||||
for (label i=0; i<p.boundaryField().size(); i++)
|
||||
{
|
||||
if (p.boundaryField()[i].fixesValue())
|
||||
{
|
||||
pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
|
||||
}
|
||||
}
|
||||
@ -0,0 +1,114 @@
|
||||
rho = thermo.rho();
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
rho.relax();
|
||||
|
||||
volScalarField rAU(1.0/UEqn().A());
|
||||
volVectorField HbyA("HbyA", U);
|
||||
HbyA = rAU*UEqn().H();
|
||||
|
||||
if (pimple.nCorrPISO() <= 1)
|
||||
{
|
||||
UEqn.clear();
|
||||
}
|
||||
|
||||
if (pimple.transonic())
|
||||
{
|
||||
surfaceScalarField phid
|
||||
(
|
||||
"phid",
|
||||
fvc::interpolate(psi)
|
||||
*(
|
||||
(fvc::interpolate(HbyA) & mesh.Sf())
|
||||
+ fvc::ddtPhiCorr(rAU, rho, U, phiAbs)
|
||||
)
|
||||
);
|
||||
|
||||
fvOptions.relativeFlux(fvc::interpolate(psi), phid);
|
||||
|
||||
volScalarField Dp("Dp", rho*rAU);
|
||||
|
||||
while (pimple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix pEqn
|
||||
(
|
||||
fvm::ddt(psi, p)
|
||||
+ fvm::div(phid, p)
|
||||
- fvm::laplacian(Dp, p)
|
||||
==
|
||||
fvOptions(psi, p, rho.name())
|
||||
);
|
||||
|
||||
fvOptions.constrain(pEqn);
|
||||
|
||||
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
|
||||
|
||||
if (pimple.finalNonOrthogonalIter())
|
||||
{
|
||||
phi == pEqn.flux();
|
||||
}
|
||||
}
|
||||
}
|
||||
else
|
||||
{
|
||||
surfaceScalarField phiHbyA
|
||||
(
|
||||
"phiHbyA",
|
||||
fvc::interpolate(rho)
|
||||
*(
|
||||
(fvc::interpolate(HbyA) & mesh.Sf())
|
||||
- fvc::meshPhi(rho, U)
|
||||
+ fvc::ddtPhiCorr(rAU, rho, U, phiAbs)
|
||||
)
|
||||
);
|
||||
|
||||
fvOptions.relativeFlux(fvc::interpolate(rho), phiHbyA);
|
||||
|
||||
volScalarField Dp("Dp", rho*rAU);
|
||||
|
||||
while (pimple.correctNonOrthogonal())
|
||||
{
|
||||
// Pressure corrector
|
||||
fvScalarMatrix pEqn
|
||||
(
|
||||
fvm::ddt(psi, p)
|
||||
+ fvc::div(phiHbyA)
|
||||
- fvm::laplacian(Dp, p)
|
||||
==
|
||||
fvOptions(psi, p, rho.name())
|
||||
);
|
||||
|
||||
fvOptions.constrain(pEqn);
|
||||
|
||||
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
|
||||
|
||||
if (pimple.finalNonOrthogonalIter())
|
||||
{
|
||||
phi = phiHbyA + pEqn.flux();
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
#include "rhoEqn.H"
|
||||
#include "compressibleContinuityErrs.H"
|
||||
|
||||
// Explicitly relax pressure for momentum corrector
|
||||
p.relax();
|
||||
|
||||
// Recalculate density from the relaxed pressure
|
||||
rho = thermo.rho();
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
rho.relax();
|
||||
Info<< "rho max/min : " << max(rho).value()
|
||||
<< " " << min(rho).value() << endl;
|
||||
|
||||
U = HbyA - rAU*fvc::grad(p);
|
||||
U.correctBoundaryConditions();
|
||||
fvOptions.correct(U);
|
||||
K = 0.5*magSqr(U);
|
||||
|
||||
if (thermo.dpdt())
|
||||
{
|
||||
dpdt = fvc::ddt(p);
|
||||
}
|
||||
@ -0,0 +1,7 @@
|
||||
#include "readTimeControls.H"
|
||||
|
||||
bool correctPhi =
|
||||
pimple.dict().lookupOrDefault<Switch>("correctPhi", true);
|
||||
|
||||
bool checkMeshCourantNo =
|
||||
pimple.dict().lookupOrDefault<Switch>("checkMeshCourantNo", false);
|
||||
@ -0,0 +1,146 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2013 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/>.
|
||||
|
||||
Application
|
||||
rhoPimpleFoam
|
||||
|
||||
Description
|
||||
Transient solver for laminar or turbulent flow of compressible fluids
|
||||
for HVAC and similar applications.
|
||||
|
||||
Uses the flexible PIMPLE (PISO-SIMPLE) solution for time-resolved and
|
||||
pseudo-transient simulations.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "fvCFD.H"
|
||||
#include "dynamicFvMesh.H"
|
||||
#include "psiThermo.H"
|
||||
#include "turbulenceModel.H"
|
||||
#include "bound.H"
|
||||
#include "pimpleControl.H"
|
||||
#include "fvIOoptionList.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
#include "setRootCase.H"
|
||||
#include "createTime.H"
|
||||
#include "createDynamicFvMesh.H"
|
||||
#include "initContinuityErrs.H"
|
||||
|
||||
pimpleControl pimple(mesh);
|
||||
|
||||
#include "readControls.H"
|
||||
#include "createFields.H"
|
||||
#include "createFvOptions.H"
|
||||
#include "createPcorrTypes.H"
|
||||
#include "CourantNo.H"
|
||||
#include "setInitialDeltaT.H"
|
||||
|
||||
// Create old-time absolute flux for ddtPhiCorr
|
||||
surfaceScalarField phiAbs("phiAbs", phi);
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
Info<< "\nStarting time loop\n" << endl;
|
||||
|
||||
while (runTime.run())
|
||||
{
|
||||
#include "readControls.H"
|
||||
#include "compressibleCourantNo.H"
|
||||
|
||||
// Make the fluxes absolute before mesh-motion
|
||||
fvc::makeAbsolute(phi, rho, U);
|
||||
|
||||
// Update absolute flux for ddtPhiCorr
|
||||
phiAbs = phi;
|
||||
|
||||
#include "setDeltaT.H"
|
||||
|
||||
runTime++;
|
||||
|
||||
Info<< "Time = " << runTime.timeName() << nl << endl;
|
||||
|
||||
{
|
||||
// Store divrhoU from the previous time-step/mesh for the correctPhi
|
||||
volScalarField divrhoU(fvc::div(phi));
|
||||
|
||||
// Do any mesh changes
|
||||
mesh.update();
|
||||
|
||||
if (mesh.changing() && correctPhi)
|
||||
{
|
||||
#include "correctPhi.H"
|
||||
}
|
||||
}
|
||||
|
||||
// Make the fluxes relative to the mesh-motion
|
||||
fvc::makeRelative(phi, rho, U);
|
||||
|
||||
if (mesh.changing() && checkMeshCourantNo)
|
||||
{
|
||||
#include "meshCourantNo.H"
|
||||
}
|
||||
|
||||
if (pimple.nCorrPIMPLE() <= 1)
|
||||
{
|
||||
#include "rhoEqn.H"
|
||||
Info<< "rhoEqn max/min : " << max(rho).value()
|
||||
<< " " << min(rho).value() << endl;
|
||||
}
|
||||
|
||||
// --- Pressure-velocity PIMPLE corrector loop
|
||||
while (pimple.loop())
|
||||
{
|
||||
#include "UEqn.H"
|
||||
#include "EEqn.H"
|
||||
|
||||
// --- Pressure corrector loop
|
||||
while (pimple.correct())
|
||||
{
|
||||
#include "pEqn.H"
|
||||
}
|
||||
|
||||
if (pimple.turbCorr())
|
||||
{
|
||||
turbulence->correct();
|
||||
}
|
||||
}
|
||||
|
||||
runTime.write();
|
||||
|
||||
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
||||
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
||||
<< nl << endl;
|
||||
}
|
||||
|
||||
Info<< "End\n" << endl;
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -1,25 +1,24 @@
|
||||
EXE_INC = \
|
||||
-I.. \
|
||||
-I$(LIB_SRC)/dynamicFvMesh/lnInclude \
|
||||
-I$(LIB_SRC)/dynamicMesh/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/turbulenceModels/incompressible/turbulenceModel \
|
||||
-I$(LIB_SRC)/transportModels \
|
||||
-I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I$(LIB_SRC)/fvOptions/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude
|
||||
|
||||
-I$(LIB_SRC)/sampling/lnInclude \
|
||||
-I$(LIB_SRC)/dynamicFvMesh/lnInclude \
|
||||
-I$(LIB_SRC)/dynamicMesh/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
|
||||
EXE_LIBS = \
|
||||
-ldynamicFvMesh \
|
||||
-ltopoChangerFvMesh \
|
||||
-ldynamicMesh \
|
||||
-lmeshTools \
|
||||
-lincompressibleTransportModels \
|
||||
-lincompressibleTurbulenceModel \
|
||||
-lincompressibleRASModels \
|
||||
-lincompressibleLESModels \
|
||||
-lfiniteVolume \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
-lsampling \
|
||||
-ldynamicFvMesh \
|
||||
-ltopoChangerFvMesh \
|
||||
-ldynamicMesh \
|
||||
-lmeshTools
|
||||
|
||||
@ -1,23 +0,0 @@
|
||||
// Solve the Momentum equation
|
||||
|
||||
tmp<fvVectorMatrix> UEqn
|
||||
(
|
||||
fvm::ddt(U)
|
||||
+ fvm::div(phi, U)
|
||||
+ turbulence->divDevReff(U)
|
||||
==
|
||||
fvOptions(U)
|
||||
);
|
||||
|
||||
UEqn().relax();
|
||||
|
||||
fvOptions.constrain(UEqn());
|
||||
|
||||
rAU = 1.0/UEqn().A();
|
||||
|
||||
if (pimple.momentumPredictor())
|
||||
{
|
||||
solve(UEqn() == -fvc::grad(p));
|
||||
|
||||
fvOptions.correct(U);
|
||||
}
|
||||
@ -22,20 +22,6 @@
|
||||
}
|
||||
}
|
||||
|
||||
wordList pcorrTypes
|
||||
(
|
||||
p.boundaryField().size(),
|
||||
zeroGradientFvPatchScalarField::typeName
|
||||
);
|
||||
|
||||
forAll(p.boundaryField(), patchI)
|
||||
{
|
||||
if (p.boundaryField()[patchI].fixesValue())
|
||||
{
|
||||
pcorrTypes[patchI] = fixedValueFvPatchScalarField::typeName;
|
||||
}
|
||||
}
|
||||
|
||||
volScalarField pcorr
|
||||
(
|
||||
IOobject
|
||||
@ -51,11 +37,13 @@
|
||||
pcorrTypes
|
||||
);
|
||||
|
||||
dimensionedScalar Dp("Dp", dimTime, 1.0);
|
||||
|
||||
while (pimple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix pcorrEqn
|
||||
(
|
||||
fvm::laplacian(rAU, pcorr) == fvc::div(phi)
|
||||
fvm::laplacian(Dp, pcorr) == fvc::div(phi)
|
||||
);
|
||||
|
||||
pcorrEqn.setReference(pRefCell, pRefValue);
|
||||
@ -68,6 +56,4 @@
|
||||
}
|
||||
}
|
||||
|
||||
phi.oldTime() = phi;
|
||||
|
||||
#include "continuityErrs.H"
|
||||
|
||||
@ -40,19 +40,3 @@
|
||||
(
|
||||
incompressible::turbulenceModel::New(U, phi, laminarTransport)
|
||||
);
|
||||
|
||||
Info<< "Reading field rAU if present\n" << endl;
|
||||
volScalarField rAU
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rAU",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
runTime.deltaT(),
|
||||
zeroGradientFvPatchScalarField::typeName
|
||||
);
|
||||
|
||||
@ -0,0 +1,13 @@
|
||||
wordList pcorrTypes
|
||||
(
|
||||
p.boundaryField().size(),
|
||||
zeroGradientFvPatchScalarField::typeName
|
||||
);
|
||||
|
||||
for (label i=0; i<p.boundaryField().size(); i++)
|
||||
{
|
||||
if (p.boundaryField()[i].fixesValue())
|
||||
{
|
||||
pcorrTypes[i] = fixedValueFvPatchScalarField::typeName;
|
||||
}
|
||||
}
|
||||
@ -10,13 +10,9 @@ surfaceScalarField phiHbyA
|
||||
(
|
||||
"phiHbyA",
|
||||
(fvc::interpolate(HbyA) & mesh.Sf())
|
||||
+ fvc::ddtPhiCorr(rAU, U, phiAbs)
|
||||
);
|
||||
|
||||
if (ddtPhiCorr)
|
||||
{
|
||||
phiHbyA += fvc::ddtPhiCorr(rAU, U, phi);
|
||||
}
|
||||
|
||||
if (p.needReference())
|
||||
{
|
||||
fvc::makeRelative(phiHbyA, U);
|
||||
|
||||
@ -33,9 +33,9 @@ Description
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "fvCFD.H"
|
||||
#include "dynamicFvMesh.H"
|
||||
#include "singlePhaseTransportModel.H"
|
||||
#include "turbulenceModel.H"
|
||||
#include "dynamicFvMesh.H"
|
||||
#include "pimpleControl.H"
|
||||
#include "fvIOoptionList.H"
|
||||
|
||||
@ -44,15 +44,21 @@ Description
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
#include "setRootCase.H"
|
||||
|
||||
#include "createTime.H"
|
||||
#include "createDynamicFvMesh.H"
|
||||
#include "initContinuityErrs.H"
|
||||
|
||||
pimpleControl pimple(mesh);
|
||||
|
||||
#include "createFields.H"
|
||||
#include "createFvOptions.H"
|
||||
#include "readTimeControls.H"
|
||||
#include "createPcorrTypes.H"
|
||||
#include "CourantNo.H"
|
||||
#include "setInitialDeltaT.H"
|
||||
|
||||
pimpleControl pimple(mesh);
|
||||
// Create old-time absolute flux for ddtPhiCorr
|
||||
surfaceScalarField phiAbs("phiAbs", phi);
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
@ -66,6 +72,9 @@ int main(int argc, char *argv[])
|
||||
// Make the fluxes absolute
|
||||
fvc::makeAbsolute(phi, U);
|
||||
|
||||
// Update absolute flux for ddtPhiCorr
|
||||
phiAbs = phi;
|
||||
|
||||
#include "setDeltaT.H"
|
||||
|
||||
runTime++;
|
||||
|
||||
@ -1,12 +1,7 @@
|
||||
#include "readTimeControls.H"
|
||||
|
||||
const dictionary& pimpleDict = pimple.dict();
|
||||
|
||||
const bool correctPhi =
|
||||
pimpleDict.lookupOrDefault("correctPhi", false);
|
||||
pimple.dict().lookupOrDefault("correctPhi", false);
|
||||
|
||||
const bool checkMeshCourantNo =
|
||||
pimpleDict.lookupOrDefault("checkMeshCourantNo", false);
|
||||
|
||||
const bool ddtPhiCorr =
|
||||
pimpleDict.lookupOrDefault("ddtPhiCorr", true);
|
||||
pimple.dict().lookupOrDefault("checkMeshCourantNo", false);
|
||||
|
||||
@ -137,9 +137,6 @@ int main(int argc, char *argv[])
|
||||
while (pimple.correct())
|
||||
{
|
||||
#include "pEqn.H"
|
||||
|
||||
// Make the fluxes relative to the mesh motion
|
||||
fvc::makeRelative(phi, U);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -38,7 +38,7 @@
|
||||
pcorrTypes
|
||||
);
|
||||
|
||||
dimensionedScalar rAUf("(1|A(U))", dimTime/rho.dimensions(), 1.0);
|
||||
dimensionedScalar Dp("Dp", dimTime/rho.dimensions(), 1.0);
|
||||
|
||||
adjustPhi(phi, U, pcorr);
|
||||
|
||||
@ -46,7 +46,7 @@
|
||||
{
|
||||
fvScalarMatrix pcorrEqn
|
||||
(
|
||||
fvm::laplacian(rAUf, pcorr) == fvc::div(phi) - divU
|
||||
fvm::laplacian(Dp, pcorr) == fvc::div(phi) - divU
|
||||
);
|
||||
|
||||
pcorrEqn.solve();
|
||||
|
||||
@ -50,11 +50,11 @@
|
||||
{
|
||||
phi -= pcorrEqn.flux();
|
||||
phiAbs = phi;
|
||||
phiAbs.oldTime() = phi;
|
||||
fvc::makeRelative(phi, U);
|
||||
phi.oldTime() = phi;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
phi.oldTime() = phi;
|
||||
|
||||
#include "continuityErrs.H"
|
||||
|
||||
@ -51,10 +51,11 @@ int main(int argc, char *argv[])
|
||||
#include "createDynamicFvMesh.H"
|
||||
#include "initContinuityErrs.H"
|
||||
#include "createFields.H"
|
||||
#include "readTimeControls.H"
|
||||
|
||||
pimpleControl pimple(mesh);
|
||||
|
||||
#include "readTimeControls.H"
|
||||
|
||||
surfaceScalarField phiAbs("phiAbs", phi);
|
||||
fvc::makeAbsolute(phiAbs, U);
|
||||
|
||||
|
||||
@ -63,7 +63,6 @@ int main(int argc, char *argv[])
|
||||
#include "CourantNo.H"
|
||||
#include "setInitialDeltaT.H"
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
Info<< "\nStarting time loop\n" << endl;
|
||||
|
||||
@ -5,7 +5,10 @@ EXE_INC = \
|
||||
-I$(LIB_SRC)/transportModels/interfaceProperties/lnInclude \
|
||||
-I$(LIB_SRC)/turbulenceModels/incompressible/turbulenceModel \
|
||||
-IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/fvOptions/lnInclude\
|
||||
-I$(LIB_SRC)/sampling/lnInclude
|
||||
|
||||
EXE_LIBS = \
|
||||
-ltwoPhaseMixture \
|
||||
@ -15,4 +18,7 @@ EXE_LIBS = \
|
||||
-lincompressibleTurbulenceModel \
|
||||
-lincompressibleRASModels \
|
||||
-lincompressibleLESModels \
|
||||
-lfiniteVolume
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
|
||||
@ -16,7 +16,7 @@
|
||||
)
|
||||
+ fvc::flux
|
||||
(
|
||||
-fvc::flux(-phir, scalar(1) - alpha1, alpharScheme),
|
||||
-fvc::flux(-phir, alpha2, alpharScheme),
|
||||
alpha1,
|
||||
alpharScheme
|
||||
)
|
||||
@ -76,6 +76,7 @@
|
||||
0
|
||||
);
|
||||
|
||||
alpha2 = 1.0 - alpha1;
|
||||
rhoPhi +=
|
||||
(runTime.deltaT()/totalDeltaT)
|
||||
*(phiAlpha*(rho1 - rho2) + phi*rho2);
|
||||
|
||||
@ -40,5 +40,5 @@ surfaceScalarField rhoPhi
|
||||
#include "alphaEqn.H"
|
||||
}
|
||||
|
||||
rho == alpha1*rho1 + (scalar(1) - alpha1)*rho2;
|
||||
rho == alpha1*rho1 + alpha2*rho2;
|
||||
}
|
||||
|
||||
@ -12,20 +12,6 @@
|
||||
mesh
|
||||
);
|
||||
|
||||
Info<< "Reading field alpha1\n" << endl;
|
||||
volScalarField alpha1
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"alpha1",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
Info<< "Reading field U\n" << endl;
|
||||
volVectorField U
|
||||
(
|
||||
@ -42,14 +28,19 @@
|
||||
|
||||
#include "createPhi.H"
|
||||
|
||||
|
||||
Info<< "Creating phaseChangeTwoPhaseMixture\n" << endl;
|
||||
autoPtr<phaseChangeTwoPhaseMixture> twoPhaseProperties =
|
||||
phaseChangeTwoPhaseMixture::New(U, phi);
|
||||
|
||||
volScalarField& alpha1(twoPhaseProperties->alpha1());
|
||||
volScalarField& alpha2(twoPhaseProperties->alpha2());
|
||||
|
||||
const dimensionedScalar& rho1 = twoPhaseProperties->rho1();
|
||||
const dimensionedScalar& rho2 = twoPhaseProperties->rho2();
|
||||
const dimensionedScalar& pSat = twoPhaseProperties->pSat();
|
||||
|
||||
|
||||
// Need to store rho for ddt(rho, U)
|
||||
volScalarField rho
|
||||
(
|
||||
@ -60,11 +51,12 @@
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT
|
||||
),
|
||||
alpha1*rho1 + (scalar(1) - alpha1)*rho2,
|
||||
alpha1*rho1 + alpha2*rho2,
|
||||
alpha1.boundaryField().types()
|
||||
);
|
||||
rho.oldTime();
|
||||
|
||||
|
||||
// Construct interface from alpha1 distribution
|
||||
interfaceProperties interface(alpha1, U, twoPhaseProperties());
|
||||
|
||||
@ -113,3 +105,5 @@
|
||||
);
|
||||
p_rgh = p - rho*gh;
|
||||
}
|
||||
|
||||
fv::IOoptionList fvOptions(mesh);
|
||||
|
||||
@ -2,7 +2,7 @@
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
|
||||
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
@ -47,6 +47,7 @@ Description
|
||||
#include "phaseChangeTwoPhaseMixture.H"
|
||||
#include "turbulenceModel.H"
|
||||
#include "pimpleControl.H"
|
||||
#include "fvIOoptionList.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
@ -80,14 +81,10 @@ int main(int argc, char *argv[])
|
||||
|
||||
Info<< "Time = " << runTime.timeName() << nl << endl;
|
||||
|
||||
twoPhaseProperties->correct();
|
||||
|
||||
#include "alphaEqnSubCycle.H"
|
||||
|
||||
if (pimple.nCorrPIMPLE() == 1)
|
||||
{
|
||||
interface.correct();
|
||||
}
|
||||
|
||||
turbulence->correct();
|
||||
interface.correct();
|
||||
|
||||
// --- Pressure-velocity PIMPLE corrector loop
|
||||
while (pimple.loop())
|
||||
@ -99,9 +96,12 @@ int main(int argc, char *argv[])
|
||||
{
|
||||
#include "pEqn.H"
|
||||
}
|
||||
}
|
||||
|
||||
twoPhaseProperties->correct();
|
||||
if (pimple.turbCorr())
|
||||
{
|
||||
turbulence->correct();
|
||||
}
|
||||
}
|
||||
|
||||
runTime.write();
|
||||
|
||||
|
||||
@ -46,6 +46,7 @@
|
||||
|
||||
U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rAUf);
|
||||
U.correctBoundaryConditions();
|
||||
fvOptions.correct(U);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
@ -0,0 +1,44 @@
|
||||
/*--------------------------------*- 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 volScalarField;
|
||||
object T;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [0 0 0 1 0 0 0];
|
||||
|
||||
internalField uniform $:outerInlet.T;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
innerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:innerInlet.T;
|
||||
}
|
||||
|
||||
outerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:outerInlet.T;
|
||||
}
|
||||
|
||||
outlet { $:outlet.T; }
|
||||
staticWalls { $:wall.T; }
|
||||
movingWalls { $staticWalls; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,45 @@
|
||||
/*--------------------------------*- 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 volVectorField;
|
||||
location "0";
|
||||
object U;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [0 1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform (0 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
innerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:innerInlet.U;
|
||||
}
|
||||
|
||||
outerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:outerInlet.U;
|
||||
}
|
||||
|
||||
outlet { $:outlet.U; }
|
||||
staticWalls { $:wall.U; }
|
||||
movingWalls { $:movingWall.U; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,35 @@
|
||||
/*--------------------------------*- 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 volScalarField;
|
||||
location "0";
|
||||
object alphat;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [1 -1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform 0;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet { $:calculated.alphat; }
|
||||
outlet { $:calculated.alphat; }
|
||||
movingWalls { $:wall.alphat; }
|
||||
staticWalls { $movingWalls; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,45 @@
|
||||
/*--------------------------------*- 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 volScalarField;
|
||||
location "0";
|
||||
object epsilon;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [0 2 -3 0 0 0 0];
|
||||
|
||||
internalField uniform $:innerInlet.epsilon;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
innerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:innerInlet.epsilon;
|
||||
}
|
||||
|
||||
outerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:outerInlet.epsilon;
|
||||
}
|
||||
|
||||
outlet { $:outlet.epsilon; }
|
||||
staticWalls { $:wall.epsilon; }
|
||||
movingWalls { $staticWalls; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,45 @@
|
||||
/*--------------------------------*- 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 volScalarField;
|
||||
location "0";
|
||||
object k;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [0 2 -2 0 0 0 0];
|
||||
|
||||
internalField uniform $:innerInlet.k;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
innerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:innerInlet.k;
|
||||
}
|
||||
|
||||
outerInlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform $:outerInlet.k;
|
||||
}
|
||||
|
||||
outlet { $:outlet.k; }
|
||||
staticWalls { $:wall.k; }
|
||||
movingWalls { $staticWalls; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,34 @@
|
||||
/*--------------------------------*- 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 volScalarField;
|
||||
location "0";
|
||||
object mut;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [1 -1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform 0;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet { $:calculated.mut; }
|
||||
outlet { $:calculated.mut; }
|
||||
movingWalls { $:wall.mut; }
|
||||
staticWalls { $movingWalls; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,34 @@
|
||||
/*--------------------------------*- 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 volScalarField;
|
||||
location "0";
|
||||
object p;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dimensions [1 -1 -2 0 0 0 0];
|
||||
|
||||
internalField uniform 1e5;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet { $:inlet.p; }
|
||||
outlet { $:outlet.p; }
|
||||
staticWalls { $:wall.p; }
|
||||
movingWalls { $staticWalls; }
|
||||
|
||||
#include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
12
tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/Allclean
Executable file
12
tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/Allclean
Executable file
@ -0,0 +1,12 @@
|
||||
#!/bin/sh
|
||||
cd ${0%/*} || exit 1 # run from this directory
|
||||
|
||||
# Source tutorial clean functions
|
||||
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
|
||||
|
||||
rm -rf 0 > /dev/null 2>&1
|
||||
rm -rf constant/extendedFeatureEdgeMesh > /dev/null 2>&1
|
||||
|
||||
cleanCase
|
||||
|
||||
# ----------------------------------------------------------------- end-of-file
|
||||
25
tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/Allrun
Executable file
25
tutorials/compressible/rhoPimpleDyMFoam/annularThermalMixer/Allrun
Executable file
@ -0,0 +1,25 @@
|
||||
#!/bin/sh
|
||||
cd ${0%/*} || exit 1 # run from this directory
|
||||
|
||||
# Source tutorial run functions
|
||||
. $WM_PROJECT_DIR/bin/tools/RunFunctions
|
||||
|
||||
runApplication surfaceFeatureExtract
|
||||
|
||||
runApplication blockMesh
|
||||
runApplication snappyHexMesh -overwrite
|
||||
sed -i 's/_solid//g' constant/polyMesh/boundary
|
||||
|
||||
if [ -d 0 ] ; then
|
||||
rm -rf 0
|
||||
fi
|
||||
|
||||
runApplication createPatch -overwrite
|
||||
runApplication createBaffles -overwrite
|
||||
runApplication mergeOrSplitBaffles -split -overwrite
|
||||
|
||||
cp -r 0.org 0
|
||||
|
||||
runApplication $(getApplication)
|
||||
|
||||
# ----------------------------------------------------------------- end-of-file
|
||||
@ -0,0 +1,24 @@
|
||||
/*--------------------------------*- 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 "constant";
|
||||
object RASProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
RASModel kEpsilon;
|
||||
|
||||
turbulence on;
|
||||
|
||||
printCoeffs on;
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,135 @@
|
||||
/*--------------------------------*- 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 "";
|
||||
object caseProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
calculated
|
||||
{
|
||||
alphat
|
||||
{
|
||||
type calculated;
|
||||
value uniform 0;
|
||||
}
|
||||
mut
|
||||
{
|
||||
type calculated;
|
||||
value uniform 0;
|
||||
}
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
p
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
p
|
||||
{
|
||||
type totalPressure;
|
||||
value uniform 1e5;
|
||||
p0 uniform 1e5;
|
||||
U U;
|
||||
phi phi;
|
||||
rho rho;
|
||||
psi none;
|
||||
gamma 1.4;
|
||||
}
|
||||
|
||||
U
|
||||
{
|
||||
type pressureInletOutletVelocity;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
|
||||
T
|
||||
{
|
||||
type inletOutlet;
|
||||
inletValue uniform $:outerInlet.T;
|
||||
value $inletValue;
|
||||
}
|
||||
|
||||
k
|
||||
{
|
||||
type inletOutlet;
|
||||
inletValue uniform $:innerInlet.k;
|
||||
value $inletValue;
|
||||
}
|
||||
|
||||
epsilon
|
||||
{
|
||||
type inletOutlet;
|
||||
inletValue uniform $:innerInlet.epsilon;
|
||||
value $inletValue;
|
||||
}
|
||||
}
|
||||
|
||||
wall
|
||||
{
|
||||
p
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
U
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
|
||||
T
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
k
|
||||
{
|
||||
type compressible::kqRWallFunction;
|
||||
value uniform $:innerInlet.k;
|
||||
}
|
||||
|
||||
epsilon
|
||||
{
|
||||
type compressible::epsilonWallFunction;
|
||||
value uniform $:innerInlet.epsilon;
|
||||
}
|
||||
|
||||
mut
|
||||
{
|
||||
type mutkWallFunction;
|
||||
value uniform 0;
|
||||
}
|
||||
|
||||
alphat
|
||||
{
|
||||
type compressible::alphatWallFunction;
|
||||
Prt 0.85;
|
||||
value uniform 0;
|
||||
}
|
||||
}
|
||||
|
||||
movingWall
|
||||
{
|
||||
U
|
||||
{
|
||||
type movingWallVelocity;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,41 @@
|
||||
/*--------------------------------*- 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 "";
|
||||
object caseSummary;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
innerInlet
|
||||
{
|
||||
U (0 0 0.2);
|
||||
epsilon 5.70e-5;
|
||||
k 2.40e-5;
|
||||
T 233;
|
||||
}
|
||||
|
||||
outerInlet
|
||||
{
|
||||
U (0 0 0.1);
|
||||
epsilon 3.98e-5;
|
||||
k 6.00e-6;
|
||||
T 293;
|
||||
}
|
||||
|
||||
meshMotionProperties
|
||||
{
|
||||
radialVelocity (0 0 1440); // deg/s
|
||||
}
|
||||
|
||||
#include "${FOAM_CASE}/constant/boundaryConditions"
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,34 @@
|
||||
/*--------------------------------*- 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 "constant";
|
||||
object dynamicMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "${FOAM_CASE}/constant/caseSettings"
|
||||
|
||||
dynamicFvMesh solidBodyMotionFvMesh;
|
||||
|
||||
solidBodyMotionFvMeshCoeffs
|
||||
{
|
||||
cellZone rotatingZone;
|
||||
|
||||
solidBodyMotionFunction rotatingMotion;
|
||||
rotatingMotionCoeffs
|
||||
{
|
||||
CofG (0 0 0);
|
||||
radialVelocity $:meshMotionProperties.radialVelocity;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,51 @@
|
||||
/*--------------------------------*- 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 0.1;
|
||||
|
||||
vertices
|
||||
(
|
||||
(-1.1 -1.1 -0.1)
|
||||
( 1.1 -1.1 -0.1)
|
||||
( 1.1 1.1 -0.1)
|
||||
(-1.1 1.1 -0.1)
|
||||
|
||||
(-1.1 -1.1 2.1)
|
||||
( 1.1 -1.1 2.1)
|
||||
( 1.1 1.1 2.1)
|
||||
(-1.1 1.1 2.1)
|
||||
);
|
||||
|
||||
blocks
|
||||
(
|
||||
hex (0 1 2 3 4 5 6 7)
|
||||
(21 21 21)
|
||||
simpleGrading (1 1 1)
|
||||
);
|
||||
|
||||
edges
|
||||
(
|
||||
);
|
||||
|
||||
boundary
|
||||
(
|
||||
);
|
||||
|
||||
mergePatchPairs
|
||||
(
|
||||
);
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,112 @@
|
||||
/*--------------------------------*- 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;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
12
|
||||
(
|
||||
innerInlet
|
||||
{
|
||||
type patch;
|
||||
inGroups 1(inlet);
|
||||
nFaces 544;
|
||||
startFace 244948;
|
||||
}
|
||||
outerInlet
|
||||
{
|
||||
type patch;
|
||||
inGroups 1(inlet);
|
||||
nFaces 1396;
|
||||
startFace 245492;
|
||||
}
|
||||
innerOutlet
|
||||
{
|
||||
type patch;
|
||||
inGroups 1(outlet);
|
||||
nFaces 544;
|
||||
startFace 246888;
|
||||
}
|
||||
outerOutlet
|
||||
{
|
||||
type patch;
|
||||
inGroups 1(outlet);
|
||||
nFaces 1396;
|
||||
startFace 247432;
|
||||
}
|
||||
rotorBlades
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(movingWalls);
|
||||
nFaces 540;
|
||||
startFace 248828;
|
||||
}
|
||||
rotorBlades_slave
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(movingWalls);
|
||||
nFaces 540;
|
||||
startFace 249368;
|
||||
}
|
||||
shaft
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(movingWalls);
|
||||
nFaces 1052;
|
||||
startFace 249908;
|
||||
}
|
||||
statorBlades
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(staticWalls);
|
||||
nFaces 2128;
|
||||
startFace 250960;
|
||||
}
|
||||
statorBlades_slave
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(staticWalls);
|
||||
nFaces 2128;
|
||||
startFace 253088;
|
||||
}
|
||||
wall
|
||||
{
|
||||
type wall;
|
||||
inGroups 1(staticWalls);
|
||||
nFaces 5808;
|
||||
startFace 255216;
|
||||
}
|
||||
AMI1
|
||||
{
|
||||
type cyclicAMI;
|
||||
inGroups 1(cyclicAMI);
|
||||
nFaces 10944;
|
||||
startFace 261024;
|
||||
matchTolerance 0.0001;
|
||||
transform noOrdering;
|
||||
neighbourPatch AMI2;
|
||||
}
|
||||
AMI2
|
||||
{
|
||||
type cyclicAMI;
|
||||
inGroups 1(cyclicAMI);
|
||||
nFaces 10944;
|
||||
startFace 271968;
|
||||
matchTolerance 0.0001;
|
||||
transform noOrdering;
|
||||
neighbourPatch AMI1;
|
||||
}
|
||||
)
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,49 @@
|
||||
/*--------------------------------*- 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 "constant";
|
||||
object thermophysicalProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
thermoType
|
||||
{
|
||||
type hePsiThermo;
|
||||
mixture pureMixture;
|
||||
transport sutherland;
|
||||
thermo hConst;
|
||||
equationOfState perfectGas;
|
||||
specie specie;
|
||||
energy sensibleEnthalpy;
|
||||
}
|
||||
|
||||
mixture
|
||||
{
|
||||
specie
|
||||
{
|
||||
nMoles 1;
|
||||
molWeight 28.9;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Cp 1007;
|
||||
Hf 0;
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.4792e-06;
|
||||
Ts 116;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,22 @@
|
||||
/*--------------------------------*- 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 "constant";
|
||||
object transportProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
transportModel Newtonian;
|
||||
|
||||
nu nu [ 0 2 -1 0 0 0 0 ] 1.0e-6;
|
||||
|
||||
// ************************************************************************* //
|
||||
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
@ -0,0 +1,243 @@
|
||||
/*--------------------------------*- 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 featureEdgeMesh;
|
||||
location "constant/triSurface";
|
||||
object rotorBlades.eMesh;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// points:
|
||||
|
||||
104
|
||||
(
|
||||
(-6.88841e-18 0.0375 0.055)
|
||||
(-6.88841e-18 0.0375 0.06)
|
||||
(0.0375 9.18455e-18 0.065)
|
||||
(0.0375 9.18455e-18 0.06)
|
||||
(-0.0125 -1.53076e-18 0.06)
|
||||
(-0.0125 -1.53076e-18 0.055)
|
||||
(-4.13305e-18 0.0225 0.04)
|
||||
(-5.0515e-18 0.0275 0.04)
|
||||
(0.0225 5.51073e-18 0.08)
|
||||
(0.0275 6.73533e-18 0.08)
|
||||
(2.29614e-18 -0.0375 0.06)
|
||||
(2.29614e-18 -0.0375 0.065)
|
||||
(-2.29614e-18 0.0125 0.06)
|
||||
(-2.29614e-18 0.0125 0.055)
|
||||
(0.0125 3.06152e-18 0.06)
|
||||
(0.0125 3.06152e-18 0.065)
|
||||
(0.0325 7.95994e-18 0.04)
|
||||
(0.0275 6.73533e-18 0.04)
|
||||
(1.68383e-18 -0.0275 0.08)
|
||||
(1.37768e-18 -0.0225 0.08)
|
||||
(-0.0375 -4.59227e-18 0.06)
|
||||
(-0.0375 -4.59227e-18 0.065)
|
||||
(7.65379e-19 -0.0125 0.065)
|
||||
(7.65379e-19 -0.0125 0.06)
|
||||
(1.68383e-18 -0.0275 0.04)
|
||||
(1.98999e-18 -0.0325 0.04)
|
||||
(-0.0275 -3.36767e-18 0.08)
|
||||
(-0.0225 -2.75536e-18 0.08)
|
||||
(-6.88841e-18 0.0375 0.065)
|
||||
(0.0375 9.18455e-18 0.07)
|
||||
(-0.0125 -1.53076e-18 0.065)
|
||||
(-0.0275 -3.36767e-18 0.04)
|
||||
(-0.0325 -3.97997e-18 0.04)
|
||||
(-5.0515e-18 0.0275 0.08)
|
||||
(-4.13305e-18 0.0225 0.08)
|
||||
(2.29614e-18 -0.0375 0.07)
|
||||
(-2.29614e-18 0.0125 0.065)
|
||||
(0.0125 3.06152e-18 0.07)
|
||||
(-5.96996e-18 0.0325 0.04)
|
||||
(0.0325 7.95994e-18 0.08)
|
||||
(-0.0375 -4.59227e-18 0.07)
|
||||
(7.65379e-19 -0.0125 0.07)
|
||||
(0.0375 9.18455e-18 0.04)
|
||||
(0.0375 9.18455e-18 0.045)
|
||||
(1.98999e-18 -0.0325 0.08)
|
||||
(-6.88841e-18 0.0375 0.07)
|
||||
(0.0375 9.18455e-18 0.075)
|
||||
(-0.0125 -1.53076e-18 0.07)
|
||||
(2.29614e-18 -0.0375 0.04)
|
||||
(2.29614e-18 -0.0375 0.045)
|
||||
(0.0125 3.06152e-18 0.04)
|
||||
(0.0125 3.06152e-18 0.045)
|
||||
(0.0175 4.28612e-18 0.04)
|
||||
(-0.0325 -3.97997e-18 0.08)
|
||||
(2.29614e-18 -0.0375 0.075)
|
||||
(-2.29614e-18 0.0125 0.07)
|
||||
(0.0125 3.06152e-18 0.075)
|
||||
(-0.0375 -4.59227e-18 0.04)
|
||||
(-0.0375 -4.59227e-18 0.045)
|
||||
(7.65379e-19 -0.0125 0.04)
|
||||
(1.07153e-18 -0.0175 0.04)
|
||||
(7.65379e-19 -0.0125 0.045)
|
||||
(-5.96996e-18 0.0325 0.08)
|
||||
(-0.0375 -4.59227e-18 0.075)
|
||||
(7.65379e-19 -0.0125 0.075)
|
||||
(-6.88841e-18 0.0375 0.04)
|
||||
(-6.88841e-18 0.0375 0.045)
|
||||
(0.0375 9.18455e-18 0.05)
|
||||
(-0.0125 -1.53076e-18 0.04)
|
||||
(-0.0175 -2.14306e-18 0.04)
|
||||
(-0.0125 -1.53076e-18 0.045)
|
||||
(-6.88841e-18 0.0375 0.075)
|
||||
(0.0375 9.18455e-18 0.08)
|
||||
(-0.0125 -1.53076e-18 0.075)
|
||||
(2.29614e-18 -0.0375 0.05)
|
||||
(-2.29614e-18 0.0125 0.04)
|
||||
(-3.21459e-18 0.0175 0.04)
|
||||
(-2.29614e-18 0.0125 0.045)
|
||||
(0.0125 3.06152e-18 0.05)
|
||||
(2.29614e-18 -0.0375 0.08)
|
||||
(-2.29614e-18 0.0125 0.075)
|
||||
(0.0125 3.06152e-18 0.08)
|
||||
(0.0175 4.28612e-18 0.08)
|
||||
(-0.0375 -4.59227e-18 0.05)
|
||||
(7.65379e-19 -0.0125 0.05)
|
||||
(0.0225 5.51073e-18 0.04)
|
||||
(-0.0375 -4.59227e-18 0.08)
|
||||
(7.65379e-19 -0.0125 0.08)
|
||||
(1.07153e-18 -0.0175 0.08)
|
||||
(-6.88841e-18 0.0375 0.05)
|
||||
(0.0375 9.18455e-18 0.055)
|
||||
(-0.0125 -1.53076e-18 0.05)
|
||||
(1.37768e-18 -0.0225 0.04)
|
||||
(-6.88841e-18 0.0375 0.08)
|
||||
(-0.0125 -1.53076e-18 0.08)
|
||||
(-0.0175 -2.14306e-18 0.08)
|
||||
(2.29614e-18 -0.0375 0.055)
|
||||
(-2.29614e-18 0.0125 0.05)
|
||||
(0.0125 3.06152e-18 0.055)
|
||||
(-0.0225 -2.75536e-18 0.04)
|
||||
(-2.29614e-18 0.0125 0.08)
|
||||
(-3.21459e-18 0.0175 0.08)
|
||||
(-0.0375 -4.59227e-18 0.055)
|
||||
(7.65379e-19 -0.0125 0.055)
|
||||
)
|
||||
|
||||
|
||||
// edges:
|
||||
|
||||
104
|
||||
(
|
||||
(0 1)
|
||||
(2 3)
|
||||
(4 5)
|
||||
(6 7)
|
||||
(8 9)
|
||||
(10 11)
|
||||
(12 13)
|
||||
(14 15)
|
||||
(16 17)
|
||||
(18 19)
|
||||
(20 21)
|
||||
(22 23)
|
||||
(24 25)
|
||||
(26 27)
|
||||
(1 28)
|
||||
(29 2)
|
||||
(30 4)
|
||||
(31 32)
|
||||
(33 34)
|
||||
(11 35)
|
||||
(36 12)
|
||||
(15 37)
|
||||
(7 38)
|
||||
(9 39)
|
||||
(21 40)
|
||||
(41 22)
|
||||
(42 16)
|
||||
(43 42)
|
||||
(44 18)
|
||||
(28 45)
|
||||
(46 29)
|
||||
(47 30)
|
||||
(25 48)
|
||||
(48 49)
|
||||
(50 51)
|
||||
(52 50)
|
||||
(53 26)
|
||||
(35 54)
|
||||
(55 36)
|
||||
(37 56)
|
||||
(32 57)
|
||||
(57 58)
|
||||
(59 60)
|
||||
(61 59)
|
||||
(62 33)
|
||||
(40 63)
|
||||
(64 41)
|
||||
(38 65)
|
||||
(65 66)
|
||||
(67 43)
|
||||
(68 69)
|
||||
(70 68)
|
||||
(45 71)
|
||||
(72 46)
|
||||
(39 72)
|
||||
(73 47)
|
||||
(49 74)
|
||||
(75 76)
|
||||
(77 75)
|
||||
(51 78)
|
||||
(79 44)
|
||||
(54 79)
|
||||
(80 55)
|
||||
(56 81)
|
||||
(81 82)
|
||||
(58 83)
|
||||
(84 61)
|
||||
(85 52)
|
||||
(86 53)
|
||||
(63 86)
|
||||
(87 64)
|
||||
(88 87)
|
||||
(66 89)
|
||||
(90 67)
|
||||
(91 70)
|
||||
(60 92)
|
||||
(93 62)
|
||||
(71 93)
|
||||
(94 73)
|
||||
(95 94)
|
||||
(74 96)
|
||||
(97 77)
|
||||
(78 98)
|
||||
(69 99)
|
||||
(100 80)
|
||||
(101 100)
|
||||
(83 102)
|
||||
(103 84)
|
||||
(76 6)
|
||||
(82 8)
|
||||
(89 0)
|
||||
(3 90)
|
||||
(5 91)
|
||||
(17 85)
|
||||
(19 88)
|
||||
(96 10)
|
||||
(13 97)
|
||||
(98 14)
|
||||
(92 24)
|
||||
(27 95)
|
||||
(102 20)
|
||||
(23 103)
|
||||
(99 31)
|
||||
(34 101)
|
||||
)
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,544 @@
|
||||
# Wavefront OBJ file
|
||||
# Regions:
|
||||
# 0 solid
|
||||
#
|
||||
# points : 216
|
||||
# triangles : 320
|
||||
#
|
||||
v -5.96996e-18 0.0325 0.055
|
||||
v -6.88841e-18 0.0375 0.055
|
||||
v -5.96996e-18 0.0325 0.06
|
||||
v -6.88841e-18 0.0375 0.06
|
||||
v 0.0325 7.95994e-18 0.06
|
||||
v 0.0325 7.95994e-18 0.065
|
||||
v 0.0375 9.18455e-18 0.06
|
||||
v 0.0375 9.18455e-18 0.065
|
||||
v -0.0125 -1.53076e-18 0.055
|
||||
v -0.0175 -2.14306e-18 0.055
|
||||
v -0.0125 -1.53076e-18 0.06
|
||||
v -0.0175 -2.14306e-18 0.06
|
||||
v 1.07153e-18 -0.0175 0.05
|
||||
v 1.37768e-18 -0.0225 0.05
|
||||
v 1.07153e-18 -0.0175 0.055
|
||||
v 1.37768e-18 -0.0225 0.055
|
||||
v -4.13305e-18 0.0225 0.04
|
||||
v -5.0515e-18 0.0275 0.04
|
||||
v -4.13305e-18 0.0225 0.045
|
||||
v -5.0515e-18 0.0275 0.045
|
||||
v 0.0225 5.51073e-18 0.045
|
||||
v 0.0225 5.51073e-18 0.05
|
||||
v 0.0275 6.73533e-18 0.045
|
||||
v 0.0275 6.73533e-18 0.05
|
||||
v -4.13305e-18 0.0225 0.07
|
||||
v -5.0515e-18 0.0275 0.07
|
||||
v -4.13305e-18 0.0225 0.075
|
||||
v -5.0515e-18 0.0275 0.075
|
||||
v 0.0225 5.51073e-18 0.075
|
||||
v 0.0225 5.51073e-18 0.08
|
||||
v 0.0275 6.73533e-18 0.075
|
||||
v 0.0275 6.73533e-18 0.08
|
||||
v -0.0275 -3.36767e-18 0.065
|
||||
v -0.0325 -3.97997e-18 0.065
|
||||
v -0.0275 -3.36767e-18 0.07
|
||||
v -0.0325 -3.97997e-18 0.07
|
||||
v 1.98999e-18 -0.0325 0.06
|
||||
v 2.29614e-18 -0.0375 0.06
|
||||
v 1.98999e-18 -0.0325 0.065
|
||||
v 2.29614e-18 -0.0375 0.065
|
||||
v -2.29614e-18 0.0125 0.055
|
||||
v -3.21459e-18 0.0175 0.055
|
||||
v -2.29614e-18 0.0125 0.06
|
||||
v -3.21459e-18 0.0175 0.06
|
||||
v 0.0125 3.06152e-18 0.06
|
||||
v 0.0125 3.06152e-18 0.065
|
||||
v 0.0175 4.28612e-18 0.06
|
||||
v 0.0175 4.28612e-18 0.065
|
||||
v -0.0175 -2.14306e-18 0.05
|
||||
v -0.0225 -2.75536e-18 0.05
|
||||
v -0.0225 -2.75536e-18 0.055
|
||||
v 1.37768e-18 -0.0225 0.045
|
||||
v 1.68383e-18 -0.0275 0.045
|
||||
v 1.68383e-18 -0.0275 0.05
|
||||
v 0.0275 6.73533e-18 0.04
|
||||
v 0.0325 7.95994e-18 0.04
|
||||
v 0.0325 7.95994e-18 0.045
|
||||
v 1.37768e-18 -0.0225 0.075
|
||||
v 1.68383e-18 -0.0275 0.075
|
||||
v 1.37768e-18 -0.0225 0.08
|
||||
v 1.68383e-18 -0.0275 0.08
|
||||
v -5.0515e-18 0.0275 0.065
|
||||
v -5.96996e-18 0.0325 0.065
|
||||
v -5.96996e-18 0.0325 0.07
|
||||
v 0.0275 6.73533e-18 0.07
|
||||
v 0.0325 7.95994e-18 0.07
|
||||
v 0.0325 7.95994e-18 0.075
|
||||
v -0.0325 -3.97997e-18 0.06
|
||||
v -0.0375 -4.59227e-18 0.06
|
||||
v -0.0375 -4.59227e-18 0.065
|
||||
v 7.65379e-19 -0.0125 0.06
|
||||
v 1.07153e-18 -0.0175 0.06
|
||||
v 7.65379e-19 -0.0125 0.065
|
||||
v 1.07153e-18 -0.0175 0.065
|
||||
v -3.21459e-18 0.0175 0.05
|
||||
v -4.13305e-18 0.0225 0.05
|
||||
v -4.13305e-18 0.0225 0.055
|
||||
v 0.0175 4.28612e-18 0.055
|
||||
v 0.0225 5.51073e-18 0.055
|
||||
v 0.0225 5.51073e-18 0.06
|
||||
v -0.0225 -2.75536e-18 0.045
|
||||
v -0.0275 -3.36767e-18 0.045
|
||||
v -0.0275 -3.36767e-18 0.05
|
||||
v 1.68383e-18 -0.0275 0.04
|
||||
v 1.98999e-18 -0.0325 0.04
|
||||
v 1.98999e-18 -0.0325 0.045
|
||||
v -0.0225 -2.75536e-18 0.075
|
||||
v -0.0275 -3.36767e-18 0.075
|
||||
v -0.0225 -2.75536e-18 0.08
|
||||
v -0.0275 -3.36767e-18 0.08
|
||||
v 1.68383e-18 -0.0275 0.07
|
||||
v 1.98999e-18 -0.0325 0.07
|
||||
v 1.98999e-18 -0.0325 0.075
|
||||
v -6.88841e-18 0.0375 0.065
|
||||
v 0.0375 9.18455e-18 0.07
|
||||
v -0.0125 -1.53076e-18 0.065
|
||||
v -0.0175 -2.14306e-18 0.065
|
||||
v 1.37768e-18 -0.0225 0.06
|
||||
v -5.0515e-18 0.0275 0.05
|
||||
v 0.0275 6.73533e-18 0.055
|
||||
v -0.0275 -3.36767e-18 0.04
|
||||
v -0.0325 -3.97997e-18 0.04
|
||||
v -0.0325 -3.97997e-18 0.045
|
||||
v -4.13305e-18 0.0225 0.08
|
||||
v -5.0515e-18 0.0275 0.08
|
||||
v -0.0325 -3.97997e-18 0.075
|
||||
v 2.29614e-18 -0.0375 0.07
|
||||
v -2.29614e-18 0.0125 0.065
|
||||
v -3.21459e-18 0.0175 0.065
|
||||
v 0.0125 3.06152e-18 0.07
|
||||
v 0.0175 4.28612e-18 0.07
|
||||
v -0.0225 -2.75536e-18 0.06
|
||||
v 1.68383e-18 -0.0275 0.055
|
||||
v -5.96996e-18 0.0325 0.04
|
||||
v -5.96996e-18 0.0325 0.045
|
||||
v 0.0325 7.95994e-18 0.05
|
||||
v -5.96996e-18 0.0325 0.075
|
||||
v 0.0325 7.95994e-18 0.08
|
||||
v -0.0375 -4.59227e-18 0.07
|
||||
v 7.65379e-19 -0.0125 0.07
|
||||
v 1.07153e-18 -0.0175 0.07
|
||||
v -4.13305e-18 0.0225 0.06
|
||||
v 0.0225 5.51073e-18 0.065
|
||||
v -0.0275 -3.36767e-18 0.055
|
||||
v 1.98999e-18 -0.0325 0.05
|
||||
v 0.0375 9.18455e-18 0.04
|
||||
v 0.0375 9.18455e-18 0.045
|
||||
v 1.98999e-18 -0.0325 0.08
|
||||
v -6.88841e-18 0.0375 0.07
|
||||
v 0.0375 9.18455e-18 0.075
|
||||
v -0.0125 -1.53076e-18 0.07
|
||||
v -0.0175 -2.14306e-18 0.07
|
||||
v 1.37768e-18 -0.0225 0.065
|
||||
v -5.0515e-18 0.0275 0.055
|
||||
v 0.0275 6.73533e-18 0.06
|
||||
v -0.0325 -3.97997e-18 0.05
|
||||
v 2.29614e-18 -0.0375 0.04
|
||||
v 2.29614e-18 -0.0375 0.045
|
||||
v 0.0125 3.06152e-18 0.04
|
||||
v 0.0125 3.06152e-18 0.045
|
||||
v 0.0175 4.28612e-18 0.04
|
||||
v 0.0175 4.28612e-18 0.045
|
||||
v -0.0325 -3.97997e-18 0.08
|
||||
v 2.29614e-18 -0.0375 0.075
|
||||
v -2.29614e-18 0.0125 0.07
|
||||
v -3.21459e-18 0.0175 0.07
|
||||
v 0.0125 3.06152e-18 0.075
|
||||
v 0.0175 4.28612e-18 0.075
|
||||
v -0.0225 -2.75536e-18 0.065
|
||||
v 1.68383e-18 -0.0275 0.06
|
||||
v -5.96996e-18 0.0325 0.05
|
||||
v 0.0325 7.95994e-18 0.055
|
||||
v -0.0375 -4.59227e-18 0.04
|
||||
v -0.0375 -4.59227e-18 0.045
|
||||
v 7.65379e-19 -0.0125 0.04
|
||||
v 1.07153e-18 -0.0175 0.04
|
||||
v 7.65379e-19 -0.0125 0.045
|
||||
v 1.07153e-18 -0.0175 0.045
|
||||
v -5.96996e-18 0.0325 0.08
|
||||
v -0.0375 -4.59227e-18 0.075
|
||||
v 7.65379e-19 -0.0125 0.075
|
||||
v 1.07153e-18 -0.0175 0.075
|
||||
v -4.13305e-18 0.0225 0.065
|
||||
v 0.0225 5.51073e-18 0.07
|
||||
v -0.0275 -3.36767e-18 0.06
|
||||
v 1.98999e-18 -0.0325 0.055
|
||||
v -6.88841e-18 0.0375 0.04
|
||||
v -6.88841e-18 0.0375 0.045
|
||||
v 0.0375 9.18455e-18 0.05
|
||||
v -0.0125 -1.53076e-18 0.04
|
||||
v -0.0175 -2.14306e-18 0.04
|
||||
v -0.0125 -1.53076e-18 0.045
|
||||
v -0.0175 -2.14306e-18 0.045
|
||||
v -6.88841e-18 0.0375 0.075
|
||||
v 0.0375 9.18455e-18 0.08
|
||||
v -0.0125 -1.53076e-18 0.075
|
||||
v -0.0175 -2.14306e-18 0.075
|
||||
v 1.37768e-18 -0.0225 0.07
|
||||
v -5.0515e-18 0.0275 0.06
|
||||
v 0.0275 6.73533e-18 0.065
|
||||
v -0.0325 -3.97997e-18 0.055
|
||||
v 2.29614e-18 -0.0375 0.05
|
||||
v -2.29614e-18 0.0125 0.04
|
||||
v -3.21459e-18 0.0175 0.04
|
||||
v -2.29614e-18 0.0125 0.045
|
||||
v -3.21459e-18 0.0175 0.045
|
||||
v 0.0125 3.06152e-18 0.05
|
||||
v 0.0175 4.28612e-18 0.05
|
||||
v 2.29614e-18 -0.0375 0.08
|
||||
v -2.29614e-18 0.0125 0.075
|
||||
v -3.21459e-18 0.0175 0.075
|
||||
v 0.0125 3.06152e-18 0.08
|
||||
v 0.0175 4.28612e-18 0.08
|
||||
v -0.0225 -2.75536e-18 0.07
|
||||
v 1.68383e-18 -0.0275 0.065
|
||||
v -0.0375 -4.59227e-18 0.05
|
||||
v 7.65379e-19 -0.0125 0.05
|
||||
v 0.0225 5.51073e-18 0.04
|
||||
v -0.0375 -4.59227e-18 0.08
|
||||
v 7.65379e-19 -0.0125 0.08
|
||||
v 1.07153e-18 -0.0175 0.08
|
||||
v -6.88841e-18 0.0375 0.05
|
||||
v 0.0375 9.18455e-18 0.055
|
||||
v -0.0125 -1.53076e-18 0.05
|
||||
v 1.37768e-18 -0.0225 0.04
|
||||
v -6.88841e-18 0.0375 0.08
|
||||
v -0.0125 -1.53076e-18 0.08
|
||||
v -0.0175 -2.14306e-18 0.08
|
||||
v 2.29614e-18 -0.0375 0.055
|
||||
v -2.29614e-18 0.0125 0.05
|
||||
v 0.0125 3.06152e-18 0.055
|
||||
v -0.0225 -2.75536e-18 0.04
|
||||
v -2.29614e-18 0.0125 0.08
|
||||
v -3.21459e-18 0.0175 0.08
|
||||
v -0.0375 -4.59227e-18 0.055
|
||||
v 7.65379e-19 -0.0125 0.055
|
||||
g solid
|
||||
f 1 2 3
|
||||
f 4 3 2
|
||||
f 5 6 7
|
||||
f 8 7 6
|
||||
f 9 10 11
|
||||
f 12 11 10
|
||||
f 13 14 15
|
||||
f 16 15 14
|
||||
f 17 18 19
|
||||
f 20 19 18
|
||||
f 21 22 23
|
||||
f 24 23 22
|
||||
f 25 26 27
|
||||
f 28 27 26
|
||||
f 29 30 31
|
||||
f 32 31 30
|
||||
f 33 34 35
|
||||
f 36 35 34
|
||||
f 37 38 39
|
||||
f 40 39 38
|
||||
f 41 42 43
|
||||
f 44 43 42
|
||||
f 45 46 47
|
||||
f 48 47 46
|
||||
f 49 50 10
|
||||
f 51 10 50
|
||||
f 52 53 14
|
||||
f 54 14 53
|
||||
f 55 23 56
|
||||
f 57 56 23
|
||||
f 58 59 60
|
||||
f 61 60 59
|
||||
f 62 63 26
|
||||
f 64 26 63
|
||||
f 65 31 66
|
||||
f 67 66 31
|
||||
f 68 69 34
|
||||
f 70 34 69
|
||||
f 71 72 73
|
||||
f 74 73 72
|
||||
f 75 76 42
|
||||
f 77 42 76
|
||||
f 78 47 79
|
||||
f 80 79 47
|
||||
f 81 82 50
|
||||
f 83 50 82
|
||||
f 84 85 53
|
||||
f 86 53 85
|
||||
f 87 88 89
|
||||
f 90 89 88
|
||||
f 91 92 59
|
||||
f 93 59 92
|
||||
f 3 4 63
|
||||
f 94 63 4
|
||||
f 6 66 8
|
||||
f 95 8 66
|
||||
f 11 12 96
|
||||
f 97 96 12
|
||||
f 15 16 72
|
||||
f 98 72 16
|
||||
f 19 20 76
|
||||
f 99 76 20
|
||||
f 22 79 24
|
||||
f 100 24 79
|
||||
f 101 102 82
|
||||
f 103 82 102
|
||||
f 27 28 104
|
||||
f 105 104 28
|
||||
f 35 36 88
|
||||
f 106 88 36
|
||||
f 39 40 92
|
||||
f 107 92 40
|
||||
f 43 44 108
|
||||
f 109 108 44
|
||||
f 46 110 48
|
||||
f 111 48 110
|
||||
f 10 51 12
|
||||
f 112 12 51
|
||||
f 14 54 16
|
||||
f 113 16 54
|
||||
f 18 114 20
|
||||
f 115 20 114
|
||||
f 23 24 57
|
||||
f 116 57 24
|
||||
f 26 64 28
|
||||
f 117 28 64
|
||||
f 31 32 67
|
||||
f 118 67 32
|
||||
f 34 70 36
|
||||
f 119 36 70
|
||||
f 73 74 120
|
||||
f 121 120 74
|
||||
f 42 77 44
|
||||
f 122 44 77
|
||||
f 47 48 80
|
||||
f 123 80 48
|
||||
f 50 83 51
|
||||
f 124 51 83
|
||||
f 53 86 54
|
||||
f 125 54 86
|
||||
f 56 57 126
|
||||
f 127 126 57
|
||||
f 59 93 61
|
||||
f 128 61 93
|
||||
f 63 94 64
|
||||
f 129 64 94
|
||||
f 66 67 95
|
||||
f 130 95 67
|
||||
f 96 97 131
|
||||
f 132 131 97
|
||||
f 72 98 74
|
||||
f 133 74 98
|
||||
f 76 99 77
|
||||
f 134 77 99
|
||||
f 79 80 100
|
||||
f 135 100 80
|
||||
f 82 103 83
|
||||
f 136 83 103
|
||||
f 85 137 86
|
||||
f 138 86 137
|
||||
f 139 140 141
|
||||
f 142 141 140
|
||||
f 88 106 90
|
||||
f 143 90 106
|
||||
f 92 107 93
|
||||
f 144 93 107
|
||||
f 108 109 145
|
||||
f 146 145 109
|
||||
f 110 147 111
|
||||
f 148 111 147
|
||||
f 12 112 97
|
||||
f 149 97 112
|
||||
f 16 113 98
|
||||
f 150 98 113
|
||||
f 20 115 99
|
||||
f 151 99 115
|
||||
f 24 100 116
|
||||
f 152 116 100
|
||||
f 102 153 103
|
||||
f 154 103 153
|
||||
f 155 156 157
|
||||
f 158 157 156
|
||||
f 28 117 105
|
||||
f 159 105 117
|
||||
f 36 119 106
|
||||
f 160 106 119
|
||||
f 120 121 161
|
||||
f 162 161 121
|
||||
f 44 122 109
|
||||
f 163 109 122
|
||||
f 48 111 123
|
||||
f 164 123 111
|
||||
f 51 124 112
|
||||
f 165 112 124
|
||||
f 54 125 113
|
||||
f 166 113 125
|
||||
f 114 167 115
|
||||
f 168 115 167
|
||||
f 57 116 127
|
||||
f 169 127 116
|
||||
f 170 171 172
|
||||
f 173 172 171
|
||||
f 64 129 117
|
||||
f 174 117 129
|
||||
f 67 118 130
|
||||
f 175 130 118
|
||||
f 131 132 176
|
||||
f 177 176 132
|
||||
f 74 133 121
|
||||
f 178 121 133
|
||||
f 77 134 122
|
||||
f 179 122 134
|
||||
f 80 123 135
|
||||
f 180 135 123
|
||||
f 83 136 124
|
||||
f 181 124 136
|
||||
f 86 138 125
|
||||
f 182 125 138
|
||||
f 183 184 185
|
||||
f 186 185 184
|
||||
f 140 187 142
|
||||
f 188 142 187
|
||||
f 93 144 128
|
||||
f 189 128 144
|
||||
f 145 146 190
|
||||
f 191 190 146
|
||||
f 147 192 148
|
||||
f 193 148 192
|
||||
f 97 149 132
|
||||
f 194 132 149
|
||||
f 98 150 133
|
||||
f 195 133 150
|
||||
f 99 151 134
|
||||
f 1 134 151
|
||||
f 100 135 152
|
||||
f 5 152 135
|
||||
f 103 154 136
|
||||
f 196 136 154
|
||||
f 157 158 197
|
||||
f 13 197 158
|
||||
f 141 142 198
|
||||
f 21 198 142
|
||||
f 106 160 143
|
||||
f 199 143 160
|
||||
f 161 162 200
|
||||
f 201 200 162
|
||||
f 109 163 146
|
||||
f 25 146 163
|
||||
f 111 148 164
|
||||
f 29 164 148
|
||||
f 112 165 149
|
||||
f 33 149 165
|
||||
f 113 166 150
|
||||
f 37 150 166
|
||||
f 115 168 151
|
||||
f 202 151 168
|
||||
f 116 152 169
|
||||
f 203 169 152
|
||||
f 172 173 204
|
||||
f 49 204 173
|
||||
f 156 205 158
|
||||
f 52 158 205
|
||||
f 117 174 159
|
||||
f 206 159 174
|
||||
f 176 177 207
|
||||
f 208 207 177
|
||||
f 121 178 162
|
||||
f 58 162 178
|
||||
f 122 179 163
|
||||
f 62 163 179
|
||||
f 123 164 180
|
||||
f 65 180 164
|
||||
f 124 181 165
|
||||
f 68 165 181
|
||||
f 125 182 166
|
||||
f 209 166 182
|
||||
f 185 186 210
|
||||
f 75 210 186
|
||||
f 187 211 188
|
||||
f 78 188 211
|
||||
f 171 212 173
|
||||
f 81 173 212
|
||||
f 190 191 213
|
||||
f 214 213 191
|
||||
f 132 194 177
|
||||
f 87 177 194
|
||||
f 133 195 178
|
||||
f 91 178 195
|
||||
f 134 1 179
|
||||
f 3 179 1
|
||||
f 135 180 5
|
||||
f 6 5 180
|
||||
f 136 196 181
|
||||
f 215 181 196
|
||||
f 197 13 216
|
||||
f 15 216 13
|
||||
f 184 17 186
|
||||
f 19 186 17
|
||||
f 142 188 21
|
||||
f 22 21 188
|
||||
f 146 25 191
|
||||
f 27 191 25
|
||||
f 148 193 29
|
||||
f 30 29 193
|
||||
f 149 33 194
|
||||
f 35 194 33
|
||||
f 150 37 195
|
||||
f 39 195 37
|
||||
f 151 202 1
|
||||
f 2 1 202
|
||||
f 152 5 203
|
||||
f 7 203 5
|
||||
f 204 49 9
|
||||
f 10 9 49
|
||||
f 158 52 13
|
||||
f 14 13 52
|
||||
f 198 21 55
|
||||
f 23 55 21
|
||||
f 162 58 201
|
||||
f 60 201 58
|
||||
f 163 62 25
|
||||
f 26 25 62
|
||||
f 164 29 65
|
||||
f 31 65 29
|
||||
f 165 68 33
|
||||
f 34 33 68
|
||||
f 166 209 37
|
||||
f 38 37 209
|
||||
f 210 75 41
|
||||
f 42 41 75
|
||||
f 211 45 78
|
||||
f 47 78 45
|
||||
f 173 81 49
|
||||
f 50 49 81
|
||||
f 205 84 52
|
||||
f 53 52 84
|
||||
f 177 87 208
|
||||
f 89 208 87
|
||||
f 178 91 58
|
||||
f 59 58 91
|
||||
f 179 3 62
|
||||
f 63 62 3
|
||||
f 180 65 6
|
||||
f 66 6 65
|
||||
f 181 215 68
|
||||
f 69 68 215
|
||||
f 216 15 71
|
||||
f 72 71 15
|
||||
f 186 19 75
|
||||
f 76 75 19
|
||||
f 188 78 22
|
||||
f 79 22 78
|
||||
f 212 101 81
|
||||
f 82 81 101
|
||||
f 191 27 214
|
||||
f 104 214 27
|
||||
f 194 35 87
|
||||
f 88 87 35
|
||||
f 195 39 91
|
||||
f 92 91 39
|
||||
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
@ -0,0 +1,739 @@
|
||||
/*--------------------------------*- 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 featureEdgeMesh;
|
||||
location "constant/triSurface";
|
||||
object statorBlades.eMesh;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// points:
|
||||
|
||||
352
|
||||
(
|
||||
(-1.83691e-17 0.1 0.13)
|
||||
(-1.83691e-17 0.1 0.135)
|
||||
(0.08 1.95937e-17 0.065)
|
||||
(0.08 1.95937e-17 0.07)
|
||||
(0.08 1.95937e-17 0.135)
|
||||
(0.08 1.95937e-17 0.14)
|
||||
(4.89843e-18 -0.08 0.07)
|
||||
(4.89843e-18 -0.08 0.065)
|
||||
(4.89843e-18 -0.08 0.14)
|
||||
(4.89843e-18 -0.08 0.135)
|
||||
(0.1 2.44921e-17 0.02)
|
||||
(0.1 2.44921e-17 0.015)
|
||||
(-0.08 -9.79685e-18 0.07)
|
||||
(-0.08 -9.79685e-18 0.065)
|
||||
(0.1 2.44921e-17 0.09)
|
||||
(0.1 2.44921e-17 0.085)
|
||||
(-0.08 -9.79685e-18 0.14)
|
||||
(-0.08 -9.79685e-18 0.135)
|
||||
(6.12303e-18 -0.1 0.015)
|
||||
(6.12303e-18 -0.1 0.02)
|
||||
(-1.46953e-17 0.08 0.07)
|
||||
(-1.46953e-17 0.08 0.065)
|
||||
(0.1 2.44921e-17 0.16)
|
||||
(0.1 2.44921e-17 0.155)
|
||||
(6.12303e-18 -0.1 0.085)
|
||||
(6.12303e-18 -0.1 0.09)
|
||||
(-0.1 -1.22461e-17 0.015)
|
||||
(-0.1 -1.22461e-17 0.02)
|
||||
(-1.46953e-17 0.08 0.14)
|
||||
(-1.46953e-17 0.08 0.135)
|
||||
(6.12303e-18 -0.1 0.155)
|
||||
(6.12303e-18 -0.1 0.16)
|
||||
(-0.1 -1.22461e-17 0.085)
|
||||
(-0.1 -1.22461e-17 0.09)
|
||||
(-1.83691e-17 0.1 0.015)
|
||||
(-1.83691e-17 0.1 0.02)
|
||||
(-0.1 -1.22461e-17 0.155)
|
||||
(-0.1 -1.22461e-17 0.16)
|
||||
(-1.83691e-17 0.1 0.085)
|
||||
(-1.83691e-17 0.1 0.09)
|
||||
(0.08 1.95937e-17 0.02)
|
||||
(0.08 1.95937e-17 0.025)
|
||||
(-1.83691e-17 0.1 0.155)
|
||||
(-1.83691e-17 0.1 0.16)
|
||||
(0.08 1.95937e-17 0.09)
|
||||
(0.08 1.95937e-17 0.095)
|
||||
(4.89843e-18 -0.08 0.025)
|
||||
(4.89843e-18 -0.08 0.02)
|
||||
(0.08 1.95937e-17 0.16)
|
||||
(0.08 1.95937e-17 0.165)
|
||||
(4.89843e-18 -0.08 0.095)
|
||||
(4.89843e-18 -0.08 0.09)
|
||||
(-0.08 -9.79685e-18 0.025)
|
||||
(-0.08 -9.79685e-18 0.02)
|
||||
(4.89843e-18 -0.08 0.165)
|
||||
(4.89843e-18 -0.08 0.16)
|
||||
(0.1 2.44921e-17 0.045)
|
||||
(0.1 2.44921e-17 0.04)
|
||||
(-0.08 -9.79685e-18 0.095)
|
||||
(-0.08 -9.79685e-18 0.09)
|
||||
(-1.46953e-17 0.08 0.025)
|
||||
(-1.46953e-17 0.08 0.02)
|
||||
(0.1 2.44921e-17 0.115)
|
||||
(0.1 2.44921e-17 0.11)
|
||||
(-0.08 -9.79685e-18 0.165)
|
||||
(-0.08 -9.79685e-18 0.16)
|
||||
(6.12303e-18 -0.1 0.04)
|
||||
(6.12303e-18 -0.1 0.045)
|
||||
(-1.46953e-17 0.08 0.095)
|
||||
(-1.46953e-17 0.08 0.09)
|
||||
(0.1 2.44921e-17 0.185)
|
||||
(0.1 2.44921e-17 0.18)
|
||||
(6.12303e-18 -0.1 0.11)
|
||||
(6.12303e-18 -0.1 0.115)
|
||||
(-0.1 -1.22461e-17 0.04)
|
||||
(-0.1 -1.22461e-17 0.045)
|
||||
(-1.46953e-17 0.08 0.165)
|
||||
(-1.46953e-17 0.08 0.16)
|
||||
(6.12303e-18 -0.1 0.18)
|
||||
(6.12303e-18 -0.1 0.185)
|
||||
(-0.1 -1.22461e-17 0.11)
|
||||
(-0.1 -1.22461e-17 0.115)
|
||||
(-1.83691e-17 0.1 0.04)
|
||||
(-1.83691e-17 0.1 0.045)
|
||||
(-0.1 -1.22461e-17 0.18)
|
||||
(-0.1 -1.22461e-17 0.185)
|
||||
(-1.83691e-17 0.1 0.11)
|
||||
(-1.83691e-17 0.1 0.115)
|
||||
(0.08 1.95937e-17 0.045)
|
||||
(0.08 1.95937e-17 0.05)
|
||||
(-1.83691e-17 0.1 0.18)
|
||||
(-1.83691e-17 0.1 0.185)
|
||||
(0.08 1.95937e-17 0.115)
|
||||
(0.08 1.95937e-17 0.12)
|
||||
(4.89843e-18 -0.08 0.05)
|
||||
(4.89843e-18 -0.08 0.045)
|
||||
(0.08 1.95937e-17 0.185)
|
||||
(0.08 1.95937e-17 0.19)
|
||||
(4.89843e-18 -0.08 0.12)
|
||||
(4.89843e-18 -0.08 0.115)
|
||||
(-0.08 -9.79685e-18 0.05)
|
||||
(-0.08 -9.79685e-18 0.045)
|
||||
(4.89843e-18 -0.08 0.19)
|
||||
(4.89843e-18 -0.08 0.185)
|
||||
(0.1 2.44921e-17 0.07)
|
||||
(0.1 2.44921e-17 0.065)
|
||||
(-0.08 -9.79685e-18 0.12)
|
||||
(-0.08 -9.79685e-18 0.115)
|
||||
(-1.46953e-17 0.08 0.05)
|
||||
(-1.46953e-17 0.08 0.045)
|
||||
(0.1 2.44921e-17 0.14)
|
||||
(0.1 2.44921e-17 0.135)
|
||||
(-0.08 -9.79685e-18 0.19)
|
||||
(-0.08 -9.79685e-18 0.185)
|
||||
(6.12303e-18 -0.1 0.065)
|
||||
(6.12303e-18 -0.1 0.07)
|
||||
(-1.46953e-17 0.08 0.12)
|
||||
(-1.46953e-17 0.08 0.115)
|
||||
(6.12303e-18 -0.1 0.135)
|
||||
(6.12303e-18 -0.1 0.14)
|
||||
(-0.1 -1.22461e-17 0.065)
|
||||
(-0.1 -1.22461e-17 0.07)
|
||||
(-1.46953e-17 0.08 0.19)
|
||||
(-1.46953e-17 0.08 0.185)
|
||||
(-0.1 -1.22461e-17 0.135)
|
||||
(-0.1 -1.22461e-17 0.14)
|
||||
(-1.83691e-17 0.1 0.065)
|
||||
(-1.83691e-17 0.1 0.07)
|
||||
(0.08 1.95937e-17 0)
|
||||
(0.08 1.95937e-17 0.005)
|
||||
(0.085 2.08183e-17 0)
|
||||
(-1.83691e-17 0.1 0.14)
|
||||
(0.08 1.95937e-17 0.075)
|
||||
(4.89843e-18 -0.08 0)
|
||||
(5.20458e-18 -0.085 0)
|
||||
(4.89843e-18 -0.08 0.005)
|
||||
(0.08 1.95937e-17 0.145)
|
||||
(4.89843e-18 -0.08 0.075)
|
||||
(-0.08 -9.79685e-18 0)
|
||||
(-0.085 -1.04092e-17 0)
|
||||
(-0.08 -9.79685e-18 0.005)
|
||||
(4.89843e-18 -0.08 0.145)
|
||||
(0.1 2.44921e-17 0.025)
|
||||
(-0.08 -9.79685e-18 0.075)
|
||||
(-1.46953e-17 0.08 0)
|
||||
(-1.56137e-17 0.085 0)
|
||||
(-1.46953e-17 0.08 0.005)
|
||||
(0.1 2.44921e-17 0.095)
|
||||
(-0.08 -9.79685e-18 0.145)
|
||||
(6.12303e-18 -0.1 0.025)
|
||||
(-1.46953e-17 0.08 0.075)
|
||||
(0.1 2.44921e-17 0.165)
|
||||
(6.12303e-18 -0.1 0.095)
|
||||
(-0.1 -1.22461e-17 0.025)
|
||||
(-1.46953e-17 0.08 0.145)
|
||||
(0.09 2.20429e-17 0)
|
||||
(6.12303e-18 -0.1 0.165)
|
||||
(-0.1 -1.22461e-17 0.095)
|
||||
(-1.83691e-17 0.1 0.025)
|
||||
(-0.1 -1.22461e-17 0.165)
|
||||
(-1.83691e-17 0.1 0.095)
|
||||
(5.51073e-18 -0.09 0)
|
||||
(0.08 1.95937e-17 0.03)
|
||||
(-1.83691e-17 0.1 0.165)
|
||||
(0.08 1.95937e-17 0.1)
|
||||
(-0.09 -1.10215e-17 0)
|
||||
(4.89843e-18 -0.08 0.03)
|
||||
(0.08 1.95937e-17 0.17)
|
||||
(4.89843e-18 -0.08 0.1)
|
||||
(-1.65322e-17 0.09 0)
|
||||
(-0.08 -9.79685e-18 0.03)
|
||||
(4.89843e-18 -0.08 0.17)
|
||||
(0.1 2.44921e-17 0.05)
|
||||
(-0.08 -9.79685e-18 0.1)
|
||||
(-1.46953e-17 0.08 0.03)
|
||||
(0.095 2.32675e-17 0)
|
||||
(0.1 2.44921e-17 0.12)
|
||||
(-0.08 -9.79685e-18 0.17)
|
||||
(6.12303e-18 -0.1 0.05)
|
||||
(-1.46953e-17 0.08 0.1)
|
||||
(0.1 2.44921e-17 0.19)
|
||||
(6.12303e-18 -0.1 0.12)
|
||||
(-0.1 -1.22461e-17 0.05)
|
||||
(5.81688e-18 -0.095 0)
|
||||
(-1.46953e-17 0.08 0.17)
|
||||
(6.12303e-18 -0.1 0.19)
|
||||
(-0.1 -1.22461e-17 0.12)
|
||||
(-1.83691e-17 0.1 0.05)
|
||||
(-0.1 -1.22461e-17 0.19)
|
||||
(-0.095 -1.16338e-17 0)
|
||||
(-1.83691e-17 0.1 0.12)
|
||||
(0.08 1.95937e-17 0.055)
|
||||
(-1.83691e-17 0.1 0.19)
|
||||
(0.08 1.95937e-17 0.125)
|
||||
(-1.74506e-17 0.095 0)
|
||||
(4.89843e-18 -0.08 0.055)
|
||||
(0.08 1.95937e-17 0.195)
|
||||
(4.89843e-18 -0.08 0.125)
|
||||
(0.1 2.44921e-17 0)
|
||||
(0.1 2.44921e-17 0.005)
|
||||
(-0.08 -9.79685e-18 0.055)
|
||||
(4.89843e-18 -0.08 0.195)
|
||||
(0.1 2.44921e-17 0.075)
|
||||
(-0.08 -9.79685e-18 0.125)
|
||||
(6.12303e-18 -0.1 0)
|
||||
(6.12303e-18 -0.1 0.005)
|
||||
(-1.46953e-17 0.08 0.055)
|
||||
(0.1 2.44921e-17 0.145)
|
||||
(-0.08 -9.79685e-18 0.195)
|
||||
(6.12303e-18 -0.1 0.075)
|
||||
(-0.1 -1.22461e-17 0)
|
||||
(-0.1 -1.22461e-17 0.005)
|
||||
(-1.46953e-17 0.08 0.125)
|
||||
(6.12303e-18 -0.1 0.145)
|
||||
(-0.1 -1.22461e-17 0.075)
|
||||
(-1.83691e-17 0.1 0)
|
||||
(-1.83691e-17 0.1 0.005)
|
||||
(-1.46953e-17 0.08 0.195)
|
||||
(-0.1 -1.22461e-17 0.145)
|
||||
(-1.83691e-17 0.1 0.075)
|
||||
(0.08 1.95937e-17 0.01)
|
||||
(-1.83691e-17 0.1 0.145)
|
||||
(0.08 1.95937e-17 0.08)
|
||||
(4.89843e-18 -0.08 0.01)
|
||||
(0.08 1.95937e-17 0.15)
|
||||
(4.89843e-18 -0.08 0.08)
|
||||
(-0.08 -9.79685e-18 0.01)
|
||||
(4.89843e-18 -0.08 0.15)
|
||||
(0.1 2.44921e-17 0.03)
|
||||
(-0.08 -9.79685e-18 0.08)
|
||||
(-1.46953e-17 0.08 0.01)
|
||||
(0.1 2.44921e-17 0.1)
|
||||
(-0.08 -9.79685e-18 0.15)
|
||||
(6.12303e-18 -0.1 0.03)
|
||||
(-1.46953e-17 0.08 0.08)
|
||||
(0.1 2.44921e-17 0.17)
|
||||
(6.12303e-18 -0.1 0.1)
|
||||
(-0.1 -1.22461e-17 0.03)
|
||||
(-1.46953e-17 0.08 0.15)
|
||||
(6.12303e-18 -0.1 0.17)
|
||||
(-0.1 -1.22461e-17 0.1)
|
||||
(-1.83691e-17 0.1 0.03)
|
||||
(-0.1 -1.22461e-17 0.17)
|
||||
(-1.83691e-17 0.1 0.1)
|
||||
(0.08 1.95937e-17 0.035)
|
||||
(-1.83691e-17 0.1 0.17)
|
||||
(0.08 1.95937e-17 0.105)
|
||||
(4.89843e-18 -0.08 0.035)
|
||||
(0.08 1.95937e-17 0.175)
|
||||
(4.89843e-18 -0.08 0.105)
|
||||
(-0.08 -9.79685e-18 0.035)
|
||||
(4.89843e-18 -0.08 0.175)
|
||||
(0.1 2.44921e-17 0.055)
|
||||
(-0.08 -9.79685e-18 0.105)
|
||||
(-1.46953e-17 0.08 0.035)
|
||||
(0.1 2.44921e-17 0.125)
|
||||
(-0.08 -9.79685e-18 0.175)
|
||||
(6.12303e-18 -0.1 0.055)
|
||||
(-1.46953e-17 0.08 0.105)
|
||||
(0.1 2.44921e-17 0.195)
|
||||
(6.12303e-18 -0.1 0.125)
|
||||
(-0.1 -1.22461e-17 0.055)
|
||||
(-1.46953e-17 0.08 0.175)
|
||||
(6.12303e-18 -0.1 0.195)
|
||||
(-0.1 -1.22461e-17 0.125)
|
||||
(-1.83691e-17 0.1 0.055)
|
||||
(-0.1 -1.22461e-17 0.195)
|
||||
(-1.83691e-17 0.1 0.125)
|
||||
(0.08 1.95937e-17 0.06)
|
||||
(-1.83691e-17 0.1 0.195)
|
||||
(0.08 1.95937e-17 0.13)
|
||||
(4.89843e-18 -0.08 0.06)
|
||||
(0.08 1.95937e-17 0.2)
|
||||
(0.085 2.08183e-17 0.2)
|
||||
(4.89843e-18 -0.08 0.13)
|
||||
(4.89843e-18 -0.08 0.2)
|
||||
(5.20458e-18 -0.085 0.2)
|
||||
(0.1 2.44921e-17 0.01)
|
||||
(-0.08 -9.79685e-18 0.06)
|
||||
(-0.08 -9.79685e-18 0.2)
|
||||
(-0.085 -1.04092e-17 0.2)
|
||||
(0.1 2.44921e-17 0.08)
|
||||
(-0.08 -9.79685e-18 0.13)
|
||||
(6.12303e-18 -0.1 0.01)
|
||||
(-1.46953e-17 0.08 0.06)
|
||||
(0.1 2.44921e-17 0.15)
|
||||
(6.12303e-18 -0.1 0.08)
|
||||
(-1.46953e-17 0.08 0.2)
|
||||
(-1.56137e-17 0.085 0.2)
|
||||
(-0.1 -1.22461e-17 0.01)
|
||||
(-1.46953e-17 0.08 0.13)
|
||||
(6.12303e-18 -0.1 0.15)
|
||||
(-0.1 -1.22461e-17 0.08)
|
||||
(0.09 2.20429e-17 0.2)
|
||||
(-1.83691e-17 0.1 0.01)
|
||||
(-0.1 -1.22461e-17 0.15)
|
||||
(5.51073e-18 -0.09 0.2)
|
||||
(-1.83691e-17 0.1 0.08)
|
||||
(0.08 1.95937e-17 0.015)
|
||||
(-1.83691e-17 0.1 0.15)
|
||||
(-0.09 -1.10215e-17 0.2)
|
||||
(0.08 1.95937e-17 0.085)
|
||||
(4.89843e-18 -0.08 0.015)
|
||||
(0.08 1.95937e-17 0.155)
|
||||
(-1.65322e-17 0.09 0.2)
|
||||
(4.89843e-18 -0.08 0.085)
|
||||
(0.095 2.32675e-17 0.2)
|
||||
(-0.08 -9.79685e-18 0.015)
|
||||
(4.89843e-18 -0.08 0.155)
|
||||
(0.1 2.44921e-17 0.035)
|
||||
(-0.08 -9.79685e-18 0.085)
|
||||
(5.81688e-18 -0.095 0.2)
|
||||
(-1.46953e-17 0.08 0.015)
|
||||
(0.1 2.44921e-17 0.105)
|
||||
(-0.08 -9.79685e-18 0.155)
|
||||
(6.12303e-18 -0.1 0.035)
|
||||
(-0.095 -1.16338e-17 0.2)
|
||||
(-1.46953e-17 0.08 0.085)
|
||||
(0.1 2.44921e-17 0.175)
|
||||
(6.12303e-18 -0.1 0.105)
|
||||
(-0.1 -1.22461e-17 0.035)
|
||||
(-1.74506e-17 0.095 0.2)
|
||||
(-1.46953e-17 0.08 0.155)
|
||||
(6.12303e-18 -0.1 0.175)
|
||||
(-0.1 -1.22461e-17 0.105)
|
||||
(-1.83691e-17 0.1 0.035)
|
||||
(0.1 2.44921e-17 0.2)
|
||||
(-0.1 -1.22461e-17 0.175)
|
||||
(6.12303e-18 -0.1 0.2)
|
||||
(-1.83691e-17 0.1 0.105)
|
||||
(0.08 1.95937e-17 0.04)
|
||||
(-0.1 -1.22461e-17 0.2)
|
||||
(-1.83691e-17 0.1 0.175)
|
||||
(-1.83691e-17 0.1 0.2)
|
||||
(0.08 1.95937e-17 0.11)
|
||||
(4.89843e-18 -0.08 0.04)
|
||||
(0.08 1.95937e-17 0.18)
|
||||
(4.89843e-18 -0.08 0.11)
|
||||
(-0.08 -9.79685e-18 0.04)
|
||||
(4.89843e-18 -0.08 0.18)
|
||||
(0.1 2.44921e-17 0.06)
|
||||
(-0.08 -9.79685e-18 0.11)
|
||||
(-1.46953e-17 0.08 0.04)
|
||||
(0.1 2.44921e-17 0.13)
|
||||
(-0.08 -9.79685e-18 0.18)
|
||||
(6.12303e-18 -0.1 0.06)
|
||||
(-1.46953e-17 0.08 0.11)
|
||||
(6.12303e-18 -0.1 0.13)
|
||||
(-0.1 -1.22461e-17 0.06)
|
||||
(-1.46953e-17 0.08 0.18)
|
||||
(-0.1 -1.22461e-17 0.13)
|
||||
(-1.83691e-17 0.1 0.06)
|
||||
)
|
||||
|
||||
|
||||
// edges:
|
||||
|
||||
352
|
||||
(
|
||||
(0 1)
|
||||
(2 3)
|
||||
(4 5)
|
||||
(6 7)
|
||||
(8 9)
|
||||
(10 11)
|
||||
(12 13)
|
||||
(14 15)
|
||||
(16 17)
|
||||
(18 19)
|
||||
(20 21)
|
||||
(22 23)
|
||||
(24 25)
|
||||
(26 27)
|
||||
(28 29)
|
||||
(30 31)
|
||||
(32 33)
|
||||
(34 35)
|
||||
(36 37)
|
||||
(38 39)
|
||||
(40 41)
|
||||
(42 43)
|
||||
(44 45)
|
||||
(46 47)
|
||||
(48 49)
|
||||
(50 51)
|
||||
(52 53)
|
||||
(54 55)
|
||||
(56 57)
|
||||
(58 59)
|
||||
(60 61)
|
||||
(62 63)
|
||||
(64 65)
|
||||
(66 67)
|
||||
(68 69)
|
||||
(70 71)
|
||||
(72 73)
|
||||
(74 75)
|
||||
(76 77)
|
||||
(78 79)
|
||||
(80 81)
|
||||
(82 83)
|
||||
(84 85)
|
||||
(86 87)
|
||||
(88 89)
|
||||
(90 91)
|
||||
(92 93)
|
||||
(94 95)
|
||||
(96 97)
|
||||
(98 99)
|
||||
(100 101)
|
||||
(102 103)
|
||||
(104 105)
|
||||
(106 107)
|
||||
(108 109)
|
||||
(110 111)
|
||||
(112 113)
|
||||
(114 115)
|
||||
(116 117)
|
||||
(118 119)
|
||||
(120 121)
|
||||
(122 123)
|
||||
(124 125)
|
||||
(126 127)
|
||||
(128 129)
|
||||
(130 128)
|
||||
(1 131)
|
||||
(3 132)
|
||||
(133 134)
|
||||
(135 133)
|
||||
(5 136)
|
||||
(137 6)
|
||||
(138 139)
|
||||
(140 138)
|
||||
(141 8)
|
||||
(142 10)
|
||||
(143 12)
|
||||
(144 145)
|
||||
(146 144)
|
||||
(147 14)
|
||||
(148 16)
|
||||
(19 149)
|
||||
(150 20)
|
||||
(151 22)
|
||||
(25 152)
|
||||
(27 153)
|
||||
(154 28)
|
||||
(155 130)
|
||||
(31 156)
|
||||
(33 157)
|
||||
(35 158)
|
||||
(37 159)
|
||||
(39 160)
|
||||
(134 161)
|
||||
(41 162)
|
||||
(43 163)
|
||||
(45 164)
|
||||
(139 165)
|
||||
(166 46)
|
||||
(49 167)
|
||||
(168 50)
|
||||
(145 169)
|
||||
(170 52)
|
||||
(171 54)
|
||||
(172 56)
|
||||
(173 58)
|
||||
(174 60)
|
||||
(175 155)
|
||||
(176 62)
|
||||
(177 64)
|
||||
(67 178)
|
||||
(179 68)
|
||||
(180 70)
|
||||
(73 181)
|
||||
(75 182)
|
||||
(161 183)
|
||||
(184 76)
|
||||
(79 185)
|
||||
(81 186)
|
||||
(83 187)
|
||||
(85 188)
|
||||
(165 189)
|
||||
(87 190)
|
||||
(89 191)
|
||||
(91 192)
|
||||
(93 193)
|
||||
(169 194)
|
||||
(195 94)
|
||||
(97 196)
|
||||
(197 98)
|
||||
(198 175)
|
||||
(199 198)
|
||||
(200 100)
|
||||
(201 102)
|
||||
(202 104)
|
||||
(203 106)
|
||||
(183 204)
|
||||
(204 205)
|
||||
(206 108)
|
||||
(207 110)
|
||||
(208 112)
|
||||
(115 209)
|
||||
(189 210)
|
||||
(210 211)
|
||||
(212 116)
|
||||
(119 213)
|
||||
(121 214)
|
||||
(194 215)
|
||||
(215 216)
|
||||
(217 122)
|
||||
(125 218)
|
||||
(127 219)
|
||||
(129 220)
|
||||
(131 221)
|
||||
(132 222)
|
||||
(223 135)
|
||||
(136 224)
|
||||
(225 137)
|
||||
(226 140)
|
||||
(227 141)
|
||||
(228 142)
|
||||
(229 143)
|
||||
(230 146)
|
||||
(231 147)
|
||||
(232 148)
|
||||
(149 233)
|
||||
(234 150)
|
||||
(235 151)
|
||||
(152 236)
|
||||
(153 237)
|
||||
(238 154)
|
||||
(156 239)
|
||||
(157 240)
|
||||
(158 241)
|
||||
(159 242)
|
||||
(160 243)
|
||||
(162 244)
|
||||
(163 245)
|
||||
(164 246)
|
||||
(247 166)
|
||||
(167 248)
|
||||
(249 168)
|
||||
(250 170)
|
||||
(251 171)
|
||||
(252 172)
|
||||
(253 173)
|
||||
(254 174)
|
||||
(255 176)
|
||||
(256 177)
|
||||
(178 257)
|
||||
(258 179)
|
||||
(259 180)
|
||||
(181 260)
|
||||
(182 261)
|
||||
(262 184)
|
||||
(185 263)
|
||||
(186 264)
|
||||
(187 265)
|
||||
(188 266)
|
||||
(190 267)
|
||||
(191 268)
|
||||
(192 269)
|
||||
(193 270)
|
||||
(271 195)
|
||||
(196 272)
|
||||
(272 273)
|
||||
(274 197)
|
||||
(275 201)
|
||||
(276 275)
|
||||
(277 199)
|
||||
(278 200)
|
||||
(279 208)
|
||||
(280 279)
|
||||
(281 202)
|
||||
(282 203)
|
||||
(205 283)
|
||||
(284 206)
|
||||
(285 207)
|
||||
(209 286)
|
||||
(287 217)
|
||||
(288 287)
|
||||
(211 289)
|
||||
(290 212)
|
||||
(213 291)
|
||||
(214 292)
|
||||
(273 293)
|
||||
(216 294)
|
||||
(218 295)
|
||||
(296 276)
|
||||
(219 297)
|
||||
(220 298)
|
||||
(221 299)
|
||||
(300 280)
|
||||
(222 301)
|
||||
(302 223)
|
||||
(224 303)
|
||||
(304 288)
|
||||
(305 225)
|
||||
(293 306)
|
||||
(307 226)
|
||||
(308 227)
|
||||
(309 228)
|
||||
(310 229)
|
||||
(311 296)
|
||||
(312 230)
|
||||
(313 231)
|
||||
(314 232)
|
||||
(233 315)
|
||||
(316 300)
|
||||
(317 234)
|
||||
(318 235)
|
||||
(236 319)
|
||||
(237 320)
|
||||
(321 304)
|
||||
(322 238)
|
||||
(239 323)
|
||||
(240 324)
|
||||
(241 325)
|
||||
(326 259)
|
||||
(306 326)
|
||||
(242 327)
|
||||
(328 311)
|
||||
(263 328)
|
||||
(243 329)
|
||||
(244 330)
|
||||
(331 316)
|
||||
(266 331)
|
||||
(245 332)
|
||||
(333 321)
|
||||
(269 333)
|
||||
(246 334)
|
||||
(335 247)
|
||||
(248 336)
|
||||
(337 249)
|
||||
(338 250)
|
||||
(339 251)
|
||||
(340 252)
|
||||
(341 253)
|
||||
(342 254)
|
||||
(343 255)
|
||||
(344 256)
|
||||
(257 345)
|
||||
(346 258)
|
||||
(260 347)
|
||||
(261 348)
|
||||
(349 262)
|
||||
(264 350)
|
||||
(265 351)
|
||||
(267 0)
|
||||
(268 2)
|
||||
(270 4)
|
||||
(7 271)
|
||||
(9 274)
|
||||
(11 277)
|
||||
(13 278)
|
||||
(15 281)
|
||||
(17 282)
|
||||
(283 18)
|
||||
(21 284)
|
||||
(23 285)
|
||||
(286 24)
|
||||
(289 26)
|
||||
(29 290)
|
||||
(291 30)
|
||||
(292 32)
|
||||
(294 34)
|
||||
(295 36)
|
||||
(297 38)
|
||||
(298 40)
|
||||
(299 42)
|
||||
(301 44)
|
||||
(47 302)
|
||||
(303 48)
|
||||
(51 305)
|
||||
(53 307)
|
||||
(55 308)
|
||||
(57 309)
|
||||
(59 310)
|
||||
(61 312)
|
||||
(63 313)
|
||||
(65 314)
|
||||
(315 66)
|
||||
(69 317)
|
||||
(71 318)
|
||||
(319 72)
|
||||
(320 74)
|
||||
(77 322)
|
||||
(323 78)
|
||||
(324 80)
|
||||
(325 82)
|
||||
(327 84)
|
||||
(329 86)
|
||||
(330 88)
|
||||
(332 90)
|
||||
(334 92)
|
||||
(95 335)
|
||||
(336 96)
|
||||
(99 337)
|
||||
(101 338)
|
||||
(103 339)
|
||||
(105 340)
|
||||
(107 341)
|
||||
(109 342)
|
||||
(111 343)
|
||||
(113 344)
|
||||
(345 114)
|
||||
(117 346)
|
||||
(347 118)
|
||||
(348 120)
|
||||
(123 349)
|
||||
(350 124)
|
||||
(351 126)
|
||||
)
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
File diff suppressed because it is too large
Load Diff
@ -0,0 +1,20 @@
|
||||
/*--------------------------------*- 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 "constant";
|
||||
object turbulenceProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
simulationType RASModel;
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -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 | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object controlDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
application rhoPimpleDyMFoam;
|
||||
|
||||
startFrom latestTime;
|
||||
|
||||
startTime 0;
|
||||
|
||||
stopAt endTime;
|
||||
|
||||
endTime 2;
|
||||
|
||||
deltaT 1e-3;
|
||||
|
||||
writeControl adjustableRunTime;
|
||||
|
||||
writeInterval 0.02;
|
||||
|
||||
purgeWrite 0;
|
||||
|
||||
writeFormat ascii;
|
||||
|
||||
writePrecision 6;
|
||||
|
||||
writeCompression off;
|
||||
|
||||
timeFormat general;
|
||||
|
||||
timePrecision 6;
|
||||
|
||||
runTimeModifiable true;
|
||||
|
||||
adjustTimeStep yes;
|
||||
|
||||
maxCo 1;
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,46 @@
|
||||
/*--------------------------------*- 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 createBafflesDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
internalFacesOnly true;
|
||||
|
||||
baffles
|
||||
{
|
||||
baffleFaces
|
||||
{
|
||||
type faceZone;
|
||||
zoneName rotatingZone;
|
||||
|
||||
patches
|
||||
{
|
||||
master
|
||||
{
|
||||
name AMI1;
|
||||
type cyclicAMI;
|
||||
matchTolerance 0.0001;
|
||||
neighbourPatch AMI2;
|
||||
transform noOrdering;
|
||||
}
|
||||
slave
|
||||
{
|
||||
$master;
|
||||
name AMI2;
|
||||
neighbourPatch AMI1;
|
||||
}
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,28 @@
|
||||
/*--------------------------------*- 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 createPatchDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// Do a synchronisation of coupled points after creation of any patches.
|
||||
// Note: this does not work with points that are on multiple coupled patches
|
||||
// with transformations (i.e. cyclics).
|
||||
pointSync false;
|
||||
|
||||
// Patches to create. An empty patch list just removes patches with zero
|
||||
// faces from $FOAM_CASE/constant/polyMesh/boundary.
|
||||
patches
|
||||
(
|
||||
);
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,63 @@
|
||||
/*--------------------------------*- 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 fvSchemes;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
ddtSchemes
|
||||
{
|
||||
default Euler;
|
||||
}
|
||||
|
||||
gradSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
}
|
||||
|
||||
divSchemes
|
||||
{
|
||||
default none;
|
||||
|
||||
div(phi,U) Gauss linearUpwind grad(U);
|
||||
div(phi,h) Gauss linearUpwind grad(h);
|
||||
div(phi,K) Gauss linear;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,epsilon) Gauss upwind;
|
||||
div((muEff*dev2(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
laplacianSchemes
|
||||
{
|
||||
default Gauss linear corrected;
|
||||
}
|
||||
|
||||
interpolationSchemes
|
||||
{
|
||||
default linear;
|
||||
}
|
||||
|
||||
snGradSchemes
|
||||
{
|
||||
default corrected;
|
||||
}
|
||||
|
||||
fluxRequired
|
||||
{
|
||||
default no;
|
||||
p ;
|
||||
pcorr ;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,86 @@
|
||||
/*--------------------------------*- 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
|
||||
{
|
||||
p
|
||||
{
|
||||
solver GAMG;
|
||||
|
||||
smoother GaussSeidel;
|
||||
cacheAgglomeration on;
|
||||
agglomerator faceAreaPair;
|
||||
nCellsInCoarsestLevel 10;
|
||||
mergeLevels 1;
|
||||
|
||||
tolerance 1e-6;
|
||||
relTol 0.01;
|
||||
}
|
||||
|
||||
pFinal
|
||||
{
|
||||
$p;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
pcorr
|
||||
{
|
||||
$p;
|
||||
tolerance 1e-2;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
"(rho|U|h|k|epsilon|omega)"
|
||||
{
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-06;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
"(rho|U|h|k|epsilon|omega)Final"
|
||||
{
|
||||
$U;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
PIMPLE
|
||||
{
|
||||
momentumPredictor yes;
|
||||
transonic no;
|
||||
nOuterCorrectors 1;
|
||||
nCorrectors 3;
|
||||
nNonOrthogonalCorrectors 0;
|
||||
rhoMin rhoMin [ 1 -3 0 0 0 ] 0.5;
|
||||
rhoMax rhoMax [ 1 -3 0 0 0 ] 2.0;
|
||||
}
|
||||
|
||||
relaxationFactors
|
||||
{
|
||||
fields
|
||||
{
|
||||
}
|
||||
equations
|
||||
{
|
||||
"(U|h|k|epsilon|omega).*" 1;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,274 @@
|
||||
/*--------------------------------*- 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 snappyHexMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
castellatedMesh true;
|
||||
snap true;
|
||||
addLayers false;
|
||||
|
||||
geometry
|
||||
{
|
||||
AMI.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name AMI;
|
||||
}
|
||||
|
||||
innerInlet.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name innerInlet;
|
||||
}
|
||||
|
||||
outerInlet.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name outerInlet;
|
||||
}
|
||||
|
||||
innerOutlet.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name innerOutlet;
|
||||
}
|
||||
|
||||
outerOutlet.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name outerOutlet;
|
||||
}
|
||||
|
||||
rotorBlades.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name rotorBlades;
|
||||
}
|
||||
|
||||
shaft.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name shaft;
|
||||
}
|
||||
|
||||
statorBlades.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name statorBlades;
|
||||
}
|
||||
|
||||
wall.obj
|
||||
{
|
||||
type triSurfaceMesh;
|
||||
name wall;
|
||||
}
|
||||
};
|
||||
|
||||
castellatedMeshControls
|
||||
{
|
||||
maxLocalCells 100000;
|
||||
maxGlobalCells 2000000;
|
||||
minRefinementCells 0;
|
||||
maxLoadUnbalance 0.10;
|
||||
nCellsBetweenLevels 2;
|
||||
|
||||
features
|
||||
(
|
||||
{
|
||||
file "AMI.eMesh";
|
||||
level 2;
|
||||
}
|
||||
{
|
||||
file "rotorBlades.eMesh";
|
||||
level 2;
|
||||
}
|
||||
{
|
||||
file "shaft.eMesh";
|
||||
level 1;
|
||||
}
|
||||
{
|
||||
file "statorBlades.eMesh";
|
||||
level 2;
|
||||
}
|
||||
{
|
||||
file "wall.eMesh";
|
||||
level 1;
|
||||
}
|
||||
);
|
||||
|
||||
refinementSurfaces
|
||||
{
|
||||
AMI
|
||||
{
|
||||
level (2 2);
|
||||
faceZone rotatingZone;
|
||||
cellZone rotatingZone;
|
||||
cellZoneInside inside;
|
||||
}
|
||||
|
||||
innerInlet
|
||||
{
|
||||
level (1 1);
|
||||
patchInfo
|
||||
{
|
||||
type patch;
|
||||
inGroups (inlet);
|
||||
}
|
||||
}
|
||||
|
||||
outerInlet
|
||||
{
|
||||
level (1 1);
|
||||
patchInfo
|
||||
{
|
||||
type patch;
|
||||
inGroups (inlet);
|
||||
}
|
||||
}
|
||||
|
||||
innerOutlet
|
||||
{
|
||||
level (1 1);
|
||||
patchInfo
|
||||
{
|
||||
type patch;
|
||||
inGroups (outlet);
|
||||
}
|
||||
}
|
||||
|
||||
outerOutlet
|
||||
{
|
||||
level (1 1);
|
||||
patchInfo
|
||||
{
|
||||
type patch;
|
||||
inGroups (outlet);
|
||||
}
|
||||
}
|
||||
|
||||
rotorBlades
|
||||
{
|
||||
level (2 2);
|
||||
patchInfo
|
||||
{
|
||||
type wall;
|
||||
inGroups (movingWalls);
|
||||
}
|
||||
faceZone rotorBlades;
|
||||
faceType baffle;
|
||||
}
|
||||
|
||||
shaft
|
||||
{
|
||||
level (1 1);
|
||||
patchInfo
|
||||
{
|
||||
type wall;
|
||||
inGroups (movingWalls);
|
||||
}
|
||||
}
|
||||
|
||||
statorBlades
|
||||
{
|
||||
level (2 2);
|
||||
patchInfo
|
||||
{
|
||||
type wall;
|
||||
inGroups (staticWalls);
|
||||
}
|
||||
faceZone statorBlades;
|
||||
faceType baffle;
|
||||
}
|
||||
|
||||
wall
|
||||
{
|
||||
level (1 1);
|
||||
patchInfo
|
||||
{
|
||||
type wall;
|
||||
inGroups (staticWalls);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
resolveFeatureAngle 30;
|
||||
|
||||
refinementRegions
|
||||
{
|
||||
}
|
||||
|
||||
locationInMesh (0.035 0.043 0.017);
|
||||
|
||||
allowFreeStandingZoneFaces true;
|
||||
}
|
||||
|
||||
snapControls
|
||||
{
|
||||
nSmoothPatch 3;
|
||||
tolerance 1.0;
|
||||
nSolveIter 100;
|
||||
nRelaxIter 5;
|
||||
nFeatureSnapIter 5;
|
||||
}
|
||||
|
||||
addLayersControls
|
||||
{
|
||||
relativeSizes true;
|
||||
layers
|
||||
{
|
||||
}
|
||||
expansionRatio 1.0;
|
||||
finalLayerThickness 0.3;
|
||||
minThickness 0.1;
|
||||
nGrow 0;
|
||||
featureAngle 30;
|
||||
nRelaxIter 3;
|
||||
nSmoothSurfaceNormals 1;
|
||||
nSmoothNormals 3;
|
||||
nSmoothThickness 10;
|
||||
maxFaceThicknessRatio 0.5;
|
||||
maxThicknessToMedialRatio 0.3;
|
||||
minMedianAxisAngle 90;
|
||||
nBufferCellsNoExtrude 0;
|
||||
nLayerIter 50;
|
||||
}
|
||||
|
||||
meshQualityControls
|
||||
{
|
||||
maxNonOrtho 65;
|
||||
|
||||
maxBoundarySkewness 20;
|
||||
maxInternalSkewness 4;
|
||||
|
||||
maxConcave 80;
|
||||
minVol 1e-13;
|
||||
minTetQuality -1;
|
||||
minArea -1;
|
||||
minTwist 0.01;
|
||||
minDeterminant 0.001;
|
||||
minFaceWeight 0.05;
|
||||
minVolRatio 0.01;
|
||||
minTriangleTwist -1;
|
||||
nSmoothScale 4;
|
||||
errorReduction 0.75;
|
||||
relaxed
|
||||
{
|
||||
maxNonOrtho 75;
|
||||
}
|
||||
}
|
||||
|
||||
debug 0;
|
||||
mergeTolerance 1e-6;
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,42 @@
|
||||
/*--------------------------------*- 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 surfaceFeatureExtractDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
AMI.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
shaft.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
wall.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
statorBlades.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
rotorBlades.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,31 @@
|
||||
/*--------------------------------*- 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 surfaceFeatureExtractDictDefaults;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
extractionMethod extractFromSurface;
|
||||
|
||||
extractFromSurfaceCoeffs
|
||||
{
|
||||
includedAngle 150; // Identifes a feature when angle
|
||||
} // between faces < includedAngle
|
||||
|
||||
trimFeatures
|
||||
{
|
||||
minElem 10; // minimum edges within a feature
|
||||
}
|
||||
|
||||
writeObj yes; // writes out .obj files to view features
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
@ -23,7 +23,6 @@ ddtSchemes
|
||||
gradSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
grad(p) Gauss linear;
|
||||
}
|
||||
|
||||
divSchemes
|
||||
@ -35,18 +34,13 @@ divSchemes
|
||||
|
||||
laplacianSchemes
|
||||
{
|
||||
default none;
|
||||
laplacian(nu,U) Gauss linear corrected;
|
||||
laplacian(rAU,pcorr) Gauss linear corrected;
|
||||
laplacian(rAU,p) Gauss linear corrected;
|
||||
default Gauss linear corrected;
|
||||
laplacian(diffusivity,cellMotionU) Gauss linear uncorrected;
|
||||
laplacian(nuEff,U) Gauss linear uncorrected;
|
||||
}
|
||||
|
||||
interpolationSchemes
|
||||
{
|
||||
default linear;
|
||||
interpolate(HbyA) linear;
|
||||
}
|
||||
|
||||
snGradSchemes
|
||||
|
||||
@ -94,8 +94,8 @@ PIMPLE
|
||||
{
|
||||
momentumPredictor no;
|
||||
nOuterCorrectors 1;
|
||||
nCorrectors 2;
|
||||
nNonOrthogonalCorrectors 1;
|
||||
nCorrectors 3;
|
||||
nNonOrthogonalCorrectors 0;
|
||||
|
||||
cAlpha 0;
|
||||
nAlphaCorr 1;
|
||||
|
||||
Reference in New Issue
Block a user