A vertically moving grid finite‐element modelling of tidal flow in the Changjiang Estuary,China

An estuarine two‐dimensional vertical finite‐element model of tidal flow has been established by laterally integrating Navier–Stokes equation. To this end, a moving grid finite‐element method has been used. An arbitrarily shaped quadrilateral element has been selected. This model has been validated by using field data from two monitoring stations at the North Passage of the Changjiang Estuary. Using this numerical model, two types of modelled results were obtained: (1) vertical distributions of tidal current velocities at the North Passage of the Changjiang Estuary; (2) longitudinal distributions of tidal current velocities at maximum flood tide, at high slack water, at maximum ebb tide and at low slack water tide at the North Passage of the Changjiang Estuary. The conclusion is that the model provides a reasonable agreement with observed data. Copyright © 2003 John Wiley & Sons, Ltd.     

Correlations for fire-wind enhancement flow characteristics based on LES simulations

Unraveling the physics of fire-wind interaction has long been a subject of interest. Among all the physics involved, enhancement of wind by fire deserves great attention due to its potential effects on building structures downstream of the fire source in bushfire attack events. Predominantly, two contributing factors determine the extent to which wind is enhanced by fire: freestream wind velocity and fire intensity. This study employs Large-Eddy Simulation (LES) to fundamentally investigate the combined effects of freestream wind velocity and fire intensity on fire-wind enhancement. An added module was implemented to an open-source transient fire solver in order to analyze the effects of freestream wind velocity and fire intensity based on the analysis of interactions between momentum and fire-induced buoyancy forces. Simulations are performed for parametric combinations of wind velocity and fire intensity. The LES results demonstrate that the normalized maximum wind enhancement increases with a reduction of freestream wind velocity and an increase in fire intensity. The non-dimensional Froude number, Fr, and normalized fire intensity, I*, were employed to quantify the effects of freestream wind velocity and fire intensity, respectively. A correlation was developed to determine the maximum wind enhancement as a function of Fr and I*. The location corresponding to maximum wind enhancement occurs further downstream of the fire source as freestream wind velocity or fire intensity increases. A correlation based on the Fr number and I* was developed for the location at which maximum wind enhancement occurs. Furthermore, the concept of wind enhancement plume line was defined as a line along which the local wind enhancement occurs at a given longitudinal location downstream of the fire source, for which a correlation was also developed. Moreover, a gradual decaying trend is observed in wind enhancement after reaching a peak along the wind enhancement plume line in all simulation scenarios for which a correlation was also developed as a function normalized longitudinal direction.     

Steam plume length diagram for direct contact condensation of steam injected into water

Direct Contact Condensation (DCC) of steam in water occurs when steam is introduced into water. It is a phenomenon of high importance in many industrial applications. An important feature of the DCC is the length of the steam plume. Correlations for a steam plume length presently available are accurate only for limited conditions. In this paper, a new two-dimensional steam plume length diagram is presented capable of predicting length accurately for a wide range of conditions. The diagram is validated against experiments. Furthermore, corrections necessary to adopt the diagram for steam injection into a water flow are discussed in the paper.     

Influences of the fish-mouth project and the groins on the flow and sediment ratio of the Yangtze River waterway

A depth-integrated two-dimensional numerical model of current, salinity andsediment transport was proposed and calibrated by the observation data in the Yangtze RiverEstuary. It was then applied to investigate the flow and sediment ratio of the navigationchannel, i. e. the North Channel of the Yangtze River Estuary, before and after the firstphase waterway project is implemented. Particularly, the influences of groin length and theorientation of the submerged dam on the flow ratio and sediment load discharging into theNorth Channel were discussed. The numerical results demonstrate that less sediment loaddischarges into the navigation channel, which unburdens the waterway dredging, but in themeantime the flow ratio is also decreased. The flow and sediment ratio can be adjusted bychanging layout and dimensions of the hydro-structures, such as the groin length, the topheight, etc. The effect of the orientation of the submerged dam is more obvious than thegroin lengh.     

Experimental study on the condensation of supersonic steam jet submerged in quiescent subcooled water: Steam plume shape and heat transfer

The condensation of supersonic steam jet submerged in the quiescent subcooled water was investigated experimentally. The results indicated that the shape of steam plume was controlled by the steam exit pressure and water temperature. Six different shapes of steam plume were observed under the present test conditions. Their distribution as a function of the steam exit pressures and water temperatures was given. As the steam mass velocity and water temperature increase, the measured maximum expansion ratio and dimensionless penetration length of steam plume were in the ranges of 1.08–1.95 and 3.05–13.15, respectively. The average heat transfer coefficient of supersonic steam jet condensation was found to be in the range of 0.63–3.44 MW/m²K. An analytical model of steam plume was found and the correlations to predict the maximum expansion ratio, dimensionless penetration length and average heat transfer coefficient were also investigated.     

Numerical Simulation of Two Phase Turbulent Flow of Nanofluids in Confined Slot Impinging Jet

In this article, a numerical investigation is performed on flow and heat transfer of confined impinging slot jet, with a mixture of water and Al₂O₃ nanoparticles as the working fluid. Two-dimensional turbulent flow is considered and a constant temperature is applied on the impingement surface. The k ? ω turbulence model is used for the turbulence computations. Two-phase mixture model is implemented to study such a flow field. The governing equations are solved using the finite volume method. In order to consider the effect of obstacle angle on temperature fields in the channel, the numerical simulations were performed for different obstacle angles of 0° ? 60°. Also different geometrical parameters, volume fractions and Reynolds numbers have been considered to study the behavior of the system in terms of stagnation point, average and local Nusselt number and stream function contours. The results showed that the intensity and size of the vortex structures depend on jet- impingement surface distance ratio (H/W) and volume fraction. The maximum Nusselt number occurs at the stagnation point with the highest values at about H/W = 1. Increasing obstacle angle, from 15° to 60°, enhances the heat transfer rate. It was also revealed that the minimum value of average Nusselt number occurs in higher H/W ratios with decreasing the channel length.     

A note on mixed convection of laminar plane water plumes at temperatures around the density extremum

In all studies concerning mixed convection in plane laminar plumes a linear relationship between fluid density and temperature has been used. However, it is known that the water density-temperature relationship is non-linear at low temperatures with a density maximum at 3.98°C for pure water. In this note the problem of plane laminar water plume in a coflowing vertical free stream has been investigated taking into account the non-linearity between density and temperature. This is the first work in the literature which treats plane mixed convection plumes with nonlinear relation between density and temperature. Both rising and descending plumes have been investigated. It was found that the ambient water temperature plays an important role on the results. When the ambient temperature is greater than maximum density temperature (T _a > T _m), the water plume behavior is similar to that of the classical plume with linear density-temperature relationship. However, when the ambient temperature is equal or lower than the maximum density temperature the water plume behavior is completely different from the classical plume with linear density-temperature relationship. The centerline velocity shows a series of maxima and minima which are produced by the combination of the nonlinear density-temperature relation and the free stream.     

Analytical and experimental investigations of the pulsed air–water jet

This paper describes a new way of generating pulsed air–water jet by entraining and mixing air into the cavity of a pulsed water jet nozzle. Based on the theory of hydro-acoustics and fluid dynamics, a theoretical model which describes the frequency characteristic of the pulsed air–water jet is outlined aimed at gaining a better understanding of this nozzle for generating pulses. The calculated result indicates that as the air hold-up increases, the jet oscillation frequency has an abrupt decrease firstly, and then reaches a minimum gradually at α (air hold-up)=0.5, finally it gets increased slightly. Furthermore, a vibration test was conducted to validate the present theoretical result. By this way, the jet oscillation frequency can be obtained by analyzing the vibration acceleration of the equal strength beam affected by the jet impinging. Thereby, it is found that the experimental result shows similar trend with the prediction of the present model. Also, the relationship between vibration acceleration and cavity length for the pulsed water jet follows a similar tendency in accord with the pulsed air–water jet, i.e. there exists a maximum for each curve and the maximum occurs at the ratio of L/d₁ (the ratio of cavity length and upstream nozzle diameter) =2.5 and 2.2, respectively. In addition, experimental results on specimens impinged by the pulsed water jet and pulsed air–water jet show that the erosion depth increases slightly with air addition within a certain range of cavity length. Further, this behavior is very close to the vibration test results. As for erosion volume, the air entrained into the cavity significantly affects the material removal rate.     

Numerical simulation of secondary vortex chamber effect on the cooling capacity enhancement of vortex tube

A vortex tube with additional chamber is investigated by computational fluid mechanics techniques to realize the effects of additional chamber in Ranque–Hilsch vortex tube and to understand optimal length for placing the second chamber in order to have maximum cooling effect. Results show that by increasing the distance between two chambers, both minimum cold and maximum hot temperatures increase and maximum cooling effect occurs at Z/L = 0.047 (dimensionless distance).     

Buoyant plume above a line source of heat in an anisotropic porous medium

Buoyancy-induced convection arising from a horizontal line heat source embedded in an anisotropic porous medium is investigated analytically. The porous medium is anisotropic is permeability with its principal axes oriented in a direction that is oblique to the gravity vector. Assuming the boundary layer approximation, closed-form exact similarity solutions for both flow and temperature fields are presented and compared with those of isotropic case. Scale analysis is applied to predict the order of magnitudes involved in the boundary layer regime for which the conditions of validity are obtained. Effects of both anisotropic parameters (K* and %) and Rayleigh number Ra_L are observed to be strongly significant. It is demonstrated that a minimum (maximum) intensity of the thermal convective plume above the line source of heat can be obtained if the porous matrix is oriented with its principal axis with higher permeability parallel (perpendicular) to the vertical direction.     

Turbulent Transition in the Inlet Region of a Circular Pipe

Calculated and experimental data on turbulent transition in a circular pipe are analyzed. The calculations were performed using the three-parameter turbulence model. The dependence of the distance from the inlet to the point of minimum friction during transition on the Reynolds number for fixed inlet conditions and the distribution of the turbulence parameters over the pipe length and radius are obtained. The dependence of the maximum (critical) Reynolds number, Re^*, for which there is no transition in the pipe, on the inlet intensity and scale of turbulence is found. It is suggested that Re^* depends on the inlet perturbation parameters up to Re^* = 1000, where the friction coefficients for laminar and turbulent flows coincide.     

Three-dimensional green water velocity on a model structure

Flow kinematics of green water due to plunging breaking waves impinging on a simplified, 3D model structure was investigated in the laboratory. Two breaking wave conditions were tested: one with waves impinging on the vertical wall of the model at still water level, and the other with waves impinging on the horizontal deck surface. The bubble image velocimetry (BIV) technique was used to measure flow velocities. Measurements were taken on both vertical and horizontal planes. Evolution of green water flow kinematics in time and space was revealed and was found to be quite different between the two wave conditions, even though the incoming waves are essentially identical. The time history of maximum velocity is demonstrated and compared. In both cases, the maximum velocity occurs near the green water front and beneath the free surface. The maximum horizontal velocity for the deck impinging case is 1.44C with C being the wave phase speed, which is greater than 1.24C for the wall impingement case. The overall turbulence level is about 0.3 of the corresponding maximum velocity in each wave condition. The results were also compared with 2D experimental results to examine the 3D effect. It was found that the magnitude of the maximum vertical velocity during the runup process is 1.7C in the 3D model study and 2.9C in the 2D model study, whereas the maximum horizontal velocity on the deck is similar, 1.2C in both 3D and 2D model studies.     

A method to compute ship exhaust plumes with waves and wind

Thermal and concentration transport models are implemented in CFDShip‐Iowa version 4.5, a semi‐coupled solver for air/water free surface flow (Int. J. Numer. Meth. Fluids 2008; 58 (6):591–624), to investigate the exhaust plume around ship superstructures. An incompressible, variable density approximation is implemented where the density can change in all governing equations due to temperature variations only. The thermal and concentration models are tested for the cases of steady and unsteady flow with thermal and solution transport in a 2D square cavity, and for a 3D thermal plume in an open environment, showing good agreement between computational results and experimental data. To test the method in an extreme motions condition, the exhaust plume of the ONR Tumblehome model DTMB 5613 is studied, showing complicated vortical structures in air including a pair of counter‐rotating vortices downstream of the stack for cross‐flow, and bended bird‐plume shape in the symmetry plane and varying arc‐shape in axial sections both for temperature and NO_x concentration fields. Effects of smoke exhaust speed and wind speed on the temperature and concentration distributions are studied. Finally, a smoke downwash computation is performed for a ship free to move in 6 degrees of freedom in a sea state 8 condition. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of light transmission and reflection techniques to determine concentrations in flow tank experiments

Transmissive and reflective intensity measurements for visual concentration determinations in 2D flow tank experiments were compared and evaluated for their applicability in the study of flow and transport phenomena. A density-dependent heterogeneous flow experiment was conducted and transmission and reflection images of the dyed saltwater plume were analyzed. A single light source and dark curtains forced the light to pass through the porous media only, thus facilitating the transmission measurements. The reflection images delivered a more homogeneous spatial illumination than the transmission images. Major perturbations of the transmission images were lens flare effects and light dispersion within the bead–water–Plexiglas system which smear the front of the plume. Based on the conducted evaluation of transmissive and reflective intensity measurements, the reflection data delivered more reliable intensity values to derive solute concentrations in intermediate scale flow tank experiments.     

Simulation of rapidly varying flow using an efficient TVD–MacCormack scheme

An efficient numerical scheme is outlined for solving the SWEs (shallow water equations) in environmental flow; this scheme includes the addition of a five‐point symmetric total variation diminishing (TVD) term to the corrector step of the standard MacCormack scheme. The paper shows that the discretization of the conservative and non‐conservative forms of the SWEs leads to the same finite difference scheme when the source term is discretized in a certain way. The non‐conservative form is used in the solution outlined herein, since this formulation is simpler and more efficient. The time step is determined adaptively, based on the maximum instantaneous Courant number across the domain. The bed friction is included either explicitly or implicitly in the computational algorithm according to the local water depth. The wetting and drying process is simulated in a manner which complements the use of operator‐splitting and two‐stage numerical schemes. The numerical model was then applied to a hypothetical dam‐break scenario, an experimental dam‐break case and an extreme flooding event over the Toce River valley physical model. The predicted results are free of spurious oscillations for both sub‐ and super‐critical flows, and the predictions compare favourably with the experimental measurements. Copyright © 2006 John Wiley & Sons, Ltd.     

Sensitivity of plume dynamics to the parameterization of vertical mixing

A three‐dimensional primitive equation, baroclinic numerical model incorporating a range of turbulence closure schemes is used to investigate the effects of vertical diffusion of momentum and density upon the spread of a freshwater plume, with particular reference to the Ebro plume. Initial calculations show that there are some differences in the horizontal spread and vertical mixing of the plume when diffusion coefficients are computed from a two‐equation turbulence energy model compared with a one‐equation model. To understand results from the turbulence energy models, the sensitivity of the plume dynamics to variations in the coefficient of vertical eddy viscosity and diffusivity is also considered, with increases in these parameters having a significant effect upon the cross‐shore and along‐shore spread of the plume. Also, increasing these parameters changes the plume characteristics from supercritical to subcritical and reduces the occurrence of meandering and baroclinic instability along the plume's off‐shore edge. However, differences in the southerly spread (the direction of Kelvin wave propagation in the model) of the plume (although not its northerly spread) produced by changes in diffusion coefficients are small compared with the influence of changes in the bottom slope, upon the along‐shore southerly spread of the plume, which moves in the direction of Kelvin wave propagation in the near coastal region. Results from the series of calculations are used as a guide in experimental design, with reference to a planned experiment in the Ebro region involving a coastal HF Radar deployment, as well as off‐shore measurements. Calculations suggest that surface current measurements from a coastal HF Radar, together with a detailed survey of the density field associated with the plume, may be an appropriate, although indirect, means of determining suitable mixing coefficients to use in plume discharge problems. Detailed measurements of water depth variation will also be required. Copyright © 1999 John Wiley & Sons, Ltd.     

Crack paths in three-dimensional elastic solids. ii: three-term expansion of the stress intensity factors—applications and perspectives

This work continues the calculation of the stress intensity factors, as a function of position s along the front of an arbitrary (kinked and curved) infinitesimal extension of some arbitrary crack on some three-dimensional body. More precisely, ε denoting a small parameter which the crack extension length is proportional to, what is studied here is the third term, proportional to εfn2 = ε and noted K ⁽¹⁾ (s) ε, of the expansion of these stress intensity factors at the point s of the crack front in powers of ε. The novelties with respect to previous works due to Gao and Rice on the one hand and Nazarov on the other hand, are that both the original crack and its extension need not necessarily be planar, and that a kink (discontinuity of the tangent plane to the crack) can occur all along the original crack front. Two expressions of K ⁽¹⁾ (s) are obtained; the difference is that the first one is more synthetic whereas the second one makes the influence of the kink angle (which can vary along the original crack front) more explicit. Application of some criterion then allows to obtain the apriori unknown geometric parameters of the small crack extension (length, kink angle, curvature parameters). The small scale segmentation of the crack front which is observed experimentally in the presence of mode III is disregarded here because a large scale point of view is adopted; this phenomenon will be discussed in a separate paper. It is shown how these results can be used to numerically predict crack paths over arbitrary distances in three dimensions. Simple applications to problems of configurational stability and bifurcation of the crack front are finally presented.     

Injection and Storage of CO<Subscript>2</Subscript> in Deep Saline Aquifers: Analytical Solution for CO<Subscript>2</Subscript> Plume Evolution During Injection

Injection of fluids into deep saline aquifers is practiced in several industrial activities, and is being considered as part of a possible mitigation strategy to reduce anthropogenic emissions of carbon dioxide into the atmosphere. Injection of CO₂ into deep saline aquifers involves CO₂ as a supercritical fluid that is less dense and less viscous than the resident formation water. These fluid properties lead to gravity override and possible viscous fingering. With relatively mild assumptions regarding fluid properties and displacement patterns, an analytical solution may be derived to describe the space–time evolution of the CO₂ plume. The solution uses arguments of energy minimization, and reduces to a simple radial form of the Buckley–Leverett solution for conditions of viscous domination. In order to test the applicability of the analytical solution to the CO₂ injection problem, we consider a wide range of subsurface conditions, characteristic of sedimentary basins around the world, that are expected to apply to possible CO₂ injection scenarios. For comparison, we run numerical simulations with an industry standard simulator, and show that the new analytical solution matches a full numerical solution for the entire range of CO₂ injection scenarios considered. The analytical solution provides a tool to estimate practical quantities associated with CO₂ injection, including maximum spatial extent of a plume and the shape of the overriding less-dense CO₂ front.     

甲烷/空气预混气体火焰的传播特征

利用高速纹影摄像等技术探讨了密闭管道内不同当量比的甲烷/空气预混气体火焰的传播特征。结果表明,当甲烷含量接近当量值时,预混气体火焰传播中会发生火焰阵面由向未燃区弯曲到向已燃区弯曲的转折过程,逐渐由层流燃烧转变成湍流燃烧,并形成Tulip火焰结构;当甲烷含量偏离当量值一定程度时,预混火焰呈现出典型的层流燃烧特征,不会发生火焰阵面由向未燃区弯曲到向已燃区弯曲的转折过程。Tulip火焰结构形成于火焰传播速度迅速降低的区间里,且只有当减速阶段的最大加速度的绝对值大于某一数值时才能形成;Tulip火焰结构是预混火焰由层流燃烧向湍流燃烧转变的一个中间过程。     

Runup and green water velocities due to breaking wave impinging and overtopping

The present study investigates, through measurements in a 2D wave tank, the velocity fields of a plunging breaking wave impinging on a structure. As the wave breaks and overtops the structure, so-called green water is generated. The flow becomes multi-phased and chaotic as a large aerated region is formed in the flow in the vicinity of the structure while water runs up onto the structure. In this study, particle image velocimetry (PIV) and its derivative, bubble image velocimetry (BIV), were employed to measure the velocity field in front and on top of the structure. Mean and turbulence properties were obtained through ensemble averaging repeated tests. The dominant and maximum velocity of the breaking wave and associated green water are discussed for the three distinct phases of the impingement–runup–overtopping sequence. Initially the flow is mainly horizontal right before the breaking wave impinges on the structure. The flow then becomes primarily vertical and rushes upward along the front wall of the structure right after the impingement. Subsequently, the flow becomes mainly horizontal on top of the structure as the remaining momentum in the wave crest carries the green water through. The distribution of the green water velocity along the top of the structure has a nonlinear profile and the maximum velocity occurs near the front of the fast moving water. Using the measured data and applying dimensional analysis, a similarity profile for the green water flow on top of the structure was obtained, and a prediction equation was formulated. The prediction equation may be used to predict the green water velocity caused by extreme waves in a hurricane.