旋转圆管湍流的大涡模拟数值研究

采用大涡模拟（LES）方法,并结合动力学亚格子尺度应力（SGS）模型,通过数值求解柱坐标系下的滤波Navier-Stokes方程,研究了绕管轴旋转圆管内的湍流流动特性.为验证计算的可靠性,以及动力学SGS模型对于旋转湍流的适用性,将大涡模拟计算所得的结果,与相应的直接模拟（DNS）结果和实验数据进行了对比验证,吻合良好.进一步对旋转圆管湍流的物理机理进行了探讨,研究了湍流特性随旋转速率的变化规律.当旋转速率增加时,湍流流动有层流化的发展趋势.基于湍动能变化的关系,分析了旋转效应对湍流脉动生成的抑制作用.     

基于大涡模拟的移动粒子半隐式法研究

将大涡模拟(LES)和无网格的移动粒子半隐式法(MPS)相结合, 以求解湍流中的自由表面问题. 对N-S方程进行滤波计算可得到大涡模拟的控制方程, 大涡模拟的控制方程相对于以往的移动粒子半隐式法而言仅多出雷诺应力项, 通过亚粒子应力(sub-particle-scale,SPS)模型并引入Smagorinsky涡黏模型将雷诺应力模型化, 可实现移动粒子半隐式法的大涡模拟. 将MPS-LES应用至具有大变形自由表面的共振晃荡中, 其模拟结果同实验及其他数值模拟结果都相当接近.      

An efficient very large eddy simulation model for simulation of turbulent flow

Among the various hybrid methodologies, Speziale's very large eddy simulation (VLES) is one that was proposed very early. It is a unified simulation approach that can change seamlessly from Reynolds Averaged Navier–Stokes (RANS) to direct numerical simulation (DNS) depending on the numerical resolution. The present study proposes a new improved variant of the original VLES model. The advantages are achieved in two ways: (i) RANS simulation can be recovered near the wall which is similar to the detached eddy simulation concept; (ii) a LES subgrid scale model can be reached by the introduction of a third length scale, that is, the integral turbulence length scale. Thus, the new model can provide a proper LES mode between the RANS and DNS limits. This new methodology is implemented in the standard k ? ? model. Applications are conducted for the turbulent channel flow at Reynolds number of Re_τ = 395, periodic hill flow at Re = 10,595, and turbulent flow past a square cylinder at Re = 22,000. In comparison with the available experimental data, DNS or LES, the new VLES model produces better predictions than the original VLES model. Furthermore, it is demonstrated that the new method is quite efficient in resolving the large flow structures and can give satisfactory predictions on a coarse mesh. Copyright © 2012 John Wiley & Sons, Ltd.     

Improving microwave brightness temperature predictions based on Bayesian model averaging ensemble approach

The choices of the parameterizations for each component in a microwave emission model have significant effects on the quality of brightness temperature (T _b) simulation. How to reduce the uncertainty in the T _b simulation is investigated by adopting a statistical post-processing procedure with the Bayesian model averaging (BMA) ensemble approach. The simulations by the community microwave emission model (CMEM) coupled with the community land model version 4.5 (CLM4.5) over mainland China are conducted by the 24 configurations from four vegetation opacity parameterizations (VOPs), three soil dielectric constant parameterizations (SDCPs), and two soil roughness parameterizations (SRPs). Compared with the simple arithmetical averaging (SAA) method, the BMA reconstructions have a higher spatial correlation coefficient (larger than 0.99) than the C-band satellite observations of the advanced microwave scanning radiometer on the Earth observing system (AMSR-E) at the vertical polarization. Moreover, the BMA product performs the best among the ensemble members for all vegetation classes, with a mean root-mean-square difference (RMSD) of 4K and a temporal correlation coefficient of 0.64.     

Utilization of nonlinear vibrations of soft pipe conveying fluid for driving underwater bio-inspired robot

Creatures with longer bodies in nature like snakes and eels moving in water commonly generate a large swaying of their bodies or tails, with the purpose of producing significant frictions and collisions between body and fluid to provide the power of consecutive forward force. This swaying can be idealized by considering oscillations of a soft beam immersed in water when waves of vibration travel down at a constant speed. The present study employs a kind of large deformations induced by nonlinear vibrations of a soft pipe conveying fluid to design an underwater bio-inspired snake robot that consists of a rigid head and a soft tail. When the head is fixed, experiments show that a second mode vibration of the tail in water occurs as the internal flow velocity is beyond a critical value. Then the corresponding theoretical model based on the absolute nodal coordinate formulation (ANCF) is established to describe nonlinear vibrations of the tail. As the head is free, the theoretical modeling is combined with the computational fluid dynamics (CFD) analysis to construct a fluid-structure interaction (FSI) simulation model. The swimming speed and swaying shape of the snake robot are obtained through the FSI simulation model. They are in good agreement with experimental results. Most importantly, it is demonstrated that the propulsion speed can be improved by 21% for the robot with vibrations of the tail compared with that without oscillations in the pure jet mode. This research provides a new thought to design driving devices by using nonlinear flow-induced vibrations.

     

Data-driven optimization study of the multi-relaxation-time lattice Boltzmann method for solid-liquid phase change

Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes. The multi-relaxation-time lattice Boltzmann method(MRT-LBM) shows great numerical performance during simulation;however, the value method of the relaxation parameters needs to be specified. Therefore,in this study, a random forest(RF) model is used to discriminate the importance of different relaxation parameters to the convergence, and a support vector machine(SV...     

SECTIONAL FINITE ELEMENT ANALYSIS OF COUPLED DEFORMATION BETWEEN ELASTOPLASTIC SHEET METAL AND VISCO-ELASTOPLASTIC BODY

The present paper is devoted to developing a new numerical simulation method for the analysis of viscous pressure forming (VPF), which is a sheet flexible-die forming (FDF) process. The pressure-carrying medium used in VPF is one kind of semisolid, flowable and viscous material and its deformation behavior can be described by the visco-elastoplastic constitutive model. A sectional finite element model for the coupled deformation analysis between the visco- elastoplastic pressure-carrying medium and the elastoplastic sheet metal is proposed. The resolution of the Updated Lagrangian (UL) formulation is based on a static explicit approach. The frictional contact between sheet metal and visco-elastoplastic pressure-carrying medium is treated by the penalty function method. Coupled deformation between sheet metal and visco-elastoplastic pressure-carrying medium with large slip is analyzed to validate the developed algorithm. Finally, the viscous pressure bulging (VPB) process of DC06 sheet metal is simulated. Good agreement between numerical simulation results and experimental measurements shows the validity of the developed algorithm.     

Numerical study of an oscillatory turbulent flow over a flat plate

Oscillatory turbulent flow over a flat plate is studied using large eddy simulation (LES) and Reynolds-average Navier-Stokes (RANS) methods. A dynamic subgrid-scale model is employed in LES and Saffman's turbulence model is used in RANS. The flow behaviors are discussed for the accelerating and decelerating phases during the oscillating cycle. The friction force on the wall and its phase shift from laminar to turbulent regime are also investigated for different Reynolds numbers. The project supported by the Youngster Funding of Academia Sinica and by the National Natural Science Foundation of China     

Parallel FEM LES with one-equation subgrid-scale model for incompressible flows

This article develops a parallel large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model based on the Galerkin finite element method and three-dimensional (3D) brick elements. The governing filtered Navier–Stokes equations were solved by a second-order accurate fractional-step method, which decomposed the implicit velocity–pressure coupling in incompressible flow and segregated the solution to the advection and diffusion terms. The transport equation for the SGS turbulent kinetic energy was solved to calculate the SGS processes. This FEM LES model was applied to study the turbulence of the benchmark open channel flow at a Reynolds number Re_τ = 180 (based on the friction velocity and channel height) using different model constants and grid resolutions. By comparing the turbulence statistics calculated by the current model with those obtained from direct numerical simulation (DNS) and experiments in literature, an optimum set of model constants for the current FEM LES model was established. The budgets of turbulent kinetic energy and vertical Reynolds stress were then analysed for the open channel flow. Finally, the flow structures were visualised to further reveal some important characteristics. It was demonstrated that the current model with the optimum model constants can predict well the organised structure near the wall and free surface, and can be further applied to other fundamental and engineering applications.     

MPS-FEM数值分析带自由面的流固耦合问题

随着计算科学的发展,研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中,耦合的移动粒子半隐式方法 MPS(Moving Particle Semi-Implicit method)和有限单元法FEM(Finite Element method)为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上,对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外,在该求解器框架内,基于FEM法拓展了针对结构场进行求解的功能。首先,对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后,采用该自研MPS-FEM耦合求解器,数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比,验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。     

Simulation of a soil loosening process by means of the modified distinct element method

We apply the Distinct Element Method (DEM) to analyze the dynamic behavior of soil. However, the conventional DEM model for calculation of contact forces between elements has some problems; for example, the movement of elements is too discrete to simulate real soil particle movement. Therefore, we modify the model to solve the difficulties. To investigate the validity of the modified model, we conduct an experiment in which soil is cut with a pendulum-typeblade, and simulate the soil loosening process with the modified DEM model. This paper presents details of the experimental apparatus and the comparison of soil behavior and energy absorption between the simulation and the experiment. Some characteristic phenomena of the experiment are reproduced in the simulation giving us confidence that the modified model is better than the conventional model for the simulation of soil behavior.     

Numerical simulation with a TVD–FVM method for circular cylinder wake control by a fairing

A water drop-shaped fairing is applied to control the wake behind a circular cylinder and to suppress the formation of Karman vortex street in this paper. The results are evaluated using high resolution CFD technique. A finite-volume total variation diminishing (TVD) approach based upon the recently proposed elemental velocity vector transformation (EVVT) method, which aims at solving the incompressible turbulent flow for irregular boundary conditions with renormalization group (RNG) turbulence model, is used to simulate the flow field around circular cylinder systems. The calculations are carried out with cylinder systems with and without fairings, while the fairings have different top shape angles within the range of $30°~90°$. The Reynolds number ranges from 1000 to 50 000. It is shown that the simulation results of present numerical method reaches good agreement with the available experimental and numerical simulation data of typical circular cylinder flow and a fixed fairing cylinder system flow. Compared with bare cylinder, the faired bluff structures can obviously reduce the lift and drag forces and alter the vortex shedding frequency. Overall, the mean drag coefficient can be reduced up to about (10–31)% and the RMS lift coefficient can be reduced up to (30–99)% for all faired systems at given Reynolds numbers. The influence of Reynolds number and attack angles on the flow field characters of bare cylinder and faired cylinders is also discussed. The faired structures with shape angles within $30°~45°$under zero-attack-angle-inflow case are considered as the optimal structures, with which the mean drag coefficient and the RMS lift coefficient can be reduced up to (26–31)% and (98–99)%, respectively. Considering the influence of attack angles on lift and drag coefficients reduction, 75° shaped faired structure may be taken as a proper option.     

Large eddy simulation of high-Reynolds-number atmospheric boundary layer flow with improved near-wall correction

It is highly attractive to develop an efficient and flexible large eddy simulation(LES)technique for high-Reynolds-number atmospheric boundary layer(ABL)simulation using the low-order numerical scheme on a relatively coarse grid,that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer.In this study,an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the strain-rate tensor is proposed to calculate the eddy viscosity coefficient in the subgrid-scale(SGS)model.The LES code is realized with a second-order finite-difference scheme,the scale-dependent dynamic SGS stress model,the equilibrium wall stress model,and the proposed correction scheme.Very-high-Reynolds-number ABL flow simulation under the neutral stratification condition is conducted to assess the performance of the method in predicting the mean and fluctuating characteristics of the rough-wall turbulence.It is found that the logarithmic profile of the mean streamwise velocity and some key features of large-scale coherent structures can be reasonably predicted by adopting the proposed correction method and the low-order numerical scheme.     

Electro-mechanical coupling analysis of MEMS structures by boundary element method

In this paper, we present the applications of Boundary Element Method (BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems (MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics. Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples. The project supported by the 973 Program (G199033108) and the national Natural Science Foundation of China (10125211)     

A large eddy simulation of the near wake of a circular cylinder

Viscous flow around a circular cylinder at a subcritical Reynolds number is investigated using a large eddy simulation (LES) coupled with the Smagorinsky subgrid-scale (SGS) model. A fractional-step method with a second-order in time and a combined finite-difference/spectral approximations are used to solve the filtered three-dimensional incompressible Navier-Stokes equations. Calculations have been performed with and without the SGS model. Turbulence statistical behaviors and flow structures in the near wake of the cylinder are studied. Some calculated results, including the lift and drag coefficients, shedding frequency, peak Reynolds stresses, and time-average velocity profile, are in good agreement with the experimental and computational data, which shows that the Smagorinsky model can reasonably predict the global features of the flow and some turbulent statistical behaviors. The project supported by the National Science Fund for Distinguished Scholars (10125210), the Special Funds for Major State Basic Research Project (G1999032801) and the National Natural Science Foundation of China (19772062)     

Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

基于混合方法的二维水力压裂数值模拟

水力压裂在页岩气开采中被广泛使用,采用数值方法研究压裂机理具有重要意义.基于连续-非连续单元法(CDEM) 和中心型有限体积法(FVM),提出解决水力压裂流固耦合问题的二维混合数值计算模型.该混合模型中,使用CDEM 求解应力场和裂缝扩展过程,使用FVM 求解裂隙渗流场.应力场裂缝扩展和渗流场均使用显式迭代求解, 并通过相互之间数据交换实现流固耦合.通过与KGD 理论模型进行对比, 验证数值模型的正确性.通过与颗粒离散元数值结果进行对比,验证数值模型的有效性.通过计算复杂缝网压裂模型,研究水力压裂机理,并说明该数值模型在水力压裂模拟中具有很好的前景.      

Modelling of thrust generated by oscillation caudal fin of underwater bionic robot

A simplified model of the thrust force is proposed based on a caudal fin oscillation of an underwater bionic robot. The caudal fin oscillation is generalized by central pattern generators (CPGs). In this model, the drag coefficient and lift coefficient are the two critical parameters which are obtained by the digital particle image velocimetry (DPIV) and the force transducer experiment. Numerical simulation and physical experiments have been performed to verify this dynamic model.     

大涡模拟及其在湍流燃烧中的应用

大涡模拟作为一种研究湍流流动和湍流燃烧的有效手段,在国际上已经得到广泛应用。本文在回顾了大涡模拟(LES)的基本思想及其实施方法的基础上着重介绍了前人在大涡模拟的亚格子湍流模式和亚格子燃烧模式中的研究成果,同时给出了采用不同亚格子模式的大涡模拟在湍流燃烧中的应用实例,指出了大涡模拟在湍流燃烧中的重要作用,为大涡模拟的进一步发展和应用提供参考。      

Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.