Fourier analysis of an equal‐order incompressible flow solver stabilized by pressure gradient projection

Fourier analysis techniques are applied to the stabilized finite element method (FEM) recently proposed by Codina and Blasco for the approximation of the incompressible Navier–Stokes equations, here denoted by pressure gradient projection (SPGP) method. The analysis is motivated by spurious waves that pollute the computed pressure in start‐up flow simulation. An example of this spurious phenomenon is reported. It is shown that Fourier techniques can predict the numerical behaviour of stabilized methods with remarkable accuracy, even though the original Navier–Stokes setting must be significantly simplified to apply them. In the steady state case, good estimates for the stabilization parameters are obtained. In the transient case, spurious long waves are shown to be persistent when the element Reynolds number is large and the Courant number is small. This can be avoided by treating the pressure gradient projection implicitly, though this implies additional computing effort. Standard extrapolation variants are unfortunately unstable. Comparisons with Galerkin–least‐squares (GLS) method and Chorin's projection method are also addressed. Copyright © 2000 John Wiley & Sons, Ltd.     

The consistency of pressure‐gradient approximations used in multi‐dimensional shock hydrodynamics

This paper presents a study of the consistency properties of the pressure‐gradient approximation used in multi‐dimensional finite‐element shock hydrodynamics codes today. In specific, consideration is given to the so‐called ‘bent‐element blues’ problem associated with the pressure‐gradient approximation when using the Q1Q0 element. On arbitrary grids comprised of distorted elements, the piecewise‐constant representation of the pressure field leads to a low‐order pressure‐gradient approximation at the global (nodal) level. This results in spurious nodal forces that are not aligned with the pressure gradient. There are several side‐effects of this behavior that include (a) incorrectly exciting physical modes in problems that exhibit unstable behavior, e.g. Rayleigh–Taylor problems (both magnetic and hydrodynamic), (b) potentially seeding hourglass modes, and (c) exhibiting non‐stationary behavior for steady‐state problems. A series of commonly used pressure‐gradient approximations are reviewed and evaluated based on linear consistency—the ability of the approximation to annihilate constant terms and exactly reproduce a linear gradient. The deeper theoretical issues associated with the proper selection of function spaces for the finite‐element hydro formulation are not discussed here. There are two gradient approximations that use piecewise‐constant data and deliver a consistent pressure‐gradient approximation on arbitrary grids. The first is the well‐known least‐squares gradient construction, and the second is a corrected gradient approximation that imposes linear consistency at the (global) nodal level. At the time of this writing, the corrected gradient approximation appears to be the most viable candidate for resolving the consistency issues associated with the Q1Q0 element technology. Copyright © 2009 John Wiley & Sons, Ltd.     

Wake-boundary layer interaction subject to convex and concave curvatures and adverse pressure gradient

Measurements of mean velocity and turbulent quantities have been carried out when the wake of a symmetrical airfoil interacts with the boundary layer on the (i) walls of a straight duct/diffuser and (ii) convex and concave walls of a curved duct/diffuser. The effects of adverse pressure gradient and of curvatures on the interaction are studied separately and in combination. Six cases are considered, viz. with (i) neither pressure gradient nor curvature, (ii) adverse pressure gradient and no curvature, (iii) and (iv) convex curvature with zero and adverse pressure gradients, respectively, (v) and (vi) concave curvature with zero and adverse pressure gradients, respectively. For the flows with curvature, the curvature parameter δ/R is 0.023, and for the flows with adverse pressure gradient, the Clauser pressure gradient parameter β is 0.62. The individual influences of adverse pressure gradient and convex and concave curvatures on the boundary layer are similar to those observed by earlier investigations. It is further observed that the combined effect of concave/convex curvature and the adverse pressure gradient causes higher turbulence intensities than the sum of the individual effects. The effect of curvature is to make the wake asymmetric, and in combination with adverse pressure gradient the asymmetry increases. It is observed that the adverse pressure gradient causes faster wake–boundary-layer interaction. Comparing measurements in a straight duct, a curved duct, a curved diffuser and a straight diffuser, it is seen that the convex curvature reduces the boundary layer thickness. The asymmetry in wake development compensates for this effect and the wake–boundary-layer interaction on a convex surface is almost the same as that on a straight surface. In the case of interaction with the boundary layer on a concave surface, the curvature increases the boundary layer thickness and causes enhanced turbulence intensities. However, the asymmetry in wake is such that the extent of wake is lower towards the boundary layer side. As a result, the wake–boundary-layer interaction on concave surface is almost the same as on a straight surface. The interaction is faster in the presence of adverse pressure gradient. Received: 16 June 2000 / Accepted: 17 May 2001     

A model for pressure loss, film thickness, and entrained fraction for gas-liquid annular flow

A two-component (air-water) annular flow model is presented requiring only flow rates, absolute pressure, temperature, and tube diameter. Film thicknesses (base film and wave height) are calculated from a critical film thickness model. Modeled pressure gradient is weighted by wave intermittency to compute average pressure gradient. Film flow rate and wave velocity are estimated using the universal velocity profile in the waves and a piecewise linear profile in the base film. For vertical flow, mean absolute errors for film thickness, wave velocity, and pressure gradient are 9%, 9%, and 19%, respectively. In horizontal flow, mean absolute errors for pressure gradient, base film thickness, and disturbance wave velocity are 17%, 10%, and 14%, respectively, on par with those from single-behavior models that require additional film thickness or other data as inputs.     

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

A parametric study of adverse pressure gradient turbulent boundary layers

There are many open questions regarding the behaviour of turbulent boundary layers subjected to pressure gradients and this is confounded by the large parameter space that may affect these flows. While there have been many valuable investigations conducted within this parameter space, there are still insufficient data to attempt to reduce this parameter space. Here, we consider a parametric study of adverse pressure gradient turbulent boundary layers where we restrict our attention to the pressure gradient parameter, β, the Reynolds number and the acceleration parameter, K. The statistics analyzed are limited to the streamwise fluctuating velocity. The data show that the mean velocity profile in strong pressure gradient boundary layers does not conform to the classical logarithmic law. Moreover, there appears to be no measurable logarithmic region in these cases. It is also found that the large-scale motions scaling with outer variables are energised by the pressure gradient. These increasingly strong large-scale motions are found to be the dominant contributor to the increase in turbulence intensity (scaled with friction velocity) with increasing pressure gradient across the boundary layer.     

Magnetothermoelastic transient response of a functionally graded thick hollow sphere subjected to magnetic and thermoelastic fields

In this article, an analytical method is developed to obtain the response of magnetothermoelastic stress and perturbation of the magnetic field vector for a thick-walled spherical functionally graded materials (FGM) vessel. The vessel, which is placed in a uniform magnetic field, is subjected to an internal pressure and transient temperature gradient. Using the Hankel and Laplace transform techniques, the dynamic equation of magnetothermoelastic is solved and the radial and circumferential stresses as well as the perturbation of the magnetic field vector for a typical material are obtained. Moreover, the effect of magnetic field vector and material inhomogeneity on the stresses is investigated.     

Transient pressure of percolation through one dimension porous media with threshold pressure gradient

IntroductionManyexperimentsshow[1,2,3]thatthepercolationinlowPermeabilityisnotfittoDarcy'sLaw.Theparticularcharacteristheexistenceofthresholdpressuregradient[1],inotherwords:thefluidcanflow,onlywhentherealpressuregradientislargerthanthresholdpressuregradient.Thedevelopmentinthisresearchingfieldhasn'tmuchprogressbecauseofquestionsineconomicsandteelmologysincethisconceptwasformedin1951.ThisfieldisbecominganewwarmresearchingPOintasChinadevelopedthelowPermeabilityreservoirsinlastfewyears.Thego…     

The initial-value problem for elastodynamics of incompressible bodies

We prove short-time well-posedness of the Cauchy problem for incompressible strongly elliptic hyperelastic materials. Our method consists in:

1. Reformulating the classical equations in order to solve for the pressure gradient (The pressure is the Lagrange multiplier corresponding to the constraint of incompressibility.) This formulation uses both spatial and material variables.
2. Solving the reformulated equations by using techniques which are common for symmetric hyperbolic systems. These are:

1. Using energy estimates to bound the growth of various Sobolev norms of solutions.
2. Finding the solution as the limit of a sequence of solutions of linearized problems.

Our equations differ from hyperbolic systems, however, in that the pressure gradient is a spatially non-local function of the position and velocity variables.     

Unsteady flow of viscoelastic fluid with fractional Maxwell model in a channel

The unsteady flow of viscoelastic fluid with the fractional derivative Maxwell model (FDMM) in a channel is studied in this note. The exact solutions are obtained for an arbitrary pressure gradient by means of the finite Fourier cosine transform and the Laplace transform. Two special cases of pressure gradient are discussed. Some results given by the classical models with integer-order are included in this note.     

Response of a skewed turbulent boundary layer to favourable pressure gradient

Experimental results are reported for the response to a favourable pressure gradient of an initially turbulent boundary layer (Re _θ?≈?1600) developing on a flat plate with its leading edge skewed at 60° to the approach flow. The pressure gradient orthogonal to the leading edge is nominally the same as that which was shown by Escudier et?al. [(1998) Exp Fluids 25: 491–502] to cause extreme thinning of a two-dimensional (2D) (i.e. unskewed) turbulent boundary layer and the intermittency in the immediate vicinity of the surface to fall to zero, i.e. an apparent laminarisation of the boundary layer. In the case of the skewed boundary layer, the responses of the turbulence and mean-flow structures are qualitatively similar to those for the 2D situation. However, the streamwise pressure gradient is much weaker than for the 2D experiment and the extent of the changes it produces is much reduced. Even so, the changes are considerably greater than would be expected from the magnitude of the streamwise pressure gradient.     

Exact Solution of the Problem on a Six‐Constant Jeffreys Model of Fluid in a Plane Channel

An exact analytic solution of the problem of a generalized viscoelastic Jeffreys fluid flow in a plane channel under the action of a pressure gradient is found. The velocity profiles are obtained in a parametric form with a velocity gradient taken as a parameter. The critical values of the pressure gradient are determined, which, when exceeded, lead to weak tangential discontinuities in the longitudinal velocity profile. When the pressure gradient changes smoothly over some range of parameters, a hysteresis loop emerges on the graph of the flow rate versus the pressure gradient.     

Analytical solutions of the boundary-value problem of nonstationary flow of viscoplastic medium between two plates

Summary Nonstationary flow of a viscoplastic medium between two parallel plates is considered for the case of a varying pressure gradient. The problem is reduced to the Stephan problem, with the condition on the boundary separating the flow domain from the quasi-rigid domain. Four multiparameter families of exact solutions are found. The first family describes the flow decelerations up to a full stop. The second family determines the development of the flow from the state of rest as the pressure gradient increases. The third family describes the development of the flow for the case where (1) the pressure gradient is constant and exceeds the threshold value related to the yield stress, (2) the upper plate does not move, and (3) the lower plate moves with a constant acceleration. Finally, the fourth family determines the flow retardation, when the pressure gradient is constant and is less than the threshold value. The decrease in the flow of the viscoplastic medium can be achieved for certain values of parameters by increasing the quasi-rigid domain, whereas the viscoplastic flow remains unchanged. Received 7 October 1998; accepted for publication 8 April 1999     

一种力-电协同驱动的细胞微流控培养腔理论模型

细胞培养液在微流控生物反应器中受到外界物理场（如压力梯度或者电场）作用流动而产生流体剪应力,并进一步刺激种子细胞调控其内部基因的表达,从而促进细胞的分化和生长,这个过程在自然生命组织内的微管中亦是如此。考虑到细胞培养微腔隙中液体流动行为很难实验量化测定,理论建模分析是目前可行的研究手段。因此建立了矩形截面的细胞微流控培养腔理论模型,将外部的物理驱动场（压力梯度与电场）与培养腔内液体的流速、切应力和流率联系起来,分别得到了压力梯度驱动（Pressure gradient driven,PGD）、电场驱动（Electric field driven,EFD）及力-电协同驱动（Pressure-electricity synergic driven,P-ESD）三种驱动方式下的液体流动理论模型。结果表明该理论模型与现有的实验结果基本一致,具体地:力-电协同作用下的解答为压力梯度驱动和电场驱动结果的叠加。细胞培养腔内的流体流速、剪应力及流率幅值均正比于外部物理场强幅值,但随着压力梯度驱动载荷频率的增大而减小,随着电场驱动频率的变化不明显。在压力梯度驱动作用下,细胞贴壁处的切应力随着腔高的增大而线性增大,流率则随着腔高的增大而非线性增大,而电场驱动下的结果不受腔高的影响。生理范围内的温度场变化对压力和电场驱动的结果影响不大。另外,在引起细胞响应的流体切应力水平,电场驱动能提供较大的切应力幅值而压力梯度驱动则能提供较大的流率幅值。该理论模型的建立为细胞微流控生物反应器实验系统的设计及参数优化提供理论参考,同时也为力-电刺激细胞生长、分化机理的研究的提供基础。      

非达西渗流拟启动压力梯度推算

岩体启动压力梯度的大小是石油开采,工程防渗处理等方面需要考虑的影响因素。非线性渗流现象已被广泛所认同,低渗透性岩渗流存在启动压力梯度,它基本代表了流体产生渗流时的压力梯度大小。试验对试样进水端水压进行长期稳定控制,出水端采用精确测定流出水体积变化量的方法进行流量测定; 当渗出端水体积的变化量与时间呈直线关系时,即认为该时段渗透特征符合达西定律。通过对具有不同渗透性的软弱岩系统的试验测试,根据实测渗透系数与压力梯度的正相关关系,提出了拟启动压力梯度推算的"三次平均法",即:首先统一调整实测渗透系数值,然后确定绝对渗透系数值,再根据该值确定V0值, 3次取平均推算拟启动压力梯度的方法。分析表明:拟启动压力梯度随围压升高而增大,随岩石的绝对渗透系数降低而增大; 在渗透压差增大过程测得的拟启动压力梯度大于降低过程的拟启动压力梯度,这些都与理论分析结果相吻合。实例说明拟启动压力梯度推算"三次平均法"较为合理,易于推广使用。这为通过室内试验确定工程岩体启动渗流压力梯度的大小,提出了一种新的思路。     

LES of the adverse-pressure gradient turbulent boundary layer

We describe large-eddy simulations (LES) of the flat-plate turbulent boundary layer in the presence of an adverse pressure gradient. The stretched-vortex subgrid-scale model is used in the domain of the flow coupled to a wall model that explicitly accounts for the presence of a finite pressure gradient. The LES are designed to match recent experiments conducted at the University of Melbourne wind tunnel where a plate section with zero pressure gradient is followed by section with constant adverse pressure gradient. First, LES are described at Reynolds numbers based on the local free-stream velocity and the local momentum thickness in the range 6560–13,900 chosen to match the experimental conditions. This is followed by a discussion of further LES at Reynolds numbers at approximately 10 times and 100 times these values, which are well out of range of present day direct numerical simulation and wall-resolved LES. For the lower Reynolds number runs, mean velocity profiles, one-point turbulent statistics of the velocity fluctuations, skin friction and the Clauser and acceleration parameters along the streamwise, adverse pressure-gradient domain are compared to the experimental measurements. For the full range of LES, the relationship of the skin-friction coefficient, in the form of the ratio of the local free-stream velocity to the local friction velocity, to both Reynolds number and the Clauser parameter is explored. At large Reynolds numbers, a region of collapse is found that is well described by a simple log-like empirical relationship over two orders of magnitude. This is expected to be useful for constant adverse-pressure gradient flows. It is concluded that the present adverse pressure gradient boundary layers are far from an equilibrium state.     

MLPG_R method for modelling 2D flows of two immiscible fluids

This is a first attempt to develop the Meshless Local Petrov–Galerkin method with Rankine source solution (MLPG_R method) to simulate multiphase flows. In this paper, we do not only further develop the MLPG_R method to model two‐phase flows but also propose two new techniques to tackle the associated challenges. The first technique is to form an equation for pressure on the explicitly identified interface between different phases by considering the continuity of the pressure and the discontinuity of the pressure gradient (i.e. the ratio of pressure gradient to fluid density), the latter reflecting the fact that the normal velocity is continuous across the interface. The second technique is about solving the algebraic equation for pressure, which gives reasonable solution not only for the cases with low density ratio but also for the cases with very high density ratio, such as more than 1000. The numerical tests show that the results of the newly developed two‐phase MLPG_R method agree well with analytical solutions and experimental data in the cases studied. The numerical results also demonstrate that the newly developed method has a second‐order convergent rate in the cases for sloshing motion with small amplitudes. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of Fluid Flow on Fractures and Implications for Reactive Transport Simulations

The diverse numerical simulation techniques employed to predict fluid flow properties of fractures yield differing results which limits their applicability for reactive transport simulations. Basically the fluid flow simulation techniques can be divided in two groups: (i) techniques that yield average fluid flow characteristics and (ii) techniques that produce space-resolved properties. These differences may have substantial impact on the reactive transport simulations but may also depend on the fracture characteristics. For this reason, a sensitivity analysis of the geometrical properties of fractures on the fluid flow properties is conducted and evaluated with respect to their impact on reactive transport modeling. Although employing space-resolved simulation techniques, the results of the tests show average values for permeability and fluid velocity that are comparable to previous studies that used other simulation techniques. Observed fluid flow channeling appears to be related to fracture surfaces matching and anisotropy. However, average flow velocities at potential sites for reactive transport differ up to a factor of five from the average ones for the entire fracture. Furthermore, extreme values at reactive transport sites may differ even more and the flow may be directed against the applied pressure gradient. For studies concerned with simulation of reactive transport, these deviations are crucial and should be explicitly considered in the calculations. Hence space-resolved fluid flow simulations should be employed for the simulation of reactive transport.     

压力梯度对湍流逆梯度输运的影响

研究Birkhoff系统的一般Lie对称性导致的非Noether守恒量. 得到非Noether守恒 量的存在定理，举例说明结果的应用.     

Nodal operator splitting adaptive finite element algorithms for the Navier–Stokes equations

Solution algorithms for solving the Navier–Stokes equations without storing equation matrices are developed. The algorithms operate on a nodal basis, where the finite element information is stored as the co-ordinates of the nodes and the nodes in each element. Temporary storage is needed, such as the search vectors, correction vectors and right hand side vectors in the conjugate gradient algorithms which are limited to one-dimensional vectors. The nodal solution algorithms consist of splitting the Navier–Stokes equations into equation systems which are solved sequencially. In the pressure split algorithm, the velocities are found from the diffusion–convection equation and the pressure is computed from these velocities. The computed velocities are then corrected with the pressure gradient. In the velocity–pressure split algorithm, a velocity approximation is first found from the diffusion equation. This velocity is corrected by solving the convection equation. The pressure is then found from these velocities. Finally, the velocities are corrected by the pressure gradient. The nodal algorithms are compared by solving the original Navier–Stokes equations. The pressure split and velocity–pressure split equation systems are solved using ILU preconditioned conjugate gradient methods where the equation matrices are stored, and by using diagonal preconditioned conjugate gradient methods without storing the equation matrices. © 1998 John Wiley & Sons, Ltd.