A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid

This paper presents a finite element method for incompressible multiphase flows with capillary interfaces based on a (formally) second‐order projection scheme. The discretization is on a fixed Eulerian grid. The fluid phases are identified and advected using a level set function. The grid is temporarily adapted around the interfaces in order to maintain optimal interpolations accounting for the pressure jump and the discontinuity of the normal velocity derivatives. The least‐squares method for computing the curvature is used, combined with piecewise linear approximation to the interface. The time integration is based on a formally second order splitting scheme. The convection substep is integrated over an Eulerian grid using an explicit scheme. The remaining generalized Stokes problem is solved by means of a formally second order pressure‐stabilized projection scheme. The pressure boundary condition on the free interface is imposed in a strong form (pointwise) at the pressure‐computation substep. This allows capturing significant pressure jumps across the interface without creating spurious instabilities. This method is simple and efficient, as demonstrated by the numerical experiments on a wide range of free‐surface problems. Copyright © 2004 John Wiley & Sons, Ltd.     

The analysis of unsteady incompressible flows by a three-step finite element method

This paper describes a three-step finite element method and its applications to unsteady incompressible fluid flows. Stability analysis of the one-dimensional pure convection equation shows that this method has third-order accuracy and an extended numerical stability domain in comparison with the Lax--Wendroff finite element method. The method is cost-effective for incompressible flows because it permits less frequent updates of the pressure field with good accuracy. In contrast with the Taylor-Galerkin method, the present method does not contain any new higher-order derivatives, which makes it suitable for solving non-linear multidimensional problems and flows with complicated boundary conditions. The three-step finite element method has been used to simulate unsteady incompressible flows. The numerical results obtained are in good agreement with those in the literature.     

半浮区液桥热毛细振荡流

采用非定常、三维直接数值模拟方法研究大Pr数半浮区液桥热毛细对流从定常流向振荡流的过渡过程．文中详细描述了热毛细振荡流的起振和振荡特征,给出了液桥横截面上振荡流的流场和温度分布．在地面引力场条件下计算的结果与地面实验的结果进行比较,得出液桥水平截面上的流场和温度分布图样以一定的速度旋转,自由表面固定点处流体的环向流速正、负交替变化的一致结论．     

A fluid solver based on vorticity–helical density equations with application to a natural convection in a cubic cavity

We study numerically a recently introduced formulation of incompressible Newtonian fluid equations in vorticity–helical density and velocity–Bernoulli pressure variables. Unlike most numerical methods based on vorticity equations, the current approach provides discrete solutions with mass conservation, divergence‐free vorticity, and accurate kinetic energy balance in a simple and natural way. The method is applied to compute buoyancy‐driven flows in a differentially heated cubic enclosure in the Boussinesq approximation for Ra ∈ {10⁴,10⁵,10⁶}. The numerical solutions on a finer grid are of benchmark quality. The computed helical density allows quantification of the three‐dimensional nature of the flow. Copyright © 2011 John Wiley & Sons, Ltd.     

Non-darcy double-diffusive mixed convection from heated vertical and horizontal plates in saturated porous media

The non-darcy mixed convection flows from heated vertical and horizontal plates in saturated porous media have been considered using boundary layer approximations. The flows are considered to be driven by multiple buoyancy forces. The similarity solutions for both vertical and horizontal plates have been obtained. The governing equations have been solved numerically using a shooting method. The heat transfer, mass transfer and skin friction are reduced due to inertial forces. Also, they increase with the buoyancy parameter for aiding flow and decrease for the opposing flow. For aiding flow, the heat and mass transfer coefficients are found to approach asymptotically the forced or free convection values as the buoyancy parameter approaches zero or infinity.     

Heat transfer in rotating cylindrical cells with partitions

A numerical analysis and an experimental study of heat transfer rate in rotating cylindrical cells with partitions are performed. The work is done mainly in the Ekman suction regime, where the Coriolis force dominates over centrifugal buoyancy. It is shown that the heat transfer rate increases substantially by placing partitions in the cell. The partitions suppress the Coriolis force so that convection induced by the centrifugal buoyancy becomes important. It is found that the Nusselt number correlates with the parameter PrβΔTEk^?1/2 with the partitions. The partitions have no effect on the heat transfer in the centrifugal buoyancy convection regime.     

Pseudostress–velocity formulation for incompressible Navier–Stokes equations

This paper presents a numerical algorithm using the pseudostress–velocity formulation to solve incompressible Newtonian flows. The pseudostress–velocity formulation is a variation of the stress–velocity formulation, which does not require symmetric tensor spaces in the finite element discretization. Hence its discretization is greatly simplified. The discrete system is further decoupled into an H ( div ) problem for the pseudostress and a post‐process resolving the velocity. This can be done conveniently by using the penalty method for steady‐state flows or by using the time discretization for nonsteady‐state flows. We apply this formulation to the 2D lid‐driven cavity problem and study its grid convergence rate. Also, computational results of the time‐dependent‐driven cavity problem and the flow past rectangular problem are reported. Copyright © 2009 John Wiley & Sons, Ltd.     

An implicit Lagrangian lattice Boltzmann method for the compressible flows

In this paper, we propose a new Lagrangian lattice Boltzmann method (LBM) for simulating the compressible flows. The new scheme simulates fluid flows based on the displacement distribution functions. The compressible flows, such as shock waves and contact discontinuities are modelled by using Lagrangian LBM. In this model, we select the element in the Lagrangian coordinate to satisfy the basic fluid laws. This model is a simpler version than the corresponding Eulerian coordinates, because the convection term of the Euler equations disappears. The numerical simulations conform to classical results. Copyright © 2006 John Wiley & Sons, Ltd.     

Simultaneous heat and mass transfer in unsteady free convection from two-dimensional and axisymmetric bodies

The unsteady laminar free convection boundary layer flows around two-dimensional and axisymmetric bodies placed in an ambient fluid of infinite extent have been studied when the flow is driven by thermal buoyancy forces and buoyancy forces from species diffusion. The unsteadiness in the flow field is caused by both temperature and concentration at the wall which vary arbitrarily with time. The coupled nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. Computations have been performed for a circular cylinder and a sphere. The skin friction, heat transfer and mass transfer are strongly dependent on the variation of the wall temperature and concentration with time. Also the skin friction and heat transfer increase or decrease as the buoyancy forces from species diffusion assist and oppose, respectively, the thermal buoyancy force, whereas the mass transfer rate is higher for small values of the ratio of the buoyancy parameters than for large values. The local heat and mass transfer rates are maximum at the stagnation point and they decrease progressively with increase of the angular position from the stagnation point.     

Numerical simulation of transient planar flow of a viscoelastic material with two moving free surfaces

In this study, we examine the numerical simulation of transient viscoelastic flows with two moving free surfaces. A modified Galerkin finite element method is implemented to the two-dimensional non-steady motion of the fluid of the Oldroyd-B type. The fluid is initially placed between two parallel plates and bounded by two straight free boundaries. In this Lagrangian finite element method, the spatial mesh deforms in time along with the moving free boundaries. The unknown shape of the free surfaces is determined with the flow field u, v, τ, p by the deformable finite element method, combined with a predictor-corrector scheme in an uncoupled fashion. The moving free surfaces and fluid motion of both Newtonian and non-Newtonian flows are investigated. The results include the influence of surface tension, fluid inertia and elasticity.     

浮力对混合对流流动及换热特性的影响

用热线和冷线相结合的技术测量垂直圆管内逆混合对流流体的平均速度、 温度以及它们的脉动. 较详细地研究了浮力对逆混合对流的流动特性和传热特性的影响. 评 估了实验中采用的冷线测量温度补偿速度探头温度敏感的影响. 逆混合对流的传热结果用无 量纲参数Ω (Ω= Grd / Red2 )来表示，其中，基于管道直 径的雷诺数Red变化范围为900~18000, 浮力参数Ω变化范围为 0.004899~0.5047. 研究结果表明，浮力对逆混合对流的换热有强化作用. 随着葛拉晓夫数Grd的增加，温度脉动，流向雷诺正应力和流向温度通量增 大，并且在靠近壁面的流体区域尤其明显. 热线与冷线相结合的技术适合于研究非绝热的流 动测量，可以用于研究浮力对流动和换热特性的影响.     

Simulation of non‐hydrostatic,density‐stratified flow in irregular domains

A numerical model has been developed for simulating density‐stratified flow in domains with irregular but simple topography. The model was designed for simulating strong interactions between internal gravity waves and topography, e.g. exchange flows in contracting channels, tidally or convectively driven flow over two‐dimensional sills or waves propagating onto a shoaling bed. The model is based on the non‐hydrostatic, Boussinesq equations of motion for a continuously stratified fluid in a rotating frame, subject to user‐configurable boundary conditions. An orthogonal boundary fitting co‐ordinate system is used for the numerical computations, which rely on a fourth‐order compact differentiation scheme, a third‐order explicit time stepping and a multi‐grid based pressure projection algorithm. The numerical techniques are described and a suite of validation studies are presented. The validation studies include a pointwise comparison of numerical simulations with both analytical solutions and laboratory measurements of non‐linear solitary wave propagation. Simulation results for flows lacking analytical or laboratory data are analysed a posteriori to demonstrate satisfaction of the potential energy balance. Computational results are compared with two‐layer hydraulic predictions in the case of exchange flow through a contracting channel. Finally, a simulation of circulation driven by spatially non‐uniform surface buoyancy flux in an irregular basin is discussed. Copyright © 2000 John Wiley & Sons, Ltd.     

Mixed convection boundary layer flow of non-Newtonian fluids along vertical wavy plates

Mixed convection boundary layer flows of non-Newtonian fluids over the wavy surfaces are studied by the coordinate transformation and the cubic spline collocation numerical method. The effects of the wavy geometry, the buoyancy parameter and the generalized Prandtl number for pseudoplastic fluids, Newtonian fluids and dilatant fluids on the skin-friction coefficient, local and mean Nusselt numbers have been graphically studied. Results show that both higher generalized Prandtl numbers and buoyancy parameters are seen to enhance the influence of wavy surfaces on the local Nusselt number, irrespective of whether the fluids are Newtonian fluids or non-Newtonian fluids. Moreover, the irregular surfaces have higher total heat flux than that of corresponding flats plate for any fluid.     

数值摄动算法及其CFD格式

作者提出的数值摄动算法把流体动力学效应耦合进NS方程组和对流扩散（CD）方程离散的数学基本格式（MBS）,特别是耦合进最简单的MBS即一阶迎风和二阶中心格式之中,由此构建成一系列新格式,称呼方便和强调耦合流体动力学起见,称它们为流体力学基本格式（FMBS）。构建FMBS的主要步骤是把MBS中的通量摄动重构为步长的幂级数,利用空间分裂和导出的高阶流体动力学关系式,把结点变量展开成Taylor级数,通过消除重构格式修正微分方程的截断误差诸项求出幂级数的待定系数,由此获得非线性FMBS。FMBS的公式是MBS与 （及 ）之简单多项式的乘积, 和 分别是网格Reynolds数和网格CFL数。FMBS和MBS使用相同结点,简单性彼此相当,但FMBS精度高稳定范围大,例如FMBS包含了许多绝对稳定和绝对正型、高阶迎风和中心有限差分（FD）格式和有限体积（FV）格式,这些格式对网格Reynolds数的任意值均为不振荡格式。可见对不振荡CFD格式的构建,数值摄动算法提供了不同于调节数值耗散等常见的人为构建方法,而利用流体力学自身关系以及把迎风机制通过上、下游摄动重构引入中心MBS的解析构建方法,FMBS除了直接应用于流体计算外;对于通过调节数值耗散、色散和数值群速度特性重构高分辨率格式的研究,最简单FMBS提供了比最简单MBS更精确、但同样简单的基础和起步格式。FMBS用于计算不可压缩流,可压缩流,液滴萃取传质,微通道两相流等,均获得良好数值结果或与已有Benchmark解一致的数值结果。已有文献称数值摄动算法为新型高精度格式和高的算法和高的格式;本文FMBS比数值摄动格式的称呼可更好反映FMBS的物理内容。文中也讨论了值得进一步研究的一些课题,该法亦可用于其它一些数学物理方程（例如,简化Boltzmann方程、磁流体方程、KdV-Burgers方程等）MBS耦合物理动力学效应的重构。      

Natural convection of non-Newtonian power law fluids in a shallow horizontal rectangular cavity uniformly heated from below

Analytical and numerical study is conducted for two-dimensional steady-state buoyancy driven flow of a non-Newtonian power law fluid confined in a shallow rectangular horizontal cavity uniformly heated from below, while its short vertical rigid sides are considered adiabatic. The effect of the non-Newtonian behaviour on the onset of convection, fluid flow, temperature field, and heat transfer is examined. A closed approximate analytical solution is developed on the basis of the parallel flow assumption and the obtained results are validated numerically by solving the full governing equations.     

Unsteady boundary layer and its separation over a heated circular cylinder

A numerical method is developed to solve the coupled unsteady laminar momentum and thermal boundary layers over a circular cylinder impulsively started from rest. The present non-iterative finite difference method, which requires relatively fewer grid points in the reversed flow region than any other method, can easily handle the separating boundary layer flows. The results indicate that the present method has accuracy comparable with the earlier methods, while consuming computer time approximately one order of magnitude less. The present numerical method allowed investigation of the effect of buoyancy parameter on the starting boundary layer. The time-dependent behaviour of the boundary layer is studied in terms of the appearance of the singularity, the distribution of skin friction and wall heat flux, and the wall position of the inflection point of the velocity profile. The transient as well as buoyancy-dependent patterns of the streamlines and isotherms are also studied.     

Similarity solutions of mixed convection heat and mass transfer in combined stagnation and rotation-induced flows over a rotating disk

The objective of the present study is to develop a novel similarity model for analysis of mixed convection heat and mas transfer in combined stagnation and rotation-induced flows over a rotating disk. Thermal and concentration (solutal) buoyancy effects stemmed from temperature and concentration gradients in rotational as well as gravitational forces fields are all taken into account. The influences of the forced flow, disk rotation, thermal buoyancy, buoyancy ratio and the fluid properties, i.e. Prandtl and Schmidt numbers, on the flow, temperature and concentration fields and the associated friction factors, heat and mass transfer rates are investigated. The present results reveal the effects of various buoyancy modes with combined forces on the transport phenomena in rotating-disk flows, and the analysis is also useful in understanding the mechanisms of mixed convection in the class of rotating fluids. Received on 30 December 1997     

Turbulence Models and Large Eddy Simulations Applied to Ascending Mixed Convection Flows

Coolant flows in the cores of current gas-cooled nuclear reactors consist of ascending vertical flows in a large number of parallel passages. Under post-trip conditions such heated turbulent flows may be significantly modified from the forced convection condition by the action of buoyancy, and the thermal-hydraulic regime is no longer one of pure forced convection. These modifications are associated primarily with changes to the turbulence structure. Flow laminarization may occur, and in that event heat transfer rates may be as low as 40% of those in the corresponding forced convection case. The present work is concerned with the modelling of such ??mixed?? convection flows in a vertical heated pipe. All fluid properties are assumed to be constant and buoyancy is accounted for within the Boussinesq approximation. Six different Eddy Viscosity Models (EVMs) are examined against experimental measurements and the direct numerical simulation (DNS) data of You et?al. (Int J Heat Mass Transfer 46:1613?C1627, 2003). The EVMs selected for study embody distinct physical refinements with respect to the parent high-Reynolds-number k-?? model. Large Eddy Simulations employing the classical Smagorinsky sub-grid-scale model are also presented. It is found that the early EVM scheme of Launder and Sharma (1974) is the turbulence model in the closest agreement with direct simulation results for the ratio of mixed-to-forced convection Nusselt number, Nu/Nu ₀; the model in the poorest accord with the DNS data is a k-??-SST formulation. However, in relation to comparisons with both numerical and experimental data for forced convection Nusselt number, Nu ₀, the present work reveals that some of the more recent models perform better than the Launder-Sharma scheme. No single scheme is in consistently close agreement with the numerical simulation flow profiles.