基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界,避免生成复杂贴体网格及动网格,极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法,将简化附加体积力以源项形式嵌入动量方程迭代中,通过用户自定义函数对CFD软件FLUENT二次开发,实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证,并探讨了插值函数对计算精度的影响。研究表明,通过引入浸入边界模型,能够提高计算效率,并实现结构网格背景下复杂边界和动边界的高效建模。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界，避免生成复杂贴体网格及动网格，极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法，将简化附加体积力以源项形式嵌入动量方程迭代中，通过用户自定义函数对CFD软件FLUENT二次开发，实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证，并探讨了插值函数对计算精度的影响。研究表明，通过引入浸入边界模型，能够提高计算效率，并实现结构网格背景下复杂边界和动边界的高效建模。     

大霍尔系数下电离气体与磁场相互作用规律数值研究

通过洛伦兹力与焦耳热耦合流场N-S方程和电势场泊松方程, 实现对低磁雷诺数磁流体流场的数值模拟, 并考虑霍尔效应和外电路对计算结果的影响. 采用无虚拟时间步的LUSGS预处理BI-CGSTAB算法解决了大霍尔系数下泊松方程病态矩阵求解问题. 对霍尔系数10$^2$的量级的高超声速磁流体圆管绕流和进气道隔离段能量沉积两种现象的数值模拟表明, 外电路、电极冷却措施、等离子体均匀度对电磁力作用效果影响很大.      

紧致插值曲线方法在流向强迫振荡圆柱绕流中的应用

利用紧致插值曲线(constrained interpolation profile method in Zhejiang University, CIP-ZJU) 数学模型, 对低科勒冈-卡朋特(Keulegan–Carpenter) 数KC 静止流体中振荡圆柱以及雷诺数Re = 200 时流向强迫振荡圆柱绕流进行了数值模拟. 模型在直角坐标系统下建立, 采用紧致插值曲线方法作为流场的基本求解器离散了纳维-斯托克斯方程, 基于多相流的理论实现流固耦合同步求解, 利用浸入边界方法处理固体边界. 模拟结果与现有文献结果进行比较, 二者吻合情况较好, 验证了此方法对于计算复杂流动问题的可靠性.      

泊松方程与双向张弦梁找形

根据双向张弦梁上、下弦的平衡荷载与索曲面及索拉力之间的关系，给出了上、下弦 所在曲面的泊松方程；给出了矩形边界和椭圆边界泊松方程的求解方法；通过算例验证了泊 松方程找形方法的正确性. 给出的泊松方程求解方法可以非常准确计算出双向张弦梁的 形状，便于在工程中应用.     

AN EFFICIENT SOLUTION FOR 2D RAYLEIGH-BÉNARD CONVECTION USING FFT

研究提高二维方腔瑞利-贝纳德对流 直接数值模拟求解方法的计算效率问题.对于非定常湍流热对流, 压力泊松方程的求解是影响整个计算效率的关键. 利用快速傅里叶变换(fast Fourier transform,FFT)解耦并结合追赶法, 可实现压力泊松方程的直接求解.通过与跳点超松弛迭代法在求解精度和计算速度对比, 可以看到, 利用FFT压力泊松方程直接方法计算热对流问题是高效的.还给出了典型状态的热对流初始羽流和大尺度环流温度场, 以及系列瑞利数(Ra)计算结果的宏观传热努塞数(Nu)变化.     

一种结构动力响应分析的快速高阶算法

为了提高基于高阶格式的结构动力响应微分求积分析方法的计算效率,发展了一种求解动力方程的快速算法．利用微分求积原理将结构动力方程转化为标准Sylvester方程的形式,通过对系数矩阵进行矩阵分解,进而将动力响应Sylvester方程化为一系列标准线性方程组,采用相关成熟算法求解这些线性方程组后即可获得结构动力时程响应的全部解答．结构动力响应微分求积分析方法为高阶数值方法,一步计算可以获得多个时点处的动力响应．基于本文快速算法,不必直接对矩阵方程进行求解．数值算例表明,本文快速算法能够准确地计算出结构动力响应,具有数值精度高、收敛性好的优点．     

高精度PPM格式在湍流燃烧中的应用

利用处理三维可压缩粘性流体流动问题中的沉浸边界法,并结合基于PPM方法的高精度TVD格式,对三维方形管道中部的圆柱火焰绕流及惰性气体绕流问题进行了数值模拟。计算湍流时采用大涡模拟(LES),化学反应速率采用EBU漩涡破碎模型。通过计算结果与实验结果的比较,发现高精度PPM格式能精确模拟两类圆柱绕流问题。计算中还发现,火焰圆柱绕流算例中,在火焰翻越圆柱前,由于燃烧的膨胀作用,使得火焰正面前的未燃气体流动并形成惰性气体绕流,这与无燃烧时的惰性气体绕流类似。但当火焰翻越圆柱过程中及完全翻越圆柱后,两种算例绕流流场出现明显变化。     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

Stress Separation of Interferometrically Measured Isopachics in a Perforated Plate

A method for the stress separation of interferometrically measured isopachics using an Airy stress function is proposed in this study. A Poisson equation that represents the relationship between the sum of principal stresses and an Airy stress function is solved using a finite element method. The Dirichlet boundary condition for solving the Poisson equation is determined by the approximation of an assumed Airy stress function along the boundary of the model. Therefore, the distribution of the Airy stress function is obtained from the measured isopachic contours. Then, the stresses are obtained from the computed Airy stress function. The effectiveness of the proposed method is validated by applying the proposed method to the isopachic contours in a perforated plate obtained by Mach-Zehnder interferometry. Results indicate that stress components around a hole in a plate can be obtained from isopachics by the proposed method.     

一种模拟动边界绕流的锐利界面浸入边界法

When the structural wall moves over a fixed grid, the structure coverage will change, resulting in many dead and emerging elements. To avoid the influence of malformation and reconstruction of body-fitted grids on the calculation efficiency and accuracy of the fluid-structure interaction problems with coupled boundary movement on the fixed grid, an improved numerical method for describing the interaction between an immersed rigid body and fluid based on a sharp-interface is proposed. In this method, both the fluid and solid are regarded as pure fluid domains in the whole computational domain, and the solid boundary is divided into several Lagrangian grid points. The flow parameter or velocity is reconstructed by interpolation at the interface element, which is then directly used as the boundary condition of the flow field, thus reflecting the influence of the wall boundary conditions. The method constructs the calculation structure of “virtual point, force point and vertical foot point”, and the velocity of the virtual point is obtained by bilinear interpolation. Then, the velocity of the force point is calculated by forcing the solid boundary to meet the no-slip condition, and the equations of the coupling system based on the immersion boundary method are finally solved to realize the numerical simulation of the flow with a complex moving boundary. The numerical program for this immersed boundary method is established using C++, then the accuracy and reliability of the proposed method are validated by comparison with the literature and experimental results of the basic numerical example of flow around a cylinder. Furthermore, the effects of the structural shape and the angle of attack on the trailing vortex structure, the vortex shedding frequency, and the lift/ coefficient characteristics of the flow around the elliptical cylinder have been analyzed. The anti-symmetric S-type, “P+S” Ⅰ-type and “P+S” Ⅱ-type trailing vortex shedding modes, as well as the variation laws of the vortex structure size, vortex shedding frequency and lift-drag coefficients ratio with axis ratio and angle of attack, are captured. The critical angle of attack (25°) corresponding to the maximum lift-drag ratio is determined as 25°.     

A study of boundary integral method for a type of nonlinear problem based on generalized quasilinearlization theory

非线性泊松问题在热传导和多孔催化粒子的扩散反应等问题中是非常常见的,为此,利用广 义拟线性化迭代理论,提出了一种非线性泊松问题的新的数值迭代方法. 该方法将非线性方 程转化成一序列线性方程的迭代,其优点是初始值的选取具有一定的理论基础,并且在一定 的初始值条件下,迭代结果将单调地收敛于非线性问题的解. 将此迭代方法与边界元和双互 易杂交边界点方法结合,并用于非线性泊松问题的求解,比较了两种方法的结果精度,收敛 速度及不同初始值下的稳定性. 结果显示,基于拟线性化的双互易杂交边界点法具有较高的 稳定性和计算效率,并且收敛速度为平方阶.     

A staggered grid‐based explicit jump immersed interface method for two‐dimensional Stokes flows

In this paper the explicit jump immersed interface method (EJIIM) is applied to stationary Stokes flows. The boundary value problem in a general, non‐grid aligned domain is reduced by the EJIIM to a sequence of problems in a rectangular domain, where staggered grid‐based finite differences for velocity and pressure variables are used. Each of these subproblems is solved by the fast Stokes solver, consisting of the pressure equation (known also as conjugate gradient Uzawa) method and a fast Fourier transform‐based Poisson solver. This results in an effective algorithm with second‐order convergence for the velocity and first order for the pressure. In contrast to the earlier versions of the EJIIM, the Dirichlét boundary value problem is solved very efficiently also in the case when the computational domain is not simply connected. Copyright © 2007 John Wiley & Sons, Ltd.     

Moving immersed boundary method

A novel implicit immersed boundary method of high accuracy and efficiency is presented for the simulation of incompressible viscous flow over complex stationary or moving solid boundaries. A boundary force is often introduced in many immersed boundary methods to mimic the presence of solid boundary, such that the overall simulation can be performed on a simple Cartesian grid. The current method inherits this idea and considers the boundary force as a Lagrange multiplier to enforce the no‐slip constraint at the solid boundary, instead of applying constitutional relations for rigid bodies. Hence excessive constraint on the time step is circumvented, and the time step only depends on the discretization of fluid Navier‐Stokes equations, like the CFL condition in present work. To determine the boundary force, an additional moving force equation is derived. The dimension of this derived system is proportional to the number of Lagrangian points describing the solid boundaries, which makes the method very suitable for moving boundary problems since the time for matrix update and system solving is not significant. The force coefficient matrix is made symmetric and positive definite so that the conjugate gradient method can solve the system quickly. The proposed immersed boundary method is incorporated into the fluid solver with a second‐order accurate projection method as a plug‐in. The overall scheme is handled under an efficient fractional step framework, namely, prediction, forcing, and projection. Various simulations are performed to validate current method, and the results compare well with previous experimental and numerical studies.     

An immersed boundary vector potential-vorticity meshless solver of the incompressible Navier–Stokes equation

We present a strong form meshless solver for numerical solution of the nonstationary, incompressible, viscous Navier–Stokes equations in two (2D) and three dimensions (3D). We solve the flow equations in their stream function-vorticity (in 2D) and vector potential-vorticity (in 3D) formulation, by extending to 3D flows the boundary condition-enforced immersed boundary method, originally introduced in the literature for 2D problems. We use a Cartesian grid, uniform or locally refined, to discretize the spatial domain. We apply an explicit time integration scheme to update the transient vorticity equations, and we solve the Poisson type equation for the stream function or vector potential field using the meshless point collocation method. Spatial derivatives of the unknown field functions are computed using the discretization-corrected particle strength exchange method. We verify the accuracy of the proposed numerical scheme through commonly used benchmark and example problems. Excellent agreement with the data from the literature was achieved. The proposed method was shown to be very efficient, having relatively large critical time steps.     

A high‐order element‐based Galerkin method for the barotropic vorticity equation

A high‐order element‐based Galerkin method is developed to solve the non‐divergent barotropic vorticity equation (BVE). The solution process involves solving a conservative transport equation for the vorticity fields and a Poisson equation for the stream function fields. The discontinuous Galerkin method is employed for solving the transport equation and a spectral element method (continuous Galerkin) is used for the Poisson equation. A third‐order strong stability preserving explicit Runge–Kutta scheme is used for time integration. A series of tests have been performed to validate the model, which include the evolution of an idealized tropical cyclone and interaction of dual vortices in close proximity. The numerical convergence study is performed by solving the BVE on the sphere where the analytic solution is known. The test results are consistent with physical observations, and the model exhibits exponential convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamic Analysis of Linear and Nonlinear Oscillations of a Beam Under Axial and Transversal Random Poisson Pulses

This paper proposes an approximate explicit probability density function for the analysis of external and parametric oscillation of a simply supported beam driven by random pulses. The impulsive loading model adopted is Poisson white noise, which is a process having Dirac delta occurrences with random intensity distributed in time according to Poisson's law. The response probability density function can be obtained by solving the related Kolmogorov–Feller integro-differential equation. An approximate solution is derived by transforming this equation to a first-order partial differential equation. The method of characteristics is then applied to obtain an explicit solution. The theory has been validated through numerical simulations.