弹性力学应力状态理论求解器开发及应用研究

论文基于吴家龙教授编著的《弹性力学》教材内容和课程"三难"特点,以应力状态理论部分常见题型求解任意微分斜面应力矢量、正应力、切应力、主应力和应力主方向及最大切应力为例,应用Matlab-GUI模块编制应力状态理论求解器,通过界面物理参数输入和结果显示过程,再现弹性力学求解思维,验证学生数理求解结论;把繁琐抽象的力学公式...     

Solution of incompressible flows with or without a free surface using the finite volume method on unstructured triangular meshes

An incompressible Navier–Stokes solver based on a cell‐centre finite volume formulation for unstructured triangular meshes is developed and tested. The solution methodology makes use of pseudocompressibility, whereby the convective terms are computed using a Godunov‐type second‐order upwind finite volume formulation. The evolution of the solution in time is obtained by subiterating the equations in pseudotime for each physical time step, with the pseudotime step set equal to infinity. For flows with a free surface the computational mesh is fitted to the free surface boundary at each time step, with the free surface elevation satisfying a kinematic boundary condition. A ‘leakage coefficient’, ε, is introduced for the calculation of flows with a free surface in order to control the leakage of flow through the free surface. This allows the assumption of stationarity of mesh points to be made during the course of pseudotime iteration. The solver is tested by comparing the output with a wide range of documented published results, both for flows with and without a free surface. The presented results show that the solver is robust. © 1999 John Wiley & Sons, Ltd.     

A CELL-CENTERED ALE METHOD WITH HLLC-2D RIEMANN SOLVER IN 2D CYLINDRICAL GEOMETRY

This paper presents a second-order direct arbitrary Lagrangian Eulerian (ALE) method for compressible flow in two-dimensional cylindrical geometry.This algorithm has half-face fluxes and a nodal velocity solver,which can ensure the compatibility between edge fluxes and the nodal flow intrinsically.In two-dimensional cylindrical geometry,the control vol-ume scheme and the area-weighted scheme are used respectively,which are distinguished by the discretizations for the source term in the momentum equation.The two-dimensional second-order extensions of these schemes are constructed by employing the monotone up-wind scheme of conservation law (MUSCL) on unstructured meshes.Numerical results are provided to assess the robustness and accuracy of these new schemes.     

NEWTON–GMRES ALGORITHM APPLIED TO COMPRESSIBLE FLOWS

This paper addresses the resolution of non-linear problems arising from an implicit time discretization in CFD problems. We study the convergence of the Newton–GMRES algorithm with a Jacobian approximated by a finite difference scheme and with restarting in GMRES. In our numerical experiments we observe, as predicted by the theory, the impact of the matrix-free approximations. A second-order scheme clearly improves the convergence in the Newton process.     

A NON–LINEAR ADAPTED TRI–TREE MULTIGRID SOLVER FOR FINITE ELEMENT FORMULATION OF THE NAVIER–STOKES EQUATIONS

An iterative adaptive equation multigrid solver for solving the implicit Navier–Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structur e which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element multigrid the Navier–Stokes equations are solved approximately. The solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. The iterative multigrid equation solver is preconditioned by incomplete LU factorization with coupled node fill-in. The non-linear Navier–Stokes equations are linearized by both the Newton method and grid adaption. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in.     

FAST ITERATIVE SOLVERS FOR THE DISCRETIZED INCOMPRESSIBLE NAVIER–STOKES EQUATIONS

In this paper some iterative solution methods of the GMRES type for the discretized Navier–Stokes equations are treated. The discretization combined with a pressure correction scheme leads to two different types of systems of linear equations: the momentum system and the pressure system. These systems may be coupled to one or more transport equations. For every system we specify a particular ILU-type preconditioner and show how to vectorize these preconditions. Finally, some numerical experiments to show the efficiency of the proposed methods are presented.     

Fast solution of three-dimensional elliptic equations with randomly generated jumping coefficients by using tensor-structured preconditioners

In this paper, we propose and analyze the numerical algorithms for fast solution of periodic elliptic problems in random media in ${ℝ}^d$, $d=2,3$. Both the two-dimensional (2D) and three-dimensional (3D) elliptic problems are considered for the jumping equation coefficients built as a checkerboard type configuration of bumps randomly distributed on a large $L\times L$, or $L\times L\times L$ lattice, respectively. The finite element method discretization procedure on a 3D $n\times n\times n$ uniform tensor grid is described in detail, and the Kronecker tensor product approach is proposed for fast generation of the stiffness matrix. We introduce tensor techniques for the construction of the low Kronecker rank spectrally equivalent preconditioner in a periodic setting to be used in the framework of the preconditioned conjugate gradient iteration. The discrete 3D periodic Laplacian pseudo-inverse is first diagonalized in the Fourier basis, and then the diagonal matrix is reshaped into a fully populated third-order tensor of size $n\times n\times n$. The latter is approximated by a low-rank canonical tensor by using the multigrid Tucker-to-canonical tensor transform. As an example, we apply the presented solver in numerical analysis of stochastic homogenization method where the 3D elliptic equation should be solved many hundred times, and where for every random sampling of the equation coefficient one has to construct the new stiffness matrix and the right-hand side. The computational characteristics of the presented solver in terms of a lattice parameter $L$ and the grid-size, $n^d$, in both 2D and 3D cases are illustrated in numerical tests. Our solver can be used in various applications where the elliptic problem should be solved for a number of different coefficients for example, in many-particle dynamics, protein docking problems or stochastic modeling.     

A conforming auxiliary space preconditioner for the mass conserving stress-yielding method

We are studying the efficient solution of the system of linear equations stemming from the mass conserving stress-yielding (MCS) discretization of the Stokes equations. We perform static condensation to arrive at a system for the pressure and velocity unknowns. An auxiliary space preconditioner for the positive definite velocity block makes use of efficient and scalable solvers for conforming Finite Element spaces of low order and is analyzed with emphasis placed on robustness in the polynomial degree of the discretization. Numerical experiments demonstrate the potential of this approach and the efficiency of the implementation.     

Tensor rank bounds and explicit QTT representations for the inverses of circulant matrices

In this paper, we are concerned with the inversion of circulant matrices and their quantized tensor-train (QTT) structure. In particular, we show that the inverse of a complex circulant matrix $A$, generated by the first column of the form ${(a_0,\dots ,a_{m-1},0,\dots ,0,a_{-n},\dots ,a_{-1})}^{\top }$ admits a QTT representation with the QTT ranks bounded by $(m+n)$. Under certain assumptions on the entries of $A$, we also derive an explicit QTT representation of $A^{-1}$. The latter can be used, for instance, to overcome stability issues arising when numerically solving differential equations with periodic boundary conditions in the QTT format.     

基于WENO重构的真正多维黎曼求解器

具有良好守恒性与网格适应性的有限体积格式在流体力学的数值计算中占有重要地位。其中，求解数值流通量是实施有限体积法的关键步骤。一维情形下，通过求解局部黎曼问题来获得数值流通量的相关理论已经比较成熟。但是在计算多维问题时，传统的维度分裂方法仅考虑沿界面法向传播的信息，这不仅影响格式的精度，还可能会造成数值不稳定性从而诱发非物理现象。本文基于对流-压力通量分裂方法来构造真正多维的黎曼求解器，通过求解网格顶点处的多维黎曼问题来实现格式的多维特性。采用五阶WENO重构方法来获得空间的高阶精度，时间离散采用三阶TVD龙格-库塔格式。一系列数值实验的结果表明，真正多维的黎曼求解器不仅具有更高的分辨率还能有效克服多维强激波模拟中的数值不稳定性。