CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   simpleDyMFoam on v1806 (https://www.cfd-online.com/Forums/openfoam-programming-development/210954-simpledymfoam-v1806.html)

gamemakerh November 8, 2018 09:15

simpleDyMFoam on v1806
 
2 Attachment(s)
Hey guys,

I'm having problems with adaptive mesh refinement integration into the OpenFOAM v1806.


The aim of the project is relatively simple - execute Steady-State, 2D thermal (buoyantBousinessqSimpleFoam) simulation with Adaptive Mesh Refinement (AMR).

I managed to compile 2D mesh refinement addon which was created by Luca Cornolti. Big Clap for him !!!

Unfortunately, his solver was working only with transient simulations (interFoam etc.)

I found simpleDyMFoam solver created by some Japanese guy named "Penguin" (unfortunately the explanation on how to develop and compile this solver is no longer available on his website). I was unable to compile this solver for OpenFOAM-v1806.

I've tried to include more header files according to the original simpleFoam solver, but the compilation was unsuccessful.

Maybe someone has working Steady-State solver which are able to work using 2D adaptive mesh refinement?

Please check attached files below.



*******
UPDATE
*******
After some adjustments I managed to to launch simpleDyMFoam using 2D Adaptive Mesh Refinement. But I encountered the problem regarding the unrefinement. maybe someone interested into deeper look at this problem ? I could upload working simpleDyMFoam solver for v1806.

I'm getting this error log

Code:

Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicRefineFvMesh2D
Detected 438 cells that are protected from refinement. Writing these to cellSet protectedCells.

SIMPLE: convergence criteria
    field U      tolerance 0.001
    field "p.*"  tolerance 0.001
    field T      tolerance 0.001
    field k      tolerance 0.001
    field omega  tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

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

No finite volume options present

Starting time loop

Reading surface description:
    yNormal

streamLine streamLines:
    Employing velocity field U
    automatic track length specified through number of sub cycles : 1

runTimePostProcessing postPro1: reading post-processing data
Time = 1

Selected 13294 cells for refinement out of 13440.
Refined from 13440 to 53322 cells.
streamLine streamLines:
    Employing velocity field U
    automatic track length specified through number of sub cycles : 1

Selected 0 split edges out of a possible 13294.
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.005390501, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.006547034, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.9997273, Final residual = 0.006281913, No Iterations 16
GAMG:  Solving for p, Initial residual = 0.01727778, Final residual = 0.0001585316, No Iterations 2
time step continuity errors : sum local = 2.528902e-05, global = 1.842507e-08, cumulative = 1.842507e-08
ExecutionTime = 7.11 s  ClockTime = 8 s

Time = 2

Selected 51713 cells for refinement out of 53322.
Refined from 53322 to 208461 cells.
streamLine streamLines:
    Employing velocity field U
    automatic track length specified through number of sub cycles : 1

Selected 0 split edges out of a possible 51954.
Setting residual field for first solver iteration for solver field: Ux
smoothSolver:  Solving for Ux, Initial residual = 0.2805062, Final residual = 0.001492275, No Iterations 2
Setting residual field for first solver iteration for solver field: Uy
smoothSolver:  Solving for Uy, Initial residual = 0.2748859, Final residual = 0.001400735, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.6047205, Final residual = 0.005748263, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.004981731, Final residual = 4.600868e-05, No Iterations 9
time step continuity errors : sum local = 6.968632e-05, global = 5.084545e-09, cumulative = 2.350962e-08
ExecutionTime = 17.19 s  ClockTime = 22 s

Time = 3

Selected 0 split edges out of a possible 51954.
Setting residual field for first solver iteration for solver field: Ux
smoothSolver:  Solving for Ux, Initial residual = 0.2000615, Final residual = 0.0007128696, No Iterations 3
Setting residual field for first solver iteration for solver field: Uy
smoothSolver:  Solving for Uy, Initial residual = 0.2088684, Final residual = 0.0005146371, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.04054855, Final residual = 0.0002980836, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.006104533, Final residual = 4.91343e-05, No Iterations 4
time step continuity errors : sum local = 0.0001698377, global = -6.749423e-10, cumulative = 2.283467e-08
ExecutionTime = 21.99 s  ClockTime = 27 s

Time = 4

Selected 3 split edges out of a possible 51954.


--> FOAM FATAL ERROR:
Problem : edge has too few face neighbours:1(271578)
edge:13241 vertices:(27584 27592) coords:(-39.15004 -14.59622 0)(-39.36746 -14.70492 0.5)

    From function void Foam::removeFaces::setRefinement(const labelList&, const labelList&, const labelList&, Foam::polyTopoChange&) const
    in file polyTopoChange/polyTopoChange/removeFaces.C at line 867.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::removeFaces::setRefinement(Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::polyTopoChange&) const at ??:?
#3  Foam::hexRef2D::setUnrefinement(Foam::List<int> const&, Foam::polyTopoChange&) at ??:?
#4  Foam::dynamicRefineFvMesh2D::unrefine(Foam::List<int> const&) at ??:?
#5  Foam::dynamicRefineFvMesh2D::update() at ??:?
#6  ? at ??:?
#7  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8  ? at ??:?
Aborted (core dumped)



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