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

Near vacuum for Euler equations

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

Like Tree2Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 26, 2023, 05:57
Default
  #21
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,180
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
I'm all but an expert in this field but, as a practictioner in need of a positivity preserving scheme, I gathered the following information:

1) A 1D, first order positively preserving scheme can be made nD and high order (EDIT: I now see this is similar to a work you already cited, so you probably know this as well)

https://deepblue.lib.umich.edu/bitst...pdf?sequence=1

2) The scheme that Einfeldt defines as HLLEM in "On Godunov-type methods near low densities" is proven by himself to be positively preserving and has been later used as a basis for a positivity preserving entropy fix for a Roe scheme, e.g.,

https://perso.ensta-paris.fr/~pelant...nti_oxford.pdf

Altough, I see there might be misalignment in nomenclature across references

3) My personal experience is with a preconditioned Roe scheme using an entropy fix based on the preconditioned HLLE+ scheme

https://arc.aiaa.org/doi/10.2514/1.12176

and it is positivity preserving at 2nd order for a shock tube on an unstructured grid (that is, an actual 3D tube with a finite cross section). Of course, a working limiter is needed here.

May I suggest you to use the very test in the Einfeldt paper to check your scheme?
sbaffini is offline   Reply With Quote

Old   April 26, 2023, 06:35
Default
  #22
New Member
 
Davide
Join Date: Apr 2023
Posts: 10
Rep Power: 3
dappoli is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
Altough, I see there might be misalignment in nomenclature across references
I agree, for instance looking at the paper you linked for the HLLE+, it seems that I have been using the HLLE+ all the time, where \bar{u} = \widehat{u} while the HLLEM scheme does not use the Roe average for the contact discontinuity velocity.

However in the paper I followed originally:

https://www.sciencedirect.com/scienc...092?via%3Dihub (page 3) it uses the Roe averages.



Anyway I will try Einfeldt test to see what happens
dappoli is offline   Reply With Quote

Old   April 26, 2023, 08:37
Default
  #23
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by dappoli View Post
First of all thank you for the reply.

I think we are talking about the same thing, for me a positive preserving scheme is a scheme that generates "physically admissible solutions from physical data" [1], or in other words it keeps positive pressure and density if the initial condition are positive. For the Euler equations the physically admissible state in this case is defined as:

G = \{ \textbf{U} \vert \rho > 0 \quad \text{and} \quad E-(m_x^2 + m_y^2 + m_z^2)/2 \rho >0 \}
being U the conserved variable vector and m the momentum along the 3 directions. For example Roe scheme is not positive preserving and neither is Osher's. Einfeldt [2] demonstrated that the Rusanov and HLLE are positive preserving, but for example HLLEM is not. All this holds true for 1st order schemes, once 2nd order schemes are taken into account, from what I found there is no proof of satisfying the positive preserving property. That is why the Muscl reconstruction near low density fails also for Rusanov. Clearly if the timestep is strongly reduced this is in general postponed, but inevitable. Having said this, I am new to the study and application of Riemann solvers so I might have said something wrong.



Going back to Noh:

You are right, but being focused on other problems, I forgot that I am using a variation of the MINMOD slope limiter described in Kurganov-Tadmor [3] (the scheme I am using which is equal to a MUSCL Rusanov) defined as:

u_x(x,t)\approx \text{minmod}\left(\theta \frac{u(x+a)-u(x)}{a}, \frac{u(x+a)-u(x-b)}{a+b}, \theta \frac{u(x)-u(x-b)}{b}\right)
which allows to reduce the diffusivity of the minmod by tuning the parameter \theta\in [1,2]. Setting the parameter to 1 recovers the original minmod. I had it at 1.4. Attached you can find the solution with the parameter = 1. Now the solution is sharp, although the HLLEM keeps behaving strangely at the centre of the domain (note that the simulation is the same, I just cut the plot). Do you have any idea of why it behaves like that at the origin?



Finally, I coded the "hurricane" test described in [4] where a uniform domain is initialized with a vortex like velocity. This creates a vacuum point in the centre of the domain and if the scheme is not positive preserving, eventually it will fail. Of course 2nd order MUSCL Rusanov, HLLE and HLLEM fail in this, but first order HLLE and Rusanov do not. What I did in the end is to use a posteriori approach:

I advance 1 timestep,

I check for negative pressure, temperature or NANs ,

If I detect errors, I reject the step and re do it with 1st order Rusanov in the troubled cell and its neighbours.



I already tried doing this with little success, but I was not considering the neighbouring cells and just using the Rusanov flux in the troubled one. This eventually does not work. Now I will have to deal with some special cases such as when the troubled cell is on the boundary between processors (I use domain decomposition with MPI).

Attached you can see a video of the hurricane like simulation with the HLLEM a polsteriori corrected scheme. Sorry for the low quality, the solution is smooth in reality but the forum has a very low size limit



Now I will have to try this on my problem, probably with a careful optimization of the minmod limiter and neighbour radius to avoid failure.


[1] https://link.springer.com/chapter/10...-011-5169-6_16
[2] https://www.sciencedirect.com/scienc...21999191902113
[3] https://www.math.umd.edu/~tadmor/pub...or.JCP-00I.pdf
[4] https://www.math.hkust.edu.hk/~makxu...D-inviscid.pdf
Yes I think so. Now it makes sense, the simulation looks similar now. The HLLEM (HLLE+) for the Noh test might be a normal behavior.

I have also testest the Hurrican example, the high speed example, with reconstruction (TVD, MINMOD) + Rusanov and even HLLC and had no problems. Your video seems unsymmetric, the solution should be symmetric.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   April 29, 2023, 13:15
Default
  #24
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
May i ask, how much is lowest density you are running into?
arjun is offline   Reply With Quote

Old   April 30, 2023, 04:11
Default
  #25
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by arjun View Post
May i ask, how much is lowest density you are running into?
Depends on the mesh resolution and the settings. With linear TVD reconstruction (MINMOD), CFL=0.5, explicit Euler O(2) and HLLC I get following values for the hurricane example (high speed):

\rho_{\text{N}=160}=9.4\times 10^{-4}
\rho_{\text{N}=320}=3.7\times 10^{-4}
\rho_{\text{N}=640}=1.6\times 10^{-4}
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   April 30, 2023, 13:42
Default
  #26
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by Eifoehn4 View Post
Depends on the mesh resolution and the settings. With linear TVD reconstruction (MINMOD), CFL=0.5, explicit Euler O(2) and HLLC I get following values for the hurricane example (high speed):

\rho_{\text{N}=160}=9.4\times 10^{-4}
\rho_{\text{N}=320}=3.7\times 10^{-4}
\rho_{\text{N}=640}=1.6\times 10^{-4}



Thank you. I was just curious. I actually have worked on cases with very low density and I decided against Riemann type solvers. I decided to go pressure based solver route. With that in my case i could go as low as 1E-5 without any problems, and i can go till 1E-7 if i push it.
arjun is offline   Reply With Quote

Old   April 30, 2023, 14:58
Default
  #27
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by arjun View Post
Thank you. I was just curious. I actually have worked on cases with very low density and I decided against Riemann type solvers. I decided to go pressure based solver route. With that in my case i could go as low as 1E-5 without any problems, and i can go till 1E-7 if i push it.
Sounds great.

We should keep in mind, that the example is designed in such a way, that it contains a physical singularity in the center (vaccum point). Every consistent, convergent and positivy preserving solver should reach a density of 1 \times 10^{-15} - 1\times 10^{-16} without much effort (under mesh refinement) until round-off errors come into play.

The validity of the Euler Equations under very high Knudsen numbers is a different topic.
arjun likes this.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   April 30, 2023, 15:31
Default
  #28
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by Eifoehn4 View Post
Sounds great.

We should keep in mind, that the example is designed in such a way, that it contains a physical singularity in the center (vaccum point). Every consistent, convergent and positivy preserving solver should reach a density of 1 \times 10^{-15} - 1\times 10^{-16} without much effort (under mesh refinement) until round-off errors come into play.

The validity of the Euler Equations under very high Knudsen numbers is a different topic.



Can you point out which test you are talking about. Sorry a lot going on in this thread so for me hard to keep track.



I can try with pressure based solver and see what happens. I expect it to break down too but things can't improve if they don't break.
arjun is offline   Reply With Quote

Old   April 30, 2023, 15:37
Default
  #29
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by arjun View Post
Can you point out which test you are talking about. Sorry a lot going on in this thread so for me hard to keep track.



I can try with pressure based solver and see what happens. I expect it to break down too but things can't improve if they don't break.
The numbers are given for the hurricane example, reference [4], see the previous answers.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   April 30, 2023, 15:46
Default
  #30
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by dappoli View Post
These are the equations:

\frac{\partial n_\alpha}{\partial t} + \nabla \cdot \left(n_\alpha\textbf{u}_\alpha\right) = S_\alpha
\frac{\partial \left(n_\alpha \textbf{u}_\alpha\right)}{\partial t} + \nabla \cdot \left(\left(n_\alpha \textbf{u}_\alpha\right)\textbf{u}_\alpha\right) = \frac{qn_\alpha}{m_\alpha}\left[\textbf{u}_\alpha\times \textbf{B} + \nabla \phi\right] - \frac{1}{m_\alpha}\nabla P_\alpha
\frac{\partial \left(\mathcal{E}_\alpha\right)}{\partial t} + \nabla \cdot \left(\left(\mathcal{E}_\alpha + p_\alpha\right)\textbf{u}_\alpha \right) = qn_\alpha \textbf{u}_\alpha \cdot \nabla \phi + {Q}_\alpha
which are integrated for electrons and ions where q is the charge (\pm e, the fundamental charge), B the magnetic field, in my case perpendicular to the domain,
\mathcal{E}_\alpha = \frac{3}{2}q n_\alpha T_\alpha + \frac{1}{2}m_\alpha n_\alpha u_\alpha^2 the total energy ( Note that T is expressed in eV here), Q_\alpha = \mathcal{E}_{inj}S_\alpha the energy injection source term and \phi the elctric potential obtained from:



\nabla^2\phi = \frac{e}{\epsilon_0}\left(n_e-n_i\right)


At the moment I am trying also with a Strang splitting, so that I solve for the homogeneous system which is identical to Euler. The eigenvalues are the usual eigenvalues. In the algorithm I advance one dt the single species and then compute the new potential for the next timestep.

If you are not familiar with plasma, just as a reference, the problem is complicated by the broad range of timescales. For instance, for Xenon m_i = 2.18e-25 kg and m_e = 9.1e-31 kg, so that there is a huge difference in the timescales of ions and electrons. Thermal velocities for electrons are in the order of 1e6 m/s and the temperature is 1-2 order of magnitude larger than the one of ions. Moreover, to solve the non-neutrality of the problem, the debye length and plasma frequency must be resolved, limiting the problem to very fine meshes and timesteps in the 0.5-10 ps order.



I know that in the pressure-density-velocity the TVD is not granted, but more in general the positivity of schemes such as the Rusanov is not granted either for multidimensional system with reconstruction > order 1. The negative temperature sometimes occurs in the face reconstruction, but that is easily solvable (or at least more easily with a dedicated limiter), but sometimes it just pops up in the cell centre. After I advance in time the conservative variables, recovering the temperature in the cell centre my give rise to a negative value. Note that in this problem it is not necessary to know the temperature at cell centre, but if one has to compute collision frequencies ecc. it is mandatory.



I know that there is a work from Linde and Roe where multidimensional schemes of order >1 are made positive preserving using a dedicated limiter, which however is somehow not straightforward for the pressure. Indeed they state that the pressure must be limited iteratively to guarantee the positivity. I would try to avoid this given the fact that the timestep is already incredibly small and time can be an issue.



https://link.springer.com/chapter/10...-011-5169-6_16


By the way, does it make sense for you that i can see a discontinuity in the solution (again in the temperature, but I guess only because I am dividing by the discontinuous density) where I use 1st order rusanov instead of 2nd order HLLEM? Should the interface be smooth? Because I might have an error in the code if that is the case



I am not an expert with plasma physics and I never worked on your problem but I would ask you how do you get the Poisson equations for phi from the second equation, what are the BC.s you prescribed and how do you solve.

If I understand correctly, the solved function is a "source" term in the energy (third equation).


What is more, did you evaluate if you have a stiff hyperbolic problem due to the coupling?
FMDenaro is offline   Reply With Quote

Old   April 30, 2023, 17:45
Default
  #31
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I am not an expert with plasma physics and I never worked on your problem but I would ask you how do you get the Poisson equations for phi from the second equation, what are the BC.s you prescribed and how do you solve.

If I understand correctly, the solved function is a "source" term in the energy (third equation).


What is more, did you evaluate if you have a stiff hyperbolic problem due to the coupling?
The model originally stems from the coupled Euler-Maxwell equations, in the limit of electro-statics. Therefore, the potential is the electro-static potenial. The gradients (electric field) exerts forces (Lorentz force) on the electron or ion species, which carry charge, which is some how similar to pressure for the incompressible Navier-Stokes equations. The model is well-known in the semi-conductor industrie.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   May 1, 2023, 03:39
Default
  #32
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Eifoehn4 View Post
The model originally stems from the coupled Euler-Maxwell equations, in the limit of electro-statics. Therefore, the potential is the electro-static potenial. The gradients (electric field) exerts forces (Lorentz force) on the electron or ion species, which carry charge, which is some how similar to pressure for the incompressible Navier-Stokes equations. The model is well-known in the semi-conductor industrie.
Ok, thanks, that confirms the reason of my question.
From application of the divergence operator on the second equation ?
What about the BCs and the discretization of the Poisson problem?
FMDenaro is offline   Reply With Quote

Old   May 1, 2023, 11:53
Default
  #33
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Ok, thanks, that confirms the reason of my question.
From application of the divergence operator on the second equation ?
What about the BCs and the discretization of the Poisson problem?
The potential rests on the assumption that the electric field is irrotational:

\vec{\nabla}\times\vec{E} = 0.

The boundary conditions and discretization may be choosen independantly.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   May 1, 2023, 11:58
Default
  #34
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,839
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Eifoehn4 View Post
The potential rests on the assumption that the electric field is irrotational:

\vec{\nabla}\times\vec{E} = 0.

The boundary conditions and discretization may be choosen independantly.





Could you please better detail? I mean, apply the divergence operator to the second equation and write down the Neumann BCs.
Is that sufficient to ensure the compatibility condition for the existance of a solution?
FMDenaro is offline   Reply With Quote

Old   May 1, 2023, 14:28
Default
  #35
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Could you please better detail? I mean, apply the divergence operator to the second equation and write down the Neumann BCs.
Is that sufficient to ensure the compatibility condition for the existance of a solution?
To provide a little clarity: The model here is an extention of the standard Euler equation with additional external forces (source terms). Here, the additional external force is the electric/(magnetic) field, either supplied by the full Maxwell equations or by the Coloumb law via the electro-static potential.

You may solve the compressible Euler-Maxwell or Euler-Poisson equation in a pressure based fashion rather than a density based fashion. In this case you would end up with an additional Poisson like equation, however for pressure, including additional source terms. Except external forces, there should be no difference to handle BCs in order to fullfill the compatibility condition for pressure.

The former electro-static Poisson equation should of course also satisfy the compatibilty condition in order to be unique.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   May 3, 2023, 02:41
Default
  #36
New Member
 
Davide
Join Date: Apr 2023
Posts: 10
Rep Power: 3
dappoli is on a distinguished road
Hi everyone, sorry for not replying for a while.

Eifoehn4 is right, the system comes from the Euler-Maxwell in the electrostatic limit. In my case boundary conditions for Poisson are simple dirichlet.


Quote:
Originally Posted by arjun View Post
May i ask, how much is lowest density you are running into?

In my specific problem I got number densities as low as 0.7 m^-3 while other regions of the domain are at 10^15, so very low. It turns out that even with a more diffusive scheme sometimes the simulation fails, in particular when I have a sudden diffusion of density that expands in the low density region.

This creates a non physical shock with extreme temperatures (to reproduce the problem try to do a free expansion in quasi-vacuum). In order to avoid this I decided to try another approach and coded a "vacuum tracking" algorithm which follows the interface between absolute vacuum (just 0 ) and solves Riemann problems with vacuum. This works fine especially in the free expansion case since the approximate/exact riemann solution is available (see Munz 1994), however once I have source terms, such as magnetic confinement ecc, it is not anymore free expansion and I think I have to modify the riemann solution to account for this. In particular I guess I have to remove the rarefaction wave since nothing can expand beyond the magnetic line (free expansion always assumes a rarefaction wave at the interface, being the contact discontinuity impossible to obtain due to the pressure gradient). I am not really sure on how to modify the flux function to account for this but I will look into it.



I never thought of using pressure based solver, I honestly know only the one like Piso or Simple, but I always skip them for this kind of problems. Might be due to the fact that I studied on LeVeque and he just goes with Riemann problems. I will look into them as well, thanks
dappoli is offline   Reply With Quote

Old   May 3, 2023, 03:26
Default
  #37
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Please check the Toro's book as that is the only book that kind of talk about Riemann solvers with Vacuum.

I thought about the appraoch you talk about detecting Vacuum and this is why i was reading the Toro's book for it.

The major reason why I gave up on Riemann solver was that i had a calculation with Vacuum pump (2 phase flow) and Riemann solver did not give me good flow profile in rotating geometry while pressure based solver did give me good profile. In this pump low pressures are generated so was exploring the idea of detecting vacuum and then going forward with modifying Riemann solver.

The unreliability of the Riemann solver made me give up on it.



I point out that the solver was validated for standard test problem so it was not inaccurate just that when rotation was involved things did not add up.

I actually wrote the Riemann solver in long double precision because of low densities involved.


Anyway keep us updated and also check the Toro's book. I do work on same problem actually and will be spending next 2 to 3 weeks on it (i resume where i left off). I have electron transport and now i will add ion transport with it.



Quote:
Originally Posted by dappoli View Post
Hi everyone, sorry for not replying for a while.

Eifoehn4 is right, the system comes from the Euler-Maxwell in the electrostatic limit. In my case boundary conditions for Poisson are simple dirichlet.





In my specific problem I got number densities as low as 0.7 m^-3 while other regions of the domain are at 10^15, so very low. It turns out that even with a more diffusive scheme sometimes the simulation fails, in particular when I have a sudden diffusion of density that expands in the low density region.

This creates a non physical shock with extreme temperatures (to reproduce the problem try to do a free expansion in quasi-vacuum). In order to avoid this I decided to try another approach and coded a "vacuum tracking" algorithm which follows the interface between absolute vacuum (just 0 ) and solves Riemann problems with vacuum. This works fine especially in the free expansion case since the approximate/exact riemann solution is available (see Munz 1994), however once I have source terms, such as magnetic confinement ecc, it is not anymore free expansion and I think I have to modify the riemann solution to account for this. In particular I guess I have to remove the rarefaction wave since nothing can expand beyond the magnetic line (free expansion always assumes a rarefaction wave at the interface, being the contact discontinuity impossible to obtain due to the pressure gradient). I am not really sure on how to modify the flux function to account for this but I will look into it.



I never thought of using pressure based solver, I honestly know only the one like Piso or Simple, but I always skip them for this kind of problems. Might be due to the fact that I studied on LeVeque and he just goes with Riemann problems. I will look into them as well, thanks
arjun is offline   Reply With Quote

Old   May 3, 2023, 03:54
Default
  #38
Senior Member
 
Eifoehn4's Avatar
 
-
Join Date: Jul 2012
Location: Germany
Posts: 184
Rep Power: 14
Eifoehn4 is on a distinguished road
Quote:
Originally Posted by arjun View Post
Please check the Toro's book as that is the only book that kind of talk about Riemann solvers with Vacuum.

I thought about the appraoch you talk about detecting Vacuum and this is why i was reading the Toro's book for it.

The major reason why I gave up on Riemann solver was that i had a calculation with Vacuum pump (2 phase flow) and Riemann solver did not give me good flow profile in rotating geometry while pressure based solver did give me good profile. In this pump low pressures are generated so was exploring the idea of detecting vacuum and then going forward with modifying Riemann solver.

The unreliability of the Riemann solver made me give up on it.



I point out that the solver was validated for standard test problem so it was not inaccurate just that when rotation was involved things did not add up.

I actually wrote the Riemann solver in long double precision because of low densities involved.


Anyway keep us updated and also check the Toro's book. I do work on same problem actually and will be spending next 2 to 3 weeks on it (i resume where i left off). I have electron transport and now i will add ion transport with it.
Did you consider All-Mach Riemann solvers? Without any fix, Riemann solvers become quite diffusive in the low Mach regime. The reason for this is that the jump terms (responsible for numerical diffusion) depend on all characteristics, also the highest, e.g. the speed of sound.
__________________
Check out my side project:

A multiphysics discontinuous Galerkin framework: Youtube, Gitlab.
Eifoehn4 is offline   Reply With Quote

Old   May 3, 2023, 11:06
Default
  #39
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,283
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by Eifoehn4 View Post
Did you consider All-Mach Riemann solvers? Without any fix, Riemann solvers become quite diffusive in the low Mach regime. The reason for this is that the jump terms (responsible for numerical diffusion) depend on all characteristics, also the highest, e.g. the speed of sound.

Thank you. I will check it once i have some time to spare. I used hllc in that sim.
arjun is offline   Reply With Quote

Reply

Tags
euler, plasma, vacuum

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
Solving coupled scalar transport equations in openFoam pavaninguva OpenFOAM Running, Solving & CFD 1 December 28, 2019 08:28
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Riemann invariants of adjoint equations of shallow water equations zqb0929 Main CFD Forum 0 March 15, 2012 00:54
CFD governing equations m.gos Main CFD Forum 0 April 30, 2011 14:21
? fluctuating equations for homogenous shear turb. ff_fan Main CFD Forum 1 September 20, 2002 07:39


All times are GMT -4. The time now is 13:46.