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

Asking help on derivation of ALE differential form of momentum conservation equation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 26, 2024, 05:02
Thumbs up Derivation of ALE differential form of momentum conservation equation
  #1
New Member
 
pengfeiguo
Join Date: Nov 2017
Posts: 10
Rep Power: 8
pfguo is on a distinguished road
I want to find the differential equation form of momentum conservation equation that holds at any point
within an arbitrary moving and deforming grid, i.e., the momentum differential equation in ALE description.

After reading the wonderful book "fluid mechanics" sixth edition by PIJUSH K. KUNDU et. al, i.e., equation 4.17,
if I understan right, the ALE form of Reynolds transport is

\frac{D}{D t} \int_{V^*(t)} \rho \mathbf{u} d V^*=\frac{\partial}{\partial t} \int_{V^*(t)} \rho \mathbf{u} d V^*+\underbrace{\int_{A^*(t)} \rho \mathbf{u}(\mathbf{u}-\mathbf{b}) \cdot d \mathbf{A}^*}_{\text {Term } *}

Let handle the volume intergral in the above equation as

\frac{\partial}{\partial t} \int_{V^*(t)} \rho \mathbf{u} d V^*=\int_{V^*(t)} \frac{\partial(\rho \mathbf{u})}{\partial t} d V^*+\underbrace{\int_{A^*(t)} \rho \mathbf{u} \mathbf{b} \cdot d \mathbf{A}^*}_{\text {Term }{ }^{**}}

You see that when I moved the differentiation into the integral,
a surface integral term (Term**) containing the boundary velocity appeared.

This term would cancel out the term containing the surface velocity (\textbf{b}) in the Reynolds transport equation,
ultimately leading to the absence of surface velocity in the transport equation.

I believe this is an incorrect result.

Can you give me some suggestions on where I am doing wrong?

Thanks!

Sincerely,

Guo Pengfei

Last edited by pfguo; April 30, 2024 at 05:51. Reason: edit for more clear question
pfguo is offline   Reply With Quote

Old   April 29, 2024, 15:50
Default
  #2
agd
Senior Member
 
Join Date: Jul 2009
Posts: 357
Rep Power: 19
agd is on a distinguished road
Try using Gauss' theorem to convert the area integrals to volume integrals. Then you have to make the argument that if the integral equals to zero then under sufficient conditions for continuity of the integrand that the integrand must equal zero over the domain. Since these conditions fail in the presence of discontinuities most people use the integral form, or figure out the singular limit conditions to apply to the jumps where the differential form is invalid.
agd is offline   Reply With Quote

Old   April 29, 2024, 15:53
Default
  #3
agd
Senior Member
 
Join Date: Jul 2009
Posts: 357
Rep Power: 19
agd is on a distinguished road
You will also need to move the temporal derivative inside the volume integral, even in the case of deforming boundaries. There is a mathematical theorem to handle that also - I'll leave that as an exercise for the reader.
agd is offline   Reply With Quote

Old   April 30, 2024, 05:56
Default
  #4
New Member
 
pengfeiguo
Join Date: Nov 2017
Posts: 10
Rep Power: 8
pfguo is on a distinguished road
Quote:
Originally Posted by agd View Post
You will also need to move the temporal derivative inside the volume integral, even in the case of deforming boundaries. There is a mathematical theorem to handle that also - I'll leave that as an exercise for the reader.
I really appreciate for your attention.

Actually, I got stuck even before applying the divergence theorem to the surface integral term.

I got stuck on the term involving the derivative of the volume integral with respect to time.

When I moved the differentiation inside the integral, a surface integral term (Term**) containing the boundary velocity appeared.

This term would cancel out the term containing the surface velocity (\textbf{b}) in the Reynolds transport equation,
ultimately leading to the absence of surface velocity in the transport equation.

I believe this is an incorrect result.

By the way, for better readability, I spent some time replacing the original image with LaTeX formulas
pfguo is offline   Reply With Quote

Old   April 30, 2024, 10:48
Default
  #5
agd
Senior Member
 
Join Date: Jul 2009
Posts: 357
Rep Power: 19
agd is on a distinguished road
Are you sure you aren't applying Liebniz's theorem twice? Consider mass conservation for a fluid element - mass conservation in a Lagrangian sense is simply the time derivative of the integral of density over the volume is zero. Then you can apply Liebniz's theorem to get the equation where the temporal derivative is inside the integral and the relative velocity appears in the area integral. Then apply divergence to that. It seems like you are double-counting the surface velocity.
agd is offline   Reply With Quote

Old   April 30, 2024, 10:55
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Just follow the idea of the derivation of the Reynolds transport theorem.
You can read a rigorous description in the textbook of Chorin & Marsden.
FMDenaro is offline   Reply With Quote

Old   May 1, 2024, 00:01
Default
  #7
New Member
 
pengfeiguo
Join Date: Nov 2017
Posts: 10
Rep Power: 8
pfguo is on a distinguished road
Quote:
Originally Posted by agd View Post
Are you sure you aren't applying Liebniz's theorem twice? Consider mass conservation for a fluid element - mass conservation in a Lagrangian sense is simply the time derivative of the integral of density over the volume is zero. Then you can apply Liebniz's theorem to get the equation where the temporal derivative is inside the integral and the relative velocity appears in the area integral. Then apply divergence to that. It seems like you are double-counting the surface velocity.
Currently, the part where I'm stuck doesn't involve Gauss's theorem, and I don't think my problem lies in the use of Leibniz's theorem either. I strongly feel that I may have misunderstood the quantity whether in Euler's description or or ALE's description and their derivatives with respect to space or time.
pfguo is offline   Reply With Quote

Old   May 1, 2024, 00:03
Default
  #8
New Member
 
pengfeiguo
Join Date: Nov 2017
Posts: 10
Rep Power: 8
pfguo is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Just follow the idea of the derivation of the Reynolds transport theorem.
You can read a rigorous description in the textbook of Chorin & Marsden.
Thank you for recommending the book. I've been reading it for a while, but I'm still unclear about where my problem lies. Nevertheless, I sincerely appreciate your attention.
pfguo is offline   Reply With Quote

Old   May 1, 2024, 12:34
Default
  #9
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by pfguo View Post
Thank you for recommending the book. I've been reading it for a while, but I'm still unclear about where my problem lies. Nevertheless, I sincerely appreciate your attention.



Follow carefully the derivation of Chorin
FMDenaro is offline   Reply With Quote

Old   May 4, 2024, 04:09
Default How the lapacian eqation for mesh motion is solved in openfoam
  #10
@ss
New Member
 
shristi
Join Date: Jun 2023
Posts: 26
Rep Power: 3
@ss is on a distinguished road
Hello,
I am working on the FSI problem for vortex-induced vibration. The ALE approach has been used to solve coupling motion for fluid and solid. The fluid equation is discretized based on FVM, and the structure equation is the newmark method used, and that is FEM. However, I am not sure about how the mesh motion equation, which is a Laplacian equation, is solved to update mesh motion velocity. and also, don't you know if this method is used for it FVM or FEM? Please help me to clear this up and also suggest some books that can give a good understanding of the ALE approach computationally. I am using OpenFOAM for my computations.


Thank you
@ss is offline   Reply With Quote

Old   May 5, 2024, 21:34
Default
  #11
New Member
 
pengfeiguo
Join Date: Nov 2017
Posts: 10
Rep Power: 8
pfguo is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Follow carefully the derivation of Chorin
Thank you for your suggestion.

I have been reading the content you mentioned in the book you recommended these past few days.

Basically, I can understand the proofs in the book, but it doesn't relate directly to the question I asked.

I kindly request you to read my question again.
pfguo is offline   Reply With Quote

Reply


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
Adding a term in the diagonal of the matrix form of the momentum equation Agavi OpenFOAM Programming & Development 1 August 3, 2020 14:12
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 06:21
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
Question on the discretization of momentum equation in icoFoam MPJ OpenFOAM 3 October 4, 2011 10:44
Constant velocity of the material Sas CFX 15 July 13, 2010 09:56


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