Direct Numerical Simulations of Interfacial Turbulence at Low Froude and Weber Numbers

Sea surface temperature accessible through use of remote sensing techniques (IR imaging, etc.) suggests abundant flow and thermal field information at the ocean surface that is closely related to subsurface turbulent activities. The suggested information includes wind stress, surface dissipation, underneath velocity and vorticity, and heat and gas transportation. Due to the constantly outgoing interfacial latent and sensible heat flux, the very surface of the ocean is often cooler than the bulk. This so called ‘cool skin layer’ below the very surface is greatly involved in the underlying interfacial turbulence and is the primary support of using sea surface temperature imaging to detect the subsurface activities. In addition, studies have shown that for this detection method the effects of ubiquitous surfactants (surface free agents) to the subsurface turbulence should also be considered. In the case when the wind stress at the surface is far less significant than the buoyancy force in the water phase, the cool skin layer accumulates and triggers free convection. A series of numerical simulations is conducted to reproduce such a free convection flow to obtain detailed statistics and structural features in order to investigate the correlation between the surface temperature and the subsurface activities of the flow. The simulations are also aimed at the quantitative evaluation of the surfactant effects on the flow. The results of the simulations demonstrate that the surface temperature is statistically and structurally correlated to the subsurface activities in various patterns, and that surfactant has a certain influence to the subsurface turbulence with an overall effect of reducing the average surface temperature. Based upon the framework of the controlled flux method, a novel approach to actively determine the interfacial gas transfer velocity at the free convection surface is proposed and numerically investigated. The proposed and simulated approach employs a temporal volumetric heating source to suppress the free convection. The heating source is defined and parameterized with respect to the physical properties of radiation absorption in water phase. Observation and interpretation of the surface temperature evolution and the flow features during and after the heating suggest the effective suppression of the free convection, the onset of the Rayleigh instability and the re-establishment of the free convection. Based on that, an analytical conduction model is formulated to obtain the heat transfer velocity at the free surface from the surface temperature. The gas transfer velocity is then inferred through similarity

Ultra-fast Imaging of Two-Phase Flow in Structured Monolith Reactors; Techniques and Data Analysis

This thesis will address the use of nuclear magnetic resonance (NMR) and magnetic resonance imaging (MRI) techniques to probe the “monolith reactor”, which consists of a structured catalyst over which reactions may occur. This reactor has emerged as a potential alternative to more traditional chemical engineering systems such as trickle bed and slurry reactors. However, being a relatively new design, its associated flow phenomena and design procedures are not rigorously understood, which is retarding its acceptance in industry. Traditional observations are unable to provide the necessary information for design since the systems are opaque and dynamic. Therefore, NMR is proposed as an ideal tool to probe these systems in detail. The theory of NMR is summarised and the development of novel NMR techniques is presented. Novel techniques are validated in simple systems, and tested in more complex systems to ascertain their quantitative nature, and to find their limitations. These techniques are improvements over existing techniques in that they either decrease the acquisition time (allowing the observation of dynamically-changing systems) or allow us to probe systems in different ways to extract useful information. The goal of this research is to better understand the flow phenomena present in such systems, and to use this information to design better, more efficient, more controllable industrial reactors. The analysis of the NMR data acquired is discussed in detail, and several novel image-processing techniques have been developed to aid in the quantification of features within the images, and also to measure quantities such as holdup and velocity. These novel techniques are validated, and then applied to the systems of interest. Various configurations of monolith reactor, ranging from low flow rate systems to more challenging (and more industrially relevant) turbulent systems, are probed using these methods, and the contrasting flow phenomena and performance of these systems are discussed, with a view to optimisation of the choice of design parameters

Study of Kinematics of Extreme Waves Impacting Offshore and Coastal Structures by Non Intrusive Measurement Techniques

Extreme wave flows associated with a large scale wave breaking during interactions with marine structures or complex coastal geography of is one of the major concerns in a design of coastal and ocean structures. In order to properly understand the impact mechanisms of breaking extreme waves, full field evaluations of impacting multiphase flow velocities should be properly conducted first. In this context, this present dissertation experimentally investigated velocity structures of turbulent, multiphase wave flow velocities during active interactions with various offshore and onshore ocean environments. First, initial inundation flow structures of tsunami-like long waves interacting with complex coastal topography are experimentally investigated. Turbulent wave surface velocities were effectively measured by introducing a non-intrusive video imagery technique, the “wave front tracing method”. Three distinctive configurations for patch layouts that vary either in characteristic patch diameter (D) or in center-to-center spacing between patches (ΔS) were employed. That is, patch layouts consisted of six (G1) and twelve (G2), “small” circular macro roughness patches of D =1.2 m and six, “large” circular macro roughness patches (G3) of D = 1.7 m were employed, respectively. A patch layout employed for G1 appears to be effective in reducing the u velocities along the centerlines of the reference patch that consistently decreased to 85% of a convergence velocity U = 2m/s and to 45% of U. However, in the channel, u velocities hardly reduced below the convergence velocity. On the other hand, the patch layout G2 is observed as rather effective in uniformly reducing the u velocities alongshore. The hand, the patch layout G3 is observed as effective in suppressing the alongshore variability in flow behind the frontal patches. This may be due to the "holding-up" effects produced by the large patches holding the flow within the patch for a longer duration. Furthermore, such a "holding-up" effect from G3 appears to induce a large inundation depth in the flow along the opening. Next, green water velocities and dynamic impacts of the extreme ocean waves on a fixed offshore deck structure are investigated. The experiments focused on the impacting waves generated in a large-scale, three-dimensional ocean wave basin. Using the BIV technique, overall flow structures and temporal and spatial distributions of the maximum velocities were successfully evaluated. The most significant spatial variability in mean velocities in the propagating direction was found from the protruding wave front near the center of the deck during early stages of the wave run-up. The maximum front speed of 1.4C was first observed in the center of the deck near y = 0 at a midpoint of the deck (x = 0.5L), where C is the wave phase speed. The flow velocities started decreasing below 1C over all fields once the wave frontal flow passed the rear edge and started leaving the deck. Pressure measurements were also conducted at four different vertical positions on vertical measurement planes at three different locations on the horizontal plane. Most of measured pressures showed impulsive impact patterns with sudden rises of pressure peaks. The highest pressure was observed as 1.56pC^(2) at x = L/2. Correlations between wave kinematic energy and dynamic pressure were examined to determine the impact coefficients ci'. ci' varied within relatively narrow ranges 0.29 ≤ ci' ≤ 1.56. In the present large scale experiments, the impact pressures on the structures are strongly affected by both variability of flow structures and impulsiveness of impacting waves containing considerable air volumes. Lastly, the study is extended for more violent sloshing wave flows. The study experimentally investigated flow kinematics and impact pressures of a partially filled liquid sloshing flow during the periodic longitudinal motion of a rectangular tank. The horizontal velocities near the free surface reached 1.6C with C being the wave phase speed calculated based on the shallow water assumption. As the tank reached its maximum displacement and about to reverse, the dominant flow changed its direction rapidly to vertical upward after the breaking wave crest impinging on the side wall and forming an up-rushing jet. The vertical velocity of the rising jet reached 3.4C before it impacted the top wall. During the flip-through event as the fast moving wave crest collided with the side wall, the steep wave crest resulted in a focused impact on the side wall at the SWL. The resulting impulsive peak pressure was recorded as about 10ghρ immediately followed by the evident pressure oscillation with a frequency approximately 500 Hz. After the wall impact, the multiphase up-rushing jet shot up and impacted the top wall. The magnitude of the pressure was again about 10ghρ, similar to that recorded by the breaking wave impact on the side wall. Correlating the dynamic impact pressures with the corresponding local maximum flow velocities in the direction normal to the walls was performed by introducing the impact coefficient ic and the modified impact coefficient c'_(i) , defined as p_(max)=c_(i)pV^(2)= c'_(i)pC^(2) with V_(max) being the magnitude of the maximum local velocities. The average values of the modified impact coefficient c′_(i) between the side wall impacts and the top wall impacts were nearly identical, with the average value of c'_(i)=5.2

Modeling, Design, Packaging and Experimental Analysis of Liquid-Phase Shear-Horizontal Surface Acoustic Wave

Recent advances in microbiology, computational capabilities, and microelectromechanical-system fabrication techniques permit modeling, design, and fabrication of low-cost, miniature, sensitive and selective liquid-phase sensors and labon- a-chip systems. Such devices are expected to replace expensive, time-consuming, and bulky laboratory-based testing equipment. Potential applications for devices include: fluid characterization for material science and industry; chemical analysis in medicine and pharmacology; study of biological processes; food analysis; chemical kinetics analysis; and environmental monitoring. When combined with liquid-phase packaging, sensors based on surface-acoustic-wave (SAW) technology are considered strong candidates. For this reason such devices are focused on in this work; emphasis placed on device modeling and packaging for liquid-phase operation. Regarding modeling, topics considered include mode excitation efficiency of transducers; mode sensitivity based on guiding structure materials/geometries; and use of new piezoelectric materials. On packaging, topics considered include package interfacing with SAW devices, and minimization of packaging effects on device performance. In this work novel numerical models are theoretically developed and implemented to study propagation and transduction characteristics of sensor designs using wave/constitutive equations, Green’s functions, and boundary/finite element methods. Using developed simulation tools that consider finite-thickness of all device electrodes, transduction efficiency for SAW transducers with neighboring uniform or periodic guiding electrodes is reported for the first time. Results indicate finite electrode thickness strongly affects efficiency. Using dense electrodes, efficiency is shown to approach 92% and 100% for uniform and periodic electrode guiding, respectively; yielding improved sensor detection limits. A numerical sensitivity analysis is presented targeting viscosity using uniform-electrode and shear-horizontal mode configurations on potassium-niobate, langasite, and quartz substrates. Optimum configurations are determined yielding maximum sensitivity. Results show mode propagation-loss and sensitivity to viscosity are correlated by a factor independent of substrate material. The analysis is useful for designing devices meeting sensitivity and signal level requirements. A novel, rapid and precise microfluidic chamber alignment/bonding method was developed for SAW platforms. The package is shown to have little effect on device performance and permits simple macrofluidic interfacing. Lastly, prototypes were designed, fabricated, and tested for viscosity and biosensor applications; results show ability to detect as low as 1% glycerol in water and surface-bound DNA crosslinking