CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Moving Mesh

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By eysteinn

Reply
 
LinkBack Thread Tools Display Modes
Old   May 24, 2012, 02:48
Default Moving Mesh
  #1
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
Dear All,

I am trying to run a moving-mesh problem when I give the moveDynamicMesh command, I get this error:

Quote:
lab@lab-laptop:~/Documenti/cases_OF/OF_case18_moving_door/moving_door_0$ moveDynamicMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.0-0bc225064152
Exec : moveDynamicMesh
Date : May 24 2012
Time : 10:00:10
Host : "lab-laptop"
PID : 3028
Case : /home/lab/Documenti/cases_OF/OF_case18_moving_door/moving_door_0
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 dynamicMotionSolverFvMesh
Selecting motion solver: displacementLaplacian
Selecting motion diffusion: inverseDistance
Time = 5e-06
GAMG: Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for cellDisplacementz, Initial residual = 1, Final residual = 8.56595e-07, No Iterations 9
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Mesh topology OK.
Boundary openness (-3.45182e-17 2.95949e-17 3.53759e-17) OK.
***High aspect ratio cells found, Max aspect ratio: 1.41654e+198, number of cells 4786
Minumum face area = 4.60908e-06. Maximum face area = 0.0398707. Face area magnitudes OK.
Min volume = 1.33333e-300. Max volume = 0.00143918. Total volume = 5.5607. Cell volumes OK.
Mesh non-orthogonality Max: 179.655 average: 62.4743
*Number of severely non-orthogonal faces: 34479.
***Number of non-orthogonality errors: 9301.
***Error in face pyramids: 19112 faces are incorrectly oriented.
***Max skewness = 2617.13, 2542 highly skew faces detected which may impair the quality of the results
Failed 4 mesh geometry checks.
Failed 1 mesh checks.
ExecutionTime = 1.6 s ClockTime = 2 s

Time = 1e-05
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 void Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vect or<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::fv::limitedSnGrad<Foam::Vector<double> >::correction(Foam::GeometricField<Foam::Vector<do uble>, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfvMotionSolvers.so"
#8 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfvMotionSolvers.so"
#9 Foam::displacementLaplacianFvMotionSolver::solve() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfvMotionSolvers.so"
#10 Foam::motionSolver::newPoints() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libdynamicMesh.so"
#11 Foam::dynamicMotionSolverFvMesh::update() in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libdynamicFvMesh.so"
#12
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/moveDynamicMesh"
#13 __libc_start_main in "/lib/libc.so.6"
#14
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/moveDynamicMesh"
Floating point exception
Is this a problem with the mesh? Which are the steps that may help to get a good result?

Thanks a lot,

Samuele

Last edited by samiam1000; May 24, 2012 at 04:04.
samiam1000 is offline   Reply With Quote

Old   May 25, 2012, 02:48
Default
  #2
Member
 
Jan
Join Date: Dec 2009
Location: Berlin
Posts: 50
Rep Power: 9
SirWombat is on a distinguished road
Send a message via Skype™ to SirWombat
Hi Samuele,

there are these three asterisks "***" which tell you, that there's something terribly wrong. So your mesh features the following bad cells:

4786 High aspect ratio cells
34479 severly non-othogonal faces
19112 incorretly oriented faces
2542 highly skewed faces

so YES, this is most probably a problem with the mesh an you will need to increase your mesh quality!

I suppose you are using snappyHexMesh, so then it's a good start to find the problematic cells by running "checkMesh" after you created your mesh. checkMesh will create "sets" of those bad cells und you can view them in paraview/paraFoam.

A good start for snappyHexMesh is to switch off all quality controls, and then have a look at the bad cells with checkMesh and paraview. After that you may slowly try to rise the quality parameters in SHM und see how that effects your mesh.
I suggest to do that on a tiny part of your mesh, e.g. by using an appropriate blockMeshDict

Greets,
Jan
__________________
~~~_/)~~~
SirWombat is offline   Reply With Quote

Old   May 28, 2012, 05:59
Default
  #3
New Member
 
Giuliano Lotta
Join Date: May 2012
Posts: 12
Rep Power: 5
Giuliano69 is on a distinguished road
May I ask where can be found "documented" examples of moving mash, to learn how to configure and deal with them ?
Giuliano69 is offline   Reply With Quote

Old   May 28, 2012, 06:34
Default
  #4
kid
Senior Member
 
cfdkid
Join Date: Mar 2009
Posts: 133
Rep Power: 8
kid is on a distinguished road
Use "moveDynamicMesh --help"
This is the best help. Also try to redo any tutorial on moveDynamicMesh.
__________________
________________________________________
Regards,
CFDkid

It never gets easier You just get Better
kid is offline   Reply With Quote

Old   May 30, 2012, 08:22
Default
  #5
Member
 
Eysteinn Helgason
Join Date: Sep 2009
Location: Gothenburg, Sweden
Posts: 52
Rep Power: 7
eysteinn is on a distinguished road
Quote:
Originally Posted by Giuliano69 View Post
May I ask where can be found "documented" examples of moving mash, to learn how to configure and deal with them ?
Hi Giuliano,

There are a few tutorials where you can read about the different settings for the mesh motion, i.e. here, here, here and here.

More tutorials can be found here by changing the year from 2007 to 2011.

Hope this helps
Eysteinn
nimasam likes this.
eysteinn is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SnappyHexMesh for internal Flow vishwa OpenFOAM Native Meshers: snappyHexMesh and Others 23 August 6, 2014 03:50
How to let the mesh motion solver just solve a small region near a moving boundary? zhajingjing OpenFOAM 8 April 7, 2014 11:10
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 06:42
2D Mesh Generation Tutorial for GMSH aeroslacker Open Source Meshers: Gmsh, Netgen, CGNS, ... 12 January 19, 2012 04:52
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43


All times are GMT -4. The time now is 11:33.