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

1D Riemann Problem (Shock Tube) with JWL EOS

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 28, 2019, 10:42
Default 1D Riemann Problem (Shock Tube) with JWL EOS
  #1
nsi
New Member
 
Join Date: Oct 2017
Posts: 7
Rep Power: 8
nsi is on a distinguished road
Dear CFD fellows!

I am wondering if anybody knows the analytic solution to 1D Riemann Problem (Shock Tube) if the Equation of State is Johns-Wilkins-Lee (JWL)? What I know is that an isentropic condition is imposed across a rarefaction while the Rankine-Hugoniot jump condition is imposed across a shock. But I have no idea what the equations are, neither how they are derived. Any type of help is much appreciated!

WB
Attached Images
File Type: png sodschematic.png (113.3 KB, 38 views)
File Type: png Screen Shot 2019-01-28 at 9.38.05 AM.png (77.3 KB, 42 views)
nsi is offline   Reply With Quote

Old   January 28, 2019, 18:32
Default
  #2
Senior Member
 
Join Date: Oct 2011
Posts: 239
Rep Power: 16
naffrancois is on a distinguished road
Hello,

Did you have a look to the book of Toro? He derives the exact solution for the Riemann problem associated to the Euler equations and closed by ideal gas eos. Once you understand how it is developed for such a simple eos you should see the additional difficulties with more sophisticated eos. Keep in mind that even with ideal gas the solution is not analytic, it involves a non linear scalar equation for the star pressure.
naffrancois is offline   Reply With Quote

Old   January 29, 2019, 07:54
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I suggest this report from which you could be able to compute the solution using your equation for the pressure
https://books.google.it/books/about/...AJ&redir_esc=y
FMDenaro is offline   Reply With Quote

Old   January 29, 2019, 13:58
Default
  #4
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13
Eifoehn4 is on a distinguished road
The exact Riemann solution for general EOS is a non-trivial task. If you can exclude non-convex and non-smooth cases you at least need two informations from your thermodynamics.
  1. The Caloric EOS
  2. A Model for speed of sound

Then you have to iterate your solution for
  • Shock via calculation of Rankine Hugeniot condition and for
  • Rarefaction Wave via path integration of Generalized Riemann Invariants

This has to be done with two nested iteration loops.

The non-nonvex and non-smooth case is much more challenging!

Regards
Eifoehn4 is offline   Reply With Quote

Old   January 30, 2019, 12:41
Default
  #5
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13
Eifoehn4 is on a distinguished road
Here is a reference of Kamm:

https://cstools.asme.org/csconnect/F...w=yes&ID=45379

The JWL EOS is mentioned there.
Eifoehn4 is offline   Reply With Quote

Old   January 30, 2019, 20:43
Default Exact soln for Convex Equation of state
  #6
New Member
 
Rakesh Raghavaraju
Join Date: Jul 2017
Location: Canada
Posts: 6
Rep Power: 8
RakeshRR is on a distinguished road
There is a paper by James R.Kamm:
'An exact compressible one dimensional Riemann solver for general, convex equation of state'.

This paper outlines the numerical procedure for solving the exact solution for Sod's problem with JWL-EoS. As it was mentioned in previous reply, analytical solution is not available and can only be solved using numerical techniques such as Secant method or Newton-Rhapson method.

Also refer to the paper by Collela and Glaz,1985 ( http://seesar.lbl.gov/anag/publicati...ella/A_1_7.pdf ) on this topic.
RakeshRR is offline   Reply With Quote

Old   July 13, 2019, 23:49
Default
  #7
Member
 
Join Date: Feb 2019
Posts: 64
Rep Power: 7
cfdnewb123 is on a distinguished road
Quote:
Originally Posted by Eifoehn4 View Post
The exact Riemann solution for general EOS is a non-trivial task. If you can exclude non-convex and non-smooth cases you at least need two informations from your thermodynamics.
  1. The Caloric EOS
  2. A Model for speed of sound

Then you have to iterate your solution for
  • Shock via calculation of Rankine Hugeniot condition and for
  • Rarefaction Wave via path integration of Generalized Riemann Invariants

This has to be done with two nested iteration loops.

The non-nonvex and non-smooth case is much more challenging!

Regards
May I know why is the Generalised Riemann Invariant used for deriving rarefaction waves property? I understand why Rankine-Hugoniot condition is applied for shock but not expansion fan as I don't understand what Generalised Riemann Invariant means physically? Is Generalised Riemann Invariant the same as Riemann invariant? As far as I know, Riemann Invariant is the property transported along the characteristics.

I am reading the book from E.F. Toro but he stop elaborating on Generalised Riemann Invariant and mentioned: "For a detailed discussion see the book by Jeffrey [259]". Can someone explain what the Generalised Riemann Invariant means?
cfdnewb123 is offline   Reply With Quote

Old   July 14, 2019, 11:48
Default
  #8
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13
Eifoehn4 is on a distinguished road
You may also derive all thermodynamic quantities in the rarefaction fan with the isentropic condition, which is equivalent to some particular Riemann invariants. Generalized means that the Riemann invariants are written in abstract mathematical notation, e.g. ratio of dw to the respective right eigenvectors. Note that the velocity can not be calculated with the isentropic relation, because it is no thermodynamic quantity.

Another reference for general EOS with NIST and multi-components is given here

https://link.springer.com/article/10...93-019-00896-1

Further informations for non-classical Riemann theory e.g. split and composite waves, can be found in papers of Menikoff and Plohr (1989) or Smith (1979).
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.

Last edited by Eifoehn4; July 14, 2019 at 14:35.
Eifoehn4 is offline   Reply With Quote

Old   July 14, 2019, 16:55
Default
  #9
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by cfdnewb123 View Post
May I know why is the Generalised Riemann Invariant used for deriving rarefaction waves property? I understand why Rankine-Hugoniot condition is applied for shock but not expansion fan as I don't understand what Generalised Riemann Invariant means physically? Is Generalised Riemann Invariant the same as Riemann invariant? As far as I know, Riemann Invariant is the property transported along the characteristics.

I am reading the book from E.F. Toro but he stop elaborating on Generalised Riemann Invariant and mentioned: "For a detailed discussion see the book by Jeffrey [259]". Can someone explain what the Generalised Riemann Invariant means?



The GRI can be deduced by introducing the integral curves in the eigenvectors space as described in the Leveque's textbook.


As far as the Riemann invariant in the physical space are concerned, they exist only in the case of omo-entropic flows, different is the case of iso-entropic flows.
FMDenaro is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Solution of SOD Shock Tube Problem by using HARTEN Modified Flux TVD Scheme xue sheng Main CFD Forum 5 August 27, 2022 08:06
Exact solution for Water-Air Shock Tube Problem agent1039 Main CFD Forum 0 January 25, 2017 06:38
Two solutions of Sod's shock tube problem orxan.shibli Main CFD Forum 14 May 20, 2016 22:41
the shock tube problem and its artifacts with sonicFoam emjay OpenFOAM Running, Solving & CFD 0 December 16, 2011 08:58
exact Riemann solver for JWL EOS abouziar Main CFD Forum 0 January 14, 2009 09:47


All times are GMT -4. The time now is 21:48.