# Mesh movement during runtime

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

May 23, 2013, 11:10
Mesh movement during runtime
#1
New Member

Matthias Stammen
Join Date: Oct 2010
Posts: 8
Rep Power: 7
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
 t925_diffusivity_uniform_1000.jpg (73.0 KB, 29 views)

 May 23, 2013, 12:00 #2 Senior Member     Marco A. Turcios Join Date: Mar 2009 Location: Vancouver, BC, Canada Posts: 727 Rep Power: 20 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.

 May 23, 2013, 15:03 #3 New Member   Matthias Stammen Join Date: Oct 2010 Posts: 8 Rep Power: 7 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?

May 23, 2013, 15:52
#4
Senior Member

Marco A. Turcios
Join Date: Mar 2009
Posts: 727
Rep Power: 20
Quote:
 Originally Posted by MatzeS 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!

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post siw ANSYS Meshing & Geometry 4 September 3, 2014 05:25 froztbear ANSYS Meshing & Geometry 4 November 10, 2011 09:52 Zymon enGrid 31 August 29, 2011 13:40 abouziar Main CFD Forum 1 May 30, 2008 04:08 SSL FLUENT 2 January 26, 2008 12:55

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