自由上浮气泡运动特性的光滑粒子流体动力学模拟

基于虚功原理,在Hu X Y等和Grenier N等的研究结果基础上推导了多相流光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)控制方程,采用精度较高的黏性力和表面张力模型,发展了一套适用于具有大密度比和大黏性比界面的多相流SPH方法.首先,通过施加人工位移修正,适当背景压力和异相界面力,使得计算全程粒子分布相对均匀,改善了界面处的失稳现象,防止了异相界面处粒子的非物理性穿透;在此基础上,利用方形流体团振荡模型对表面张力模型进行了验证,数值结果与解析解甚为吻合;然后采用上浮气泡经典数值算例对比研究了不同黏性力计算方法、不同核函数的适用性以及人工位移修正的效果;最后,对单个气泡的上浮、变形、撕裂以及垂向两个气泡的追赶、融合等现象进行了模拟,初步揭示了气泡上浮过程中各种有趣物理现象的细节过程和动力学机理.

Numerical simulation of the self-propulsive motion of a fishlike swimming foil using the δ^+-SPH model

The present work is dedicated to the application of the recently developed （δ＋ -SPH） scheme to the self-propulsive fishlike swimming hydrodynamics. In the numerical method, a particle shifting technique （PST） is implemented in the framework of δ-SPH, combining with an adaptive particle refinement （APR） which is a numerical technique adopted to refine the particle resolution in the local region and de-refine particles outside that region. This comes into being the so-called δ＋- SPH scheme which contributes to higher numerical accuracy and efficiency. In the fishlike swimming modeling, a NACA0012 profile is controlled to perform a wavy motion mimicking the fish swimming in water. Thanks to the mesh-free characteristic of SPH method, the NACA0012 profile can undergo a wavy motion with large amplitude and move forward freely, avoiding the problem of mesh distortion. A parallel staggered algorithm is adopted to perform the fluidstructure interaction between the foil and the surrounding fluid. Two different approaches are adopted for the fishlike swimming problem. In Approach 1, the foil is fixed and flaps in a free stream and in Approach 2, the wavy foil can move forward under the self-driving force. The numerical results clearly demonstrate the capability of the δ＋ -SPH scheme in modeling such kind of self-propulsive fishlike swimming problems.

Smoothed particle hydrodynamics and its applications in fluid-structure interactions

In ocean engineering, the applications are usually related to a free surface which brings so many interesting physical phenomena （e.g. water waves, impacts, splashing jets, etc.）. To model these complex free surface flows is a tough and challenging task for most computational fluid dynamics （CFD） solvers which work in the Eulerian framework. As a Lagrangian and meshless method, smoothed particle hydrodynamics （SPH） offers a convenient tracking for different complex boundaries and a straightforward satisfaction for different boundary conditions. Therefore SPH is robust in modeling complex hydrodynamic problems characterized by free surface boundaries, multiphase interfaces or material discontinuities. Along with the rapid development of the SPH theory, related numerical techniques and high-performance computing technologies, SPH has not only attracted much attention in the academic community, but also gradually gained wide applications in industrial circles. This paper is dedicated to a review of the recent developments of SPH method and its typical applications in fluid-structure interactions in ocean engineering. Different numerical techniques for improving numerical accuracy, satisfying different boundary conditions, improving computational efficie- ncy, suppressing pressure fluctuations and preventing the tensile instability, etc., are introduced. In the numerical results, various typical fluid-structure interaction problems or multiphase problems in ocean engineering are described, modeled and validated. The prospective developments of SPH in ocean engineering are also discussed.

Recent progress in CFD for naval architecture and ocean engineering

An overview is provided of CFDShip-Iowa modeling, numerical methods and high performance computing （HPC）, including both current V4.5 and V5.5 and next generation V6. Examples for naval architecture highlight capability and needs. High fidelity V6 simulations for ocean engineering and fundamental physics describe increased resolution for analysis of physics of fluids. Uncertainty quantification research is overviewed as the first step towards development stochastic optimization.

A one-dimensional polynomial chaos method in CFD–Based uncertainty quantification for ship hydrodynamic performance

A one-dimensional non-intrusive Polynomial Chaos （PC） method is applied in Uncertainty Quantification （UQ） studies for CFD-based ship performances simulations. The uncertainty properties of Expected Value （EV） and Standard Deviation （SD） are evaluated by solving the PC coefficients from a linear system of algebraic equations. The one-dimensional PC with the Legendre polynomials is applied to： （1） stochastic input domain and （2） Cumulative Distribution Function （CDF） image domain, allowing for more flexibility. The PC method is validated with the Monte-Carlo benchmark results in several high-fidelity, CFD-based, ship UQ problems, evaluating the geometrical, operational and environmental uncertainties for the Delft Catamaran 372. Convergence is studied versus PC order P for both EV and SD, showing that high order PC is not necessary for present applications. Comparison is carried out for PC with/without the least square minimization when solving the PC coefficients. The least square minimization, using larger number of CFD samples, is recommended for current test cases. The study shows the potentials of PC method in Robust Design Optimization （RDO） and Reliability-Based Design Optimization （RBDO） of ship hydrodynamic performances.

High-speed water impacts of flat plates in different ditching configuration through a Riemann--ALE SPH model

The violent water entry of flat plates is investigated using a Riemann-arbitrary Eulerian-Lagrangian（ALE） smoothed particle hydrodynamics（SPH） model. The test conditions are of interest for problems related to aircraft and helicopter emergency landing in water. Three main parameters are considered： the horizontal velocity, the approach angle（i.e., vertical to horizontal velocity ratio） and the pitch angle, a. Regarding the latter, small angles are considered in this study. As described in the theoretical work by Zhao and Faltinsen（1993）, for small a a very thin, high-speed jet of water is formed, and the time-spatial gradients of the pressure field are extremely high. These test conditions are very challenging for numerical solvers. In the present study an enhanced SPH model is firstly tested on a purely vertical impact with deadrise angle α=4°. An in-depth validation against analytical solutions and experimental results is carried out, highlighting the several critical aspects of the numerical modelling of this kind of flow, especially when pressure peaks are to be captured. A discussion on the main difficulties when comparing to model scale experiments is also provided. Then, the more realistic case of a plate with both horizontal and vertical velocity components is discussed and compared to ditching experiments recently carried out at CNR-INSEAN. In the latter case both 2-D and 3-D simulations are considered and the importance of 3-D effects on the pressure peak is discussed for α=4° and α=10°.

Marine propellers performance and flow-field prediction by a free-wake panel method

A Boundary Element Method （BEM） hydrodynamics combined with a flow-alignment technique to evaluate blades shed vorticity is presented and applied to a marine propeller in open water. Potentialities and drawbacks of this approach in capturing propeller performance, slipstream velocities, blade pressure distribution and pressure disturbance in the flow-field are highlighted by comparisons with available experiments and RANSE results. In particular, correlations between the shape of the convected vortex- sheet and the accuracy of BEM results are discussed throughout the paper. To this aim, the analysis of propeller thrust and torque is the starting point towards a detailed discussion on the capability of a 3-D free-wake BEM hydrodynamic approach to describe the local features of the flow-field behind the propeller disk, in view of applications to propulsive configurations where the shed wake plays a dominant role.

Simplified hydrodynamic models for the analysis of marine propellers in a wake-field

This paper presents a comparison among different hydrodynamic models for the analysis of the unsteady loads delivered by a marine propeller working in an axial, non-uniform inflow. Specifically, for a propeller subjected to a wake-field dominated by local high-frequency changes in space, the unsteady hydroloads predicted by the Nakatake formulation are compared with those given by the Theodorsen and Sears theories, respectively. Drawbacks and potentialities of these approaches are highlighted to assess a computationally efficient hydrodynamic solver for the analysis of operating conditions where propeller blades are significantly perturbed by a multi-harmonic onset-flow. Guidelines coming from this investigation may drive the choice of a fast and reliable unsteady propeller modeling that represents a good trade-off between accuracy of simulation and cost of computation within implementation in Computational Fluid Dynamics （CFD） solvers. The hydrodynamic formulations herein proposed are validated through numerical comparisons with the （accurate but computationally expensive） propeller loads predicted by a fully 3-D panel-method Boundary Element Method （BEM） solver, suited for the analysis of propellers operating in a complex hydrodynamic environment.

A 3-D SPH model for simulating water flooding of a damaged floating structure

With the quasi-static analysis method, the terminal floating state of a damaged ship is usually evaluated for the risk assessment. But this is not enough since the ship has the possibility to lose its stability during the transient flooding process. Therefore, an enhanced smoothed particle hydrodynamics（SPH） model is applied in this paper to investigate the response of a simplified cabin model under the condition of the transient water flooding. The enhanced SPH model is presented firstly including the governing equations, the diffusive terms, the boundary implementations and then an algorithm regarding the coupling motions of six degrees of freedom（6-DOF） between the structure and the fluid is described. In the numerical results, a non-damaged cabin floating under the rest condition is simulated. It is shown that a stable floating state can be reached and maintained by using the present SPH scheme. After that, three-dimensional（3-D） test cases of the damaged cabin with a hole at different locations are simulated. A series of model tests are also carried out for the validation. Fairly good agreements are achieved between the numerical results and the experimental data. Relevant conclusions are drawn with respect to the mechanism of the responses of the damaged cabin model under water flooding conditions.