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

InterDyMFoam dynamic refinement

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

Like Tree1Likes
  • 1 Post By olwi

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2009, 09:49
Default hallo all, i am student and n
  #1
ala
New Member
 
Alicja M
Join Date: Mar 2009
Location: Erlangen, DE
Posts: 26
Rep Power: 8
ala is on a distinguished road
hallo all,
i am student and new at/in (?) OpenFOAM

i am just trying my first case in OpenFOAM
(surface tension of a water , in a square, with contact angle with and
without gravity)
solver: interFoam

so, it's all really pretty
at yet i just want to deal with interDyMFoam and dynamic refinement.
in 3D its all ok,
but i have no idea how can i dynamic refine in 2D

i get following error:


Starting time loop

Courant Number mean: 0 max: 0
deltaT = 0.001
Time = 0.001

Selected 56 cells for refinement out of 5400.
Refined from 5400 to 5792 cells.
Selected 0 split points out of a possible 56.
Execution time for mesh.update() = 0.18 s


This mesh contains patches of type empty but is not 1D or 2D
by virtue of the fact that the number of faces of this
empty patch is not divisible by the number of cells.

From function emptyFvPatchField<type>::updateCoeffs()
in file fields/fvPatchFields/constraint/empty/emptyFvPatchField.C at
line 148.

FOAM exiting

Thanks in advance. Greetings, Ala
ala is offline   Reply With Quote

Old   January 30, 2009, 10:39
Default Hello Ala and welcome, The ab
  #2
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
Hello Ala and welcome,
The above message says that you have few boundary surfaces set as "empty", but your domain is not one cell thick. You will have to change the "empty" to something else like "wall".
Also, try to use the search feature of this forum, since there are many people that encountered the same problem.

I hope this is helpful,
Dragos
dmoroian is offline   Reply With Quote

Old   January 31, 2009, 09:56
Default Hello Dragos, first of all
  #3
ala
New Member
 
Alicja M
Join Date: Mar 2009
Location: Erlangen, DE
Posts: 26
Rep Power: 8
ala is on a distinguished road
Hello Dragos,

first of all thank you very much for your reply.

I followed your approach and changed the empty boundary surfaces to wall.
Unfortunately it seems that the changes didn't have an effect.
InterDyMFoam still refined my mesh in 3D.
I have a really simple mesh. See my blockMeshDict:

vertices
(
( 0 0 0 )
( 2 0 0 )
( 2 4 0 )
( 0 4 0 )
( 0 0 1 )
( 2 0 1 )
( 2 4 1 )
( 0 4 1 )
);

blocks
(
hex ( 0 1 2 3 4 5 6 7 ) ( 16 32 1 ) simpleGrading ( 1 1 1 )
);

edges
(
);

patches
(
wall leftWall
(
( 0 3 7 4 )
)

wall rightWall
(
( 1 2 6 5 )
)

wall lowerWall
(
( 0 4 5 1 )
)

wall emptyWalls
(
( 0 1 2 3 )
( 4 5 6 7 )
)

patch atmosphere
(
( 3 2 6 7 )
)
);

mergePatchPairs
(
);


At the 0 directory I have changed the patch emptyWall as follows :

gamma -> zeroGradient, pd -> zeroGradient, U -> fixedValue, value ( 0 0 0 ).

After a few seconds of simulation InterDyMFoam still refined the mesh in 3D, also at walls with one cell thickness. How can I do it, that InterDyMFoam just refines the mesh in 2D? And not 3D! (Cmp. the pictures below).

I also used the search of the OpenFOAM forum but I did not find something that helped on! Any hints are welcome!

Thanks,
Ala

Mesh before: and the mesh afterwards:
ala is offline   Reply With Quote

Old   January 31, 2009, 10:08
Default Sorry, there was an error whil
  #4
ala
New Member
 
Alicja M
Join Date: Mar 2009
Location: Erlangen, DE
Posts: 26
Rep Power: 8
ala is on a distinguished road
Sorry, there was an error while uploading the pictures from my last post. Here are the two pictures:
Mesh before: and the mesh afterwards:
ala is offline   Reply With Quote

Old   February 1, 2009, 17:25
Default Hi Ala, The dynamic mesh cl
  #5
Member
 
Ola Widlund
Join Date: Mar 2009
Location: Sweden
Posts: 87
Rep Power: 8
olwi is on a distinguished road
Hi Ala,

The dynamic mesh class used in interDyMFoam handles only refinement in 3D. When a hexahedron is refined, it is split in eight. Always. Note that the same hex splitter is used in the snappyHexMesh mesh generator, which consequently will always give you a 3D mesh.

To make the dynamic mesh class refine/coarsen in 2D, someone would have to develop a new algorithm which does the cell and edge splitting only in two directions. That is certainly possible, but I have not heard of anyone spending time on it. I think it takes quite some experience and knowledge to do it. And time. Or finance Icon or OpenCFD to do it.

Depending on what you want to do, you might look at the (older?) approach of the refineMesh utility. It works in 2D as well, but it is not meant to run during the simulation. Maybe you could hack a solver of your own to do the refinement action between time steps, but I don't think this is a trivial task, either.

Good luck!

Ola
banji likes this.
olwi is offline   Reply With Quote

Old   June 23, 2009, 05:11
Default 2D mesh refinement
  #6
New Member
 
Luther Terblanche
Join Date: May 2009
Posts: 15
Rep Power: 8
luther is on a distinguished road
Hi All,

Has this issue been resolved by someone? I would also like to be able to run Dynamic mesh refinement in a 2D case.

Please let me know if anyone has been able to accomplish this.

Regards
Luther
luther is offline   Reply With Quote

Old   June 23, 2009, 05:36
Default empty to symmetry
  #7
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 9
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
Hi

i also face the same problem .
so i set the empty to symmetry plane.

wayne
waynezw0618 is offline   Reply With Quote

Old   June 23, 2009, 05:40
Default
  #8
New Member
 
Luther Terblanche
Join Date: May 2009
Posts: 15
Rep Power: 8
luther is on a distinguished road
hi Wayne,

Thank you. Will try it.

Do you perhaps know what the implications is on computation time changing empty to symmetry?

Regards
Luther
luther is offline   Reply With Quote

Old   June 23, 2009, 05:55
Default
  #9
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 307
Rep Power: 9
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to waynezw0618
I am sorry i don`t..i try this for we do 2D case in CFX is setting this face to symmetry and i think it will restrict something.but it also split the mesh in these direction


Quote:
Originally Posted by luther View Post
hi Wayne,

Thank you. Will try it.

Do you perhaps know what the implications is on computation time changing empty to symmetry?

Regards
Luther
waynezw0618 is offline   Reply With Quote

Old   April 11, 2010, 13:40
Default Problem with dynamicRefineFvMesh
  #10
New Member
 
Join Date: Nov 2009
Posts: 4
Rep Power: 7
Brunno is on a distinguished road
Hello,
I am also trying to employ the dynamic mesh refinement tool. It works fine for any mesh generated by blockMesh; whenever I try to run a case whose mesh was created by snappyHM I get the following error:

Only call if constructed with history capability#0 Foam::error:rintStack(Foam::Ostream&) in "/home/brunno/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/brunno/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::hexRef8::getSplitPoints() const in "/home/brunno/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicMesh.so"
#3 Foam::dynamicRefineFvMesh::selectUnrefinePoints(do uble, Foam::PackedList<1u> const&, Foam::Field<double> const&) const in "/home/brunno/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicFvMesh.so"
#4 Foam::dynamicRefineFvMesh::update() in "/home/brunno/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libdynamicFvMesh.so"
#5 main in "/home/brunno/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/interDyMFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116


From function hexRef8::getSplitPoints()
in file polyTopoChange/polyTopoChange/hexRef8.C at line 4873.

FOAM aborting

Aborted

Has anyone faced this problem before ?
Brunno is offline   Reply With Quote

Old   June 28, 2010, 06:17
Default
  #11
New Member
 
Victor Fleischer
Join Date: Nov 2009
Posts: 21
Rep Power: 7
Victor is on a distinguished road
Hi at all,
I'm also looking for dynamicrefineMesh for 2D!
Does anybody have a solution for this? I tried to understand the code but its dufficult for me....
Thanks a lot in advance!
Victor
Victor is offline   Reply With Quote

Reply

Tags
112

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
InterDyMFoam SDA SKA coefficients markc OpenFOAM Running, Solving & CFD 12 June 20, 2015 15:53
InterDyMFoam dynamic messing in parallel fails under nonquiescent conditions adona058 OpenFOAM Running, Solving & CFD 5 August 19, 2010 11:47
Choise of ddtSchemes regarding Dynamic Mesh Refinement TaylorTayleighInstability sega OpenFOAM Running, Solving & CFD 0 March 9, 2009 10:34
Dynamic mesh refinement in a 2D or axissymmetric case christian_kunkelmann OpenFOAM Running, Solving & CFD 9 January 19, 2009 13:15
InterDyMFoam dynamic meshing in parallel fails under nonquiescent conditions adona058 OpenFOAM Bugs 7 November 18, 2008 15:58


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