This paper proposes and analyzes arbitrarily high-order discontinuous Galerkin (DG) and finite volume methods which provably preserve the positivity of density and pressure for the ideal magnetohydrodynamics (MHD) on general meshes. Unified auxiliary theories are built for rigorously analyzing the positivity-preserving (PP) property of numerical MHD schemes with a Harten–Lax–van Leer (HLL) type flux on polytopal meshes in any space dimension. The main challenges overcome here include establishing certain relation between the PP property and a discrete divergence of magnetic field on general meshes, and estimating proper wave speeds in the HLL flux to ensure the PP property. In the 1D case, we prove that the standard DG and finite volume methods with the proposed HLL flux are PP, under a condition accessible by a PP limiter. For the multidimensional conservative MHD system, the standard DG methods with a PP limiter are not PP in general, due to the effect of unavoidable divergence error in the magnetic field. We construct provably PP high-order DG and finite volume schemes by proper discretization of the symmetrizable MHD system, with two divergence-controlling techniques: the locally divergence-free elements and suitably discretized Godunov–Powell source term. The former technique leads to zero divergence within each cell, while the latter controls the divergence error across cell interfaces. Our analysis reveals in theory that a coupling of these two techniques is very important for positivity preservation, as they exactly contribute the discrete divergence terms which are absent in standard multidimensional DG schemes but crucial for ensuring the PP property. Several numerical tests further confirm the PP property and the effectiveness of the proposed PP schemes. Unlike the conservative MHD system, the exact smooth solutions of the symmetrizable MHD system are proved to retain the positivity even if the divergence-free condition is not satisfied. Our analysis and findings further the understanding, at both discrete and continuous levels, of the relation between the PP property and the divergence-free constraint.

The gas dynamics equations, coupled with a static gravitational field, admit the hydrostatic balance where the flux produced by the pressure is exactly canceled by the gravitational source term. Many astrophysical problems involve the hydrodynamical evolution in a gravitational field, therefore it is essential to correctly capture the effect of gravitational force in the simulations. Improper treatment of the gravitational force can lead to a solution which either oscillates around the equilibrium, or deviates from the equilibrium after a long time run. In this paper we design high order well-balanced finite difference WENO schemes to this system, which can preserve the hydrostatic balance state exactly and at the same time can maintain genuine high order accuracy. Numerical tests are performed to verify high order accuracy, well-balanced property, and good resolution for smooth and discontinuous solutions. The main purpose of the well-balanced schemes designed in this paper is to well resolve small perturbations of the hydrostatic balance state on coarse meshes. The more difficult issue of convergence towards such hydrostatic balance state from an arbitrary initial condition is not addressed in this paper.

No abstract uploaded!

For cable bridges, the cable tension force plays a crucial role in their construction, assessment and long-term structural health monitoring. Cable tension forces vary in real time with the change of the moving vehicle loads and environmental effects, and this continual variation in tension force may cause fatigue damage of a cable. Traditional vibration-based cable tension force estimation methods can only obtain the time-averaged cable tension force and not the instantaneous force. This paper proposes a new approach to identify the time-varying cable tension forces of bridges based on an adaptive sparse time-frequency analysis method. This is a recently developed method to estimate the instantaneous frequency by looking for the sparsest time-frequency representation of the signal within the largest possible time-frequency dictionary (i.e. set of expansion functions). In the proposed approach, first, the time-varying modal frequencies are identified from acceleration measurements on the cable, then, the time-varying cable tension is obtained from the relation between this force and the identified frequencies. By considering the integer ratios of the different modal frequencies to the fundamental frequency of the cable, the proposed algorithm is further improved to increase its robustness to measurement noise. A cable experiment is implemented to illustrate the validity of the proposed method. For comparison, the Hilbert–Huang transform is also employed to identify the time-varying frequencies, which are then used to calculate the time-varying cable-tension force. The results show that the adaptive sparse time-frequency analysis method produces more accurate estimates of the time-varying cable tension forces than the Hilbert–Huang transform method.

In this paper, we propose a novel and simple technique to construct effective difference schemes for solving systems of singularly perturbed convection-diffusionreaction equations, whose solutions may display boundary or interior layers. We illustrate the technique by taking the Il’in-Allen-Southwell scheme for 1-D scalar equations as a basis to derive a formally second-order scheme for 1-D coupled systems and then extend the scheme to 2-D case by employing an alternating direction approach. Numerical examples are given to demonstrate the high performance of the obtained scheme on uniform meshes as well as piecewise-uniform Shishkin meshes.