CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

simpleDyMFoam on v1806

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 8, 2018, 08:15
Post simpleDyMFoam on v1806
  #1
New Member
 
Join Date: Sep 2014
Posts: 9
Rep Power: 8
gamemakerh is on a distinguished road
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)
Attached Files
File Type: zip simpleDyMFoam.zip (3.4 KB, 13 views)
File Type: zip simpleDyMFoam_1806_Test.zip (3.1 KB, 9 views)

Last edited by gamemakerh; November 9, 2018 at 00:59.
gamemakerh is offline   Reply With Quote

Reply

Tags
adaptive mesh refinement, amr, dynamic mesh adaption, solver, steady state

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
[OpenFOAM.com] Problem installing v1806 jormollo OpenFOAM Installation 2 October 10, 2018 09:48
[OpenFOAM.com] compile v1806 on Ubuntu 1804? phsieh2005 OpenFOAM Installation 6 August 23, 2018 12:38
v1806 soloshingTank3D6DoF error when running in parallel phsieh2005 OpenFOAM 0 July 15, 2018 12:31
[Commercial meshers] v1806 ccm utility shixiangyu OpenFOAM Meshing & Mesh Conversion 2 July 11, 2018 20:49
[OpenFOAM.com] v1806: Problems compiling on Ubuntu via WSL dela OpenFOAM Installation 1 July 10, 2018 02:45


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