A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics

We present a highly robust second order accurate scheme for the Euler equations and the ideal MHD equations. The scheme is of predictor–corrector type, with a MUSCL scheme following as a special case. The crucial ingredients are an entropy stable approximate Riemann solver and a new spatial reconstruction that ensures positivity of mass density and pressure. For multidimensional MHD, a new discrete form of the Powell source terms is vital to ensure the stability properties. The numerical examples show that the scheme has superior stability compared to standard schemes, while maintaining accuracy. In particular, the method can handle very low values of pressure (i.e. low plasma β$\beta$ or high Mach numbers) and low mass densities.     

An unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations

Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [J.A. Rossmanith, An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM J. Sci. Comput. 28 (2006) 1766], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J.O. Langseth, R.J. LeVeque, A wave propagation method for threedimensional hyperbolic conservation laws, J. Comput. Phys. 165 (2000) 126]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The resulting scheme is applied to several numerical test cases.     

An adaptive version of ghost-cell immersed boundary method for incompressible flows with complex stationary and moving boundaries

An adaptive version of immersed boundary method for simulating flows with complex stationary and moving boundaries is presented.The method employs a ghost-cell methodology which allows for a sharp representation of the immersed boundary.To simplify the implementation of the methodology,a volume-of-fluid method is introduced to identify the immersed boundary.In addition,the domain is spatially discretized using a tree-based discretization which is relatively simple to implement a fully flexible adaptive refi...     

An interaction potential based lattice Boltzmann method with adaptive mesh refinement (AMR) for two-phase flow simulation

The lattice Boltzmann method (LBM) for two-phase flow simulation is often hindered by insufficient resolution at the interface. As a result, the LBM simulation of bubbles in bubbling flows is commonly limited to spherical or slightly deformed bubble shapes. In this study, the adaptive mesh refinement method for the LBM is developed to overcome such a problem. The approach for this new method is based on the improved interaction potential model, which is able to maintain grid-independent fluid properties in the two-fluid phases and at the interface. The LBM–AMR algorithm is described, especially concerning the LBM operation on a non-uniform mesh and the improved interaction potential model. Numerical simulations have been performed to validate the method in both single phase and multiphase flows. The 2D and 3D simulations of the buoyant rise of bubbles are conducted under various conditions. The agreement between the simulated bubble shape and velocity with experiments illustrates the capability of the LBM–AMR approach in predicting bubble dynamics even under the large bubble deformation conditions. Further, the LBM–AMR technique is capable of simulating a complex topology change of the interface. Integration of LBM with AMR can significantly improve the accuracy and reduce computation cost. The method developed in this study may appreciably enhance the capability of LBM in the simulation of complex multiphase flows under realistic conditions.     

Simulations of multiphase flows with multiple length scales using moving mesh interface tracking with adaptive meshing

In multiphase flows, the length scales of thin regions, such as thin films between nearly touching drops and thin threads formed during the interface pinch-off, are usually several orders of magnitude smaller than the size of the drops. In this paper, a number of extra length criteria for adaptive meshes are developed and implemented in the moving mesh interface tracking method to solve these multiple-length-scale problems with high fidelity. A nominal length scale based on the solutions of Laplace’s equations with the unit normal vectors of surfaces as the boundary conditions is proposed for the adaptive mesh refinement in the thin regions. For almost flat interfaces/boundaries which are near to the thin regions, the averaged length of the interior edges sharing the two nodes with the boundary edge is introduced for the mesh adaptation. The averaged length of the interfacial edges is used for the interior elements near the interfaces but outside of the thin regions. For the interior mesh away from the interfaces/boundaries, different averaged length scales based on the initial mesh are employed for the adaptive mesh refining and coarsening. Numerous cases are simulated to demonstrate the capability of the proposed schemes in handling multiple length scales, which include the relaxation and necking of an elongated droplet, droplet–droplet head-on approaching, droplet-wall interactions, and a droplet pair in a shear flow. The smallest length resolved for the thin regions is three orders of magnitude smaller than the largest characteristic length of the problem.     

Feature-driven Cartesian adaptive mesh refinement for vortex-dominated flows

We develop locally normalized feature-detection methods to guide the adaptive mesh refinement (AMR) process for Cartesian grid systems to improve the resolution of vortical features in aerodynamic wakes. The methods include: the Q-criterion [1], the λ₂ method [2], the λ_ci method [3], and the λ₊ method [4]. Specific attention is given to automate the feature identification process by applying a local normalization based upon the shear-strain rate so that they can be applied to a wide range of flow-fields without the need for user intervention. To validate the methods, we assess tagging efficiency and accuracy using a series of static vortex-dominated flow-fields, and use the methods to drive the AMR process for several theoretical and practical simulations. We demonstrate that the adaptive solutions provide comparable accuracy to solutions obtained on uniformly refined meshes at a fraction of the computational cost. Overall, the normalized feature detection methods are shown to be effective in driving the AMR process in an automated and efficient manner.     

Compact integration factor methods for complex domains and adaptive mesh refinement

Implicit integration factor (IIF) method, a class of efficient semi-implicit temporal scheme, was introduced recently for stiff reaction–diffusion equations. To reduce cost of IIF, compact implicit integration factor (cIIF) method was later developed for efficient storage and calculation of exponential matrices associated with the diffusion operators in two and three spatial dimensions for Cartesian coordinates with regular meshes. Unlike IIF, cIIF cannot be directly extended to other curvilinear coordinates, such as polar and spherical coordinates, due to the compact representation for the diffusion terms in cIIF. In this paper, we present a method to generalize cIIF for other curvilinear coordinates through examples of polar and spherical coordinates. The new cIIF method in polar and spherical coordinates has similar computational efficiency and stability properties as the cIIF in Cartesian coordinate. In addition, we present a method for integrating cIIF with adaptive mesh refinement (AMR) to take advantage of the excellent stability condition for cIIF. Because the second order cIIF is unconditionally stable, it allows large time steps for AMR, unlike a typical explicit temporal scheme whose time step is severely restricted by the smallest mesh size in the entire spatial domain. Finally, we apply those methods to simulating a cell signaling system described by a system of stiff reaction–diffusion equations in both two and three spatial dimensions using AMR, curvilinear and Cartesian coordinates. Excellent performance of the new methods is observed.     

A moving Kriging interpolation-based boundary node method for two-dimensional potential problems

In this paper,a meshfree boundary integral equation(BIE) method,called the moving Kriging interpolationbased boundary node method(MKIBNM),is developed for solving two-dimensional potential problems.This study combines the BIE method with the moving Kriging interpolation to present a boundary-type meshfree method,and the corresponding formulae of the MKIBNM are derived.In the present method,the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker’s delta property,then the boundary conditions can be imposed directly and easily.To verify the accuracy and stability of the present formulation,three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.     

A meshless method based on moving Kriging interpolation for a two-dimensional time-fractional diffusion equation

Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.     

An improved boundary element-free method （IBEFM） for two-dimensional potential problems

The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker δ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker δ function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.     

An Exner-based coupled model for two-dimensional transient flow over erodible bed

Transient flow over erodible bed is solved in this work assuming that the dynamics of the bed load problem is described by two mathematical models: the hydrodynamic model, assumed to be well formulated by means of the depth averaged shallow water equations, and the Exner equation. The Exner equation is written assuming that bed load transport is governed by a power law of the flow velocity and by a flow/sediment interaction parameter variable in time and space. The complete system is formed by four coupled partial differential equations and a genuinely Roe-type first order scheme has been used to solve it on triangular unstructured meshes. Exact solutions have been derived for the particular case of initial value Riemann problems with variable bed level and depending on particular forms of the solid discharge formula. The model, supplied with the corresponding solid transport formulae, is tested by comparing with the exact solutions. The model is validated against laboratory experimental data of different unsteady problems over erodible bed.     

An efficient three-dimensional adaptive quasicontinuum method using variable-node elements

A new quasicontinuum (QC) implementation using the so-called “variable-node finite elements” is reported in this work. Tetrahedral elements, which have been exclusively utilized for the conventional QC are replaced by hexahedral elements in conjunction with the so-called variable-node elements. This enables an effective adaptive mesh refinement in QC, leading to fast and efficient simulations compared with the conventional QC. To confirm the solution accuracy, comparison is made for a nanoindentation problem with a molecular dynamics simulation as well as a molecular mechanics solution. Further examples of nanoindentation are shown and discussed to demonstrate the effectiveness of the present scheme.     

一种混响背景下的自适应动目标检测方法

本文提出了一种混响和噪声干扰下动目标检测的自适应方法,其结构为自适应混响抵消器（ARC）加推广的自适应相干累积器（GACI）。用自适应混响抵消实现干扰背景的预白化,用推广的自适应相干累积实现匹配滤波的功能。该方法不需要混响干扰的先验知识,也不需要知道信号参量,就可以实现混响干扰背景下动目标的检测。仿真和试验表明,其具有良好的动目标检测性能,且具有很好的海洋信道适应能力。     

Novel energy dissipative method on the adaptive spatial discretization for the Allen–Cahn equation

We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids. For spatial discretization, the classical central difference method is utilized, while the average vector field method is applied for time discretization. Compared with the average vector field method on the uniform mesh, the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation. This is due to the moving mesh method, which can concentrate the grid points more densely where the solution changes drastically. Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.     

An adaptive regression method for infrared blind-pixel compensation

Blind pixel compensation is an ill-posed inverse problem of infrared imaging systems and image restoration. The performance of a blind pixel compensation algorithm depends on the accuracy of estimation for the underlying true infrared images. We propose an adaptive regression method (ARM) for blind pixel compensation that integrates the multi-scale framework with a regression model. A blind-pixel is restored by exploiting the intra-scale properties through the nonparametric regressive estimation and the inter-scale characteristics via parametric regression for continuous learning. Combining the respective strengths of a parametric model and a nonparametric model, ARM establishes a set of multi-scale blind-pixel compensation method to correct the non-uniformity based on key frame extraction. Therefore, it is essentially different from the traditional frameworks for blind pixel compensation which are based on filtering and interpolation. Experimental results on some challenging cases of blind compensation show that the proposed algorithm outperforms existing methods by a significant margin in both isolated blind restoration and clustered blind restoration.     

An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems

Heterogeneous anisotropic diffusion problems arise in the various areas of science and engineering including plasma physics, petroleum engineering, and image processing. Standard numerical methods can produce spurious oscillations when they are used to solve those problems. A common approach to avoid this difficulty is to design a proper numerical scheme and/or a proper mesh so that the numerical solution validates the discrete counterpart (DMP) of the maximum principle satisfied by the continuous solution. A well known mesh condition for the DMP satisfaction by the linear finite element solution of isotropic diffusion problems is the non-obtuse angle condition that requires the dihedral angles of mesh elements to be non-obtuse. In this paper, a generalization of the condition, the so-called anisotropic non-obtuse angle condition, is developed for the finite element solution of heterogeneous anisotropic diffusion problems. The new condition is essentially the same as the existing one except that the dihedral angles are now measured in a metric depending on the diffusion matrix of the underlying problem. Several variants of the new condition are obtained. Based on one of them, two metric tensors for use in anisotropic mesh generation are developed to account for DMP satisfaction and the combination of DMP satisfaction and mesh adaptivity. Numerical examples are given to demonstrate the features of the linear finite element method for anisotropic meshes generated with the metric tensors.     

A fourth-order divergence-free method for MHD flows

This paper extends our previous third-order method [S. Li, High order central scheme on overlapping cells for magneto-hydrodynamic flows with and without constrained transport method, J. Comput. Phys. 227 (2008) 7368–7393] to the fourth-order. Central finite-volume schemes on overlapping grid are used for both the volume-averaged variables and the face-averaged magnetic field. The magnetic field at the cell boundaries falls within the dual grid and is naturally continuous so that our method eliminates the instability triggered by the discontinuity in the normal component of the magnetic field. Our fourth-order scheme has much smaller numerical dissipation than the third-order scheme. The divergence-free condition of the magnetic field is preserved by our fourth-order divergence-free reconstruction and the constrained transport method. Numerical examples show that the divergence-free condition is essential to the accuracy of the method when a limiter is used in the reconstruction. The high-order, low-dissipation, and divergence-free properties of this method make it an ideal tool for direct magneto-hydrodynamic turbulence simulations.     

A new electromagnetic particle-in-cell model with adaptive mesh refinement for high-performance parallel computation

A new electromagnetic particle-in-cell (EMPIC) model with adaptive mesh refinement (AMR) has been developed to achieve high-performance parallel computation in distributed memory system. For minimizing the amount and frequency of inter-processor communications, the present study uses the staggering grid scheme with the charge conservation method, which consists only of the local operations. However, the scheme provides no numerical damping for electromagnetic waves regardless of the wavenumber, which results in significant noise in the refinement region that eventually covers over physical signals. In order to suppress the electromagnetic noise, the present study introduces a smoothing method which gives numerical damping preferentially for short wavelength modes. The test simulations show that only a weak smoothing results in drastic reduction in the noise, so that the implementation of the AMR is possible in the staggering grid scheme. The computational load balance among the processors is maintained by a new method termed the adaptive block technique for the domain decomposition parallelization. The adaptive block technique controls the subdomain (block) structure dynamically associated with the system evolution, such that all the blocks have almost the same number of particles. The performance of the present code is evaluated for the simulations of the current sheet evolution. The test simulations demonstrate that the usage of the adaptive block technique as well as the staggering grid scheme enhances significantly the parallel efficiency of the AMR-EMPIC model.     

An improved complex variable element-free Galerkin method for two-dimensional elasticity problems

In this paper, the improved complex variable moving least-squares (ICVMLS) approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares (CVMLS) approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin (EFG) method and the ICVMLS approximation, the improved complex variable element-free Galerkin (ICVEFG) method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFG method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.