Efficient Simulation of Wave Propagation with Implicit Finite Difference Schemes

Finite difference method is an important methodology in the approximation of waves.In this paper,we will study two implicit finite difference schemes for the simulation of waves.They are the weighted alternating direction implicit(ADI)scheme and the locally one-dimensional(LOD)scheme.The approximation errors,stability conditions,and dispersion relations for both schemes are investigated.Our analysis shows that the LOD implicit scheme has less dispersion error than that of the ADI scheme.Moreover,the unconditional stability for both schemes with arbitrary spatial accuracy is established for the first time.In order to improve computational efficiency,numerical algorithms based on message passing interface(MPI)are implemented.Numerical examples of wave propagation in a three-layer model and a standard complex model are presented.Our analysis and comparisons show that both ADI and LOD schemes are able to efficiently and accurately simulate wave propagation in complex media.

Homogenization with the Quasistatic Tresca Friction Law:Qualitative and Quantitative Results

Modeling of frictional contacts is crucial for investigating mechanical performances of composite materials under varying service environments.The paper considers a linear elasticity system with strongly heterogeneous coefficients and quasistatic Tresca friction law,and studies the homogenization theories under the frameworks of H-convergence and small ε-periodicity.The qualitative result is based on H-convergence,which shows the original oscillating solutions will converge weakly to the homogenized solution,while the author’s quantitative result provides an estimate of asymptotic errors in H^(1)-norm for the periodic homogenization.This paper also designs several numerical experiments to validate the convergence rates in the quantitative analysis.

A Sub-Grid Structure Enhanced Discontinuous Galerkin Method for Multiscale Diffusion and Convection-Diffusion Problems

In this paper,we present an efficient computational methodology for diffusion and convection-diffusion problems in highly heterogeneous media as well as convection-dominated diffusion problem.It is well known that the numerical computation for these problems requires a significant amount of computermemory and time.Nevertheless,the solutions to these problems typically contain a coarse component,which is usually the quantity of interest and can be represented with a small number of degrees of freedom.There are many methods that aim at the computation of the coarse component without resolving the full details of the solution.Our proposed method falls into the framework of interior penalty discontinuous Galerkin method,which is proved to be an effective and accurate class of methods for numerical solutions of partial differential equations.A distinctive feature of our method is that the solution space contains two components,namely a coarse space that gives a polynomial approximation to the coarse component in the traditional way and a multiscale space which contains sub-grid structures of the solution and is essential to the computation of the coarse component.In addition,stability of the method is proved.The numerical results indicate that the method can accurately capture the coarse behavior of the solution for problems in highly heterogeneous media as well as boundary and internal layers for convection-dominated problems.

A TV-Based Iterative Regularization Method for the Solutions of Thermal Convection Problems

Linear/nonlinear and Stokes based-stabilizations for the filter equations for damping out primitive variable(PV)solutions corrupted by uniformly distributed randomnoises are numerically studied through the natural convection(NC)aswell as the mixed convection(MC)environment.The most recognizable filter-scheme is based on a combination of the negative Laplace equation multiplied with the selection of the spatial scale and a linear function in order to preserve the uniqueness of the filtered solution.A more complicated filter-scheme,based on a Stokes problem which couples a filtered velocity and a filtered(artificial)pressure(or Lagrange multiplier)in order to enforce the incompressibility constraint,is also studied.Linear and Stokes basedfilters via nested iterative(NI)filters and the consistent splitting scheme(CSS)are proposed for the NC/MC problems.Inspired by the total-variation(TV)model of image diffusion,well preserved feature flow patterns from the corrupted NC/MC environment are obtained by TV-Stokes based-filters together with the CSS.Our experimental results show that our proposed algorithms are effective and efficient in eliminating the unwanted spurious oscillations and preserving the accuracy of thermal convective fluid flows.