
[Sponsors] 
Update of the variables after dynamic mesh motion. 

LinkBack  Thread Tools  Search this Thread  Display Modes 
December 21, 2012, 12:19 
Update of the variables after dynamic mesh motion.

#1 
New Member
Join Date: Jun 2011
Posts: 13
Rep Power: 11 
Dear all,
I have a question regarding the update of the different variables or fields after a mesh deformation. I solve a solid problem, and for a while was using the standard procedure for the deformation of the solid phases in the solid solvers offered in openFOAM (the latest I used is moveMesh from elasticNonLinSolidFoam). From my understanding, this deformation is Lagrangian, meaning that each of the mesh points is moved by the displacement U obtained at that location. Such procedure works fine for close to solid body deformations. In my case, I have quite large compressions in some regions and not in others. As a result, a purely Lagrangian formulation leads to a deterioration of the mesh quality. I thus adopted the procedure that is used for fluid phases and am using the dynamic mesh motion capabilities (more precisely the displacement Laplacian with inverse distance diffusivity with respect to the moving boundary). This yields a much better mesh, but .. I am not sure that my fields are transferred properly: The solution reached differs depending on my choice of mesh motion procedure. As an example, Fig 1 (compareOldNewMotion_T50.png) shows the deformation of a cylinder subjected to a pressure load on the inner surface, with two different mesh motion procedures (moveMesh vs. dynamic mesh). Dynamic mesh definitely yields a better mesh quality, but the displacement and pressure profiles differ … The procedure I follow for the mesh deformation is as follows: Code:
pointVectorField pointsDispl = cpi.interpolate(U); displacementSolidPatch == pointsDispl.boundaryField([solidPatchID].patchInternalField(); mesh.update(); I think I am missing something about the transfer of variables but haven't been able to find the answer in the different tutorials. Should I do an interpolation from the Lagrangian mesh to the smoother one obtained using mesh.update()? Would you have an example solver in mind that I could follow? Thank you for any help, Diane 

December 21, 2012, 15:59 

#2  
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 829
Rep Power: 28 
Quote:
You are right that the elasticNonLinULSolidFoam solver adopts a Lagrangian approach, so it assumes that there is no mass transfer in or out of each cell. So when the mesh is moved, it is moved using the solvedfor displacements meaning there is no convection (mass transfer) between the cells. But in your case when you move the mesh using the fluid mesh motion classes then there will be a flow of mass across the cell surfaces and the momentum carried by this convection needs to be accounted for. So essentially you need a solver using an Arbitrary Lagrangian Eulerian (ALE) approach instead of an/a updated/total Lagrangian approach. An ALE approach is where the mesh coincides with the material deformations only at the boundary and independent of each other elsewhere. ALE methods are often employed for modelling solids which undergo very large deformation e.g. extrusion and metal forming processes. The bad news is that there is no OpenFOAM ALE for solids (although I think it is used for interface tracking in interTrackFoam). The good news is that OpenFOAM is perfectly suited for ALE for finite volume (FV) solids as it has all the convection and mapping stuff. Also there is considerable literature on ALE methods for solids. If you do intend to develop an ALE FV solid solver then I can help you if you like. Best regards, Philip 

January 7, 2013, 02:27 

#3 
New Member
Join Date: Jun 2011
Posts: 13
Rep Power: 11 
Hi Philip,
First of all, happy new year!! Sorry for the late reply, I went on a internetbreak over the Christmas vacations! I have been reading about ALE already, and believe that this is indeed what I will have to do as my solids deform to much for me not to loose mesh quality over time. I would be delighted to discuss the matter with you when you have time. Especially mapping as it is something I am not familiar with and I guess that talking about it with you could make things go much faster! Let me know when and how you would prefer to do that. Best regards, Diane 

January 12, 2013, 16:10 

#5 
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 829
Rep Power: 28 
Hi Hisham,
I will post here when I start getting some interesting results. A nice description of ALE methods is given in the online Encyclopedia of Computational Mechanics  Chapter 14  Arbitrary LagrangianEulerian Methods. Best regards, Philip 

February 14, 2013, 14:49 

#7 
Senior Member
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 144
Rep Power: 13 
Hi,
I am working currently on a solver with a moving mesh. Actually one boundary of a domain is changing according to the solution of some other equation. I implemented the mesh changes similar to icoSructFoam solver. But to update my mesh I need to add following line: HTML Code:
mesh.update(); I am confused because my mesh is not get updated without this command, but mesh in icoSructFoam is updated... Anyway more important/confusing is that my solution, after mesh is updated, behaves somehow strange. Energy is added to the system, what is not a case if mesh is not deformed... It happend for a advectiondiffusion equation, which should just diffuse with time. It behaves like that if mesh is not changed. But if it is changed, energy is somehow added to the system and field gets stronger and stronger with time... I hope anyone will have any suggestions. ZMM 

February 14, 2013, 17:24 

#8 
Senior Member
Mieszko Młody
Join Date: Mar 2009
Location: POLAND, USA
Posts: 144
Rep Power: 13 
Hi,
HTML Code:
mesh.update(); But I still have the problem with energy being added to the system after mesh is changed. It happens for simple diffusion equation: HTML Code:
fvScalarMatrix omEqn ( dimDt*fvm::ddt(om)  fvm::laplacian(dimensionedScalar("1",dimensionSet(0, 2, 0, 0, 0),1), om) ); omEqn.solve(); But if mesh is changed "little more than a little" then "om" gets stronger and stronger with each time step... Do you have any clue why it is like that ?? thanks ZMM 

January 16, 2014, 17:28 

#9  
Member
yijin Mao
Join Date: May 2010
Location: Columbia, MO
Posts: 49
Rep Power: 12 
Quote:
I am also stick to this problem. Do you have any clue? 

January 18, 2014, 10:08 

#10  
Member
yijin Mao
Join Date: May 2010
Location: Columbia, MO
Posts: 49
Rep Power: 12 
Quote:
fvScalarMatrix omEqn ( dimDt*fvm::ddt(om)  fvm::laplacian(dimensionedScalar("1",dimensionSet( 0, 2, 0, 0, 0),1)  fvm::div(mesh.phi(), om), om) ); 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
potentialFreeSurfaceFoam with Dynamic Mesh  GuilhermeMP  OpenFOAM Programming & Development  1  October 6, 2013 03:05 
Dynamic Mesh moving interface help  akash.iitb  FLUENT  0  August 23, 2010 23:53 
Dynamic Mesh Problem.  Tom Clark  FLUENT  9  July 7, 2010 07:56 
dynamic mesh motion (angularOscillatingDisplacement) in OF1.6  tamas.baksa.hu  OpenFOAM Running, Solving & CFD  4  May 13, 2010 12:22 
Convergence moving mesh  lr103476  OpenFOAM Running, Solving & CFD  30  November 19, 2007 14:09 