gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity
2,174 Commits 5 Branches 0 Tags
Commit Graph

48 Commits

Author SHA1 Message Date
Grzegorz Kowal
cac043487c SCHEMES: Slightly rewrite HLLD solver for adiabatic MHD. 
This fixes some weird behavior when Bx is small too.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-08-20 17:39:03 -03:00
Grzegorz Kowal
40619aa119 SCHEMES: Slightly rewrite HLLD solver for isothermal MHD. 
This fixes some weird behavior when Bx is small.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-08-20 14:52:43 -03:00
Grzegorz Kowal
8634f89ce0 SCHEMES: Properly average total pressure in MHD HLLC solver. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-08-20 11:27:00 -03:00
Grzegorz Kowal
fda340b21b EQUATIONS: Rename gamma to adiabatic_index. 
Keyword 'gamma' is reserved in Fortran for calculating Gamma function.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-08-16 16:41:28 -03:00
Grzegorz Kowal
8b905114b6 SCHEMES: Remove unused variables. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-08-06 18:56:57 -03:00
Grzegorz Kowal
88138269b7 SCHEMES: Remove 'hlld-m' from the available Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-07-31 21:12:58 -03:00
Grzegorz Kowal
0e7f1a3f15 SCHEMES: Use Mignone's isothermal Riemann solver as HLLD. 
Remove the modified version which allows for the density jump in the
intermediate state.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-07-31 21:10:20 -03:00
Grzegorz Kowal
e6e73dd4f1 SCHEMES: Unify numerical_flux() subroutines. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:46:02 -03:00
Grzegorz Kowal
1ad1ca8890 SCHEMES: Introduce pointers to Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:58 -03:00
Grzegorz Kowal
733dca391d SCHEMES: Separate ROE Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:52 -03:00
Grzegorz Kowal
45222ea63b SCHEMES: Separate ROE Riemann solver for isothermal MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:47 -03:00
Grzegorz Kowal
20a48c2543 SCHEMES: Separate ROE Riemann solver for adiabatic HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:42 -03:00
Grzegorz Kowal
ea4b69050b SCHEMES: Separate ROE Riemann solver for isothermal HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:36 -03:00
Grzegorz Kowal
a9483f320a SCHEMES: Separate HLLC Riemann solver for adiabatic SRMHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:31 -03:00
Grzegorz Kowal
8f27c838df SCHEMES: Separate HLLC Riemann solver for adiabatic SRHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:25 -03:00
Grzegorz Kowal
bae69d1a95 SCHEMES: Separate HLLD Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:20 -03:00
Grzegorz Kowal
366ad64f43 SCHEMES: Separate Mignone HLLD Riemann solver for isothermal MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:15 -03:00
Grzegorz Kowal
1bddacc2d8 SCHEMES: Separate HLLD Riemann solver for isothermal MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:10 -03:00
Grzegorz Kowal
5f366ab2eb EQUATIONS: Add indicator of positive variables. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:06 -03:00
Grzegorz Kowal
c2c2a65dbf SCHEMES: Unify subroutine update_flux(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:45:01 -03:00
Grzegorz Kowal
60539e0558 SCHEMES: Move subroutines update_flux() up. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:52 -03:00
Grzegorz Kowal
be548fde57 SCHEMES: Separate HLLC Riemann solver for adiabatic MHD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:45 -03:00
Grzegorz Kowal
dc44a3ba41 SCHEMES: Separate HLLC Riemann solver for adiabatic HD. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:40 -03:00
Grzegorz Kowal
f3aa475f63 SCHEMES: Unify HLL Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:33 -03:00
Grzegorz Kowal
bf9ef52fb9 SCHEMES: Rename riemann pointer to numerical_flux. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:28 -03:00
Grzegorz Kowal
c54505917a SCHEMES: Use rather subroutines instead of pointers from equations. 
Use compatible (for the same set of equations and equation of state) in
Riemann solvers, not the procedure pointers.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:23 -03:00
Grzegorz Kowal
9801ac6d8b SCHEMES: Fix description of riemann_hll(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:19 -03:00
Grzegorz Kowal
fed3900810 SCHEMES: Separate universal HLL solver. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:14 -03:00
Grzegorz Kowal
1154cb5ec0 SCHEMES: Use higher_order_flux_correction(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:10 -03:00
Grzegorz Kowal
1acb7cca04 SCHEMES: Add a separate subroutine for higher order flux corrections. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:44:04 -03:00
Grzegorz Kowal
ad37de3a99 EQUATIONS, SCHEMES: Merge speeds cm and cp into one variable. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-03 10:43:55 -03:00
Grzegorz Kowal
ad0303d3ba Update copyright. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-03-02 13:18:53 -03:00
Grzegorz Kowal
f36979d74a SCHEMES: Let user control the high order flux correction - corrected. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-04 11:12:05 -03:00
Grzegorz Kowal
7578f1e9f2 Revert "SCHEMES: Let user control the high order flux correction." 
This reverts commit 0dcc3cc95c384c8ed63d6b62aa1cdbcaa58e6672.
 2019-10-04 11:11:05 -03:00
Grzegorz Kowal
0dcc3cc95c SCHEMES: Let user control the high order flux correction. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-04 11:02:05 -03:00
Grzegorz Kowal
92b905cecd SCHEMES: Pass only flux variables to Riemann solvers. 
There is no need to interpolate or pass all variables to the Riemann
solvers. Just pass those related to flux calculation.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-10-03 21:48:04 -03:00
Grzegorz Kowal
5847a4f086 SCHEMES: Correct numerical fluxes for higher-orders. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-03-30 22:50:08 -03:00
Grzegorz Kowal
db1d3c6d4f SCHEMES: Make initialize_schemes() resistant to wrong parameters. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-11 09:51:02 -02:00
Grzegorz Kowal
5193ff3250 SCHEMES: Get rid of argument n in Riemann solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-05 22:08:30 -02:00
Grzegorz Kowal
e9d8bcc6cc EQUATIONS: Remove the vector length argument from prim2cons(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-05 18:07:59 -02:00
Grzegorz Kowal
218ca06190 EQUATIONS: Remove the vector length argument from fluxspeed(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-05 17:26:29 -02:00
Grzegorz Kowal
2c66f83d80 SCHEMES: Use qi and fi for directional states and fluxes. 
Since block dimensions are the same in all directions, only one state
vector and flux vector is required in update_flux_*() subroutines.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-05 15:34:35 -02:00
Grzegorz Kowal
e65d9137d6 SCHEMES: Use nb, ne instead of ib, ie, jb, je, kb, ke. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-05 14:19:49 -02:00
Grzegorz Kowal
2f41bdabf2 SCHEMES: Use bcells instead of im, jm, km. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-02-05 10:55:04 -02:00
Grzegorz Kowal
8bfe3adb7c SCHEMES: Use helper subroutines in print_schemes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-01-30 22:32:53 -02:00
Grzegorz Kowal
122d228248 SCHEMES, DRIVER: Add print_schemes(). 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-01-30 14:58:40 -02:00
Grzegorz Kowal
54f5a6084f SCHEMES: Replace explicit calls to get_parameter_*() with unified one. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-01-28 21:31:30 -02:00
Grzegorz Kowal
f8b65842af Move source files from src/ to sources/. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2019-01-28 09:14:20 -02:00