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

Dynamic moving mesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 6, 2005, 10:12
Default My automatic mesh motion solve
  #21
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
My automatic mesh motion solver does not need the crank angle - have a look at the icoFoamAutoMotion or moveMesh applications on how to set up a moving mesh problem. Once you specify the motion of the boundary, the rest will be done automatically.

As for the flow solver, you do not have any topological changes (right?), so any compressible code with mesh motion should be a good starting point.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   June 6, 2005, 12:25
Default Can I apply moveMesh utilities
  #22
albi80
Guest
 
Posts: n/a
Can I apply moveMesh utilities at the oodles solver too??
  Reply With Quote

Old   June 6, 2005, 12:27
Default Is your flow compressible or i
  #23
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
Is your flow compressible or incompressible?
henry is offline   Reply With Quote

Old   June 6, 2005, 12:33
Default I can consider incompressible
  #24
albi80
Guest
 
Posts: n/a
I can consider incompressible ...
  Reply With Quote

Old   June 6, 2005, 12:38
Default In which case you can either a
  #25
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
In which case you can either add the LES model construction and calls from oodels to icoFoamAutoMotion or add the mesh-motion calls from icoFoamAutoMotion to oodles. You might also want to consider how you want to move the mesh; if you only want to move in by extension and contraction you might find the mesh-motion solver over-kill and it might be better for you to move the points directly and call mesh.movePoints with the new point locations you compute yourself.
henry is offline   Reply With Quote

Old   June 8, 2005, 10:23
Default Hi, I'm always me .... I add
  #26
albi80
Guest
 
Posts: n/a
Hi,
I'm always me .... I add the mesh-motions calls and costructions of icoFoamAutoMotion to oodles and I simulate a test "in a box" with a moving wall. For the moment it's ok, the cell move in extension, but is it possible to create other cell or destroy it (rezone ...???) ?
Thanks
  Reply With Quote

Old   June 8, 2005, 10:32
Default Take a look at icoTopoFoam.
  #27
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
Take a look at icoTopoFoam.
henry is offline   Reply With Quote

Old   June 8, 2005, 10:59
Default Now I've tried with the tutori
  #28
albi80
Guest
 
Posts: n/a
Now I've tried with the tutorial for icoTopoFoam, blockMesh ok, but when I launch the simulation (icoTopoFoam . mixer2D), it stops and give me Segmentation fault... For what is it icoTopoFoam ? What are and how do I set the file cellToRegion, mixerProperties ecc.. ??
  Reply With Quote

Old   June 8, 2005, 11:06
Default In order to do cell addition/r
  #29
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
In order to do cell addition/removal, sliding interface, you need my dynamic mesh stuff. The details have been described in one of the threads in this discussion group in detail and Pei has promised to make a tutorial from the examples I've sent him a while back (any news?).

As for the functionality, you can take a look at the pictures on my private web site:

http://www.h.jasak.dial.pipex.com under FOAM development or

http://www.foamcfd.org (same thing, I am migrating the contents).

This is not trivial and for a better version with new stuff, some bug fixes and more pre-implemented functionality you'll need a new version of FOAM.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   June 8, 2005, 11:06
Default There have already been lots o
  #30
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
There have already been lots of discussions about running with topology change, you should be able to find the threads using the search facility.
henry is offline   Reply With Quote

Old   June 8, 2005, 11:33
Default Thank Henry and Hrv, probabl
  #31
albi80
Guest
 
Posts: n/a
Thank Henry and Hrv,
probably for the moment it's enough what I've done (I have spoken five minutes ago with the professor......) so now I try oodles with the call of mesh-motion..... and probably I disturb you later....
Bye
Alberto
  Reply With Quote

Old   June 9, 2005, 05:37
Default Hi , Anyone can help me ? Ho
  #32
albi80
Guest
 
Posts: n/a
Hi ,
Anyone can help me ?
How can I set for my movingWall (in motionU file) a motion law (for example a piston motion, simple harmonic) ?

Thanks
Alberto
  Reply With Quote

Old   June 9, 2005, 05:43
Default Take a look at engineFoam and
  #33
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
Take a look at engineFoam and the mesh-motion used by it in OpenFOAM-1.1/src/engine/include/movePiston.H

You will probably find the simple approach for layered meshes at the top of this file adequate but if your mesh is not layered then you will have to use the more complex and expensive mesh-motion solver.
henry is offline   Reply With Quote

Old   June 9, 2005, 07:22
Default But isn't it possible to set d
  #34
albi80
Guest
 
Posts: n/a
But isn't it possible to set directly in the Umotion file, instead "fixedValue" for movingWall, something like ... "variableValue" and I can specify the exact law motion ?

Thanks Alberto
  Reply With Quote

Old   June 9, 2005, 13:09
Default In order to do that, you will
  #35
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
In order to do that, you will need your own boundary condition type - the function that described the motion needs to exist is C++ and only the parameters can be read in. Search for my oscillating fixed value boundary condition, posted on this group a while back.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   June 10, 2005, 05:20
Default So -I should add #include
  #36
albi80
Guest
 
Posts: n/a
So
-I should add
#include "setBoundaryMotion.H"

in the moveMesh.C file.
-I create the setBoundaryMotion.H file, I compile it ... and... it's ok?

The amplitude and frequency I can set directly in the setBoundaryMotion.H file?

What should I set in the motionU file in the /0 directory?
  Reply With Quote

Old   June 10, 2005, 08:27
Default @Hrvs posting from the 8th:
  #37
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
@Hrvs posting from the 8th:

Some days ago I copied the descriptions Hrv was refering to to the Wiki (with some very minor editing). It's the first How-To on the Wiki that really deserves that name:

http://openfoamwiki.net/index.php/HowTo_setting_u p_dynamic_mesh_cases

(I know it's a bit off topic here, but I wanted to say it)
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   June 10, 2005, 10:30
Default Thanks a lot Bernhard, but in
  #38
albi80
Guest
 
Posts: n/a
Thanks a lot Bernhard,
but in the page that you give me I haven't found any example of "oscillating fixed value boundary condition" that Hrv has posted February 19. Do you know how to "create" step by step the boundary oscillating of which Hrv told me??
Thanks Alberto.
  Reply With Quote

Old   June 14, 2005, 06:45
Default Sorry Alberto! Personally I
  #39
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Sorry Alberto!

Personally I don't have any experience with dynamic meshes. I just copied Hrvs explanation to the Wiki because I thought it was an excellent start for a How-To (The original thread is referenced at the top of the page, there are some examples there - don't know if they apply to your problem - which I havn't edited to the Wiki yet; I plan to do that when I have some time, but I wouldn't mind if somebody who actually tried them did that)
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   December 1, 2006, 17:50
Default I have tried to implement the
  #40
Member
 
stefan
Join Date: Mar 2009
Posts: 96
Rep Power: 8
stefanke is on a distinguished road
I have tried to implement the layer addition/removal things into an existing engine solver. So far in the compression stroke all works fine but as soon as there is a layer addition in the expansion stroke some strange things happen:

...
Added all cells. Final cell count = 600

Added zone points: untouched = 0() modified = 0() added = 0(). Points per zone = 0()
Added zone faces: untouched = 1(0) modified = 1(100) added = 1(0). Faces per zone = 1(100)
Added zone cells: untouched = 0() modified = 0() added = 0(). Cells per zone = 0()
void polyMesh::clearGeom() : clearing geometric data
void polyMesh::clearAddressing() : clearing topology
void polyMesh::setInstance(const fileName& inst) : Resetting file instance to "69"
void polyMesh::initMesh() : calculating faceCells
Foam::polyMesh::morph
(
const polyTopoChange& ref
) : completed topological change.


From function void polyTopoChanger::update(const mapPolyMesh& m)
in file polyTopoChange/polyTopoChanger/polyTopoChanger.C at line 247
Updating mesh modifiers
layerAdditionRemoval::updateMesh(const mapPolyMesh&) for object pistonLayer : Clearing addressing on external request. Pointers not set.
void polyMesh::setInstance(const fileName& inst) : Resetting file instance to "69"
Mapping field meshPhi
Mapping field meshPhi_0
Piston layering mode
tmp<scalarfield> polyMesh::movePoints(const pointField&) : Moving points for time 0.00575 index 249
bool primitiveMesh::checkMeshMotion(const pointField& newPoints, const bool report) const: checking mesh motion
Min volume = 3.90451e-08. Total volume = 2.63342e-05. Cell volumes OK.
Zero or negative face area detected for internal face 1400 between cells 500 and 501. Face area magnitude = 0
--> FOAM Warning :
From function primitiveMesh::checkMeshMotion(const pointField& newPoints, const bool report) const
in file meshes/primitiveMesh/primitiveMeshCheckMotion.C at line 200
Severe non-orthogonality in mesh motion for face 1400 between cells 500 and 501: Angle = 90 deg.
Zero or negative face area detected for internal face 1401 between cells 500 and 510. Face area magnitude = 0
...

Has anyone an idea what is going wrong here. The strange thing is that the layer removement in the compression stroke works perfectly !?!
stefanke is offline   Reply With Quote

Reply

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
moving (dynamic) mesh gives fluctuating pressure daniel FLUENT 1 July 31, 2010 14:32
dynamic meshing, rotating moving mesh...... madhan FLUENT 1 January 24, 2007 00:55
Moving Dynamic Mesh in Fluent 6.2 R.M.Bharath Somayaji FLUENT 0 January 9, 2006 21:36
Moving(dynamic) mesh in Fluent yonghyun FLUENT 1 January 25, 2005 15:50
moving solid in dynamic mesh James FLUENT 1 October 4, 2002 07:22


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