Solving wave propagation problems within heterogeneous media has been of great interest and has a wide range of applications in physics and engineering. The design of numerical methods for such general wave propagation problems is challenging because the energy conserving property has to be incorporated in the numerical algorithms in order to minimize the phase or shape errors after long time integration. In this paper, we focus on multi-dimensional wave problems and consider linear second-order wave equation in heterogeneous media. We develop and analyze an LDG method, in which numerical fluxes are carefully designed to maintain the energy preserving property and accuracy. Compatible high order energy conserving time integrators are also proposed. The optimal error estimates and the energy conserving property are proved for the semi-discrete methods. Our numerical experiments demonstrate optimal rates of convergence, and show that the errors of the numerical solutions do not grow significantly in time due to the energy conserving property.

A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.

We consider an inverse time-dependent source problem governed by a distributed time-fractional diffusion equation using interior measurement data. Such a problem arises in some ultra-slowdiffusion phenomena in many applied areas. Based on the regularity result of the solution to the direct problem, we establish the solvability of this inverse problem as well as the conditional stability in suitable function space with a weak norm. By a variational identity connecting the unknown time-dependent source and the interior measurement data, the conjugate gradient method is also introduced to construct the inversion algorithm under the framework of regularizing scheme. We show the validity of the proposed scheme by several numerical examples.

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.

No abstract uploaded!