标准线性固体土模型条件下桩的动力响应分析与应用

在桩侧土为标准线性固体模型条件下，建立了成层土中考虑桩端和桩侧土作用的有限长桩，桩顶受纵向激励情况下的定解问题，用拉氏变化(Laplace)求得桩顶阻抗函数和速度频率响应的函数，然后利用卷积定理和傅立叶(Fourier)逆交换求得任意激振力作用下桩顶速度响应的半解析解，并研究了模型中各个主要参数对桩顶响应及速度响应曲线的影响，并分析了缺陷桩和不同土参数对曲线的影响，得到了许多新的有意义的结论。     

Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.     

Analysis of a turbulent buoyant confined jet modeled using realizable k–ɛ model

Through this paper, analyses of components of the unheated/heated turbulent confined jet are introduced and some models to describe them are developed. Turbulence realizable k–? model is used to model the turbulence of this problem. Numerical simulations of 2D axisymmetric vertical hot water confined jet into a cylindrical tank have been done. Solutions are obtained for unsteady flow while velocity, pressure, temperature and turbulence distributions inside the water tank are analyzed. For seeking verification, an experiment was conducted for measuring of the temperature of the same system, and comparison between the measured and simulated temperature shows a good agreement. Using the simulated results, some models are developed to describe axial velocity, centerline velocity, radial velocity, dynamic pressure, mass flux, momentum flux and buoyancy flux for both unheated (non-buoyant) and heated (buoyant) jet. Finally, the dynamics of the heated jet in terms of the plume function which is a universal quantity and the source parameter are studied and therefore the maximum velocity can be predicted theoretically.     

A structural model for the vortex tubes of isotropic turbulence

A class of exact solutions of the Navier–Stokes equations is introduced to model the fine-scale, tubular structures of isotropic turbulence. The model vortices exhibit slow algebraic fall-off of the induced velocity, and accurately reproduce the velocity signatures observed in DNS and experiments. The proposed model has interesting implications for the theoretical analysis of turbulence, supporting the view that the inertial range energy scaling may have a link with the near-singular velocity field induced by vortex tubes produced by the roll-up of vortex sheets.      

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.     

Parameter estimation of engineering turbulence model

A parameter estimation algorithm is introduced and used to determine the parameters in the standardk-∈ two equation turbulence model (SKE). It can be found from the estimation results that although the parameter estimation method is an effective method to determine model parameters, it is difficult to obtain a set of parameters for SKE to suit all kinds of separated flow and a modification of the turbulence model structure should be considered. So, a new nonlineark-∈ two-equation model (NNKE) is put forward in this paper and the corresponding parameter estimation technique is applied to determine the model parameters. By implementing the NNKE to solve some engineering turbulent flows, it is shown that NNKE is more accurate and versatile than SKE. Thus, the success of NNKE implies that the parameter estimation technique may have a bright prospect in engineering turbulence model research.     

THE RENORMALIZATION-GROUP ANALYSIS FOR THE STATISTICAL PROPERTIES OF RAPIDLY ROTATING TURBULENCE

利用重正化群方法对强旋转湍流场统计性质予以研究, 通过重正化微扰展开, 对高波数速度分量进行逐 阶平均.计算结果显示当旋转角速度Ω → ∞时, 用以表征高波数速度分量对低波数速度分量影响的重正化黏性将趋于0, 这表明在强旋转条件下科氏力将抑制湍流速度分量之间的非线性相互作用, 从而阻碍湍流的能量级串效应, 当Ω → ∞时湍流的能量级串效应消失, 导致湍流脉动消失, 流动将层流化.理论计算结果还显示对于强旋转湍流, 时域-空域联立Fourier的湍流速度分量存在二维化趋势, 球面平均能谱函数有标度关系E(k) ∝ k^-3.     

Spatial structures and scaling in the Convective Boundary Layer

We performed an investigation on spatial features of the Convective Boundary Layer (CBL) of the atmosphere, which was simulated in a laboratory model and analyzed by means of image analysis techniques. This flow is dominated by large, anisotropic vortical structures, whose spatial organization affects the scalar transport and therefore the fluxes across the boundary layer. With the aim of investigating the spatial structure and scaling in the Convective Boundary Layer, two-dimensional velocity fields were measured, on a vertical plane, by means of a pyramidal Lucas–Kanade algorithm. The coherent structures characterizing the turbulent convection were educed by analyzing the Finite-Time Lyapunov Exponent fields, which also revealed interesting phenomenological features linked to the mixing processes occurring in the Convective Boundary Layer. Both velocity and vorticity fields were analyzed in a scale-invariance framework. Data analysis showed that normalized probability distribution functions for velocity differences are dependent on the scale and tend to become Gaussian for large separations. Extended Self Similarity holds true for velocity structure functions computed within the mixing layer, and their scaling exponents are interpreted well in the phenomenological framework of the Hierarchical Structure Model. Specifically, β parameter, which is related to the similarity between weak and strong vortices, reveals a higher degree of intermittency for the vertical velocity component with respect to the horizontal one. On the other hand, the analysis of circulation structure functions shows that scaling exponents are fairly constant in the lowest part of the mixed layer, and their values are in agreement with those reported in Benzi et al. (Phys Rev E 55:3739–3742, 1997) for shear turbulence. Moreover, the relationship between circulation and velocity scaling exponents is analyzed, and it is found to be linear in the bottom part of the mixing layer. The investigation of the CBL spatial features, which has seldom been studied experimentally, has important implications for the comprehension of the mixing dynamics, as well as in turbulence closure models.     

Measurements of ensemble averaged flame dynamics using spatially resolved analysis

When applying flame sheet models to predict the dynamics of turbulent flames, it is common to model turbulence using ensemble averaging of the velocity. Measurements of the flame dynamics were made to support use this type of methodology, by measuring the dynamic volume of the flame using phase averaged images of the CH^∗ chemiluminescence. The dynamics agreed with the common behavior described in the literature, namely frequency scaling according to Strouhal number based on flow convective timescales. However, slightly different timescales were observed for the response magnitude and phase, indicating the possibility of different scaling mechanisms at work between these phenomena. The flame heat release rate dynamics were found to be identical to the dynamic response of the flame volume to inlet velocity perturbations, suggesting a simple proportionality between heat release rate and the flame volume. This result supports the use of ensemble averaging for modeling of the turbulent velocity for predictions of flame dynamics.     

Effect of Lorentz forces on forced-convection nanofluid flow over a stretched surface

Magnetic nanofluid hydrothermal analysis over a plate is studied that includes consideration of thermal radiation. The Runge–Kutta (RK4) method is utilized to get solution of ODEs which are obtained from similarity solution. In considering the impacts of Brownian motion, we applied Koo–Kleinstreuer–Li correlation to simulate the properties of CuO–water. The influence is discussed of important parameters such as the temperature index, magnetic, radiation, and velocity ratio parameters and volume fraction of nanoparticle on hydrothermal behavior. Results illustrate that the coefficient of skin friction enhances with enhancing magnetic parameter while reduces with enhancing velocity ratio parameter. Also the Nusselt number was found to directly depend on the velocity ratio and temperature index parameters but has an inverse dependence on the magnetic and radiation parameters.     

用子波变换检测壁湍流信号局部标度指数

本文用子波变换检测了刻画壁湍流脉动信号自相似性的局部标度指数，研究了不同尺度的湍流结构的自相似性，发现在湍流边界层猝发过程中，喷射和扫掠发生时刻小尺度脉动速度信号的局部标度指数为负值，说明在大尺度猝发事件发生的时刻小尺度结构具有奇异的自相似性，在猝发过程中其作用不仅仅是对湍能的耗散．     

Numerical Simulation of Evolution of Helical-Vortex Instability in a Fluid Layer Heated from Below

Large-scale helical-vortex instability is simulated numerically for laminar convection regimes. The simulation is based on the nonlinear Boussinesq equations supplemented by an external force of special form, whose structure is identical to the tensor structure of the generative term in the equation for the hydrodynamic alpha-effect. Introducing the external force makes it possible to model the average effect produced by small-scale helical turbulence: generation of a positive feedback between the toroidal and poloidal components of the velocity vector field. The structure and energetics of helical-vortex flows developing in a convectively unstable medium are analyzed. The parameters of helical and non-helical convection regimes are compared. Certain novel effects related to spiral feedback are discussed: generation of a toroidal velocity field in structures with the poloidal circulation typical of ordinary convection; enlargement of the horizontal scale of supercritical motions as a result of the merging of vortex cells, accompanied by a significant increase in the circulation strength in the larger helical vortices formed; qualitative changes in the energetics.     

The analysis of the channel flow sensitivity to the parameters of the k–ε method

This paper deals with the problem of using sensitivity analysis for fluid mechanics solutions to the constants of the standard k–ε method for 2D, incompressible and steady flows. The problem is described and analysed on the basis of a channel flow. Sensitivity coefficients of the following properties were determined: a pressure, two components of a velocity, a turbulence kinetic energy, a dissipation rate of turbulence kinetic energy and a turbulence dynamic viscosity. The calculated property values depend on five model constants that are parameters of the sensitivity analysis in this paper. Sensitivity coefficients are derivatives of the above properties, for individual parameters. In this paper these coefficients are determined using a finite difference approximation to the sensitivities coefficients. The author of this paper compares three models of the boundary layer with regard to the sensitivity of properties to the parameters. Irrespective of the boundary layer model used here, the analysis of sensitivity coefficients for the channel flow properties shows that the most sensitive property is the turbulence dissipation rate. Next properties of consequence, although of significantly smaller values of sensitivity coefficients, are the turbulence viscosity and the turbulence kinetic energy. All flow properties are mostly sensitive to the C_µ parameter. One of the final conclusions in this paper is that the analysis of sensitivity coefficient fields allows the reliable checking of results and indicates those areas most prone to calculation difficulties. Copyright © 2008 John Wiley & Sons, Ltd.     

Bifurcation response and Melnikov chaos in the dynamic of a Bose–Einstein condensate loaded into a moving optical lattice

We investigate global bifurcation of a Bose–Einstein condensate with both repulsive two-body interaction between atoms and attractive three-body interaction loaded into a traveling optical lattice. Slow-flow equations of the traveling wave function are the first to derive and the reduced amplitude equation is obtained. The Melnikov method is applied on the reduced parametrically driven system and the Melnikov function is subsequently established. Effects of different physical parameters on the global bifurcation are studied analytically and numerically, and different chaotic regions of the parameter space are found. The results suggest that optical intensity may help to enhance chaos while the strength of the effective three-body interaction, the velocity of the optical lattice, and the damping coefficients annihilate or reduce chaotic behavior of the steady-state traveling wave solution of the particle number density of a Bose–Einstein condensate.     

Simulation of the Effect of the Freestream Turbulence Parameters on Heat Transfer in an Unsteady Boundary Layer

A variant of the two-parameter turbulence model which makes it possible continuously to calculate a flow region with laminar, transition and turbulent regimes is proposed for investigating the flow under conditions of high freestream turbulence intensity. It is shown that the properties of the thermal transition can be theoretically described using the quasi-steady turbulence model in the case of periodic freestream velocity distribution. The numerical results are compared with theoretical and experimental data. The approach proposed is developed for determining the combined effect of the parameters of harmonic fluctuations of the external velocity and freestream turbulence on the heat transfer characteristics on a flat plate with different boundary conditions for the enthalpy.     

MHD mixed convection stagnation-point flow of a nanofluid over a vertical permeable surface: a comprehensive report of dual solutions

The steady laminar magnetohydrodynamic mixed convection boundary layer flow of a nanofluid near the stagnation-point on a vertical permeable plate with prescribed external flow and surface temperature is investigated in this study. Here, both assisting and opposing flows are considered and studied. Using appropriate similarity variables, the governing equations are transformed into nonlinear ordinary differential equations in the dimensionless stream function, which is solved numerically using the Runge–Kutta scheme coupled with a conventional shooting procedure. Three different types of nanoparticles, namely copper Cu, alumina Al₂O₃ and titania TiO₂ with water as the base fluid are considered. Numerical results are obtained for the skin-friction coefficient and Nusselt number as well as for the velocity and temperature profiles for some values of the governing parameters, namely, the volume fraction of nanoparticles ?, permeability parameter f _o , magnetic parameter M and mixed convection parameter λ. It is found that dual solutions exist for both assisting and opposing flows, and the range of the mixed convection parameter for which the solution exists, increases with suction, magnetic field and volume fraction of nanoparticles.     

Spontaneous rotation in the exact solution of magnetohydrodynamic equations for flow between two stationary impermeable disks

Magnetohydrodynamic (MHD) flow of a viscous electrically conducting incompressible fluid between two stationary impermeable disks is considered. A homogeneous electric current density vector normal to the surface is specified on the upper disk, and the lower disk is nonconducting. The exact von Karman solution of the complete system of MHD equations is studied in which the axial velocity and the magnetic field depend only on the axial coordinate. The problem contains two dimensionless parameters: the electric current density on the upper plate Y and the Batchelor number (magnetic Prandtl number). It is assumed that there is no external source that produces an axial magnetic field. The problem is solved for a Batchelor number of 0–2. Fluid flow is caused by the electric current. It is shown that for small values of Y, the fluid velocity vector has only axial and radial components. The velocity of motion increases with increasing Y, and at a critical value of Y, there is a bifurcation of the new steady flow regime with fluid rotation, while the flow without rotation becomes unstable. A feature of the obtained new exact solution is the absence of an axial magnetic field necessary for the occurrence of an azimuthal component of the ponderomotive force, as is the case in the MHD dynamo. A new mechanism for the bifurcation of rotation in MHD flow is found.     

Investigation of the porous drag and permeability at the porous-fluid interface: Influence of the filtering parameters on Darcy closure

The porous-fluid interface encompasses a region bridging the flow inside a porous medium and a free-flowing fluid. In the context of volume-averaged simulations, it can be described by a set of gradually changing parameters defining the porous medium, mainly porosity and permeability. In this paper, both the permeability and the porous-induced drag force are evaluated a-priori, by explicitly filtering a set of Particle-Resolved Simulations (PRS) of the flow in the channel partially occupied by the porous medium. Different porous matrices are considered and the influence of the geometry and filtering parameters on the macroscopic quantities is studied. Especially, the focus is placed on the requirements for the kernel type and size to perform filtering accurately, and their impact on the distribution of permeability at the interface. The performance of the typically used models for the permeability is compared to the explicitly filtered results. Lastly, a new model for permeability and the drag force is introduced, taking into account the information about the filtering size and non-uniformity of the velocity field. The model greatly improves the prediction of velocity at the porous-fluid interface and serves as a proof of concept that a successful porous drag model should strive to include information about both parameters.     

Turbulence enhancement of stagnation point heat transfer on a body of revolution

Enhancement, by free stream turbulence, of convective heat transfer to the stagnation region of a hemispherical-nosed cylinder has been studied. Increases in heat transfer were found to depend primarily on the Reynolds number and turbulence intensity of the free stream, experimental results being most successfully correlated on a NuRe^?0.5 versus TuRe^0.5 basis. Flow visualization studies have demonstrated the validity of a phenomenological model of the enhancement process, predictions of this theory showing good agreement with experimental results. The effect of free stream turbulence on the stagnation point velocity gradient has also been evaluated.