Eddy diffusivity: A single dispersion analysis of high resolution drifters in a tidal shallow estuary

In an estuary, mixing and dispersion resulting from turbulence and small scale fluctuation has strong spatio-temporal variability which cannot be resolved in conventional hydrodynamic models while some models employs parameterizations large water bodies. This paper presents small scale diffusivity estimates from high resolution drifters sampled at 10 Hz for periods of about 4 hours to resolve turbulence and shear diffusivity within a tidal shallow estuary (depth < 3 m). Taylor's diffusion theorem forms the basis of a first order estimate for the diffusivity scale. Diffusivity varied between 0.001 – 0.02 m2/s during the flood tide experiment. The diffusivity showed strong dependence (R2 > 0.9) on the horizontal mean velocity within the channel. Enhanced diffusivity caused by shear dispersion resulting from the interaction of large scale flow with the boundary geometries was observed. Turbulence within the shallow channel showed some similarities with the boundary layer flow which include consistency with slope of 5/3 predicted by Kolmogorov's similarity hypothesis within the inertial subrange. The diffusivities scale locally by 4/3 power law following Okubo's scaling and the length scale scales as 3/2 power law of the time scale. The diffusivity scaling herein suggests that the modelling of small scale mixing within tidal shallow estuaries can be approached from classical turbulence scaling upon identifying pertinent parameters.

Analysis of gradually and spatially varied flow in sand-bed channels

Seepage through a sand bed affects the channel hydrodynamics, which in turn alters channel stability. Thus, the effect of seepage on its hydrodynamic parameters needs to be ascertained. The present work analyses spatially varied flow of a sand-bed channel subjected to seepage in the downward direction through a sand bed. Numerically calculated flow profiles affected by seepage have been verified using experimental observations. The present work also analyses the friction slope, velocity and bed shear stress variations along the channel for both seepage and no-seepage conditions. It was found that the downward seepage-induced channel flow has larger friction slope and bed shear stress than that of no-seepage.

Weakly non-linear stability of a hydrodynamic accretion disc

When the cold accretion disc coupling between neutral gas and a magnetic field is so weak that the magnetorotational instability is less effective or even stops working, it is of prime interest to investigate the pure hydrodynamic origin of turbulence and transport phenomena. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases. In an accretion disc where the molecular viscosity is too small, the Reynolds number is large enough for the non-linear term to have new effects. We investigate the scenario of the `weakly non-linear' evolution of the amplitude of the linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to the plane Couette flow, but with the Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, because of the self-interaction, the least stable eigenmode will grow in an intermediate phase. Later, this will lead to higher-order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy-conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy (alpha A(c)(2), where A(c) is the threshold amplitude) needed for the flow to transform to the turbulent phase. Such an unstable phase is possible only if the Reynolds number >= 10(3-4). The numerical difficulties in obtaining such a large Reynolds number might be the reason for the negative result of numerical simulations on a pure hydrodynamic Keplerian accretion disc.

Virus-sized colloid transport in a single pore: Model development and sensitivity analysis

A mathematical model is developed to simulate the transport and deposition of virus-sized colloids in a cylindrical pore throat considering various processes such as advection, diffusion, colloid-collector surface interactions and hydrodynamic wall effects. The pore space is divided into three different regions, namely, bulk, diffusion and potential regions, based on the dominant processes acting in each of these regions. In the bulk region, colloid transport is governed by advection and diffusion whereas in the diffusion region, colloid mobility due to diffusion is retarded by hydrodynamic wall effects. Colloid-collector interaction forces dominate the transport in the potential region where colloid deposition occurs. The governing equations are non-dimensionalized and solved numerically. A sensitivity analysis indicates that the virus-sized colloid transport and deposition is significantly affected by various pore-scale parameters such as the surface potentials on colloid and collector, ionic strength of the solution, flow velocity, pore size and colloid size. The adsorbed concentration and hence, the favorability of the surface for adsorption increases with: (i) decreasing magnitude and ratio of surface potentials on colloid and collector, (ii) increasing ionic strength and (iii) increasing pore radius. The adsorbed concentration increases with increasing Pe, reaching a maximum value at Pe = 0.1 and then decreases thereafter. Also, the colloid size significantly affects particle deposition with the adsorbed concentration increasing with increasing particle radius, reaching a maximum value at a particle radius of 100 nm and then decreasing with increasing radius. System hydrodynamics is found to have a greater effect on larger particles than on smaller ones. The secondary minimum contribution to particle deposition has been found to increase as the favorability of the surface for adsorption decreases. The sensitivity of the model to a given parameter will be high if the conditions are favorable for adsorption. The results agree qualitatively with the column-scale experimental observations available in the literature. The current model forms the building block in upscaling colloid transport from pore scale to Darcy scale using Pore-Network Modeling. (C) 2014 Elsevier By. All rights reserved.

Hydrodynamic instability and shear layer effect on the response of an acoustically excited laminar premixed flame

Combustion instabilities can cause serious problems which limit the operating envelope of low-emission lean premixed combustion systems. Predicting the onset of combustion instability requires a description of the unsteady heat release driving the instability, i.e., the heat release response transfer function of the system. This study focuses on the analysis of fully coupled two-way interactions between a disturbance field and a laminar premixed flame that incorporates gas expansion effects by solving the conservation equations of a compressible fluid. Results of the minimum and maximum flame front deflections are presented to underline the impact of the hydrodynamic instability on the flame and the shear layer effect on the initial flame front wrinkling which is increased at decreasing gas expansion. These phenomena influence the magnitude of the burning area and burning area rate response of the flame at lower frequency excitation more drastically than reduced-order model (ROM) predictions even for low temperature ratios. It is shown that the general trend of the flame response magnitudes can be well captured at higher frequency excitation, where stretch effects are dominant. The phase response is influenced by the DL mechanism, which cannot be captured by the ROM, and by the resulting discrepancy in the flame pocket formation and annihilation process at the flame tip. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved,

Improved Analysis Method for Wave-Induced Pipeline Stability on Sandy Seabed

The existing Det Norske Veritas DNV Recommended Practice RP E305 for pipeline on-bottom stability is mainly based on the pipe–soil interaction model reported by Wagner et al. in 1987, and the wake model reported by Lambrakos et al. in 1987, to calculate the soil resistance and the hydrodynamic forces upon pipeline, respectively. Unlike the methods in the DNV Practice, in this paper, an improved analysis method is proposed for the on-bottom stability of a submarine pipeline, which is based on the relationships between Um/ gD 0.5 and Ws / D2 for various restraint conditions obtained by the hydrodynamic loading experiments, taking into account the coupling effects between wave, pipeline, and sandy seabed. The analysis procedure is illustrated with a detailed flow chart. A comparison is made between the submerged weights of pipeline predicted with the DNV Practice and those with the new method. The proposed analysis method may provide a helpful tool for the engineering practice of pipeline on-bottom stability design.

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the earths mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.

Effects of fundamental structure parameters on dynamic responses of submerged floating tunnel under hydrodynamic loads

This paper investigates the effects of structure parameters on dynamic responses of submerged floating tunnel (SFT) under hydrodynamic loads. The structure parameters includes buoyancy-weight ratio (BWR), stiffness coefficients of the cable systems, tunnel net buoyancy and tunnel length. First, the importance of structural damp in relation to the dynamic responses of SFT is demonstrated and the mechanism of structural damp effect is discussed. Thereafter, the fundamental structure parameters are investigated through the analysis of SFT dynamic responses under hydrodynamic loads. The results indicate that the BWR of SFT is a key structure parameter. When BWR is 1.2, there is a remarkable trend change in the vertical dynamic response of SFT under hydrodynamic loads. The results also indicate that the ratio of the tunnel net buoyancy to the cable stiffness coefficient is not a characteristic factor affecting the dynamic responses of SFT under hydrodynamic loads.

Hydrodynamic interaction of two unsteady model microorganisms.

The study of pair-wise interactions between swimming microorganisms is fundamental to the understanding of the rheological and transport properties of semi-dilute suspensions. In this paper, the hydrodynamic interaction of two ciliated microorganisms is investigated numerically using a boundary-element method, and the microorganisms are modeled as spherical squirmers that swim by time-dependent surface deformations. The results show that the inclusion of the unsteady terms in the ciliary propulsion model has a large impact on the trajectories of the interacting cells, and causes a significant change in scattering angles with potential important consequences on the diffusion properties of semi-dilute suspensions. Furthermore, the analysis of the shear stress acting on the surface of the microorganisms revealed that the duration and the intensity of the near-field interaction are significantly modified by the presence of unsteadiness. This observation may account for the hydrodynamic nature of randomness in some biological reactions, and supersedes the distinction between intrinsic randomness and hydrodynamic interactions, adding a further element to the understanding and modeling of interacting microorganisms.

Hydrodynamic entrapment of bacteria swimming near a solid surface.

The near-surface motility of bacteria is important in the initial formation of biofilms and in many biomedical applications. The swimming motion of Escherichia coli near a solid surface is investigated both numerically and experimentally. A boundary element method is used to predict the hydrodynamic entrapment of E. coli bacteria, their trajectories, and the minimum separation of the cell from the surface. The numerical results show the existence of a stable swimming distance from the boundary that depends only on the shape of the cell body and the flagellum. The experimental validation of the numerical approach allows one to use the numerical method as a predictive tool to estimate with reasonable accuracy the near-wall motility of swimming bacteria of known geometry. The analysis of the numerical database demonstrated the existence of a correlation between the radius of curvature of the near-wall circular trajectory and the separation gap. Such correlation allows an indirect estimation of either of the two quantities by a direct measure of the other without prior knowledge of the cell geometry. This result may prove extremely important in those biomedical and technical applications in which the near-wall behavior of bacteria is of fundamental importance.

Analysis of the piston ring/liner oil film development during warm-up for an SI-engine

A one-dimensional ring-pack lubrication model developed at MIT is applied to simulate the oil film behavior during the warm-up period of a Kohler spark ignition engine [1]. This is done by making assumptions for the evolution of the oil temperatures during warm-up and that the oil control ring during downstrokes is fully flooded. The ring-pack lubrication model includes features such as three different lubrication regimes, i.e. pure hydrodynamic lubrication, boundary lubrication and pure asperity contact, non-steady wetting of both inlet and outlet of the piston ring, capability to use all ring face profiles that can be approximated by piece-wise polynomials and, finally, the ability to model the rheology of multi-grade oils. Not surprisingly, the simulations show that by far the most important parameter is the temperature dependence of the oil viscosity. This dependence is subsequently examined further by choosing different oils. The baseline oil is SAE 10W30 and results are compared to those using the SAE 30 and the SAE 10W50 oils.

Simplifying mooring analysis for deepwater systems using truncation

Over recent years academia and industry have engaged with the challenge of model testing deepwater structures at conventional scales. One approach to the limited depth problem has been to truncate the lines. This concept will be introduced, highlighting the need to better understand line dynamic processes. The type of line truncation developed here models the upper sections of each line in detail, capturing wave action and all coupling effects with the vessel, terminating to an approximate analytical model that aims to simulate the remainder of the line. A rationale for this is that in deep water transverse elastic waves of a line are likely to decay before they are reflected at the seabed because of nonlinear hydrodynamic drag forces. The first part of this paper is centered on verification of this rationale. A simplified model of a mooring line that describes the transverse dynamics in wave frequency is used, adopting the equation of motion of an inextensible taut string. The line is submerged in still water, one end fixed at the bottom the other assumed to follow the vessel response, which can be harmonic or random. A dimensional analysis, supported by exact benchmark numerical solutions, has shown that it is possible to produce a universal curve for the decay of transverse vibrations along the line, which is suitable for any kind of line with any top motion. This has a significant engineering benefit, allowing for a rapid assessment of line dynamics - it can be useful in deciding whether a truncated line model is appropriate, and if so, at which point truncation might be applied. This is followed by developing a truncation mechanism, formulating an end approximation that can reproduce the correct impedance, had the line been continuous to full depth. It has been found that below a certain length criterion, which is also universal, the transverse vibrational characteristics for each line are inertia driven. As such the truncated model can assume a linear damper whose coefficient depends on the line properties and frequency of vibration. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).

Analysis on suspended sediment deposition rate for muddy coast of reclaimed land

A new expression for calculating suspended fine-sediment deposition rate is developed based on theoretic analysis and experiments. The resulting equation is applied to simulation of fine sediment deposition in the reclaimed land in the Hangzhou Bay, China. The hydrodynamic environment in this area is solved by use of a long wave model, which gives the 2D-velocity field and considers bathymetric changes due to fine sediment deposition. The expression is proved convenient to use in engineering practice, and the predicted deposition rate agrees with the annual data available from field measurements from the first year to the third year after the construction of the long groin as a reclaiming method.

Radiative hydrodynamic modelling and observations of the X-class solar flare on 2011 March 9

Aims. We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters.

Methods. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 s, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300 s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy.

Results. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density (ne ≤ 1015 cm-3) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

Calibration and validation of hydrodynamic, salt and heat transport models for Ria de Aveiro lagoon (Portugal)

Ria deAveiro is a very complex shallow water coastal lagoon located on the northwest of Portugal. Important issues would be left unanswered without a good understanding of hydrodynamic and transport processes occurring in the lagoon. Calibration and validation of hydrodynamic, salt and heat transport models for Ria de Aveiro lagoon are presented. The calibration of the hydrodynamic model was performed adjusting the bottom friction coefficient, through the comparison between measured and predicted time series of sea surface elevation for 22 stations. Harmonic analysis was performed in order to evaluate the model's accuracy. To validate the hydrodynamic model measured and predicted SSE values were compared for 11 stations, as well as main flow direction velocities for 10 stations. The salt and heat transport models were calibrated comparing measured and predicted time series of salinity and water temperature for 7 stations, and the RMS of the difference between the series was determined. These models were validated comparing the model results with an independent field data set. The hydrodynamic and the salt and heat transport models for Ria de Aveiro were successfully calibrated and validated. They reproduce accurately the barotropic flows and can therefore adequately represent the salt and heat transport and the heat transfer processes occurring in Ria deAveiro.