基于二重网格的动网格松弛法改进

相比传统的弹簧法等方法,基于球松弛算法的动网格松弛法在复杂边界大变形条件下可以得到质量更高的边界网格以及更大的极限变形量,但该方法在时间效率上还有提升的空间.引入二重网格,采用动网格松弛法进行稀疏网格的网格变形,将边界位移传递到整个网格计算域;再利用二重网格映射,将稀疏网格位移映射到原有计算网格的节点上.算例表明,改进...     

A first-order system approach for diffusion equation. II: Unification of advection and diffusion

In this paper, we unify advection and diffusion into a single hyperbolic system by extending the first-order system approach introduced for the diffusion equation [J. Comput. Phys., 227 (2007) 315–352] to the advection–diffusion equation. Specifically, we construct a unified hyperbolic advection–diffusion system by expressing the diffusion term as a first-order hyperbolic system and simply adding the advection term to it. Naturally then, we develop upwind schemes for this entire   system; there is thus no need to develop two different schemes, i.e., advection and diffusion schemes. We show that numerical schemes constructed in this way can be automatically uniformly accurate, allow O(h)$O(h)$ time step, and compute the solution gradients (viscous stresses/heat fluxes for the Navier–Stokes equations) simultaneously to the same order of accuracy as the main variable, for all Reynolds numbers. We present numerical results for boundary-layer type problems on non-uniform grids in one dimension and irregular triangular grids in two dimensions to demonstrate various remarkable advantages of the proposed approach. In particular, we show that the schemes solving the first-order advection–diffusion system give a tremendous speed-up in CPU time over traditional scalar schemes despite the additional cost of carrying extra variables and solving equations for them. We conclude the paper with discussions on further developments to come.     

格林公式的教学探究

首先通过一个热点问题和两个计算题引出学习格林公式的目的和必要性,迅速抓住学生的兴趣点,从第二型曲线积分入手,通过逐步分析操作给出前提条件,并最终得到公式,打破了书上的传统推导模式.     

Moving least squares reconstruction for sharp interface immersed boundary methods

We propose a new approach for reconstructing velocity boundary conditions in sharp-inerface immersed boundary (IB) methods based on the moving least squares (MLS) interpolation method. The MLS is employed to not only reconstruct velocity boundary conditions but also to calculate the pressure and velocity gradients in the vicinity of the immersed body, which are required in fluid structure interaction problems to obtain the force exerted by the fluid on the structure. To extend the method to arbitrarily complex geometries with nonconvex shaped boundaries, the visibility method is combined with the MLS method. The performance of the proposed curvilinear IB MLS (CURVIB-MLS) is demonstrated by systematic grid-refinement studies for two- and three-dimensional tests and compared with the standard CURVIB method employing standard wall-normal interpolation for reconstructing boundary conditions. The test problems are flow in a lid-driven cavity with a sphere, uniform flow over a sphere, flow on a NACA0018 airfoil at incidence, and vortex-induced vibration of an elastically-mounted cylinder. We show that the CURVIB-MLS formulation yields a method that is easier to implement in complex geometries and exhibits higher accuracy and rate of convergence relative to the standard CURVIB method. The MLS approach is also shown to dramatically improve the accuracy of calculating the pressure and viscous forces imparted by the flow on the body and improve the overall accuracy of FSI simulations. Finally, the CURVIB-MLS approach is able to qualitatively capture on relatively coarse grids important features of complex separated flows that the standard CURVIB method is able to capture only on finer grids.     

A study of properties of adaptive quadratic grids for three-dimensional near-surface field computations

Curved geometries and the corresponding near-surface fields typically require a large number of linear computational elements. High-order numerical solvers have been primarily used with low-order meshes. There is a need for curved, high-order computational elements. Typical near-surface meshes consist of hexahedral and/or prismatic elements. The present work studies the employment of quadratic meshes that are relatively coarse for field simulations. Directionally quadratic high-order elements are proposed for the near-surface field regions. The quadratic meshes are compared with the conventional low-order ones in terms of accuracy and efficiency. The cases considered include closed surface volume calculations, as well as computation of gradients of several analytic fields. A special method of adaptive local quadratic meshes is proposed and evaluated. Truncation error analysis for quadratic grids yields comparison with the conventional linear hexahedral/prismatic meshes, which are subject to typical distortions such as stretching, skewness, and torsion.     

High-order curvilinear mesh generation technique based on an improved radius basic function approach

A high-order curvilinear hybrid mesh generation technique is developed for high-order numerical method (eg, discontinuous Galerkin method) applications to improve the accuracy for problems with curve boundary. The grid generation technique is based on an improved radius basic function (RBF) approach by which the straight-edge mesh is converted into high-order curve mesh. Firstly, an initial straight-edge mesh is prepared by traditional grid generation software. Then, high-order interpolation points are inserted into the mesh entities such as edges, faces, and cells according to the final demand of mesh order. To preserve the original geometry, the inserted points on solid wall are then projected onto the CAD model using an open source tool “Open Cascade.” Finally, other inserted points in the field near the solid wall are moved to appropriate positions by the improved RBF approach to avoid tangled cells. If we use the original RBF approach, then the inserted points on the edge and face entities normal to the solid boundary in the region of boundary layer will move to improper positions. To overcome this problem, a weighting based on the local grid aspect ratio between normal direction and tangential direction is introduced into the baseline RBF approach. Three typical configurations are tested to validate the mesh generator. Meanwhile, a third-order solution of subsonic flow over an analytical 3D body of revolution in the second International Workshop on High-Order CFD Methods is supplied by a discontinuous Galerkin solver. These numerical tests demonstrate the potential capability of present technique for high-order simulations of complex geometries.     

Numerical Experiences with New Truncated Newton Methods in Large Scale Unconstrained Optimization

Recently, in [12] a very general class oftruncated Newton methods has been proposed for solving large scale unconstrained optimization problems. In this work we present the results of an extensive numericalexperience obtained by different algorithms which belong to the preceding class. This numerical study, besides investigating which arethe best algorithmic choices of the proposed approach, clarifies some significant points which underlies every truncated Newton based algorithm.     

Numerical convergence on adaptive grids for control volume methods

An adaptive technique for control‐volume methods applied to second order elliptic equations in two dimensions is presented. The discretization method applies to initially Cartesian grids aligned with the principal directions of the conductivity tensor. The convergence behavior of this method is investigated numerically. For solutions with low Sobolev regularity, the found L² convergence order is two for the potential and one for the flow density. The system of linear equations is better conditioned for the adaptive grids than for uniform grids. The test runs indicate that a pure flux‐based refinement criterion is preferable.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

在非交错网格上求解非正交曲线坐标系下N-S方程

介绍一种在非交错网格上,求解非正交曲线坐标系下的流动控制方程的方法。为解决在三维非正交曲线坐标系下采用非交错网格时所遇到的压力波动问题,将一种网格界面速度插值方法进行推广,并导出了相应的计算公式。运用所建立的方法对90°方截面弯管流动问题进行了数值模拟。结果是令人满意的,从而表明计算方法及所开发程序的可靠性。