CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Troubleshooting errors with pimpleDyMFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 27, 2017, 11:46
Default Troubleshooting errors with pimpleDyMFoam
New Member
Join Date: May 2016
Posts: 17
Rep Power: 9
j_moulton is on a distinguished road

I am attempting to run a pimpleDyMFoam using displacementLaplacian. I have had success using velocityLaplacian, but it is very sensitive, and either prone to blowing up or the change is too small to be noticed. My case requires a specific % change, so displacementLaplacian seems natural.

Here is my mesh:

My goal is to expand outward the outside walls of the circular sac on the right

I am stuck with a float point exception, which occurs in the first pressure update (I think). The log file is here:
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: displacementLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present

Reading/calculating face velocity Uf

No finite volume options present

Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Time = 1e-05

DICPCG:  Solving for cellDisplacementx, Initial residual = 1, Final residual = 0.000870402, No Iterations 10
DICPCG:  Solving for cellDisplacementy, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 9.83755e-05, No Iterations 63
GAMG:  Solving for pcorr, Initial residual = 0.999308, Final residual = 9.47909e-05, No Iterations 55
GAMG:  Solving for pcorr, Initial residual = 0.999823, Final residual = 9.44465e-05, No Iterations 38
time step continuity errors : sum local = -23977.2, global = 49.658, cumulative = 49.658
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0937035, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0848355, No Iterations 5
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/"
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib/x86_64-linux-gnu/"
#10  ? at ??:?
Does anyone out there know where to start with tackling this issue? Is there a better choice of solver?

My full case is attached. Please help! This is for my thesis and it is the last step. The results are acceptable with velocityLaplacian, but as I mentioned, it either blows up or the change is hardly visible (though flow is clearly occuring).
Attached Files
File Type: zip (57.2 KB, 2 views)
j_moulton is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Building OpenFOAM1.7.0 from source ata OpenFOAM Installation 46 March 6, 2022 13:21
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 18:07

All times are GMT -4. The time now is 00:07.