The diagnostic of the performance of numerical methods for physical models, like those in computational fluid mechanics and other fields of continuum mechanics, rely on the preservation of statistical moments of extensive quantities. Dynamic and adaptive meshing often use interpolations to represent fields over a new set of elements and require to be conservative and moment-preserving. Denoising algorithms should not affect moment distributions of data. And numerical deltas are described using the number of moments preserved. Therefore, all these methodologies benefit from the use of moment-preserving interpolations. In this article, we review the presentation of the piecewise polynomial basis functions that provide moment-preserving interpolations, better described as the collocation basis of compact finite differences, or Z-splines. We present different applications of these basis functions that show the improvement of numerical algorithms for fluid mechanics, discrete delta functions and denoising. We also provide theorems of the extension of the properties of the basis, previously known as the Strang and Fix theory, to the case of arbitrary knot partitions.

This article deals with the derivation of an adaptive numerical method for mono-dimensional kinetic equations for gas mixtures. For classical deterministic kinetic methods, the velocity domain is chosen accordingly to the initial condition. In such methods, this velocity domain is the same for all time, all space points and all species. The idea developed in this article relies on defining velocity domains that depend on space, time and species. This allows the method to locally adapt to the support of the distribution functions. The method consists in computing macroscopic quantities by the use of conservation laws, which enables the definition of such local grids. Then, an interpolation procedure along with a upwind scheme is performed in order to treat the advection term, and an implicit treatment of the BGK operator allows for the derivation of an AP scheme, where the stability condition is independent of the relaxation rate. The method is then applied to a series of test cases and compared to the classical DVM method.

In this paper, the RSEL (Random Subfeature Ensemble Learning) algorithm is proposed to improve the forecast results of weather forecasting. Based on the classical machine learning algorithms, RSEL algorithm integrates random subfeature selection and ensemble learning combination strategy to enhance the diversity of the features and avoid the influence of a small number of unstable outliers generated randomly. Furthermore, the feature engineering schemes are designed for the weather forecast data to make full use of spatial or temporal context. RSEL algorithm is tested by forecasting the wind speed and direction, and it improves the forecast accuracy of traditional methods and has good robustness.

The high-order gas-kinetic scheme (HGKS) has achieved success in simulating compressible flows with Cartesian meshes. To study the flow problems in general geometries, such as the flow over a wing-body, the development of HGKS in general curvilinear coordinates becomes necessary. In this paper, a two-stage fourth-order gas-kinetic scheme is developed for the Euler and Navier-Stokes solutions in the curvilinear coordinates from one-dimensional to three-dimensional computations. Based on the coordinate transformation, the kinetic equation is transformed first to the computational space, and the flux function in the gas-kinetic scheme is obtained there and is transformed back to the physical domain for the update of flow variables inside each control volume. To achieve the expected order of accuracy, the dimension-by-dimension reconstruction based on the WENO scheme is adopted in the computational domain, where the reconstructed variables are the cell averaged Jacobian and the Jacobian-weighted conservative variables. In the two-stage fourth-order gas-kinetic scheme, the point values as well as the spatial derivatives of conservative variables at Gaussian quadrature points have to be used in the evaluation of the time dependent flux function. The point-wise conservative variables are obtained by ratio of the above reconstructed data, and the spatial derivatives are reconstructed through orthogonalization in physical space and chain rule. A variety of numerical examples from the accuracy tests to the solutions with strong discontinuities are presented to validate the accuracy and robustness of the current scheme for both inviscid and viscous flows. The precise satisfaction of the geometrical conservation law in non-orthogonal mesh is also demonstrated through the numerical example.

We have extended the Helfrich's spontaneous curvature model [M. Iwamoto and Z. C. Ou-Yang. Chem. Phys. Lett. 590(2013)183; Y. X. Deng, et al., EPL. 123(2018)68002] of the equilibrium vesicle shapes by adding the interaction between magnetic field and the constituent molecules to explain the phenomena of the reversibly deformation of artificial stomatocyte [P. G. van Rhee, et al., Nat. Commun. Sep 24;5:5010(2014), doi: 10.1038/ncomms6010] and the anharmonic deformation of a self-assembled nanocapsules of bola-amphiphilic molecules and the linear birefringence [O.V. Manyuhina, et al., Phys. Rev. Lett. 98(2007)146101]. However, the sophisticated mathematics in differential geometry is still covered. Here, we present the derivations of formulas in detail to reveal the perturbation of deformation $ψ$ under two cases. New features such as the influence of temperature on the bend modulus of vesicle membrane have been revealed.

The elastic transmission eigenvalue problem has important applications in the inverse elastic scattering theory. Recently, the numerical computation for this problem has attracted the attention of the researchers. In this paper, we propose the $H^2$-conforming methods including the classical $H^2$-conforming finite element method and the spectral element method, and establish the two-grid discretization scheme. Theoretical analysis and numerical experiments show that the methods presented in this paper can efficiently compute real and complex elastic transmission eigenvalues.

Here, we develop a first and a second order time stepping schemes for a binary fluid-surfactant phase field model by using the scalar auxiliary variable approach. The free energy contains a double-well potential, a nonlinear coupling entropy and a Flory-Huggins potential. The resulting coupled system consists of a Cahn-Hilliard type equation and a Wasserstein type equation which leads to a degenerate problem. By introducing only one scalar auxiliary variable, the system is transformed into an equivalent form so that the nonlinear terms can be treated semi-explicitly. Both the schemes are linear and decoupled, thus they can be solved efficiently. We further prove that these semi-discretized schemes in time are unconditionally energy stable. Some numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.

A kernel-independent treecode (KITC) is presented for fast summation of particle interactions. The method employs barycentric Lagrange interpolation at Chebyshev points to approximate well-separated particle-cluster interactions. The KITC requires only kernel evaluations, is suitable for non-oscillatory kernels, and relies on the scale-invariance property of barycentric Lagrange interpolation. For a given level of accuracy, the treecode reduces the operation count for pairwise interactions from $\mathcal{O}$($N^2$) to $\mathcal{O}$($N$log$N$), where $N$ is the number of particles in the system. The algorithm is demonstrated for systems of regularized Stokeslets and rotlets in 3D, and numerical results show the treecode performance in terms of error, CPU time, and memory consumption. The KITC is a relatively simple algorithm with low memory consumption, and this enables a straightforward OpenMP parallelization.

In this paper, we focus on the numerical simulation of the two-layer shallow water equations over variable bottom topography. Although the existing numerical schemes for the single-layer shallow water equations can be extended to two-layer shallow water equations, it is not a trivial work due to the complexity of the equations. To achieve the well-balanced property of the numerical scheme easily, the two-layer shallow water equations are reformulated into a new form by introducing two auxiliary variables. Since the new equations are only conditionally hyperbolic and their eigenstructure cannot be easily obtained, we consider the utilization of the central discontinuous Galerkin method which is free of Riemann solvers. By choosing the values of the auxiliary variables suitably, we can prove that the scheme can exactly preserve the still-water solution, and thus it is a truly well-balanced scheme. To ensure the non-negativity of the water depth, a positivity-preserving limiter and a special approximation to the bottom topography are employed. The accuracy and validity of the numerical method will be illustrated through some numerical tests.

In this paper, a high-order cell-centered discontinuous Galerkin (DG) multi-material arbitrary Lagrangian-Eulerian (MMALE) method is developed for compressible fluid dynamics. The MMALE method utilizes moment-of-fluid (MOF) interface reconstruction technology to simulate multi-materials of immiscible fluids. It is an explicit time-marching Lagrangian plus remap type. In the Lagrangian phase, an updated high-order discontinuous Galerkin Lagrangian method is applied for the discretization of hydrodynamic equations, and Tipton's pressure relaxation closure model is used in the mixed cells. A robust moment-of-fluid interface reconstruction algorithm is used to provide the information of the material interfaces for remapping. In the rezoning phase, Knupp's algorithm is used for mesh smoothing. For the remapping phase, a high-order accurate remapping method of the cell-intersection-based type is proposed. It can be divided into four stages: polynomial reconstruction, polygon intersection, integration, and detection of problematic cells and limiting. Polygon intersection is based on the "clipping and projecting" algorithm, and detection of problematic cells depends on a troubled cell marker, and a posteriori multi-dimensional optimal order detection (MOOD) limiting strategy is used for limiting. Numerical tests are given to demonstrate the robustness and accuracy of our method.

The mechanism of jet flow expanding into vacuum environment (or extremely low density environment) is important for the propulsion unit of micro-electro-mechanical systems (MEMS), the thruster of spacecraft, the attitude control system of satellite, etc.. Since its flow field is often composed of local continuum region and local rarefied region, the jet flow into vacuum has noteworthy multi-scale transportation behaviors. Therefore, the numerical study of such flows needs the multi-scale schemes which are valid for both continuum and rarefied flows. In the past few years, a series of unified methods for whole flow regime (from continuum regime to rarefied regime) have been developed from the perspective of the direct modeling, and have been verified by sufficient test cases. In this paper, the compressible conserved discrete unified gas-kinetic scheme is further developed and is utilized for predicting the jet flows into vacuum environment. In order to cover the working conditions of both aerospace and MEMS applications, the jet flows with a wide range of inlet Knudsen (Kn) numbers (from 1E-4 to 100) are considered. The evolution of flow field during the entire startup and shutdown process with Kn number 100 is predicted by the present method, and it matches well with the result of analytical collisionless Boltzmann equation. For Kn numbers from 1E-4 to 10, the flow field properties such as density, momentum, and pressure are investigated, and the results are provided in details, since the published results are not sufficient at the present stage. The extent and intensity of the jet flow influence are especially investigated, because they are strongly related to the plume contamination and momentum impact on objects facing the jet, such as the solar paddles which face the attitude control thruster during the docking process.

In this paper we consider the numerical solution of the Allen-Cahn type diffuse interface model in a polygonal domain. The intersection of the interface with the re-entrant corners of the polygon causes strong corner singularities in the solution. To overcome the effect of these singularities on the accuracy of the approximate solution, for the spatial discretization we develop an efficient finite element method with exponential mesh refinement in the vicinity of the singular corners, that is based on ($k$−1)-th order Lagrange elements, $k$≥2 an integer. The problem is fully discretized by employing a first-order, semi-implicit time stepping scheme with the Invariant Energy Quadratization approach in time, which is an unconditionally energy stable method. It is shown that for the error between the exact and the approximate solution, an accuracy of $\mathcal{O}$($h^k$+$τ$) is attained in the $L^2$-norm for the number of $\mathcal{O}$($h^{−2}$ln$h^{−1}$) spatial elements, where $h$ and $τ$ are the mesh and time steps, respectively. The numerical results obtained support the analysis made.

In this paper, we propose some effective one- and two-level domain decomposition preconditioners for elastic crack problems modeled by extended finite element method. To construct the preconditioners, the physical domain is decomposed into the "crack tip" subdomain, which contains all the degrees of freedom (dofs) of the branch enrichment functions, and the "regular" subdomains, which contain the standard dofs and the dofs of the Heaviside enrichment function. In the one-level additive Schwarz and restricted additive Schwarz preconditioners, the "crack tip" subproblem is solved directly and the "regular" subproblems are solved by some inexact solvers, such as ILU. In the two-level domain decomposition preconditioners, traditional interpolations between the coarse and the fine meshes destroy the good convergence property. Therefore, we propose an unconventional approach in which the coarse mesh is exactly the same as the fine mesh along the crack line, and adopt the technique of a non-matching grid interpolation between the fine and the coarse meshes. Numerical experiments demonstrate the effectiveness of the two-level domain decomposition preconditioners applied to elastic crack problems.

We study an identification problem which estimates the parameters of the underlying random distribution for uncertain scalar conservation laws. The hyperbolic equations are discretized with the so-called discontinuous stochastic Galerkin method, i.e., using a spatial discontinuous Galerkin scheme and a Multielement stochastic Galerkin ansatz in the random space. We assume an uncertain flux or uncertain initial conditions and that a data set of an observed solution is given. The uncertainty is assumed to be uniformly distributed on an unknown interval and we focus on identifying the correct endpoints of this interval. The first-order optimality conditions from the discontinuous stochastic Galerkin discretization are computed on the time-continuous level. Then, we solve the resulting semi-discrete forward and backward schemes with the Runge-Kutta method. To illustrate the feasibility of the approach, we apply the method to a stochastic advection and a stochastic equation of Burgers' type. The results show that the method is able to identify the distribution parameters of the random variable in the uncertain differential equation even if discontinuities are present.

A higher order interpolation scheme based on a multi-stage BVD (Boundary Variation Diminishing) algorithm is developed for the FV (Finite Volume) method on non-uniform, curvilinear structured grids to simulate the compressible turbulent flows. The designed scheme utilizes two types of candidate interpolants including a higher order linear-weight polynomial as high as eleven and a THINC (Tangent of Hyperbola for INterface Capturing) function with the adaptive steepness. We investigate not only the accuracy but also the efficiency of the methodology through the cost efficiency analysis in comparison with well-designed mapped WENO (Weighted Essentially Non-Oscillatory) scheme. Numerical experimentation including benchmark broadband turbulence problem as well as real-life wall-bounded turbulent flows has been carried out to demonstrate the potential implementation of the present higher order interpolation scheme especially in the ILES (Implicit Large Eddy Simulation) of compressible turbulence.