Ralinus |
February 11, 2014 07:30 |
1 Attachment(s)
Hello!
Im trying do simulate a moving mesh with a sliding interface with pimpleDymFoam. So far your setup and case was very helpful.
So far I got my geometry defined and also got your Valve class compiled into my openfoam 2.2.1 version. Then i succeeded in running pimpleDymFoam with my case.
The actual problem i have is, that the interface ist actually not detached and attached, which leads to enourmous cell deformation.
The picture shows my mesh after 5 steps, the lower part is moving left, the upper one is at rest.
Here a solver log:
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.1-57f3c3617a2d
Exec : pimpleDyMFoam
Date : Feb 11 2014
Time : 13:27:17
Host : "HEIPC85541-ubuntu"
PID : 1777
Case : /home/pizzard/OpenFOAM/pizzard-2.2.1/haschur/pimpleDyMFoam/cellMotionFirstTry
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Selecting dynamicFvMesh linearValveFvMesh
Selecting motion solver: velocityComponentLaplacian
Selecting motion diffusion: directional
void linearValveFvMesh::addZonesAndModifiers() : Zones and modifiers already present. Skipping.
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading field rAU if present
No finite volume options present
PIMPLE: no residual control data found. Calculations will employ 2 corrector loops
Starting time loop
Courant Number mean: 2.35719e-20 max: 9.31866e-18
deltaT = 2e-06
Time = 2e-06
linearValveFvMesh: detached!
Sliding interfaces decoupled
Sliding Interfaces sucessfully decoupled!
GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 9.91454e-09, No Iterations 321
Coupling sliding interfaces
Sliding Interfaces sucessfully decoupled!
Moving points post slider attach
linearValveFvMesh: detached!
Sliding interfaces coupled: 0
GAMG: Solving for pcorr, Initial residual = 1.50744e-05, Final residual = 5.8296e-08, No Iterations 8
GAMG: Solving for pcorr, Initial residual = 9.90667e-06, Final residual = 4.92653e-08, No Iterations 4
GAMG: Solving for pcorr, Initial residual = 3.34591e-06, Final residual = 4.46853e-08, No Iterations 3
time step continuity errors : sum local = 1.5378e-22, global = 1.7368e-23, cumulative = 1.7368e-23
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.0107594, Final residual = 0.000167086, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0945104, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00304966, Final residual = 4.73578e-05, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.00610149, No Iterations 6
GAMG: Solving for p, Initial residual = 0.0175591, Final residual = 0.000150969, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00556229, Final residual = 3.78537e-05, No Iterations 5
time step continuity errors : sum local = 8.60855e-08, global = 2.5772e-08, cumulative = 2.5772e-08
GAMG: Solving for p, Initial residual = 0.0300467, Final residual = 0.000239986, No Iterations 6
GAMG: Solving for p, Initial residual = 0.0241693, Final residual = 0.000186804, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00680939, Final residual = 5.70799e-05, No Iterations 4
time step continuity errors : sum local = 5.68957e-08, global = -1.38687e-08, cumulative = 1.19033e-08
GAMG: Solving for p, Initial residual = 0.0250689, Final residual = 0.000146277, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00795179, Final residual = 6.37941e-05, No Iterations 5
GAMG: Solving for p, Initial residual = 0.0032912, Final residual = 3.15121e-05, No Iterations 4
time step continuity errors : sum local = 2.83412e-08, global = 2.01249e-08, cumulative = 3.20282e-08
GAMG: Solving for p, Initial residual = 0.00836049, Final residual = 4.95432e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00224231, Final residual = 1.21857e-05, No Iterations 4
GAMG: Solving for p, Initial residual = 0.000868547, Final residual = 5.3291e-06, No Iterations 3
time step continuity errors : sum local = 4.6474e-09, global = -6.49243e-10, cumulative = 3.1379e-08
GAMG: Solving for p, Initial residual = 0.00334386, Final residual = 2.45829e-05, No Iterations 3
GAMG: Solving for p, Initial residual = 0.000802673, Final residual = 5.40721e-06, No Iterations 4
GAMG: Solving for p, Initial residual = 0.000319812, Final residual = 2.14775e-06, No Iterations 3
time step continuity errors : sum local = 1.87786e-09, global = 3.85007e-10, cumulative = 3.1764e-08
PIMPLE: iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.26793, Final residual = 3.48295e-08, No Iterations 8
DILUPBiCG: Solving for Uy, Initial residual = 0.140926, Final residual = 1.70975e-08, No Iterations 8
DILUPBiCG: Solving for Uz, Initial residual = 0.390905, Final residual = 4.5078e-08, No Iterations 8
GAMG: Solving for p, Initial residual = 0.155118, Final residual = 0.00104918, No Iterations 5
GAMG: Solving for p, Initial residual = 0.0158902, Final residual = 0.000133865, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00473208, Final residual = 2.87043e-05, No Iterations 4
time step continuity errors : sum local = 2.69914e-08, global = 1.26524e-08, cumulative = 4.44164e-08
GAMG: Solving for p, Initial residual = 0.0940748, Final residual = 0.00086161, No Iterations 5
GAMG: Solving for p, Initial residual = 0.0217233, Final residual = 0.000166681, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00587747, Final residual = 3.12268e-05, No Iterations 4
time step continuity errors : sum local = 1.56548e-08, global = -7.91986e-09, cumulative = 3.64965e-08
GAMG: Solving for p, Initial residual = 0.0300275, Final residual = 0.000122957, No Iterations 4
GAMG: Solving for p, Initial residual = 0.00550682, Final residual = 3.29272e-05, No Iterations 5
GAMG: Solving for p, Initial residual = 0.00201069, Final residual = 1.80117e-05, No Iterations 3
time step continuity errors : sum local = 8.98568e-09, global = -7.76606e-10, cumulative = 3.57199e-08
GAMG: Solving for p, Initial residual = 0.00832041, Final residual = 7.13671e-05, No Iterations 3
GAMG: Solving for p, Initial residual = 0.00163519, Final residual = 5.93839e-06, No Iterations 4
GAMG: Solving for p, Initial residual = 0.000627401, Final residual = 2.99847e-06, No Iterations 4
time step continuity errors : sum local = 1.51558e-09, global = -3.25407e-10, cumulative = 3.53945e-08
GAMG: Solving for p, Initial residual = 0.00338704, Final residual = 2.54999e-05, No Iterations 3
GAMG: Solving for p, Initial residual = 0.000560385, Final residual = 5.51933e-06, No Iterations 3
GAMG: Solving for p, Initial residual = 0.000209328, Final residual = 8.23839e-08, No Iterations 8
time step continuity errors : sum local = 4.17036e-11, global = 1.08622e-11, cumulative = 3.54054e-08
ExecutionTime = 6.47 s ClockTime = 6 s
I have 2 questions: in your control dict, you include libs ("libmyTopochangeFvMesh.so"); What is this and where do I get it from? Do i need it?
So far I dont include this.
What i want is that the solver splits up the cells at that interface and does a new faces merging in every step. How do I tell this in the dynamic mesh dict?
If i try to run run your case i get the follwing error:
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.1-57f3c3617a2d
Exec : moveDynamicMesh
Date : Feb 11 2014
Time : 14:14:25
Host : "HEIPC85541-ubuntu"
PID : 2824
Case : /home/pizzard/Arbeitsfläche/slideTest
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Selecting dynamicFvMesh linearValveFvMesh
--> FOAM FATAL ERROR:
Unknown dynamicFvMesh type linearValveFvMesh
Valid dynamicFvMesh types are :
6
(
dynamicInkJetFvMesh
dynamicMotionSolverFvMesh
dynamicRefineFvMesh
multiSolidBodyMotionFvMesh
solidBodyMotionFvMesh
staticFvMesh
)
From function dynamicFvMesh::New(const IOobject&)
in file dynamicFvMesh/dynamicFvMeshNew.C at line 81.
FOAM exiting
If i use moveMesh instead everything works.
And i also get a warning that libmyTopochangeFvMesh.so could not be found.
|