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

[Technical] Mesh movement during runtime

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

Like Tree1Likes
  • 1 Post By Hgholami

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 23, 2013, 11:10
Default Mesh movement during runtime
  #1
New Member
 
Matthias Stammen
Join Date: Oct 2010
Posts: 8
Rep Power: 12
MatzeS is on a distinguished road
Hello there,

although I already was searching here and there, I am still not sure how to use the diffusivity of the velocityLaplacian solver in the dynamicMeshDict.

I want to move my mesh during runtime depending on the pressure on the patch boundary "airfoil" (optimization of a 2D-airfoil).
But somehow the diffusivity (uniform, quadratic, motionDirectional etc.) does not change much the results and I always get some overlapping cells (see attachment). If I choose different values for uniform I also don't see any change in the result.

Here is the dynamicMeshDict I am using:
HTML Code:
dynamicFvMesh           dynamicMotionSolverFvMesh;

motionSolverLibs       ("libDynamicFvMesh.so");

solver                  velocityLaplacian;

diffusivity             quadratic uniform; //uniform ; //(100);

distancePatches         (airfoil);
Could anybody tell me what I am doing wrong?

Kind regards,
Matthias
Attached Images
File Type: jpg t925_diffusivity_uniform_1000.jpg (73.0 KB, 75 views)
MatzeS is offline   Reply With Quote

Old   May 23, 2013, 12:00
Default
  #2
Senior Member
 
mturcios777's Avatar
 
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 728
Rep Power: 25
mturcios777 will become famous soon enough
I'm not sure you can get the behaviour you want with the standard motion solvers. I'm not sure how you are coupling the pressure on the patch with the motionSolver. velocityLaplacian means that the displacement of the points is computed with velocity specified as initial and boundary conditions of the laplacian equation.

The diffusivity for the points mainly affects how "stiff" the mesh is. The different options are there to modify how diffusivity is computed based on several criteria. I know there is one called inverseDistance that has the diffusivity be proportional to the inverse distance from certain patches (so its highest near the patches and then falls off as you move further). I think you will need some kind of distance dependant diffusivity to help things move along.
mturcios777 is offline   Reply With Quote

Old   May 23, 2013, 15:03
Default
  #3
New Member
 
Matthias Stammen
Join Date: Oct 2010
Posts: 8
Rep Power: 12
MatzeS is on a distinguished road
Thanks, mturcios.
The coupling of pressure and motionSolver I am doing via the pointMotionU.

The diffusivity inverseDistance I also tried, but the results are just slightly different and I still get overlapping cells.
But is it normal that different values, e.g. in motionDirectional, give no difference in the result? So, is my dynamicMeshDict correct in principle?
MatzeS is offline   Reply With Quote

Old   May 23, 2013, 15:52
Default
  #4
Senior Member
 
mturcios777's Avatar
 
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 728
Rep Power: 25
mturcios777 will become famous soon enough
Quote:
Originally Posted by MatzeS View Post
Thanks, mturcios.
The coupling of pressure and motionSolver I am doing via the pointMotionU.

The diffusivity inverseDistance I also tried, but the results are just slightly different and I still get overlapping cells.
But is it normal that different values, e.g. in motionDirectional, give no difference in the result? So, is my dynamicMeshDict correct in principle?
The diffusivity is mainly affecting the stiffness of your system of points and how long/much the points are moved in the motionSolver step. You may want to try write the diffusivity to a pointField see what the absolute numbers are. There is nothing in the motionSolver about mesh quality, so you may want to check the mesh every few iterations and remesh when the quality is below a certain level. Good luck!
mturcios777 is offline   Reply With Quote

Old   August 16, 2020, 02:07
Default
  #5
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 142
Rep Power: 4
Hgholami is on a distinguished road
Dear foamer
I use dynamic mesh for "Flow-induced deformation of a flexible thin structure as manifestation of heat transfer enhancement" article. the geometry is same as turek-hron FSI benchmark.
the quadratic inverseDistance uses 1/L^2 relation of laplacian equation of mesh motion, but when deflection amplitude increases the mesh deform near the moving boundary and make it distorted as below.

I looking for a better diffusion parameter for this case. I think some mesh layers of the moving boundary should move with boundary (for example 3layers) and then diffusion spread in mesh domain to prevent distortion of mesh in first mesh layer of moving boundary. do you have any suggestion for any formulation to change the code base it?
Thanks.
salachnaj likes this.
Hgholami is offline   Reply With Quote

Old   February 3, 2021, 04:23
Default
  #6
New Member
 
Conor
Join Date: Oct 2016
Posts: 14
Rep Power: 6
ConorMD is on a distinguished road
Hi Hgholami,

I am experiencing the same issue as you when using the displacementLaplacian solver with a cantilever foil.

Have you been able to find a fix for this?

Regards,
Conor
ConorMD is offline   Reply With Quote

Old   February 3, 2021, 07:14
Default
  #7
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 142
Rep Power: 4
Hgholami is on a distinguished road
Hi Conor
I used velocity solver prefer than displacement solver. At now, the refVelocityLaplacian is prefer than other. Also "diffusivity quadratic inverseDistance" is good. But still it can't support large deflection. I also modified "quadratic" in motionDiffusivity for compatibility with geometry but still it is not sufficient. May you can use overset in foam-extend4.1. the solids4Foam for FSI problems use it.

Quote:
Originally Posted by ConorMD View Post
Hi Hgholami,

I am experiencing the same issue as you when using the displacementLaplacian solver with a cantilever foil.

Have you been able to find a fix for this?

Regards,
Conor
Hgholami is offline   Reply With Quote

Reply

Tags
displacementlaplacian, mesh movement, moving meshes, velocitylaplacian

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
sliding mesh problem in CFX Saima CFX 46 September 11, 2021 07:38
[snappyHexMesh] snappyHexMesh does not create any mesh except one for the reference cell Arman_N OpenFOAM Meshing & Mesh Conversion 1 May 20, 2019 17:16
[snappyHexMesh] Snappyhex mesh: poor inlet mesh Swagga5aur OpenFOAM Meshing & Mesh Conversion 1 December 3, 2016 16:59
[snappyHexMesh] SnappyHexMesh no layers and no decent mesh for complex geometry pizzaspinate OpenFOAM Meshing & Mesh Conversion 1 February 25, 2015 07:05
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 11:55


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