|
[Sponsors] |
January 28, 2019, 10:42 |
1D Riemann Problem (Shock Tube) with JWL EOS
|
#1 |
New Member
Join Date: Oct 2017
Posts: 7
Rep Power: 8 |
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 |
|
January 28, 2019, 18:32 |
|
#2 |
Senior Member
Join Date: Oct 2011
Posts: 239
Rep Power: 16 |
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. |
|
January 29, 2019, 07:54 |
|
#3 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71 |
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 |
|
January 29, 2019, 13:58 |
|
#4 |
Senior Member
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13 |
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.
Then you have to iterate your solution for
This has to be done with two nested iteration loops. The non-nonvex and non-smooth case is much more challenging! Regards |
|
January 30, 2019, 12:41 |
|
#5 |
Senior Member
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13 |
Here is a reference of Kamm:
https://cstools.asme.org/csconnect/F...w=yes&ID=45379 The JWL EOS is mentioned there. |
|
January 30, 2019, 20:43 |
Exact soln for Convex Equation of state
|
#6 |
New Member
Rakesh Raghavaraju
Join Date: Jul 2017
Location: Canada
Posts: 6
Rep Power: 8 |
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. |
|
July 13, 2019, 23:49 |
|
#7 | |
Member
Join Date: Feb 2019
Posts: 64
Rep Power: 7 |
Quote:
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? |
||
July 14, 2019, 11:48 |
|
#8 |
Senior Member
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 13 |
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. |
|
July 14, 2019, 16:55 |
|
#9 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,760
Rep Power: 71 |
Quote:
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. |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |