A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

基于多矩的有限体积浸入边界法研究

浸入边界法（Immersed Boundary Method）是计算流体力学中求解具有复杂、移动边界流动问题的一类有效途径，该方法在笛卡尔坐标系上离散求解流体控制方程，并通过在控制方程中添加相应源相来代表浸入边界。尽管浸入边界法借助其简单、高效的显著特点在计算流体力学应用中显示出极强的生命力，特别是针对复杂的实际流动及动边界流动问题有着无可比拟的优势，但仍有许多问题需要进一步的研究。 本论文基于浸入边界方法及多矩VSIAM3（Volume/Surface Integrated Average Multi-Moment Method）格式提出了一种不可压缩流体求解数值格式。不可压N-S方程使用VSIAM3格式进行法进行离散，引入浸入边界法处理复杂、移动流动边界条件，使用虚拟网格方法计算动量方程修正项，同时还考虑了对连续方程的修正。VSIAM3格式是一种基于多矩的有限体积法，在方程的离散中总是使用两种或两种以上的矩，如：VIA（Volume Integrated Average）和SIA（Surface Integrated Average）。而不同的矩在求解过程中依据不同形式的控制方程使用不同的离散方法进行更新。VSIAM3格式更多的局地自由度及同时使用交错网格和同位网格的特点使浸入边界法的实施更加便利、高效。研究中，浸入边界法不仅应用于处理动力边界条件，同样可以处理热动力边界条件。 研究中对大量经典算例进行了数值实验，包括一维线性初始问题、方腔流问题、二维绕静止及振荡圆柱流动、三维绕球流动及热对流问题等。数值结果同实验值及其它计算结果保持一致，该算法可准确、高效处理具有复杂、移动边界及存在热对流的不可压流动问题，为实际应用打下了基础。

Thrust generation and wake structure of wiggling hydrofoil

Marine animals and micro-machines often use wiggling motion to generate thrust. The wiggling motion can be modeled by a progressive wave where its wavelength describes the flexibility of wiggling animals. In the present study, an immersed boundary method is used to simulate the flows around the wiggling hydrofoil NACA 65-010 at low Reynolds numbers. One can find from the numerical simulations that the thrust generation is largely determined by the wavelength. The thrust coefficients decrease with the increasing wavelength while the propulsive efficiency reaches a maximum at a certain wavelength due to the viscous effects. The thrust generation is associated with two different flow patterns in the wake: the well-known reversed Karman vortex streets and the vortex dipoles. Both are jet-type flows where the thrust coefficients associated with the reversed Karman vortex streets are larger than the ones associated with the vortex diploes.

一种基于多矩的有限体积浸入边界法

基于多矩VSIAM3格式及浸入边界法, 提出一套在复杂计算区域内求解不可压缩流动的数值格式.不可压N-S方程使用VSIAM3格式进行离散, 引入浸入边界法处理复杂、移动边界, 使用虚拟网格方法计算动量方程修正项, 同时还考虑了对连续方程的修正.使用标准算例对数值模式进行验证

Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Compressible laminar boundary-layer flows of a dusty gas over a semi-infinite flat-plate

The compressible laminar boundary-layer flows of a dilute gas-particle mixture over a semi-infinite flat plate are investigated analytically. The governing equations are presented in a general form where more reasonable relations for the two-phase interaction and the gas viscosity are included. The detailed flow structures of the gas and particle phases are given in three distinct regions : the large-slip region near the leading edge, the moderate-slip region and the small-slip region far downstream. The asymptotic solutions for the two limiting regions are obtained by using a seriesexpansion method. The finite-difference solutions along the whole length of the plate are obtained by using implicit four-point and six-point schemes. The results from these two methods are compared and very good agreement is achieved. The characteristic quantities of the boundary layer are calculated and the effects on the flow produced by the particles are discussed. It is found that in the case of laminar boundary-layer flows, the skin friction and wall heat-transfer are higher and the displacement thickness is lower than in the pure-gas case alone. The results indicate that the Stokes-interaction relation is reasonable qualitatively but not correct quantitatively and a relevant non-Stokes relation of the interaction between the two phases should be specified when the particle Reynolds number is higher than unity.

Short-surface waves radiated by a partially-immersed 3-dimensional oscillating body

The short-surface waves generated by a 3-D arbitrarily oscillating body floating onwater are discussed. In the far-field off the body, the phase and the amplitude functions ofthe radiated waves are determined by the ray method. An undetermined constant is includ-ed in the amplitude function. From the result of Ref. [1], the near-field boundary layersolution near the body waterline is obtained. The amplitude of this solution depends on thewhole wall shape of the body and the slope at the body waterline on the cross-sections per-pendicular to the waterline. By matching the far-field solution with the near-field bound-ary layer solution, the undetermined constant in the amplitude function of the far-fieldradiated waves is determined. For the special case of a half-submerged sphere which per-forms vertical oscillating motion, the result obtained in this paper is in agreement withthat of Ref. [ 2 ].

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Experimental Study Of Vortex-Induced Vibrations Of A Cylinder Near A Rigid Plane Boundary In Steady Flow

In this study, the vortex-induced vibrations of a cylinder near a rigid plane boundary in a steady flow are studied experimentally. The phenomenon of vortex-induced vibrations of the cylinder near the rigid plane boundary is reproduced in the flume. The vortex shedding frequency and mode are also measured by the methods of hot film velocimeter and hydrogen bubbles. A parametric study is carried out to investigate the influences of reduced velocity, gap-to-diameter ratio, stability parameter and mass ratio on the amplitude and frequency responses of the cylinder. Experimental results indicate: (1) the Strouhal number (St) is around 0.2 for the stationary cylinder near a plane boundary in the sub-critical flow regime; (2) with increasing gap-to-diameter ratio (e (0)/D), the amplitude ratio (A/D) gets larger but frequency ratio (f/f (n) ) has a slight variation for the case of larger values of e (0)/D (e (0)/D > 0.66 in this study); (3) there is a clear difference of amplitude and frequency responses of the cylinder between the larger gap-to-diameter ratios (e (0)/D > 0.66) and the smaller ones (e (0)/D < 0.3); (4) the vibration of the cylinder is easier to occur and the range of vibration in terms of V (r) number becomes more extensive with decrease of the stability parameter, but the frequency response is affected slightly by the stability parameter; (5) with decreasing mass ratio, the width of the lock-in ranges in terms of V (r) and the frequency ratio (f/f (n) ) become larger.

Fast and effective analysis methods for fiber Raman amplifiers

A series of new single-step methods and their corresponding algorithms with automatic step size adjustment for model equations of fiber Raman amplifiers are proposed and compared in this paper. On the basis of the Newton-Raphson method, multiple shooting algorithms for the two-point boundary value problems involved in solving Raman amplifier propagation equations are constructed. A verified example shows that, compared with the traditional Runge-Kutta methods, the proposed methods can increase the accuracy by more than two orders of magnitude under the same conditions. The simulations for Raman amplifier propagation equations demonstrate that our methods can increase the computing speed by more than 5 times, extend the step size significantly, and improve the stability in comparison with the Dormand-Prince method. The numerical results show that the combination of the multiple shooting algorithms and the proposed methods has the capacity to rapidly and effectively solve the model equations of multipump Raman amplifiers under various conditions such as co-, counter- and bi-directionally pumped schemes, as well as dual-order pumped schemes.

An augmented method for free boundary problems with moving contact lines

An augmented immersed interface method (IIM) is proposed for simulating one-phase moving contact line problems in which a liquid drop spreads or recoils on a solid substrate. While the present two-dimensional mathematical model is a free boundary problem, in our new numerical method, the fluid domain enclosed by the free boundary is embedded into a rectangular one so that the problem can be solved by a regular Cartesian grid method. We introduce an augmented variable along the free boundary so that the stress balancing boundary condition is satisfied. A hybrid time discretization is used in the projection method for better stability. The resultant Helmholtz/Poisson equations with interfaces then are solved by the IIM in an efficient way. Several numerical tests including an accuracy check, and the spreading and recoiling processes of a liquid drop are presented in detail. (C) 2010 Elsevier Ltd. All rights reserved.