Towards understanding the mixing characteristics of turbulent buoyant flows

This work proposes a new simulation methodology in which variable density turbulent flows can be studied in the context of a mixing layer with or without the presence of gravity. Specifically, this methodology is developed to probe the nature of non-buoyantly-driven (i.e. isotropically-driven) or buoyantly-driven mixing deep inside a mixing layer. Numerical forcing methods are incorporated into both the velocity and scalar fields, which extends the length of time over which mixing physics can be studied. The simulation framework is designed to allow for independent variation of four non-dimensional parameters, including the Reynolds, Richardson, Atwood, and Schmidt numbers. Additionally, the governing equations are integrated in such a way to allow for the relative magnitude of buoyant energy production and non-buoyant energy production to be varied.

The computational requirements needed to implement the proposed configuration are presented. They are justified in terms of grid resolution, order of accuracy, and transport scheme. Canonical features of turbulent buoyant flows are reproduced as validation of the proposed methodology. These features include the recovery of isotropic Kolmogorov scales under buoyant and non-buoyant conditions, the recovery of anisotropic one-dimensional energy spectra under buoyant conditions, and the preservation of known statistical distributions in the scalar field, as found in other DNS studies.

This simulation methodology is used to perform a parametric study of turbulent buoyant flows to discern the effects of varying the Reynolds, Richardson, and Atwood numbers on the resulting state of mixing. The effects of the Reynolds and Atwood numbers are isolated by looking at two energy dissipation rate conditions under non-buoyant (variable density) and constant density conditions. The effects of Richardson number are isolated by varying the ratio of buoyant energy production to total energy production from zero (non-buoyant) to one (entirely buoyant) under constant Atwood number, Schmidt number, and energy dissipation rate conditions. It is found that the major differences between non-buoyant and buoyant turbulent flows are contained in the transfer spectrum and longitudinal structure functions, while all other metrics are largely similar (e.g. energy spectra, alignment characteristics of the strain-rate tensor). Also, despite the differences noted between fully buoyant and non-buoyant turbulent fields, the scalar field, in all cases, is unchanged by these. The mixing dynamics in the scalar field are found to be insensitive to the source of turbulent kinetic energy production (non-buoyant vs. buoyant).

Wave induced oscillations in harbors of arbitrary shape

Theoretical and experimental studies were conducted to investigate the wave induced oscillations in an arbitrary shaped harbor with constant depth which is connected to the open-sea.

A theory termed the “arbitrary shaped harbor” theory is developed. The solution of the Helmholtz equation, ∇²f + k²f = 0, is formulated as an integral equation; an approximate method is employed to solve the integral equation by converting it to a matrix equation. The final solution is obtained by equating, at the harbor entrance, the wave amplitude and its normal derivative obtained from the solutions for the regions outside and inside the harbor.

Two special theories called the circular harbor theory and the rectangular harbor theory are also developed. The coordinates inside a circular and a rectangular harbor are separable; therefore, the solution for the region inside these harbors is obtained by the method of separation of variables. For the solution in the open-sea region, the same method is used as that employed for the arbitrary shaped harbor theory. The final solution is also obtained by a matching procedure similar to that used for the arbitrary shaped harbor theory. These two special theories provide a useful analytical check on the arbitrary shaped harbor theory.

Experiments were conducted to verify the theories in a wave basin 15 ft wide by 31 ft long with an effective system of wave energy dissipators mounted along the boundary to simulate the open-sea condition.

Four harbors were investigated theoretically and experimentally: circular harbors with a 10° opening and a 60° opening, a rectangular harbor, and a model of the East and West Basins of Long Beach Harbor located in Long Beach, California.

Theoretical solutions for these four harbors using the arbitrary shaped harbor theory were obtained. In addition, the theoretical solutions for the circular harbors and the rectangular harbor using the two special theories were also obtained. In each case, the theories have proven to agree well with the experimental data.

It is found that: (1) the resonant frequencies for a specific harbor are predicted correctly by the theory, although the amplification factors at resonance are somewhat larger than those found experimentally,(2) for the circular harbors, as the width of the harbor entrance increases, the amplification at resonance decreases, but the wave number bandwidth at resonance increases, (3) each peak in the curve of entrance velocity vs incident wave period corresponds to a distinct mode of resonant oscillation inside the harbor, thus the velocity at the harbor entrance appears to be a good indicator for resonance in harbors of complicated shape, (4) the results show that the present theory can be applied with confidence to prototype harbors with relatively uniform depth and reflective interior boundaries.

Analysis of four decades of high elevation climate data

Four decades of instrumented climate records at D1 on Niwot Ridge suggest that high elevation data are an important - and even unique - part of the full climate picture. High elevation data provide information on changing lapse rates as well as model verification for global warming, which is predicted to occur earliest in high latitudes and at high elevations. The D1 records show climatic trends that arguably support global warming, assuming that greater planetary wave amplitude is verification of warming. Lapse rates reflect conditions of air mass stability, atmospheric moisture, and could [sic] cover, which contribute to feedback processes involving temperature, precipitation, and snowpack. The D1 record show a period, 1981-1985, when the lapse rate increased significantly, and this change was not detected by other data.

Micro-ramp control for oblique shock wave / boundary layer interactions

Supersonic engine intakes operating supercritically feature shock wave / boundary layer interactions (SBLIs), which are conventionally controlled using boundary layer bleed. The momentum loss of bleed flow causes high drag, compromising intake performance. Micro-ramp sub-boundary layer vortex generators (SBVGs) have been proposed as an alternative form of flow control for oblique SBLIs in order to reduce the bleed requirement. Experiments have been conducted at Mach 2.5 to characterise the flow details on such devices and investigate their ability to control the interaction between an oblique shock wave and the naturally grown turbulent boundary layer on the tunnel floor. Micro-ramps of four sizes with heights ranging from 25% to 75% of the uncontrolled boundary layer thickness were tested. The flow over all sizes of microramp was found to be similar, featuring streamwise counter-rotating vortices which entrain high momentum fluid, locally reducing the boundary layer displacement thickness. When installed ahead of the shock interaction it was found that the positioning of the micro-ramps is of limited importance. Micro-ramps did not eliminate flow separation. However, the previously two-dimensional separation was broken up into periodic three-dimensional separation zones. The interaction length was reduced and the pressure gradient across the interaction was increased.

Design and simulation of SMES system using YBCO tapes for direct drive wave energy converters

The ocean represents a huge energy reservoir since waves can be exploited to generate clean and renewable electricity; however, a hybrid energy storage system is needed to smooth the fluctuation. In this paper a hybrid energy storage system using a superconducting magnetic energy system (SMES) and Li-ion battery is proposed. The SMES is designed using Yttrium Barium Copper Oxide (YBCO) tapes, which store 60 kJ electrical energy. The magnet component of the SMES is designed using global optimization algorithm. Mechanical stress, coupled with electromagnetic field, is calculated using COMSOL and Matlab. A cooling system is presented and a suitable refrigerator is chosen to maintain a cold working temperature taking into account four heat sources. Then a microgrid system of direct drive linear wave energy converters is designed. The interface circuit connecting the generator and storage system is given. The result reveals that the fluctuated power from direct drive linear wave energy converters is smoothed by the hybrid energy storage system. The maximum power of the wave energy converter is 10 kW. © 2012 IEEE.

Vortex migration caused by travelling magnetic wave in a 2 in. diameter YBa2Cu3O7-δ thin film

We provide experimental evidence for a vortex migration phenomenon in YBa2Cu3O7-δ (YBCO) thin film caused by travelling magnetic wave. The experiment is carried out on a 2 in. diameter YBCO thin film with a circular-type magnetic flux pump. We found that the travelling wave helps the vortices migrate into the centre of the sample: after the zero-field cooling process, the increase of the flux density in the centre is four times larger than the amplitude of the travelling wave. The reason for this massive vortex migration is probably due to the magnetic stress variation caused by the travelling wave: the magnetic stress increases locally in the crest region while decreases locally in the trough region, which could help the vortices to move locally. A comparison shows that the magnetization by standing wave can be easily predicted by Bean's model while travelling wave causes vortex migration generally much larger than the prediction of Bean's model. It is possible that travelling magnetic wave can be an effective way to magnetize a type II superconductor in considering this unusual vortex dynamics. © 2013 AIP Publishing LLC.

Integrated ZnO surface acoustic wave microfluidic components for biosensor

This paper describes the novel nanocrystalline film ZnO surface acoustic wave devices, which demonstrate their great potential for the portable disease diagnostic system with integrated functions of microfluidic transport, mixing and biosensing. The devices can be easily integrated with electronic control circuitry and fabricated with low temperature process on Si, glass or even polymer substrates. The liquid convection and internal streaming patterns was easily induced by acoustic wave at signal voltages. With further increase in applied voltage to above 20V, the liquid droplet was pushed forward. Immunoreaction-based bio-detection PSA/ACT, all based on SAW devices on thin film piezoelectric ZnO on Si substrate was demonstrated. © 2009 CBMS.

Resonant Enhanced Wave Filter and Waveguide via Surface Plasmons

The authors present an analysis of plasmonic wave filter and curved waveguide, simulated using a 2-D finite-difference time-domain technique. With different dielectric materials or surface structures located on the interface of the metal/dielectric, the resonant enhanced wave filter can divide light waves of different wavelengths and guide them with low losses. And the straight or curved waveguide can confine and guide light waves in a subwavelength scale. Within the 20 mu m simulation region, it is found that the intensity of the guided light at the interface is roughly four times the peak intensity of the incident light.

TUNNELING ESCAPE TIME OF ELECTRONS FROM A QUANTUM-WELL WITH GAMMA-X MIXING EFFECT

By using the envelope function method we calculated the tunneling escape time of electrons from a quantum well. We adopted a simplified interface matrix to describe the GAMMA-X mixing effect, and employed a wave packet method to determine the tunneling escape time. When the GAMMA state in the well was in resonance with the X state in the barrier, the escape time reduced remarkably. However, it was possible that the wave functions in two different channels, i.e., GAMMA-GAMMA-GAMMA and GAMMA-X-GAMMA, could interfere destructively, leading the escape time greater than that of pure GAMMA-GAMMA-GAMMA tunneling.

Quantum confinement effect in thin quantum wires

The electronic states and optical transition properties of three semiconductor wires Si? GaAs, and ZnSe are studied by the empirical pseudopotential homojunction model. The energy levels, wave functions, optical transition matrix elements, and lifetimes are obtained for wires of square cross section with width from 2 to 5 (root 2a/2), where a is the lattice constant. It is found that these three kinds of wires have different quantum confinement properties. For Si wires, the energy gap is pseudodirect, and the wave function of the electronic ground state consists mainly of four bulk Delta states. The optical transition matrix elements are much smaller than that of a direct transition, and increase with decreasing wire width. Where the width of wire is 7.7 Angstrom, the Si wire changes from an indirect energy gap to a direct energy gap due to mixing of the bulk Gamma(15) state. For GaAs wires. the energy gap is also pseudodirect in the width range considered, but the optical transition matrix elements are larger than those of Si wires by two orders of magnitude for the same width. However, there is no transfer to a direct energy gap as the wire width decreases. For ZnSe wires, the energy gap is always direct, and the optical transition matrix elements are comparable to those of the direct energy gap bulk semiconductors. They decrease with decreasing wire width due to mixing of the bulk Gamma(1) state with other states. All quantum confinement properties are discussed and explained by our theoretical model and the semiconductor energy band structures derived. The calculated lifetimes of the Si wire, and the positions of photoluminescence peaks, are in good agreement with experimental results.

Development and transition of spiral wave in the coupled Hindmarsh-Rose neurons in two-dimensional space

The dynamics and the transition of spiral waves in the coupled Hindmarsh-Rose (H-R) neurons in two-dimensional space are investigated in the paper. It is found that the spiral wave can be induced and developed in the coupled HR neurons in two-dimensional space, with appropriate initial values and a parameter region given. However, the spiral wave could encounter instability when the intensity of the external current reaches a threshold value of 1.945. The transition of spiral wave is found to be affected by coupling intensity D and bifurcation parameter r. The spiral wave becomes sparse as the coupling intensity increases, while the spiral wave is eliminated and the whole neuronal system becomes homogeneous as the bifurcation parameter increases to a certain threshold value. Then the coupling action of the four sub-adjacent neurons, which is described by coupling coefficient D', is also considered, and it is found that the spiral wave begins to breakup due to the introduced coupling action from the sub-adjacent neurons (or sites) and together with the coupling action of the nearest-neighbour neurons, which is described by the coupling intensity D.

Effect of low-frequency Rossby wave on thermal structure of the upper southwestern tropical Indian Ocean

We investigate the influence of low-frequency Rossby waves on the thermal structure of the upper southwestern tropical Indian Ocean (SWTIO) using Argo profiles, satellite altimetric data, sea surface temperature, wind field data and the theory of linear vertical normal mode decomposition. Our results show that the SWTIO is generally dominated by the first baroclinic mode motion. As strong downwelling Rossby waves reach the SWTIO, the contribution of the second baroclinic mode motion in this region can be increased mainly because of the reduction in the vertical stratification of the upper layer above thermocline, and the enhancement in the vertical stratification of the lower layer under thermocline also contributes to it. The vertical displacement of each isothermal is enlarged and the thermal structure of the upper level is modulated, which is indicative of strong vertical mixing. However, the cold Rossby waves increase the vertical stratification of the upper level, restricting the variability related to the second baroclinic mode. On the other hand, during decaying phase of warm Rossby waves, Ekman upwelling and advection processes associated with the surface cyclonic wind circulation can restrain the downwelling processes, carrying the relatively colder water to the near-surface, which results in an out-of-phase phenomenon between sea surface temperature anomaly (SSTA) and sea surface height anomaly (SSHA) in the SWTIO.

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2.5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, influence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.

Tidally induced upwelling off Yangtze River estuary and in Zhejiang coastal waters in summer

MASNUM wave-tide-circulation coupled numerical model (MASNUM coupled model, hereinafter) is developed based on the Princeton Ocean Model (POM). Both POM and MASNUM coupled model are applied in the numerical simulation of the upwelling off Yangtze River estuary and in Zhejiang coastal waters in summer. The upwelling mechanisms are analyzed from the viewpoint of tide, and a new mechanism is proposed. The study suggests that the tidally inducing mechanism of the upwelling includes two dynamic aspects: the barotropic and the baroclinic process. On the one hand, the residual currents induced by barotropic tides converge near the seabed, and upwelling is generated to maintain mass conservation. The climbing of the residual currents along the sea bottom slope also contributes to the upwelling. On the other hand, tidal mixing plays a very important role in inducing the upwelling in the baroclinic sea circumstances. Strong tidal mixing leads to conspicuous front in the coastal waters. The considerable horizontal density gradient across the front elicits a secondary circulation clinging to the tidal front, and the upwelling branch appears near the frontal zone. Numerical experiments are designed to determine the importance of tide in inducing the upwelling. The results indicate that tide is a key and dominant inducement of the upwelling. Experiments also show that coupling calculation of the four main tidal constituents(M-2, S-2, K-1, and O-1), rather than dealing with the single M-2 constituent, improves the modeling precision of the barotropic tide-induced upwelling.