CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Computing the new mesh in Lagrangian co-ordinates

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 21, 2023, 13:23
Default Computing the new mesh in Lagrangian co-ordinates
  #1
Senior Member
 
Matthew
Join Date: Mar 2022
Location: United Kingdom
Posts: 184
Rep Power: 4
hunt_mat is on a distinguished road
Suppose I am working in a 2D compressible elasticity. I am using Lagrangian co-ordinates (X,Y), and I want to compute the new shape. How would I go about it? Naively I would solve the equations:
\frac{\partial x}{\partial t}=u
\frac{\partial y}{\partial t}=v

I did that in 1D but I got errors, but I'm not sure about the parameters I used.
hunt_mat is offline   Reply With Quote

Old   November 21, 2023, 15:15
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by hunt_mat View Post
Suppose I am working in a 2D compressible elasticity. I am using Lagrangian co-ordinates (X,Y), and I want to compute the new shape. How would I go about it? Naively I would solve the equations:
\frac{\partial x}{\partial t}=u
\frac{\partial y}{\partial t}=v

I did that in 1D but I got errors, but I'm not sure about the parameters I used.



what do you mean? in 1D you have only one equation but you should be aware about the global character of your problem.

For example, the Burgers equation in lagrangian form is Du/Dt=0 but you can create the evolution of the regular lagrangian grid only if the characteristic lines dx/dt=u do not intersect.

In other words until the u is regular.
FMDenaro is offline   Reply With Quote

Old   November 22, 2023, 06:26
Default
  #3
Senior Member
 
Matthew
Join Date: Mar 2022
Location: United Kingdom
Posts: 184
Rep Power: 4
hunt_mat is on a distinguished road
I use the amount of mass as my "length" co-ordinate, and as such I have this magical property that the mass in the cell is constant for all time. I compute the density at each time-step, and mass=density*length, to find the new length, you just need to compute length=mass of cell/new density.

I don't have a hyperbolic system, so I don't have characteristics to compute with.
hunt_mat is offline   Reply With Quote

Old   November 22, 2023, 08:34
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by hunt_mat View Post
I use the amount of mass as my "length" co-ordinate, and as such I have this magical property that the mass in the cell is constant for all time. I compute the density at each time-step, and mass=density*length, to find the new length, you just need to compute length=mass of cell/new density.

I don't have a hyperbolic system, so I don't have characteristics to compute with.
Could you better detail the problem you see during the computation?
FMDenaro is offline   Reply With Quote

Old   November 22, 2023, 08:53
Default
  #5
Senior Member
 
Matthew
Join Date: Mar 2022
Location: United Kingdom
Posts: 184
Rep Power: 4
hunt_mat is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Could you better detail the problem you see during the computation?
I have a 2D elastic system with diffusive terms that cause the shape to be deformed because I have an external heat source that is causing it to shrink. I believe I have the correct equations describing the process but I don't know how to compute the shape as it deforms as this changes the grid.
hunt_mat is offline   Reply With Quote

Old   November 22, 2023, 09:03
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by hunt_mat View Post
I have a 2D elastic system with diffusive terms that cause the shape to be deformed because I have an external heat source that is causing it to shrink. I believe I have the correct equations describing the process but I don't know how to compute the shape as it deforms as this changes the grid.
I think you can not use a fully Lagrangian grid. At some time you need to remapping it and restart.
Deformation will be high to have a stable computation. I don’use lagariangian grid for several issues.
FMDenaro is offline   Reply With Quote

Old   November 22, 2023, 09:23
Default
  #7
Senior Member
 
Matthew
Join Date: Mar 2022
Location: United Kingdom
Posts: 184
Rep Power: 4
hunt_mat is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I think you can not use a fully Lagrangian grid. At some time you need to remapping it and restart.
Deformation will be high to have a stable computation. I don’use lagariangian grid for several issues.
It's the remapping aspect that I need to understand. I haven't seen any references for it.

Are you referring to ALE at all?
hunt_mat is offline   Reply With Quote

Old   November 22, 2023, 10:32
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by hunt_mat View Post
It's the remapping aspect that I need to understand. I haven't seen any references for it.

Are you referring to ALE at all?
ALE is a possible way. But you could use also the way presented by Peskin in his immersed boundary method.
FMDenaro is offline   Reply With Quote

Old   November 22, 2023, 10:39
Default
  #9
Senior Member
 
Matthew
Join Date: Mar 2022
Location: United Kingdom
Posts: 184
Rep Power: 4
hunt_mat is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
ALE is a possible way. But you could use also the way presented by Peskin in his immersed boundary method.
Which book/paper is this in?
hunt_mat is offline   Reply With Quote

Old   November 22, 2023, 10:49
Default
  #10
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,896
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by hunt_mat View Post
Which book/paper is this in?
Charles Peskin is a famous scientist, you can find his home page and all his papers.
However, have a look to my references here

https://www.researchgate.net/publica...ructured_grids
FMDenaro is offline   Reply With Quote

Reply

Tags
finite volume method, grid analysis


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
[snappyHexMesh] snappyHexMesh stuck when snap is turned on yukuns OpenFOAM Meshing & Mesh Conversion 3 February 2, 2021 14:05
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 06:38
Mesh motion with Translation & Rotation Doginal CFX 2 January 12, 2014 07:21
[Gmsh] 2D Mesh Generation Tutorial for GMSH aeroslacker OpenFOAM Meshing & Mesh Conversion 12 January 19, 2012 04:52
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 19:10


All times are GMT -4. The time now is 12:20.