[j_moulton]

Hi,

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/libc.so.6"
#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/libc.so.6"
#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).