Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAchievable Efficiency of Numerical Methods for Simulations of Solar Surface Convection
106
0
0.0
(
0
)
Added by
Hannes Grimm-Strele
Publication date
2014
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We investigate the achievable efficiency of both the time and the space discretisation methods used in Antares for mixed parabolic-hyperbolic problems. We show that the fifth order variant of WENO combined with a second order Runge-Kutta scheme is not only more accurate than standard first and second order schemes, but also more efficient taking the computation time into account. Then, we calculate the error decay rates of WENO with several explicit Runge-Kutta schemes for advective and diffusive problems with smooth and non-smooth initial conditions. With this data, we estimate the computational costs of three-dimensional simulations of stellar surface convection and show that SSP RK(3,2) is the most efficient scheme considered in this comparison.
rate research
Read More
We investigate the applicability of curvilinear grids in the context of astrophysical simulations and WENO schemes. With the non-smooth mapping functions from Calhoun et al. (2008), we can tackle many astrophysical problems which were out of scope with the standard grids in numerical astrophysics. We describe the difficulties occurring when implementing curvilinear coordinates into our WENO code, and how we overcome them. We illustrate the theoretical results with numerical data. The WENO finite difference scheme works only for high Mach number flows and smooth mapping functions whereas the finite volume scheme gives accurate results even for low Mach number flows and on non-smooth grids.
We investigate the computational performance of various numerical methods for the integration of the equations of motion and the variational equations for some typical classical many-body models of condensed matter physics: the Fermi-Pasta-Ulam-Tsingou (FPUT) chain and the one- and two-dimensional disordered, discrete nonlinear Schrodinger equations (DDNLS). In our analysis we consider methods based on Taylor series expansion, Runge-Kutta discretization and symplectic transformations. The latter have the ability to exactly preserve the symplectic structure of Hamiltonian systems, which results in keeping bounded the error of the systems computed total energy. We perform extensive numerical simulations for several initial conditions of the studied models and compare the numerical efficiency of the used integrators by testing their ability to accurately reproduce characteristics of the systems dynamics and quantify their chaoticity through the computation of the maximum Lyapunov exponent. We also report the expressions of the implemented symplectic schemes and provide the explicit forms of the used differential operators. Among the tested numerical schemes the symplectic integrators $ABA864$ and $SRKN^a_{14}$ exhibit the best performance, respectively for moderate and high accuracy levels in the case of the FPUT chain, while for the DDNLS models $s9mathcal{ABC}6$ and $s11mathcal{ABC}6$ (moderate accuracy), along with $s17mathcal{ABC}8$ and $s19mathcal{ABC}8$ (high accuracy) proved to be the most efficient schemes.
WavePacket is an open-source program package for numerical simulations in quantum dynamics. Building on the previous Part I [Comp. Phys. Comm. 213, 223-234 (2017)] and Part II [Comp. Phys. Comm. 228, 229-244 (2018)] which dealt with quantum dynamics of closed and open systems, respectively, the present Part III adds fully classical and mixed quantum-classical propagations to WavePacket. In those simulations classical phase-space densities are sampled by trajectories which follow (diabatic or adiabatic) potential energy surfaces. In the vicinity of (genuine or avoided) intersections of those surfaces trajectories may switch between surfaces. To model these transitions, two classes of stochastic algorithms have been implemented: (1) J. C. Tullys fewest switches surface hopping and (2) Landau-Zener based single switch surface hopping. The latter one offers the advantage of being based on adiabatic energy gaps only, thus not requiring non-adiabatic coupling information any more. The present work describes the MATLAB version of WavePacket 6.0.2 which is essentially an object-oriented rewrite of previo
One essential component of operational space weather forecasting is the prediction of solar flares. With a multitude of flare forecasting methods now available online it is still unclear which of these methods performs best, and none are substantially better than climatological forecasts. Space weather researchers are increasingly looking towards methods used by the terrestrial weather community to improve current forecasting techniques. Ensemble forecasting has been used in numerical weather prediction for many years as a way to combine different predictions in order to obtain a more accurate result. Here we construct ensemble forecasts for major solar flares by linearly combining the full-disk probabilistic forecasts from a group of operational forecasting methods (ASAP, ASSA, MAG4, MOSWOC, NOAA, and MCSTAT). Forecasts from each method are weighted by a factor that accounts for the methods ability to predict previous events, and several performance metrics (both probabilistic and categorical) are considered. It is found that most ensembles achieve a better skill metric (between 5% and 15%) than any of the members alone. Moreover, over 90% of ensembles perform better (as measured by forecast attributes) than a simple equal-weights average. Finally, ensemble uncertainties are highly dependent on the internal metric being optimized and they are estimated to be less than 20% for probabilities greater than 0.2. This simple multi-model, linear ensemble technique can provide operational space weather centres with the basis for constructing a versatile ensemble forecasting system -- an improved starting point to their forecasts that can be tailored to different end-user needs.
We investigate the anisotropy of Alfvenic turbulence in the inertial range of slow solar wind and in both driven and decaying reduced magnetohydrodynamic simulations. A direct comparison is made by measuring the anisotropic second-order structure functions in both data sets. In the solar wind, the perpendicular spectral index of the magnetic field is close to -5/3. In the forced simulation, it is close to -5/3 for the velocity and -3/2 for the magnetic field. In the decaying simulation, it is -5/3 for both fields. The spectral index becomes steeper at small angles to the local magnetic field direction in all cases. We also show that when using the global rather than local mean field, the anisotropic scaling of the simulations cannot always be properly measured.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
