基于格子-波尔兹曼法的流体能耗求解

格子-波尔兹曼法是近年来新兴的一种计算流体力学数值方法。随着这种方法的不断发展,人们将它用于流体的仿真、优化等不同场合。与此同时,一些与流场流速和压强相关的物理量(如能耗)的求解也成为关注的焦点。本文介绍了能耗这一流体宏观量的格子-波尔兹曼法求解及其实现。与传统的有限差分法不同,本文在求解有关的速度梯度时使用了格子-波尔兹曼-矩法,这种方法不但能够避免有限差分法在边界处失效的缺点,而且计算简单,算法局部性好,适合大规模并行计算。本文在分析其数值解精度的基础上,使用这种方法进行了以能耗极小为目标的直通道内椭圆挡块的参数优化。这些分析和算例分别定量和定性地说明了本文算法的准确性。     

The designing approach of difference schemes by controlling the remainder‐effect

In this paper, based on the idea of the ‘modified partial differential equation’, a new designing approach to explicit finite difference schemes for the Burgers equation and PDE is proposed. The approach differs from other constructured methods in such a way that it considers the requests of the numerical dissipation and dispersion coefficients first. This method is much more constructional and directional. The results of numerical tests indicate that the method is quite successful. Copyright © 1999 John Wiley & Sons, Ltd.     

Nonlinear Multigrid Methods for Numerical Solution of the Unsaturated Flow Equation in Two Space Dimensions

Picard and Newton iterations are widely used to solve numerically the nonlinear Richards’ equation (RE) governing water flow in unsaturated porous media. When solving RE in two space dimensions, direct methods applied to the linearized problem in the Newton/Picard iterations are inefficient. The numerical solving of RE in 2D with a nonlinear multigrid (MG) method that avoids Picard/Newton iterations is the focus of this work. The numerical approach is based on an implicit, second-order accurate time discretization combined with a second-order accurate finite difference spatial discretization. The test problems simulate infiltration of water in 2D unsaturated soils with hydraulic properties described by Broadbridge–White and van Genuchten–Mualem models. The numerical results show that nonlinear MG deserves to be taken into consideration for numerical solving of RE.     

Alternating segment explicit-implicit scheme for nonlinear third-order KdV equation

A group of asymmetric difference schemes to approach the Korteweg-de Vries(KdV)equation is given here.According to such schemes,the full explicit difference scheme and the fun implicit one,an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed.The scheme is linear unconditionally stable by the analysis of linearization procedure,and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.     

Direct Integration Methods with Integral Model for Dynamic Systems

ＩｎｔｒｏｄｕｃｔｉｏｎＴｒａｎｓｉｅｎｔｓｔａｔｅａｎａｌｙｓｉｓｈａｓｂｅｅｎａｎａｃｔｉｖｅｒｅｓｅａｒｃｈａｒｅａｉｎｍａｎｙｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ.Ｉｎｐｒａｃｔｉｃａｌｓｉｔｕａｔｉｏｎ ,ｌｉｋｅｓｔｒｕｃｔｕｒａｌｍｅｃｈａｎｉｃｓ,ｔｈｅｓｙｓｔｅｍｓｂｅｉｎｇｓｔｕｄｉｅｄａｒｅｕｓｕａｌｌｙｎｏｎｌｉｎｅａｒａｎｄｔｉｍｅ＿ｄｅｐｅｎｄｅｎｔ.Ｔｈｅａｎａｌｙｔｉｃａｌｍｅｔｈｏｄｓｔｏｓｏｌｖｅｔｈｅｓｅｐｒｏｂｌｅｍｓａｒｅｉｍｐｏｓｓｉｂｌｅ .Ｔｈｅｎｕｍ…     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 John Wiley & Sons, Ltd.     

A problem for diffusion with relaxation in polymers

A nonlinear problem for penetrant diffusion with relaxation in polymers is considered. A numerical approach to solving this type of problems is developed. The proposed numerical scheme based on a finite element domain approximation and a time difference method can be used for numerical simulation of the considered penetrant diffusion in 2-D and 3-D domains. A numerical procedure and a corresponding computer code are created and tested for some examples in 1-D and 2-D domains.     

各种网格上统一的数值离散方法

提出一种在任意网格上计算数值微分的方法，这种方法利用各种不同网格所具有的共同性质， 基于Taylor展开和加权最小二乘法，得到了各种网格下都可以使用的数值微分格式. 有了这一技术， 可以极大地丰富已经发展起来的各种数值方法，原来只能用在结构网格上的格式，可以直接推广到 其他各种网格上，从而可以用于各种复杂区域内微分方程的数值求解. 初步的应用表明这种技术是 简单而有效的.     

Analysis of thermal responses in a two-dimensional porous medium caused by pulse heat flux

In this article, the generalized model for thermoelastic waves with two relaxation times is utilized to compute the increment of temperature, the displacement components, the stress components, and the changes in the volume fraction field in a two-dimensional porous medium. By using the Fourier-Laplace transform and the eigenvalue method, the considered variables are obtained analytically. The derived approach is estimated with numerical outcomes which are applied to the porous media with a geometrical simplification. The numerical results for the considered variables are performed and presented graphically. Finally, the outcomes are represented graphically to display the difference among the classical dynamical(CD) coupled, the Lord-Shulman(LS), and the Green-Lindsay(GL) models.     

双相各向异性介质弹性波场有限差分正演模拟

从双相各向异性介质模型出发,以Boit理论为基础,推导了斜方晶系各向异性介质-阶弹性波动方程,引入固、流体密度比和孔隙几何参数,将Biot方程系数简化为测量简单、物理意义明确的物理量,采用交错网格技术建立了各向异性孔隙介质波动方程的高精度差分格式,并首次对这类差分格式的频散特性和稳定性作了详细分析讨论,解决了计算稳定性和边界反射问题,与解析解的对比以及理论模型的数值模拟都表明,该方法不仅大大降低了计算量,提高了正演速度,并且具有良好的稳定性和精确性。     

一种基于Associated Hermite正交函数求解对流扩散方程的算法

提出了一种基于AH(Associated Hermite)正交基函数求解对流扩散方程的无条件稳定算法。该算法将方程的时间项通过Hermite多项式作为正交基函数进行展开,利用Galerkin方法消除时间变量项,从而导出有限维AH域隐式差分方程,突破了传统显式差分格式稳定性条件的限制,最后通过对AH域展开系数的求解得到该对流扩散方程的数值解。在数值算例中,将该算法与传统显示差分法和交替方向隐式差分法进行对比分析,数值计算结果表明,算法无条件稳定且其计算精度与时间步长无关,对于具有精细结构的对流换热问题,该算法具有明显的效率优势,且保持了较高的精度。     

Convergence and exact solutions of spline finite strip method using unitary transformation approach

The spline finite strip method(PSM) is one of the most popular numerical methods for analyzing prismatic structures.Efficacy and convergence of the method have been demonstrated in previous studies by comparing only numerical results with analytical results of some benchmark problems.To date,no exact solutions of the method or its explicit forms of error terms have been derived to show its convergence analytically. As such,in this paper,the mathematical exact solutions of spline finite strips in the plat...     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

Numerical investigation of two-phase rotating flow

Finite difference solutions of the two-fluid equations of motion for a particle (droplet)-fluid mixture in a rotating finite axisymmetric cylinder are presented. The numerical method, which can be regarded as an extension of the Harlow & Amsden approach, employs forward time and centred space discretization and treats implicitly the pressure, Coriolis and volume flux terms. The computed flow fields are examined via a detailed comparison to previous analytic approximations, which illuminates both the physical and numerical aspects and the validity of these approximations.     

Operator Splitting Multiscale Finite Volume Element Method for Two-Phase Flow with Capillary Pressure

A numerical method used for solving a two-phase flow problem as found in typical oil recovery is investigated in the setting of physics-based two-level operator splitting. The governing equations involve an elliptic differential equation coupled with a parabolic convection-dominated equation which poses a severe restriction for obtaining fully implicit numerical solutions. Furthermore, strong heterogeneity of the porous medium over many length scales adds to the complications for effectively solving the system. One viable approach is to split the system into three sub-systems: the elliptic, the hyperbolic, and the parabolic equation, respectively. In doing so, we allow for the use of appropriate numerical discretization for each type of equation and the careful exchange of information between them. We propose to use the multiscale finite volume element method (MsFVEM) for the elliptic and parabolic equations, and a nonoscillatory difference scheme for the hyperbolic equation. Performance of this procedure is confirmed through several numerical experiments.     

Connection machine simulation of ultrasonic wave propagation in materials. II: The two-dimensional case

A local interaction simulation approach (LISA) for the wave propagation in inhomogeneous 2D media is presented. The method is designed to take full advantage of massively parallel computing, such as provided by the Connection Machine. Crosspoints at the intersection of orthogonal interfaces separating media of different physical properties are treated in the framework of a sharp interface model. A comparison with finite difference techniques shows that the proposed method avoids the ambiguities due to the smoothing of the physical quantities, which is necessary in order to transform differential equations into finite difference equations. The smoothing procedure may cause severe numerical errors, when the variations of the physical properties across the interfaces are large.

In order to demonstrate the efficiency and reliability of the approach several examples of simulation of pulse propagation in different media are reported.     

Modeling of wrinkling of thin circular sheets

A version of the imperfection method is used to investigate the wrinkling (tension buckling) of thin, elastic, homogeneous, isotropic circular plates of uniform thickness undergoing small deflections without inplane twisting. A numerical procedure based on the finite difference approach is employed to quantitatively predict wrinkling loads and wrinkling patterns for three sets of support conditions. Representative numerical results are presented in tabular and graphical formats and used to illustrate interesting aspects of the predictions.     

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids.     

功能梯度材料结构分析的半解析数值方法研究

一种典型的半解析数值方法——线法被引入功能梯度材料的结构分析。首先推导了功能梯度材料位移形式的平衡方程和边界条件,然后阐述了线法功能梯度材料结构分析的基本步骤和数值原理。该方法的基本思想是通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条函数Gauss配点法求解该常微分方程组得到问题的解答。为演示线法在功能梯度材料结构分析中的应用,给出了线性梯度和指数梯度功能梯度材料板分别受恒定位移、均匀拉伸载荷和弯曲载荷作用的数值算例。与相应问题解析解和其他数值方法的比较表明,线法的计算结果具有很高的精度,而且不需要任何特殊的考虑就能够有效模拟材料内部物性参数的连续变化,也无需事先选取满足特定条件的待定场函数,是一种非常适合功能梯度材料结构形式和材料特点的半解析数值方法。     

A priori tests on numerical errors in large eddy simulation using finite differences and explicit filtering

When low‐order finite‐difference methods are applied in large eddy simulation (LES), the magnitude of the numerical error may be larger than that of the subgrid‐scale (SGS) term. In this paper, the effect of explicit filtering on the numerical error related to the spatial discretization of the convection term and the exact SGS term is studied a priori in the turbulent fully developed channel flow. As the filter width is increased the grid resolution is kept constant. Also filtering in the inhomogeneous wall‐normal direction is discussed. The main conclusions are related to two approaches to explicit filtering. In the traditional approach, the whole velocity field is filtered explicitly while in the alternative approach, only the non‐linear convection term of the Navier–Stokes equations is filtered explicitly. Based on the results presented in the paper it seems that the first approach leads to an unphysical situation. However, the later approach works in the desired way, and the numerical error becomes clearly smaller than the SGS term. The main difference between the two approaches seems to be the interpretation of the resolved non‐linear term in the filtered Navier–Stokes equations. Copyright © 2005 John Wiley & Sons, Ltd.