High-Order Flux Correction for Viscous Flows on Arbitrary Unstructured Grids

A novel high-order method, termed flux correction, previously formulated for inviscid flows, is extended to viscous flows on arbitrary triangular grids. The correction method involves the addition of truncation error-canceling terms to the second-order linear Galerkin (node-centered finite volume) scheme to produce a third-order inviscid and fourth-order viscous scheme. The correction requires minimal modification of the underlying second-order scheme. As such, the method retains many of the advantages of traditional finite volume schemes, including robust shock capturing, low algorithmic complexity, and solver efficiency. In addition, we extend the scheme to unsteady flows. Verification and validation studies in two dimensions are presented. Significant improvement in accuracy is observed in all cases, with between 30–70 % increase in computational cost over a second-order finite volume method.     

A FV-TD electromagnetic solver using adaptive Cartesian grids

A second-order finite-volume (FV) method has been developed to solve the time-domain (TD) Maxwell equations, which govern the dynamics of electromagnetic waves. The computational electromagnetic (CEM) solver is capable of handling arbitrary grids, including structured, unstructured, and adaptive Cartesian grids, which are topologically arbitrary. It is argued in this paper that the adaptive Cartesian grid is better than a tetrahedral grid for complex geometries considering both efficiency and accuracy. A cell-wise linear reconstruction scheme is employed to achieve second-order spatial accuracy. Second-order time accuracy is obtained through a two-step Runge-Kutta scheme. Issues on automatic adaptive Cartesian grid generation such as cell-cutting and cell-merging are discussed. A multi-dimensional characteristic absorbing boundary condition (MDC-ABC) is developed at the truncated far-field boundary to reduce reflected waves from this artificial boundary. The CEM solver is demonstrated with several test cases with analytical solutions.     

Conservative compact difference scheme for the Zakharov–Rubenchik equations

In this article, a compact difference scheme is investigated to solve the Zakharov–Rubenchik equations in one dimension. The new scheme is proved to conserve the total mass and energy in the discrete sense. Rigorous error estimates are established for the new method with the help of an induction argument in energy space which show that the new scheme has second-order accuracy in time and fourth-order accuracy in space. Extensive numerical results are provided to verify our theoretical analysis, and show the accuracy and efficiency of the new scheme.     

An implicit compact scheme solver for two-dimensional multicomponent flows

A 2D implicit compact scheme solver has been implemented for the vorticity-velocity formulation in the case of nonreacting, multicomponent, axisymmetric, low Mach number flows. To stabilize the discrete boundary value problem, two sets of boundary closures are introduced to couple the velocity and vorticity fields. A Newton solver is used for solving steady-state and time-dependent equations. In this solver, the Jacobian matrix is formulated and stored in component form. To solve the system of linearized equations within each iteration of Newton’s method, preconditioned Bi-CGSTAB is used in combination with a matrix-vector product computed in component form. The almost dense Jacobian matrix is approximated by a partial Jacobian. For the preconditioner equation, the partial Jacobian is approximately factored using several methods. In a detailed study of several preconditioning techniques, a promising method based on ILUT preconditioning in combination with reordering and double scaling using the MC64 algorithm by Duff and Koster is selected. To validate the implicit compact scheme solver, several nonreacting model problems have been considered. At least third order accuracy in space is recovered on nonuniform grids. A comparison of the results of the implicit compact scheme solver with the results of a traditional implicit low order solver shows an order of magnitude reduction of computer memory and time using the compact scheme solver in the case of time-dependent mixing problems.     

A new ADI technique for two-dimensional parabolic equation with an integral condition

Nonclassical parabolic initial-boundary value problems arise in the study of several important physical phenomena. This paper presents a new approach to treat complicated boundary conditions appearing in the parabolic partial differential equations with nonclassical boundary conditions. A new fourth-order finite difference technique, based upon the Noye and Hayman (N-H) alternating direction implicit (ADI) scheme, is used as the basis to solve the two-dimensional time dependent diffusion equation with an integral condition replacing one boundary condition. This scheme uses less central processor time (CPU) than a second-order fully implicit scheme based on the classical backward time centered space (BTCS) method for two-dimensional diffusion. It also has a larger range of stability than a second-order fully explicit scheme based on the classical forward time centered space (FTCS) method. The basis of the analysis of the finite difference equations considered here is the modified equivalent partial differential equation approach, developed from the 1974 work of Warming and Hyeet. This allows direct and simple comparison of the errors associated with the equations as well as providing a means to develop more accurate finite difference methods. The results of numerical experiments for the new method are presented. The central processor times needed are also reported. Error estimates derived in the maximum norm are tabulated.     

Evaluation of multigrid acceleration for preconditioned time-accurate Navier-Stokes algorithms

The development of a two-dimensional time-accurate dual time step Navier-Stokes flow solver with time-derivative preconditioning and multigrid acceleration is described. The governing equations are integrated in time with both an explicit Runge-Kutta scheme and an implicit lower-upper symmetric-Gauss-Seidel scheme in a finite volume framework, yielding second-order accuracy in space and time. Issues concerning the implementation of multigrid for preconditioned, dual time step algorithms are discussed. Steady and unsteady computations were made of lid driven cavity flow, thermally driven cavity flow and pulsatile channel flow for a variety of conditions to validate the schemes and evaluate the effectiveness of multigrid for time-accurate simulations. Significant speedups were observed for steady and unsteady simulations. The speedups for unsteady simulations were problem dependent, a function of how rapidly the flow varied in time and the size of the allowable time step.     

Time-accurate Navier-Stokes simulation of vortex convection using an unstructured dynamic mesh procedure

A two-dimensional Navier-Stokes flow solver is developed for the simulation of unsteady flows on unstructured adaptive meshes. The solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with Roe’s flux-difference splitting. The Spalart-Allmaras one equation model is employed for the simulation of turbulence. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of spatial and temporal variations of the flowfield. Unsteady viscous flow for a traveling vortex in a free stream is simulated to validate the accuracy of the dynamic mesh adaptation procedure. Flow around a circular cylinder and two blade-vortex interaction problems are investigated for demonstration of the present method. Computed results show good agreement with existing experimental and computational results. It was found that unsteady time-accurate viscous flows can be accurately simulated using the present unstructured dynamic mesh adaptation procedure.     

A new fourth-order difference scheme for solving an N-carrier system with Neumann boundary conditions

In this paper, a numerical method is developed to solve an N-carrier system with Neumann boundary conditions. First, we apply the compact finite difference scheme of fourth order for discretizing spatial derivatives at the interior points. Then, we develop a new combined compact finite difference scheme for the boundary, which also has fourth-order accuracy. Lastly, by using a Padé approximation method for the resulting linear system of ordinary differential equations, a new compact finite difference scheme is obtained. The present scheme has second-order accuracy in time direction and fourth-order accuracy in space direction. It is shown that the scheme is unconditionally stable. The present scheme is tested by two numerical examples, which show that the convergence rate with respect to the spatial variable from the new scheme is higher and the solution is much more accurate when compared with those obtained by using other previous methods.     

A Galerkin finite element scheme for time–space fractional diffusion equation

In this paper, a Galerkin finite element scheme to approximate the time–space fractional diffusion equation is studied. Firstly, the fractional diffusion equation is transformed into a fractional Volterra integro-differential equation. And a second-order fractional trapezoidal formula is used to approach the time fractional integral. Then a Galerkin finite element method is introduced in space direction, where the semi-discretization scheme and fully discrete scheme are given separately. The stability analysis of semi-discretization scheme is discussed in detail. Furthermore, convergence analysis of semi-discretization scheme and fully discrete scheme are given in details. Finally, two numerical examples are displayed to demonstrate the effectiveness of the proposed method.     

Simple and accurate scheme for fluid velocity interpolation for Eulerian-Lagrangian computation of dispersed flows in 3D curvilinear grids

Particle tracking in turbulent flows in complex domains requires accurate interpolation of the fluid velocity field. If grids are non-orthogonal and curvilinear, the most accurate available interpolation methods fail. We propose an accurate interpolation scheme based on Taylor series expansion of the local fluid velocity about the grid point nearest to the desired location. The scheme is best suited for curvilinear grids with non-orthogonal computational cells. We present the scheme with second-order accuracy, yet the order of accuracy of the method can be adapted to that of the Navier-Stokes solver.An application to particle dispersion in a turbulent wavy channel is presented, for which the scheme is tested against standard linear interpolation. Results show that significant discrepancies can arise in the particle displacement produced by the two schemes, particularly in the near-wall region which is often discretized with highly-distorted computational cells.     

非结构化环境下的单目视觉可通行区域检测

目前三维避障主要采用三维激光雷达或者基于深度学习的障碍物识别,但前者价格昂贵,后者训练成本高且不稳定。为了稳定、鲁棒地实现低成本三维空间避障方案,提出了一种基于单目相机的可通行区域检测方法。该方法利用特征点标识障碍物,通过对相机高度和旋转平面加以约束,解决单目SLAM中的尺度不一致问题,并设计了障碍物距离求解器和代价求解器对小车前方区域特征点进行处理,计算出视觉代价地图,划分可通行区域。该方法优势在于仅需要低成本的单目相机便可完成可通行区域检测任务,可方便地移植到轻量级的移动设备,计算得到的视觉代价地图亦有利于小车后续的路径规划任务。在KITTI数据集上进行的实验表明,该方法的平均运算速度能达到20?frame/s,能满足小车实时避障的要求,对于单目SLAM的尺度恢复误差为2.5%~4.9%。     

Locally Limited and Fully Conserved RKDG2 Shallow Water Solutions with Wetting and Drying

This work extends a well-balanced second-order Runge-Kutta discontinuous Galerkin (RKDG2) scheme to provide conservative simulations for shallow flows involving wetting and drying over irregular topographies with friction effects. For this purpose, a wetting and drying technique designed originally for a finite volume (FV) scheme is improved and implemented, which includes the discretization of friction source terms via a splitting implicit integration approach. Another focus of this work is to design a fully conserved RKDG2 scheme to provide conservative solutions for both mass and momentum through a local slope limiting process. Several steady and transient benchmark tests with/without friction effects are simulated to validate the new solver and demonstrate the effects of different slope limiting processes, i.e. globally and locally slope limiting processes.     

Fourth-order runge-kutta schemes for fluid mechanics applications

Multiple high-order time-integration schemes are used to solve stiff test problems related to the Navier-Stokes (NS) equations. The primary objective is to determine whether high-order schemes can displace currently used second-order schemes on stiff NS and Reynolds averaged NS (RANS) problems, for a meaningful portion of the work-precision spectrum. Implicit-Explicit (IMEX) schemes are used on separable problems that naturally partition into stiff and nonstiff components. Non-separable problems are solved with fully implicit schemes, oftentimes the implicit portion of an IMEX scheme. The convection-diffusion-reaction (CDR) equations allow a term by term stiff/nonstiff partition that is often well suited for IMEX methods. Major variables in CDR converge at near design-order rates with all formulations, including the fourth-order IMEX additive Runge-Kutta (ARK₂) schemes that are susceptible to order reduction. The semi-implicit backward differentiation formulae and IMEX ARK₂ schemes are of comparable efficiency. Laminar and turbulent aerodynamic applications require fully implicit schemes, as they are not profitably partitioned. All schemes achieve design-order convergence rates on the laminar problem. The fourth-order explicit singly diagonally implicit Runge-Kutta (ESDIRK4) scheme is more efficient than the popular second-order backward differentiation formulae (BDF2) method. The BDF2 and fourth-order modified extended backward differentiation formulae (MEBDF4) schemes are of comparable efficiency on the turbulent problem. High precision requirements slightly favor the MEBDF4 scheme (greater than three significant digits). Significant order reduction plagues the ESDIRK4 scheme in the turbulent case. The magnitude of the order reduction varies with Reynolds number. Poor performance of the high-order methods can partially be attributed to poor solver performance. Huge time steps allowed by high-order formulations challenge the capabilities of algebraic solver technology.     

Fast algorithms for high-order numerical methods for space-fractional diffusion equations

In this paper, fast numerical methods for solving space-fractional diffusion equations are studied in two stages. Firstly, a fast direct solver for an implicit finite difference scheme proposed by Hao et al. [A fourth-order approximation of fractional derivatives with its applications, J. Comput. Phys. 281 (2015), pp. 787–805], which is fourth-order accurate in space and second-order accurate in time, is developed based on a circulant-and-skew-circulant (CS) representation of Toeplitz matrix inversion. Secondly, boundary value method with spatial discretization of Hao et al. [A fourth-order approximation of fractional derivatives with its applications, J. Comput. Phys. 281 (2015), pp. 787–805] is adopted to produce a numerical solution with higher order accuracy in time. Particularly, a method with fourth-order accuracy in both space and time can be achieved. GMRES method is employed for solving the discretized linear system with two preconditioners. Based on the CS representation of Toeplitz matrix inversion, the two preconditioners can be applied efficiently, and the convergence rate of the preconditioned GMRES method is proven to be fast. Numerical examples are given to support the theoretical analysis.     

A multigrid scheme for 3D Monge–Ampère equations*

The elliptic Monge–Ampère equation is a fully nonlinear partial differential equation which has been the focus of increasing attention from the scientific computing community. Fast three-dimensional solvers are needed, for example in medical image registration but are not yet available. We build fast solvers for smooth solutions in three dimensions using a nonlinear full-approximation storage multigrid method. Starting from a second-order accurate centred finite difference approximation, we present a nonlinear Gauss–Seidel iterative method which has a mechanism for selecting the convex solution of the equation. The iterative method is used as an effective smoother, combined with the full-approximation storage multigrid method. Numerical experiments are provided to validate the accuracy of the finite difference scheme and illustrate the computational efficiency of the proposed multigrid solver.     

Long time stability and convergence rate of MacCormack rapid solver method for nonstationary Stokes–Darcy problem

We propose and study a combination of two second-order implicit–explicit (IMEX) methods for the coupled Stokes–Darcy system that governs flows in karst aquifers. The first is a second-order explicit two-step MacCormack scheme and the second is a second-order implicit Crank–Nicolson method. Both algorithms only require the solution of two decoupled problems at each time step, one Stokes and the other Darcy. This combination so called the MacCormack rapid solver method is very efficient (faster, at least of second order accuracy in time and space) and can be easily implemented using legacy codes. Under time step limitation of the form $\Delta t\le Ch$ (where $h,$$\Delta t$ are mesh size and time step, respectively, and $C$ is a physical parameter) we prove both long time stability and the rate of convergence of the method. Some numerical experiments are presented and discussed.     

Fourth-Order Runge–Kutta Schemes for Fluid Mechanics Applications

Multiple high-order time-integration schemes are used to solve stiff test problems related to the Navier–Stokes (NS) equations. The primary objective is to determine whether high-order schemes can displace currently used second-order schemes on stiff NS and Reynolds averaged NS (RANS) problems, for a meaningful portion of the work-precision spectrum. Implicit–Explicit (IMEX) schemes are used on separable problems that naturally partition into stiff and nonstiff components. Non-separable problems are solved with fully implicit schemes, oftentimes the implicit portion of an IMEX scheme. The convection–diffusion-reaction (CDR) equations allow a term by term stiff/nonstiff partition that is often well suited for IMEX methods. Major variables in CDR converge at near design-order rates with all formulations, including the fourth-order IMEX additive Runge–Kutta (ARK₂) schemes that are susceptible to order reduction. The semi-implicit backward differentiation formulae and IMEX ARK₂ schemes are of comparable efficiency. Laminar and turbulent aerodynamic applications require fully implicit schemes, as they are not profitably partitioned. All schemes achieve design-order convergence rates on the laminar problem. The fourth-order explicit singly diagonally implicit Runge–Kutta (ESDIRK4) scheme is more efficient than the popular second-order backward differentiation formulae (BDF2) method. The BDF2 and fourth-order modified extended backward differentiation formulae (MEBDF4) schemes are of comparable efficiency on the turbulent problem. High precision requirements slightly favor the MEBDF4 scheme (greater than three significant digits). Significant order reduction plagues the ESDIRK4 scheme in the turbulent case. The magnitude of the order reduction varies with Reynolds number. Poor performance of the high-order methods can partially be attributed to poor solver performance. Huge time steps allowed by high-order formulations challenge the capabilities of algebraic solver technology.     

Exponential integrator methods for systems of non-linear space-fractional models with super-diffusion processes in pattern formation

Nonlocality and spatial heterogeneity of many practical systems have made fractional differential equations very useful tools in Science and Engineering. However, solving these type of models is computationally demanding. In this paper, we propose an exponential integrator method for space fractional models as an attractive and easy-to-code alternative for other existing second-order exponential integrator methods. This scheme is based on using a real distinct poles discretization for the underlying matrix exponentials. One of the major benefits of the proposed scheme is that the algorithm could be easily implemented in parallel to take advantage of multiple processors for increased computational efficiency. The scheme is established to be second-order convergent; and proven to be robust for nonlinear space fractional reaction–diffusion problems involving non-smooth initial data. Our approach is exhibited by solving a system of two-dimensional problems which exhibits pattern formation and has applications in cell-division. Empirically, super-diffusion processes are displayed by investigating the effect of the fractional power of the underlying Laplacian operator on the pattern formation found in these models. Furthermore, the superiority of our method over competing second order ETD schemes, BDF2 scheme, and IMEX schemes is demonstrated.     

An unconditional stable compact fourth-order finite difference scheme for three dimensional Allen–Cahn equation

In this paper, we present an unconditional stable linear high-order finite difference scheme for three dimensional Allen–Cahn equation. This scheme, which is based on a backward differentiation scheme combined with a fourth-order compact finite difference formula, is second order accurate in time and fourth order accurate in space. A linearly stabilized splitting scheme is used to remove the restriction of time step. We prove the unconditional stability of our proposed method in analysis. A fast and efficient linear multigrid solver is employed to solve the resulting discrete system. We perform various numerical experiments to confirm the high-order accuracy, unconditional stability and efficiency of our proposed method. In particular, we show two applications of our proposed method: triply-periodic minimal surface and volume inpainting.     

三维非定常不可压涡量—速度Navier—Stokes方程组的有限差分法

提出了一种数值求解三维非定常涡量—速度形式的不可压Navier-Stokes方程组的有限差分方法,该方法在空间方向上具有二阶精度,并且系数矩阵具有对角占优性,因此适合高雷诺数问题的数值求解．同时,给出了适合的二阶涡量边界条件．通过对有精确解的狄利克雷边值问题和典型的驱动方腔流问题的数值实验,验证了本文格式的精确性、稳定性和有效性．