Search

Numerical schemes provably preserving the positivity of density and pressure are highly desirable for ideal magnetohydrodynamics (MHD), but the rigorous positivity-preserving (PP) analysis remains challenging. The difficulties mainly arise from the intrinsic complexity of the MHD equations as well as the indeterminate relation between the PP property and the divergence-free condition on the magnetic field. This paper presents the first rigorous PP analysis of conservative schemes with the Lax--Friedrichs (LF) flux for 1D and multidimensional ideal MHD. The significant innovation is the discovery of the theoretical connection between the PP property and a discrete divergence-free (DDF) condition. This connection is established through the generalized LF splitting properties, which are alternatives to the usually expected LF splitting property that does not hold for ideal MHD. The generalized LF splitting properties involve a number of admissible states strongly coupled by the DDF condition, making their derivation very difficult. We derive these properties via a novel equivalent form of the admissible state set and an important inequality, which is skillfully constructed by technical estimates. Rigorous PP analysis is then presented for finite volume and discontinuous Galerkin schemes with the LF flux on uniform Cartesian meshes. In the 1D case, the PP property is proved for the first-order scheme with proper numerical viscosity, and also for arbitrarily high-order schemes under conditions accessible by a PP limiter. In the 2D case, we show that the DDF condition is necessary and crucial for achieving the PP property. It is observed that even slightly violating the proposed DDF condition may cause failure to preserve the positivity of pressure. We prove that the 2D LF type scheme with proper numerical viscosity preserves both the positivity and the DDF condition. Sufficient conditions are derived for 2D PP high-order schemes, and extension to 3D is discussed. Numerical examples provided in the supplementary material further confirm the theoretical findings. [ABSTRACT FROM AUTHOR]

Many complex real-world systems can be accurately modeled by simulations. However, high-fidelity simulations may take hours or even days to compute. Because this can be impractical, a surrogate model is often used to approximate the dynamic behavior of the original simulator. This model can then be used as a cheap, drop-in replacement for the simulator. Because simulations can be very expensive, the data points, which are required to build the model, must be chosen as optimally as possible. Sequential design strategies offer a huge advantage over one-shot experimental designs because they can use information gathered from previous data points in order to determine I lie location of new data points. Each sequential design strategy must perform a trade-off between exploration and exploitation, where the former involves selecting data points in unexplored regions of the design space, while the latter suggests adding data points in regions which were previously identified to be interesting (for example, highly nonlinear regions). In this paper, a novel hybrid sequential design strategy is proposed which uses a Monte Carlo-based approximation of a Voronoi tessellation for exploration and local linear approximations of the simulator for exploitation. The advantage of this method over other sequential design methods is that it is independent of the model type, and can therefore be used in heterogeneous modeling environments, where multiple model types are used at the same time. The new method is demonstrated on a number of test problems, showing that it is a robust, competitive, and efficient sequential design strategy. [ABSTRACT FROM AUTHOR]

This paper is concerned with the numerical simulations for the dynamics of the molecular beam epitaxy (MBE) model. The numerical simulations of the MBE models require long time computations, and therefore large time-stepping methods become necessary. In this work, we consider some unconditionally energy stable finite difference schemes, which will be used in the time adaptivity strategies. It is found that the use of the time adaptivity cannot only resolve the steady-state solutions but also the dynamical changes of the solution accurately and efficiently. The adaptive time step is selected based on the energy variation or the change of the roughness of the solution. The numerical experiments demonstrated that the CPU time is significantly saved for long time simulations. [ABSTRACT FROM AUTHOR]

In this paper, we present a rigorous numerical analysis for a bounded artificial viscosity model (τ = μδσa(δ∥∇su∥F )∇su) for the numerical simulation of turbulent flows. In practice, the commonly used Smagorinsky model (τ = (csδ)2∥∇su∥F ∇su) is overly dissipative and yields unphysical results. To date, several methods for "clipping" the Smagorinsky viscosity have proven useful in improving the physical characteristics of the simulated flow. However, such heuristic strategies strongly rely upon a priori knowledge of the flow regime. The bounded artificial viscosity model relies on a highly nonlinear, but monotone and smooth, semilinear elliptic form for the artificial viscosity. For this model, we have introduced a variational computational strategy, provided finite element error convergence estimates, and included several computational examples indicating its improvement on the overly diffusive Smagorinsky model. [ABSTRACT FROM AUTHOR]