exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Viscoelastic models and tutorials update. Jovani Favero and HJ

Browse source
This commit is contained in:
Hrvoje Jasak 2010-10-12 12:44:15 +01:00
parent f129e6c013
commit 98647fcdb3
13 changed files with 641 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

47
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/0/U Normal file
View file

@ -0,0 +1,47 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0.1 0 0);
}
fixedWalls
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type zeroGradient;
}
simetry
{
type symmetryPlane;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

46
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/0/p Normal file
View file

@ -0,0 +1,46 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
fixedWalls
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
simetry
{
type symmetryPlane;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

46
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/0/tau Normal file
View file

@ -0,0 +1,46 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volSymmTensorField;
object tau;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform (0 0 0 0 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 0 0 0 0);
}
fixedWalls
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
simetry
{
type symmetryPlane;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

8
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/Allclean Executable file
View file

@ -0,0 +1,8 @@
#!/bin/sh
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
cleanCase
echo "done!"

16
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/Allrun Executable file
View file

@ -0,0 +1,16 @@
#!/bin/sh
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
# Get application name
application=`getApplication`
echo "Running blockMesh application"
runApplication blockMesh
echo "Running viscoelasticFluidFoam solver"
echo "White-Metzner_Carreau-Yasuda model:"
runApplication $application
echo "done!"

30
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/README Normal file
View file

@ -0,0 +1,30 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
#
#+TITLE: *White-Metzner_Carreau-Yasuda viscoelastic constitutive model README file for version 1.6-ext*
#+AUTHOR: Jovani
#+DATE: 05 October 2010
* This is a simulation case for a viscoelastic fluid described by the
White-Metzner with Carreau-Yasuda viscosity constitutive model (White-Metzner_Carreau-Yasuda).
The geometry was the 4:1 plane contraction.
*** This is only an illustrative case for the White-Metzner with Carreau-Yasuda constitutive model.
Therefore, there are not more details published until now for this model using this geometry
and viscoelasticFluidFoam solver.

98
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/constant/polyMesh/blockMeshDict Normal file
View file

@ -0,0 +1,98 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.0032;
vertices
(
(0 0 0)
(80 0 0)
(0 1 0)
(80 1 0)
(0 2.5 0)
(80 2.5 0)
(0 4 0)
(80 4 0)
(130 0 0)
(130 1 0)
(0 0 0.1)
(80 0 0.1)
(0 1 0.1)
(80 1 0.1)
(0 2.5 0.1)
(80 2.5 0.1)
(0 4 0.1)
(80 4 0.1)
(130 0 0.1)
(130 1 0.1)
);
blocks
(
hex (0 1 3 2 10 11 13 12) (150 30 1) simpleGrading (0.002 0.3 1)
hex (2 3 5 4 12 13 15 14) (150 45 1) simpleGrading (0.002 3.333 1)
hex (4 5 7 6 14 15 17 16) (150 45 1) simpleGrading (0.002 0.3 1)
hex (1 8 9 3 11 18 19 13) (90 30 1) simpleGrading (500 0.3 1)
);
edges
(
);
patches
(
patch inlet
(
(0 10 12 2)
(2 12 14 4)
(4 14 16 6)
)
wall fixedWalls
(
(6 16 17 7)
(5 7 17 15)
(3 5 15 13)
(3 13 19 9)
)
patch outlet
(
(8 9 19 18)
)
symmetryPlane simetry
(
(0 1 11 10)
(1 8 18 11)
)
empty frontAndBack
(
(0 2 3 1)
(2 4 5 3)
(4 6 7 5)
(1 3 9 8)
(10 11 13 12)
(12 13 15 14)
(14 15 17 16)
(11 18 19 13)
)
);
mergePatchPairs
(
);
// ************************************************************************* //

52
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/constant/polyMesh/boundary Normal file
View file

@ -0,0 +1,52 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
5
(
inlet
{
type patch;
nFaces 120;
startFace 41040;
}
fixedWalls
{
type wall;
nFaces 330;
startFace 41160;
}
outlet
{
type patch;
nFaces 30;
startFace 41490;
}
simetry
{
type symmetryPlane;
nFaces 240;
startFace 41520;
}
frontAndBack
{
type empty;
nFaces 41400;
startFace 41760;
}
)
// ************************************************************************* //

34
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/constant/viscoelasticProperties Normal file
View file

@ -0,0 +1,34 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object viscoelasticProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
rheology
{
type WhiteMetznerCarreauYasuda;
rho rho [1 -3 0 0 0 0 0] 850;
etaS etaS [1 -1 -1 0 0 0 0] 0.0;
etaP etaP [1 -1 -1 0 0 0 0] 3.14e2;
lambda lambda [0 0 1 0 0 0 0] 8.7487e-1;
K K [0 0 1 0 0 0 0] 6.0362e-1;
L L [0 0 1 0 0 0 0] 4.0056e1;
m m [0 0 0 0 0 0 0] 1.7314e-1;
n n [0 0 0 0 0 0 0] 6.4293e-2;
a a [0 0 0 0 0 0 0] 5.5061e-1 ;
b b [0 0 0 0 0 0 0] 2.3229e-1;
}
// ************************************************************************* //

55
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/system/controlDict Normal file
View file

@ -0,0 +1,55 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application viscoelasticFluidFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 2;
deltaT 1e-8;
writeControl adjustableRunTime;
writeInterval 1;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
graphFormat raw;
runTimeModifiable yes;
adjustTimeStep on;
maxCo 0.8;
maxDeltaT 0.01;
// ************************************************************************* //

54
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/system/decomposeParDict Normal file
View file

@ -0,0 +1,54 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
numberOfSubdomains 2;
method metis;
simpleCoeffs
{
n (2 2 1);
delta 0.001;
}
hierarchicalCoeffs
{
n (1 1 1);
delta 0.001;
order xyz;
}
metisCoeffs
{
processorWeights
(
2
2
);
}
manualCoeffs
{
dataFile "";
}
distributed no;
roots
(
);
// ************************************************************************* //

63
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/system/fvSchemes Normal file
View file

@ -0,0 +1,63 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,tau) Gauss upwind;
div(tau) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(etaPEff,U) Gauss linear corrected;
laplacian(etaPEff+etaS,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(HbyA) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //

92
tutorials/viscoelastic/viscoelasticFluidFoam/WhiteMetzner/CarreauYasuda/system/fvSolution Normal file
View file

@ -0,0 +1,92 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver PCG;
preconditioner
{
// preconditioner Cholesky;
preconditioner AMG;
cycle W-cycle;
policy AAMG;
nPreSweeps 0;
nPostSweeps 2;
groupSize 4;
minCoarseEqns 20;
nMaxLevels 100;
scale off;
smoother ILU;
}
tolerance 1e-07;
relTol 0;
minIter 0;
maxIter 800;
}
U
{
solver BiCGStab;
preconditioner
{
preconditioner Cholesky;
}
tolerance 1e-6;
relTol 0;
minIter 0;
maxIter 1000;
}
tau
{
solver BiCGStab;
preconditioner
{
preconditioner Cholesky;
}
tolerance 1e-6;
relTol 0;
minIter 0;
maxIter 1000;
};
}
PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
p 0.3;
U 0.5;
tau 0.3;
}
// ************************************************************************* //