Uniform Convergent Tailored Finite Point Method for Advection–Diffusion Equation with Discontinuous,Anisotropic and Vanishing Diffusivity

We propose two tailored finite point methods for the advection–diffusion equation with anisotropic tensor diffusivity. The diffusion coefficient can be very small in one direction in some part of the domain and be discontinuous at the interfaces. When flows advect from the vanishing-diffusivity region towards the non-vanishing diffusivity region, standard numerical schemes tend to cause spurious oscillations or negative values. Our proposed schemes have uniform convergence in the vanishing diffusivity limit, even when the solution exhibits interface and boundary layers. When the diffusivity is along the coordinates, the positivity and maximum principle can be proved. We use the value as well as their derivatives at the grid points to construct the scheme for nonaligned case, which makes it can achieve good accuracy and convergence for the derivatives as well, even when there exhibit boundary or interface layers. Numerical experiments are presented to show the performance of the proposed scheme.     

求解二维对流扩散方程的格子Boltzmann方法

针对二维对流扩散方程,基于D2Q4格子速度,用Chapman-Enskog多尺度分析技术,将时间尺度取为二阶,空间尺度取为一阶,推导了各个速度方向上的平衡态分布函数所满足的条件,给出了简单且对称的平衡态分布函数表达式,所得到的平衡态分布函数能正确地恢复出二维对流扩散方程,从而构建了一种新的求解二维对流扩散方程的D2Q4格子Boltzmann（LB）模型。用所给LB模型对扩散方程和两个不同初边界条件的对流扩散方程进行了数值求解,数值实验结果表明数值解与精确解吻合较好,与相关文献结果比较边界误差要小得多,验证了模型的有效性。     

基于格子Boltzmann方法的一维Burgers方程的数值模拟

基于格子Boltzmann方法(LBM)的一维Burgers方程的数值解法,已有2-bit和4-bit模型。文中通过选择合适的离散速度模型构造出恰当的平衡态分布函数; 然后, 利用单松弛的格子Bhatnagar-Gross-Krook模型、Chapman-Enskog展开和多尺度技术, 提出了用于求解一维Burgers方程的3-bit的格子Boltzmann模型,即D1Q3模型,并进行了数值实验。实验结果表明,该方法的数值解与解析解吻合的程度很好,且误差比现有文献中的误差更小,从而验证了格子Boltzamnn模型的有效性。     

多重网格格子Boltzmann方法的并行算法

针对复杂流动数值模拟中的格子Boltzmann方法存在计算网格量大、收敛速度慢的缺点,提出了基于三维几何边界的多重笛卡儿网格并行生成算法,并基于该网格生成方法提出了多重网格并行格子Boltzmann方法（LBM）。该方法结合不同尺度网格间的耦合计算,有效减少了计算网格量,提高了收敛速度;而且测试结果也表明该并行算法具有良好的可扩展性。     

Spectral element spatial discretization error in solving highly anisotropic heat conduction equation

This paper describes a study of the effects of the overall spatial resolution, polynomial degree and computational grid directionality on the accuracy of numerical solutions of a highly anisotropic thermal diffusion equation using the spectral element spatial discretization method. The high-order spectral element macroscopic modeling code SEL/HiFi has been used to explore the parameter space. It is shown that for a given number of spatial degrees of freedom, increasing polynomial degree while reducing the number of elements results in exponential reduction of the numerical error. The alignment of the grid with the direction of anisotropy is shown to further improve the accuracy of the solution. These effects are qualitatively explained and numerically quantified in 2- and 3-dimensional calculations with straight and curved anisotropy.     

Lattice Boltzmann simulation of some nonlinear convection–diffusion equations

In this work we proposed a lattice Boltzmann model for the nonlinear convection–diffusion equation (NCDE) with anisotropic diffusion. The constraints on the model for correctly recovering macroscopic equation are also carefully analyzed, which are ignored in some existing work. Detailed simulations of some 1D/2D NCDEs, including the nonlinear Schrödinger equation (NLSE), Buckley–Leverett equation with discontinuous initial data, NCDE with anisotropic diffusion, and generalized Zakharov system, are performed. The numerical results obtained by the proposed model agree well with the analytical solutions and/or the numerical solutions reported in previous studies. It is also found that, for complex-valued NLSE, the model using a complex distribution function is superior to that using two real distribution functions for the real and imaginary parts of the NLSE separately.     

Accelerated Lattice Boltzmann Schemes for Steady-State Flow Simulations

A matrix formulation of the steady Lattice Boltzmann equation is presented. It is shown that the strict steady-state formulation, combined with preconditioned iterative solvers, leads to significant computational savings as compared to the standard explicit LBE scheme.     

Domain decomposition based coupling between the lattice Boltzmann method and traditional CFD methods—Part I: Formulation and application to the 2-D Burgers’ equation

The lattice Boltzmann method is being increasingly employed in the field of computational fluid dynamics due to its computational efficiency. Floating-point operations in the lattice Boltzmann method involve local data and therefore allow easy cache optimization and parallelization. Due to this, the cache-optimized lattice Boltzmann method has superior computational performance over traditional finite difference methods for solving unsteady flow problems. When solving steady flow problems, the explicit nature of the lattice Boltzmann discretization limits the time step size and therefore the efficiency of the lattice Boltzmann method for steady flows. To quantify the computational performance of the lattice Boltzmann method for steady flows, a comparison study between the lattice Boltzmann method (LBM) and the alternating direction implicit (ADI) method was performed using the 2-D steady Burgers’ equation. The comparison study showed that the LBM performs comparatively poor on high-resolution meshes due to smaller time step sizes, while on coarser meshes where the time step size is similar for both methods, the cache-optimized LBM performance is superior. Because flow domains can be discretized with multiblock grids consisting of coarse and fine grid blocks, the cache-optimized LBM can be applied on the coarse grid block while the traditional implicit methods are applied on the fine grid blocks. This paper finds the coupled cache-optimized lattice Boltzmann-ADI method to be faster by a factor of 4.5 over the traditional methods while maintaining similar accuracy.     

Nitsche Type Mortaring for Singularly Perturbed Reaction-diffusion Problems

The paper is concerned with the Nitsche mortaring in the framework of domain decomposition where non-matching meshes and weak continuity of the finite element approximation at the interface are admitted. The approach is applied to singularly perturbed reaction-diffusion problems in 2D. Non-matching meshes of triangles being anisotropic in the boundary layers are applied. Some properties as well as error estimates of the Nitsche mortar finite element schemes are proved. In particular, using a suitable degree of anisotropy of triangles in the boundary layers of a rectangle, convergence rates as known for the conforming finite element method are derived. Numerical examples illustrate the approach and the results.     

An adaptive scheme using hierarchical grids for lattice Boltzmann multi-phase flow simulations

The lattice Boltzmann (LB) method is extended and adapted to simulate multi-phase flows on non-uniform tree-type grids. Our model is an extension of the model developed by Gunstensen [Gunstensen AK, Rothman D. Lattice Boltzmann model of immiscible fluids. Phys Rev A 1991;43(8):4320–4327], which is based on the Rothman–Keller model [Rothman DH, Keller JM. Immiscible cellular automaton fluids. J Stat Phys 1988;52:1119–1127]. A first approach we use an a priori grid refinement. We find that the maximum number of possible grid levels for problems with dominant capillary forces is very restricted, if the physical interface is allowed to pass over grid interfaces. Thus a second approach based on adaptive grids was developed, where the physical interface is always discretized on the finest grid level. Efficient and flexible data structures have been developed to manage the remeshing. The application of the scheme for a rising bubble in three dimensions shows very good agreement with the semi-analytical solution and demonstrates the efficiency of our approach.     

Continuous elliptic and multi-dimensional hyperbolic Darcy-flux finite-volume methods

Elliptic and hyperbolic Darcy-flux approximations are presented. Families of flux-continuous finite-volume methods are investigated for the elliptic full-tensor pressure equation with general discontinuous coefficients. Full pressure continuity across control-volume interfaces is built into the methods leading to an important distinction from the earlier pointwise continuous methods. The families of quasi-positive methods significantly reduce spurious oscillations (induced by earlier schemes) in discrete pressure solutions for strongly anisotropic full-tensor fields. Anisotropy favoring triangulation and non-linear flux splitting are also shown to be effective for computing solutions free of spurious oscillations.Multi-dimensional upwind schemes that reduce cross-wind numerical diffusion induced by the standard upwind scheme are also presented for hyperbolic Darcy-flux approximation.     

Recursive algorithm and accurate computation of dyadic Green’s functions for stratified uniaxial anisotropic media

A recursive algorithm is adopted for the computation of dyadic Green＇s functions in three-dimensional stratified uniaxial anisotropic media with arbitrary number of layers. Three linear equation groups for computing the coefficients of the Sommerfeld integrals are obtained according to the continuity condition of electric and magnetic fields across the interface between different layers, which are in correspondence with the TM wave produced by a vertical unit electric dipole and the TE or TM wave produced by a horizontal unit electric dipole, respectively. All the linear equation groups can be solved via the recursive algorithm. The dyadic Green＇s functions with source point and field point being in any layer can be conveniently obtained by merely changing the position of the elements within the source term of the linear equation groups. The problem of singularities occurring in the Sommerfeld integrals is efficiently solved by deforming the integration path in the complex plane. The expression of the dyadic Green＇s functions provided by this paper is terse in form and is easy to be programmed, and it does not overflow. Theoretical analysis and numerical examples show the accuracy and effectivity of the algorithm.     

Lattice Boltzmann Models for Anisotropic Diffusion of Images

The lattice Boltzmann method has attracted more and more attention as an alternative numerical scheme to traditional numerical methods for solving partial differential equations and modeling physical systems. The idea of the lattice Boltzmann method is to construct a simplified discrete microscopic dynamics to simulate the macroscopic model described by the partial differential equations. The use of the lattice Boltzmann method has allowed the study of a broad class of systems that would have been difficult by other means. The advantage of the lattice Boltzmann method is that it provides easily implemented fully parallel algorithms and the capability of handling complicated boundaries. In this paper, we present two lattice Boltzmann models for nonlinear anisotropic diffusion of images. We show that image feature selective diffusion (smoothing) can be achieved by making the relaxation parameter in the lattice Boltzmann equation be image feature and direction dependent. The models naturally lead to the numerical algorithms that are easy to implement. Experimental results on both synthetic and real images are described.     

基于窄带的自适应Level Set方法

Level Set模型将运动界面表示为高维场函数的零等值面,自然而鲁棒地解决了界面演化中拓扑结构改变的问题,但计算效率不高.文中提出了基于窄带的自适应Level Set方法.自适应方法首先构建粗网格满足界面演化的整体需求,同时估算粗网格点的曲率值,使用快速扩散法聚类高曲率点,通过主元分析估算点集朝向,构建细网格捕捉演化中的细节区域.粗、细网格均为独立的计算单元,定义为存储网格中的有向包围盒.这种数据结构可以有效避免频繁的坐标变换和插值操作,同时保证了数值解的精度.实验结果与误差分析表明,自适应方法能有效减少计算量,达到更好的界面跟踪效果.     

Diffusion network architectures for implementation of Gibbssamplers with applications to assignment problems

In this paper, analog circuit designs for implementations of Gibbs samplers are presented, which offer fully parallel computation. The Gibbs sampler for a discrete solution space (or Boltzmann machine) can be used to solve both deterministic and probabilistic assignment (association) problems. The primary drawback to the use of a Boltzmann machine for optimization is its computational complexity, since updating of the neurons is typically performed sequentially. We first consider the diffusion equation emulation of a Boltzmann machine introduced by Roysam and Miller (1991), which employs a parallel network of nonlinear amplifiers. It is shown that an analog circuit implementation of the diffusion equation requires a complex neural structure incorporating matched nonlinear feedback amplifiers and current multipliers. We introduce a simpler implementation of the Boltzmann machine, using a "constant gradient" diffusion equation, which eliminates the need for a matched feedback amplifier. The performance of the Roysam and Miller network and the new constant gradient (CG) network is compared using simulations for the multiple-neuron case, and integration of the Chapman-Kolmogorov equation for a single neuron. Based on the simulation results, heuristic criteria for establishing the diffusion equation boundaries, and neuron sigmoidal gain are obtained. The final CG analog circuit is suitable for VLSI implementation, and hence may offer rapid convergence.     

A direct-forcing immersed boundary method for the thermal lattice Boltzmann method

In this study, a direct-forcing immersed boundary method (IBM) for thermal lattice Boltzmann method (TLBM) is proposed to simulate the non-isothermal flows. The direct-forcing IBM formulas for thermal equations are derived based on two TLBM models: a double-population model with a simplified thermal lattice Boltzmann equation (Model 1) and a hybrid model with an advection–diffusion equation of temperature (Model 2). As an interface scheme, which is required due to a mismatch between boundary and computational grids in the IBM, the sharp interface scheme based on second-order bilinear and linear interpolations (instead of the diffuse interface scheme, which uses discrete delta functions) is adopted to obtain the more accurate results. The proposed methods are validated through convective heat transfer problems with not only stationary but also moving boundaries – the natural convection in a square cavity with an eccentrically located cylinder and a cold particle sedimentation in an infinite channel. In terms of accuracy, the results from the IBM based on both models are comparable and show a good agreement with those from other numerical methods. In contrast, the IBM based on Model 2 is more numerically efficient than the IBM based on Model 1.     

An Adaptive Staggered Discontinuous Galerkin Method for the Steady State Convection–Diffusion Equation

Staggered grid techniques have been applied successfully to many problems. A distinctive advantage is that physical laws arising from the corresponding partial differential equations are automatically preserved. Recently, a staggered discontinuous Galerkin (SDG) method was developed for the convection–diffusion equation. In this paper, we are interested in solving the steady state convection–diffusion equation with a small diffusion coefficient \(\epsilon \). It is known that the exact solution may have large gradient in some regions and thus a very fine mesh is needed. For convection dominated problems, that is, when \(\epsilon \) is small, exact solutions may contain sharp layers. In these cases, adaptive mesh refinement is crucial in order to reduce the computational cost. In this paper, a new SDG method is proposed and the proof of its stability is provided. In order to construct an adaptive mesh refinement strategy for this new SDG method, we derive an a-posteriori error estimator and prove its efficiency and reliability under a boundedness assumption on \(h/\epsilon \), where h is the mesh size. Moreover, we will present some numerical results with singularities and sharp layers to show the good performance of the proposed error estimator as well as the adaptive mesh refinement strategy.     

High Order Accurate Finite Difference Modeling of Seismo-Acoustic Wave Propagation in a Moving Atmosphere and a Heterogeneous Earth Model Coupled Across a Realistic Topography

We develop a numerical method for simultaneously simulating acoustic waves in a realistic moving atmosphere and seismic waves in a heterogeneous earth model, where the motions are coupled across a realistic topography. We model acoustic wave propagation by solving the linearized Euler equations of compressible fluid mechanics. The seismic waves are modeled by the elastic wave equation in a heterogeneous anisotropic material. The motion is coupled by imposing continuity of normal velocity and normal stresses across the topographic interface. Realistic topography is resolved on a curvilinear grid that follows the interface. The governing equations are discretized using high order accurate finite difference methods that satisfy the principle of summation by parts. We apply the energy method to derive the discrete interface conditions and to show that the coupled discretization is stable. The implementation is verified by numerical experiments, and we demonstrate a simulation of coupled wave propagation in a windy atmosphere and a realistic earth model with non-planar topography.