|
[Sponsors] |
March 17, 2023, 08:47 |
Van Leer flux for 3D
|
#1 |
Senior Member
Sayan Bhattacharjee
Join Date: Mar 2020
Posts: 495
Rep Power: 8 |
Hello everyone,
I'm trying to implement Van Leer flux for 3D, and Toro's book on Riemann solvers only mention the x-split equations for the flux. I think it's fairly understandable that they expected the reader to derive the y and z split equations.I think it's straight forward, but wanted to verify if I was correct or wrong before committing to building the solver. I tried to form the positive flux equation for the y and z direction fluxes. Are the equations correct? |
|
March 17, 2023, 14:11 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,776
Rep Power: 71 |
Quote:
Well, if you strictly rewrite the 1D counterparts in the generic (i,j,k) directions, you apply an index permutation. Formally, that should imply also a directional Mach number Mi=ui/a. Differently you have to consider the velocity magnitude. Anyway, a real multi-dimensional implementation is somehow more complicate. |
||
March 17, 2023, 15:03 |
|
#3 | |
Senior Member
Sayan Bhattacharjee
Join Date: Mar 2020
Posts: 495
Rep Power: 8 |
Quote:
Yes you're correct. The Mach number will be different for each direction and of form Mi = ui/a. So, the formulas I derived previously, were wrong, and possibly unnecessary. Because, in every reference paper/book, they're not giving the full flux equations for the x,y,z directions in 3D. Even in this paper, only the Mach number difference is shown: https://www.abcm.org.br/anais/encit/...gos/1_4669.pdf At section 4.2 it says similar expressions are obtained by using the Mach numbers My = v/a and Mz = w/a. I'm starting to believe that we only need to calculate the Mach number for each direction (Mx, My, Mz), and substitute those in the equation, and that'll actually be the flux contribution in the x,y,z directions. That is, Fx = Fx(rho,a,Mx) is given to us, so Fy = Fy(rho,a,My), and Fz = Fz(rho,a,Mz). (This kind of makes sense to me, as neither rho, or a, changes as in the different directions. They remain constant in the control volume, but the Mach number can change in different directions, as we're splinting the velocity it into parts, and using that to calculate the Mach number) That is, we use the same equation, and we just substitute different values of Mach number to get the fluxes in x,y,z direction. Am I missing anything else? |
||
March 17, 2023, 15:35 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,776
Rep Power: 71 |
||
March 18, 2023, 04:19 |
|
#6 |
Senior Member
Sayan Bhattacharjee
Join Date: Mar 2020
Posts: 495
Rep Power: 8 |
I found a paper which does multidimensional van leer flux splitting for cartesian and curvilinear coordinates: https://arc.aiaa.org/doi/10.2514/3.9858
The paper, and even Toro's book, once mentioned that we can use a rotation matrix that will allow us to split the three dimensional flux vector, by treating it similar to a one dimensional case.. Seems possible, but my trust is low. I might just use Roe flux for tetrahedral elements in the end, if this doesn't work. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Issue with my Van Leer flux vector splitting. | CFD_is_Hard | Main CFD Forum | 1 | December 14, 2022 10:27 |
Radiation in semi-transparent media with surface-to-surface model? | mpeppels | CFX | 11 | August 22, 2019 07:30 |
Udf for moving heat flux in 2D cylindrical geometry | devia21 | Fluent UDF and Scheme Programming | 0 | April 20, 2015 00:27 |
van leer second order method | k_tafazoli | Main CFD Forum | 0 | March 12, 2013 11:20 |
Jacobian Matrix of Van Leer Flux Splitting Vector | Jackson | Main CFD Forum | 2 | July 29, 2010 10:57 |