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

Version of surfactantFoam that is case-independent and a simple tutorial for it

Browse source 
--HG--
branch : proposal/famTutorialReorganizationBranch
rename : tutorials/surfaceTracking/surfactantFoam/sphereTransport/sphereSurfactantFoam/Make/files => applications/solvers/surfaceTracking/surfactantFoam/Make/files
rename : tutorials/surfaceTracking/surfactantFoam/sphereTransport/sphereSurfactantFoam/Make/options => applications/solvers/surfaceTracking/surfactantFoam/Make/options
rename : tutorials/surfaceTracking/surfactantFoam/sphereTransport/sphereSurfactantFoam/createFaFields.H => applications/solvers/surfaceTracking/surfactantFoam/createFaFields.H
rename : tutorials/surfaceTracking/surfactantFoam/sphereTransport/sphereSurfactantFoam/createFaMesh.H => applications/solvers/surfaceTracking/surfactantFoam/createFaMesh.H
rename : tutorials/surfaceTracking/surfactantFoam/sphereTransport/sphereSurfactantFoam/createVolFields.H => applications/solvers/surfaceTracking/surfactantFoam/createVolFields.H
rename : tutorials/surfaceTracking/surfactantFoam/sphereTransport/sphereSurfactantFoam/surfactantFoam.C => applications/solvers/surfaceTracking/surfactantFoam/surfactantFoam.C
This commit is contained in:
Bernhard F.W. Gschaider 2011-01-29 15:17:55 +01:00
parent 2b6b5ac90b
commit c613c505d3
21 changed files with 802 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

3
applications/solvers/surfaceTracking/surfactantFoam/Make/files Normal file
View file

@ -0,0 +1,3 @@
surfactantFoam.C
EXE = $(FOAM_APPBIN)/surfactantFoam

9
applications/solvers/surfaceTracking/surfactantFoam/Make/options Normal file
View file

@ -0,0 +1,9 @@
EXE_INC = \
-I$(LIB_SRC)/finiteArea/lnInclude \
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/cfdTools/general/lnInclude
EXE_LIBS = \
-lfiniteArea \
-lfiniteVolume \
-llduSolvers

63
applications/solvers/surfaceTracking/surfactantFoam/createFaFields.H Normal file
View file

@ -0,0 +1,63 @@
Info << "Reading field Cs" << endl;
areaScalarField Cs
(
IOobject
(
"Cs",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
aMesh
);
Info<< "Reading transportProperties\n" << endl;
IOdictionary transportProperties
(
IOobject
(
"transportProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);
Info<< "Reading diffusivity D\n" << endl;
dimensionedScalar Ds
(
transportProperties.lookup("Ds")
);
areaVectorField Us
(
IOobject
(
"Us",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
),
aMesh
);
edgeScalarField phis
(
IOobject
(
"phis",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
linearEdgeInterpolate(Us) & aMesh.Le()
);

2
applications/solvers/surfaceTracking/surfactantFoam/createFaMesh.H Normal file
View file

@ -0,0 +1,2 @@
// Create Finite Area mesh
faMesh aMesh(mesh);

36
applications/solvers/surfaceTracking/surfactantFoam/createVolFields.H Normal file
View file

@ -0,0 +1,36 @@
// Create volume-to surface mapping object
volSurfaceMapping vsm(aMesh);
volScalarField Cvf
(
IOobject
(
"Cvf",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("0", dimless/dimLength, 0)
);
vsm.mapToVolume(Cs, Cvf.boundaryField());
Cvf.write();
volVectorField U
(
IOobject
(
"U",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedVector("zero", dimVelocity, vector::zero)
);
vsm.mapToVolume(Us, U.boundaryField());
U.write();

83
applications/solvers/surfaceTracking/surfactantFoam/surfactantFoam.C Normal file
View file

@ -0,0 +1,83 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2004-2007 Z. Tukovic and H. Jasak
\\/ 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 2 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, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
Application
finiteAreaFoam
Description
\*---------------------------------------------------------------------------*/
#include "fvCFD.H"
#include "faCFD.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int main(int argc, char *argv[])
{
# include "setRootCase.H"
# include "createTime.H"
# include "createMesh.H"
# include "createFaMesh.H"
# include "createFaFields.H"
# include "createVolFields.H"
Info<< "Total mass of surfactant: "
<< sum(Cs.internalField()*aMesh.S()) << endl;
Info << "\nStarting time loop\n" << endl;
for (runTime++; !runTime.end(); runTime++)
{
Info << "Time = " << runTime.value() << endl;
faScalarMatrix CsEqn
(
fam::ddt(Cs)
+ fam::div(phis, Cs)
- fam::laplacian(Ds, Cs)
);
CsEqn.solve();
if (runTime.outputTime())
{
vsm.mapToVolume(Cs, Cvf.boundaryField());
runTime.write();
}
Info<< "Total mass of surfactant: "
<< sum(Cs.internalField()*aMesh.S()) << endl;
Info << "ExecutionTime = "
<< scalar(runTime.elapsedCpuTime())
<< " s\n" << endl << endl;
}
return(0);
}
// ************************************************************************* //

41
tutorials/surfaceTracking/surfactantFoam/planeTransport/0/Cs Normal file
View file

@ -0,0 +1,41 @@
/*--------------------------------*- 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 areaScalarField;
location "0";
object Cs;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -2 0 0 0 0 0];
internalField uniform 0;
referenceLevel 0;
boundaryField
{
inlet {
type fixedValue;
value uniform 1;
}
outlet {
type inletOutlet;
value uniform 1;
phi phis;
inletValue uniform 0;
}
bound {
type symmetryPlane;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

52
tutorials/surfaceTracking/surfactantFoam/planeTransport/0/Cvf 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 volScalarField;
location "0";
object Cvf;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 -1 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
bound
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
bottom
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
}
// ************************************************************************* //

52
tutorials/surfaceTracking/surfactantFoam/planeTransport/0/U 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 volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
inlet
{
type calculated;
value uniform (0 0 0);
}
bound
{
type calculated;
value uniform (0 0 0);
}
outlet
{
type calculated;
value uniform (0 0 0);
}
bottom
{
type calculated;
value uniform (0 0 0);
}
top
{
type calculated;
value uniform (0.05 0 0);
}
}
// ************************************************************************* //

36
tutorials/surfaceTracking/surfactantFoam/planeTransport/0/Us Normal file
View file

@ -0,0 +1,36 @@
/*--------------------------------*- 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 areaVectorField;
location "0";
object Us;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0.05 0 0);
boundaryField
{
inlet {
type fixedValue;
value $internalField;
}
outlet {
type zeroGradient;
}
bound {
type symmetryPlane;
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

7
tutorials/surfaceTracking/surfactantFoam/planeTransport/Allclean Executable file
View file

@ -0,0 +1,7 @@
#!/bin/sh
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
cleanTimeDirectories
cleanFaMesh

9
tutorials/surfaceTracking/surfactantFoam/planeTransport/Allrun Executable file
View file

@ -0,0 +1,9 @@
#!/bin/sh
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
application="surfactantFoam"
runApplication blockMesh
runApplication makeFaMesh
runApplication $application

40
tutorials/surfaceTracking/surfactantFoam/planeTransport/constant/faMesh/faMeshDefinition Normal file
View file

@ -0,0 +1,40 @@
/*--------------------------------*- 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 faMeshDefinition;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
polyMeshPatches 1( top );
boundary
{
inlet {
type patch;
ownerPolyPatch top;
neighbourPolyPatch inlet;
}
outlet {
type patch;
ownerPolyPatch top;
neighbourPolyPatch outlet;
}
bound {
type symmetryPlane;
ownerPolyPatch top;
neighbourPolyPatch bound;
}
}
// ************************************************************************** //

69
tutorials/surfaceTracking/surfactantFoam/planeTransport/constant/polyMesh/blockMeshDict Normal file
View file

@ -0,0 +1,69 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.1;
vertices
(
(0 0 0)
(3 0 0)
(3 1 0)
(0 1 0)
(0 0 0.1)
(3 0 0.1)
(3 1 0.1)
(0 1 0.1)
);
blocks
(
hex (0 1 2 3 4 5 6 7) (60 20 1) simpleGrading (1 1 1)
);
edges
(
);
patches
(
patch inlet
(
(0 4 7 3)
)
wall bound
(
(3 7 6 2)
(1 5 4 0)
)
patch outlet
(
(2 6 5 1)
)
patch bottom
(
(0 3 2 1)
)
patch top
(
(4 5 6 7)
)
);
mergePatchPairs
(
);
// ************************************************************************* //

52
tutorials/surfaceTracking/surfactantFoam/planeTransport/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 20;
startFace 2320;
}
bound
{
type wall;
nFaces 120;
startFace 2340;
}
outlet
{
type patch;
nFaces 20;
startFace 2460;
}
bottom
{
type patch;
nFaces 1200;
startFace 2480;
}
top
{
type patch;
nFaces 1200;
startFace 3680;
}
)
// ************************************************************************* //

21
tutorials/surfaceTracking/surfactantFoam/planeTransport/constant/transportProperties Normal file
View file

@ -0,0 +1,21 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Ds Ds [ 0 2 -1 0 0 0 0 ] 0.00001;
// ************************************************************************* //

48
tutorials/surfaceTracking/surfactantFoam/planeTransport/system/controlDict Normal file
View file

@ -0,0 +1,48 @@
/*--------------------------------*- 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;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application surfactantFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 6;
deltaT 0.1;
writeControl runTime;
writeInterval 0.2;
purgeWrite 0;
writeFormat ascii;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
// ************************************************************************* //

57
tutorials/surfaceTracking/surfactantFoam/planeTransport/system/faSchemes Normal file
View file

@ -0,0 +1,57 @@
/*--------------------------------*- 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;
location "system";
object faSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
}
divSchemes
{
default none;
div(phis,Cs) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(Ds,Cs) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
p;
}
// ************************************************************************* //

52
tutorials/surfaceTracking/surfactantFoam/planeTransport/system/faSolution 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 dictionary;
location "system";
object faSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
Cs
{
solver PBiCG;
preconditioner DILU;
minIter 0;
maxIter 2000;
tolerance 1e-06;
relTol 0;
}
}
PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
}
SIMPLE
{
nTimeCorrectors 6;
nNonOrthogonalCorrectors 1;
}
relaxationFactors
{
p 0.7;
U 0.7;
k 0.7;
epsilon 0.7;
R 0.7;
}
// ************************************************************************* //

47
tutorials/surfaceTracking/surfactantFoam/planeTransport/system/fvSchemes 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 dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
}
gradSchemes
{
}
divSchemes
{
}
laplacianSchemes
{
}
interpolationSchemes
{
}
snGradSchemes
{
}
fluxRequired
{
}
// ************************************************************************* //

23
tutorials/surfaceTracking/surfactantFoam/planeTransport/system/fvSolution Normal file
View file

@ -0,0 +1,23 @@
/*--------------------------------*- 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;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
}
// ************************************************************************* //