Implement Colella and Woodward's Piecewise Parabolic Method (PPM) for
block prolongation. This method is also used for the prolongation of
boundaries from blocks at lower levels to higher levels.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The -B.div(B) source term in the momentum equation was dropped for
the isothermal MHD case. Apparently, its lack causes some numerical
instabilities related to the accumulation of the divergence of B.
Therefore, take it into account for the isothermal case too.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
- Update riemann_mhd_iso_kepes() subroutine to include missing symmetry
preserving optimizations.
- Ensure consistency with the adiabatic version by adding and utilizing
the ch variable.
- Correct the initialization and usage of transformation matrices (rm and tm).
- Adjust calculations to ensure symmetric properties, particularly in the handling
of variables cs, ca, and cglm.
- Include appropriate sign functions for consistent behavior.
- Make changes to the calculation of intermediate variables to maintain symmetry
and consistency.
- Adjusted flux calculation to correctly apply symmetry preserving
terms.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This change introduces the logical flag 'noflush' to control the flushing
and synchronization of I/O buffers. Subsequently, subroutines
'enable_flush_and_sync' and 'disable_flush_and_sync' were added to enable
and disable the flushing and synchronization, respectively.
The 'flush_and_sync' subroutine was modified to respect the noflush flag.
Additionally, this change updates 'statistics.F90' to read
the 'enable_io_flush' parameter and call the appropriate subroutine
('enable_flush_and_sync' or 'disable_flush_and_sync') based on its value.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
To prevent solution instability caused by sudden increases in the
timestep, such as when the maximum refinement level decreases, a new
parameter 'dt_jump_factor' has been introduced. This parameter controls
the allowable increase in the timestep.
By default, 'dt_jump_factor' is set to 1.05, permitting a maximum
increase of 5% from the previous timestep.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
Restricted initial time steps to values smaller than the CFL time step
to ensure numerical stability. This prevents instability issues that
may arise from excessively large initial steps.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit separates the smoothness indicator and weight calculations
for 5th order WENO methods into different subroutines. This allows
modification of only the weight factor calculations for different WENO methods.
It extends the 5th order methods to include WENO5-JS, WENO5-Z+, WENO5-ZC,
WENO5-ZC+, WENO5-Z+M, and WENO5-AT.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The injected energy and rate were added to the global array within the
loop over the blocks, essentially multiplying the energy by the number
of blocks for each process. This fix moves the injected energy
accounting out of the loop over the blocks.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The kappa parameter is calculated from CFL value. It is a decreasing
value with increasing CFL. However, to maintain the precision of MP
limiter, it should be kept above 1. Therefore, instead of taking the min
values between the user defined kappa and the one calculated from the
CFL value, take its maximum.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The tearing test problem was numerically unstable in the adiabatic case.
Setting bty to 1 in the case of br beeing zero solved the issue. Make it
consistent for the isothermal ROE solver as well. Additionaly, one typo
was fixed in the calculation of the left eigenvectors.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
