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

Derivation of source-term of axisymmetric form of Navier-Stokes?

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 1 Post By FMDenaro
  • 1 Post By FMDenaro
  • 1 Post By LuckyTran

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 20, 2024, 16:31
Default Derivation of source-term of axisymmetric form of Navier-Stokes?
  #1
Member
 
Michael Jensen
Join Date: May 2022
Posts: 36
Rep Power: 4
mikethe1wheelnut is on a distinguished road
Greetings all :-)


First, on this subject: Guide: How to ask a question on the forums, if I followed those guidelines to the letter, I would likely never post anything at all because I would never be desperate enough for the answer. I don't -need- the answer, it would just be nice to have, and maybe there is a reference somewhere I've missed.


My question is closely related to this other one, in this same part of the forum, concerns exactly the same pages of the same book: 2D axisymmetric Navier-stokes equation


My focus is instead the source-term. I have not succeeded, yet, in comprehending how switching to cylindrical coordinates introduces an extra component to the source term.



I have searched for answers to this online, but while I have found several derivations of navier stokes in cylindrical, none of them present their answers in the same format as this text.


The next thing I would try, if I really needed to know, would be to attempt to re-derive the equation, with the new symmetry, following material covered earlier in the book. I could also consult derivations of the equations presented as differential equations in a few places, and attempt to convert those to integral equations. From that, it might become obvious where the source term has come from..



Enjoy your coding!
mikethe1wheelnut is offline   Reply With Quote

Old   August 20, 2024, 22:06
Default
  #2
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
There is no source term in the Navier-Stokes equation so there is nothing to derive.
LuckyTran is offline   Reply With Quote

Old   August 21, 2024, 04:41
Default
  #3
Senior Member
 
Hamid Zoka
Join Date: Nov 2009
Posts: 293
Rep Power: 19
Hamidzoka is on a distinguished road
Quote:
Originally Posted by mikethe1wheelnut View Post
Greetings all :-)


First, on this subject: Guide: How to ask a question on the forums, if I followed those guidelines to the letter, I would likely never post anything at all because I would never be desperate enough for the answer. I don't -need- the answer, it would just be nice to have, and maybe there is a reference somewhere I've missed.


My question is closely related to this other one, in this same part of the forum, concerns exactly the same pages of the same book: 2D axisymmetric Navier-stokes equation


My focus is instead the source-term. I have not succeeded, yet, in comprehending how switching to cylindrical coordinates introduces an extra component to the source term.



I have searched for answers to this online, but while I have found several derivations of navier stokes in cylindrical, none of them present their answers in the same format as this text.


The next thing I would try, if I really needed to know, would be to attempt to re-derive the equation, with the new symmetry, following material covered earlier in the book. I could also consult derivations of the equations presented as differential equations in a few places, and attempt to convert those to integral equations. From that, it might become obvious where the source term has come from..



Enjoy your coding!

Hi,
If I am right, you are talking about terms like "viscosity*Ur*Ur/r**2". These are not source terms in its classical defintion. If you see the Laplacian and divergence definitions in cylindrical reference frame, you will find some coefficiencts like r and 1/r which are not the case in cartesian ones. This terms, when derivated, generate some terms that looks like source terms.

Best
Hamidzoka is offline   Reply With Quote

Old   August 21, 2024, 10:37
Default
  #4
Member
 
Michael Jensen
Join Date: May 2022
Posts: 36
Rep Power: 4
mikethe1wheelnut is on a distinguished road
Quote:
Originally Posted by Hamidzoka View Post
Hi,
If I am right, you are talking about terms like "viscosity*Ur*Ur/r**2". These are not source terms in its classical defintion. If you see the Laplacian and divergence definitions in cylindrical reference frame, you will find some coefficiencts like r and 1/r which are not the case in cartesian ones. This terms, when derivated, generate some terms that looks like source terms.

Best

I am referring to the second attached thumb-nail image here: 2D axisymmetric Navier-stokes equation, which, when I click on it expands into a view of the text-book in question. There, you will find equation A.61, which is referred to again at the end of page 449, claiming that the conversion to cylindrical effectively introduces an extra source term, even if there wasn't one there originally. this source-term could be written as [0,0,(p-tau_theta_theta)/r,0], so only applicable to the radial component..
mikethe1wheelnut is offline   Reply With Quote

Old   August 21, 2024, 11:38
Default
  #5
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 mikethe1wheelnut View Post
I am referring to the second attached thumb-nail image here: 2D axisymmetric Navier-stokes equation, which, when I click on it expands into a view of the text-book in question. There, you will find equation A.61, which is referred to again at the end of page 449, claiming that the conversion to cylindrical effectively introduces an extra source term, even if there wasn't one there originally. this source-term could be written as [0,0,(p-tau_theta_theta)/r,0], so only applicable to the radial component..



Unfortunately, the terminology used by Blazek is misleading.

First, the general integral form of the momentum equation has only the surface integral of the normal component of the fluxes. There is no volume integral of a source Q (even in case of gravity we could introduce a potential and then convert in the surface integral).
The textbook denoted as "source" a term that appears only in the specific reference system, it has not physical meaning of source.
I supposed Blazek just organized that for the specific case of cylindrical system and the word "source" has not any physical meaning but only collect that terms outside the surface integral in the initial equation.
mikethe1wheelnut likes this.
FMDenaro is offline   Reply With Quote

Old   August 21, 2024, 13:41
Default
  #6
Member
 
Michael Jensen
Join Date: May 2022
Posts: 36
Rep Power: 4
mikethe1wheelnut is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Unfortunately, the terminology used by Blazek is misleading.

First, the general integral form of the momentum equation has only the surface integral of the normal component of the fluxes. There is no volume integral of a source Q (even in case of gravity we could introduce a potential and then convert in the surface integral).
The textbook denoted as "source" a term that appears only in the specific reference system, it has not physical meaning of source.
I supposed Blazek just organized that for the specific case of cylindrical system and the word "source" has not any physical meaning but only collect that terms outside the surface integral in the initial equation.

..this basically matches my understanding. Blazek has re-organized the terms to try to clarify what has actually changed, leading to -other- confusions. So I would have to, or I will have to, go through the derivation myself, possibly following online sources I've collected, to see for myself where this term actually comes from. Thanks!
mikethe1wheelnut is offline   Reply With Quote

Old   August 21, 2024, 14:07
Default
  #7
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 mikethe1wheelnut View Post
..this basically matches my understanding. Blazek has re-organized the terms to try to clarify what has actually changed, leading to -other- confusions. So I would have to, or I will have to, go through the derivation myself, possibly following online sources I've collected, to see for myself where this term actually comes from. Thanks!

My suggestion is to start from the differential divergence form of the momentum equation, then simply express the nabla operator and vectors in the specific reference system and apply the rule of inner product and derivation of unit vector.

Then you can analyse the integral form.
mikethe1wheelnut likes this.
FMDenaro is offline   Reply With Quote

Old   August 21, 2024, 22:49
Default
  #8
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,761
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Just to further clarify, the "source" term that you are specifically asking about is just the centripetal force. It already has a very clear and intuitive meaning to it.
mikethe1wheelnut likes this.
LuckyTran is offline   Reply With Quote

Old   August 24, 2024, 12:42
Default
  #9
Member
 
Michael Jensen
Join Date: May 2022
Posts: 36
Rep Power: 4
mikethe1wheelnut is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
Just to further clarify, the "source" term that you are specifically asking about is just the centripetal force. It already has a very clear and intuitive meaning to it.

That -does- make a lot of sense. Particularly since I've gone back and looked at the ordinary source term again, and wrapped my head, more-or-less, around gravity itself being a source term. Chemical reactions generating heat.. those are obvious.. :-)
mikethe1wheelnut is offline   Reply With Quote

Old   August 24, 2024, 13:07
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 mikethe1wheelnut View Post
That -does- make a lot of sense. Particularly since I've gone back and looked at the ordinary source term again, and wrapped my head, more-or-less, around gravity itself being a source term. Chemical reactions generating heat.. those are obvious.. :-)



Even gravity, under incompressible flow assumption, is not a source term since can be written as surface integral exactly as the pressure.
FMDenaro is offline   Reply With Quote

Reply

Tags
cylindrical coordinates, navier stokes, source terms


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
Using PengRobinsonGas EoS with sprayFoam Jabo OpenFOAM Running, Solving & CFD 36 July 16, 2024 04:52
[swak4Foam] swak4foam building problem GGerber OpenFOAM Community Contributions 54 April 24, 2015 17:02
Trouble compiling utilities using source-built OpenFOAM Artur OpenFOAM Programming & Development 14 October 29, 2013 11:59
centOS 5.6 : paraFoam not working yossi OpenFOAM Installation 2 October 9, 2013 02:41
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 18:51


All times are GMT -4. The time now is 15:18.