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

59 Commits

Author SHA1 Message Date
Grzegorz Kowal
8c46b83626 EVOLUTION: Call update_flux() in update_increment(). 
Remove update_fluxes().

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2021-07-10 18:26:52 -03:00
Grzegorz Kowal
3d16d5c13d Update copyrights. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2021-02-04 17:35:04 -03:00
Grzegorz Kowal
4833067bea SCHEMES: Improve detection of numerical degeneracies in HLLD solvers. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-12-27 11:37:07 -03:00
Grzegorz Kowal
7fdac2f584 SCHEMES: Make HD ROE solvers numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-17 17:52:22 -03:00
Grzegorz Kowal
56fb625e59 SCHEMES: Make adiabatic MHD ROE solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-17 17:49:05 -03:00
Grzegorz Kowal
52747633f8 SCHEMES: Make isothermal MHD ROE solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-17 17:27:37 -03:00
Grzegorz Kowal
16a7842fc0 SCHEMES: Make adiabatic MHD HLLD solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-10 12:04:03 -03:00
Grzegorz Kowal
ff7ba2e015 SCHEMES: Make isothermal MHD HLLD solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-10 12:03:03 -03:00
Grzegorz Kowal
44b990644d SCHEMES: Make MHD HLLC solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-10 12:02:16 -03:00
Grzegorz Kowal
87e748feb7 SCHEMES: Make HD HLLC solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-10 12:01:28 -03:00
Grzegorz Kowal
4a4c363660 SCHEMES: Make HLL solver numerically symmetric. 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
 2020-11-10 11:58:17 -03:00
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