On the feasibility of time-lapse superconducting gravimetry for reservoir monitoring

The feasibility of monitoring fluid flow subsurface processes that result in density changes, using the iGrav superconducting gravimeter, is investigated. Practical targets include steam-assisted gravity drainage (SAGD) bitumen depletion and water pumping from aquifers, for which there is currently a void in low-impact, inexpensive monitoring techniques. This study demonstrates that the iGrav has the potential to be applied to multi-scale and diverse reservoirs. Gravity and gravity gradient signals are forward modeled for a real SAGD reservoir at two time steps, and for surface-fed and groundwater-fed aquifer pumping models, to estimate signal strength and directional dependency of water flow. Time-lapse gravimetry on small-scale reservoirs exhibits two obstacles, namely, a µgal sensitivity requirement and high noise levels in the vicinity of the reservoir. In this study, both limitations are overcome by proposing (i) a portable superconducting gravimeter, and (ii) a pair of instruments under various baseline geometries. This results in improved spatial resolution for locating depletion zones, as well as the cancellation of noise common in both instruments. Results indicate that a pair of iGrav superconducting gravimeters meet the sensitivity requirements and the spatial focusing desired to monitor SAGD bitumen migration at the reservoir scales. For SAGD reservoirs, the well pair separation, reservoir depth, and survey sampling determine the resolvability of individual well pair depletion patterns during the steam chamber rising phase, and general reservoir depletion patterns during the steam chamber spreading phase. Results show that monitoring water table elevation changes due to pumping and tracking whether groundwater or surface water is being extracted are feasible.

Parametric Analysis of Concrete Solar Collectors

Concrete solar collectors offer a type of solar collector with structural, aesthetic and economic advantages over current populartechnologies. This study examines the influential parameters of concrete solar collectors. In addition to the external conditions,the performance of a concrete solar collector is influenced by the thermal properties of the concrete matrix, piping network andfluid. Geometric and fluid flow parameters also influence the performance of the concrete solar collector. A literature review ofconcrete solar collectors is conducted in order to define the benchmark parameters from which individual parameters are thencompared. The numerical model consists of a 1D pipe flow network coupled with the heat transfer in a 3D concrete domain. Thispaper is concerned with the physical parameters that define the concrete solar collector, thus a constant surface temperature isused as the exposed surface boundary condition with all other surfaces being insulated. Results show that, of the parametersinvestigated, the pipe spacing, ps, concrete conductivity, kc, and the pipe embedment depth, demb, are among those parameterswhich have greatest effect on the collector’s performance. The optimum balance between these parameters is presented withrespect to the thermal performance and discussed with reference to practical development issues.

PHYSICA: A multiphysics computational framework and its application to casting simulations

Metal casting is a process governed by the interaction of a range of physical phenomena. Most computational models of this process address only what are conventionally regarded as the primary phenomena – heat conduction and solidification. However, to predict other phenomena, such as porosity formation, requires modelling the interaction of the fluid flow, heat transfer, solidification and the development of stressdeformation in the solidified part of the casting. This paper will describe a modelling framework called PHYSICA[1] which has the capability to stimulate such multiphysical phenomena.

Dynamics of an axisymmetric electromagnetic 'crucible' melting

Different industrial induction melting processes involve free surface and melt-solid interface of the liquid metal subject to dynamic change during the technological operation. Simulation of the liquid metal dynamics requires to solve the non-linear, coupled hydrodynamic-electromagnetic-heat transfer problem accounting for the time development of the liquid metal free boundary with a suitable turbulent viscosity model. The present paper describes a numerical solution method applicable for various axisymmetric induction melting processes, such as, crucible with free top surface, levitation, semi-levitation, cold crucible and similar melting techniques. The presented results in the cases of semi-levitation and crucible with free top surface meltings demonstrate oscillating transient behaviour of the free metal surface indicating the presence of gravity-inertial-electromagnetic waves which are coupled to the internal fluid flow generated by both the rotational and potential parts of the electromagnetic force.

Modelling the mixing of lead bullion during the refining process

Abstract not available

Modelling interacting phenomena during the solidification of metals in the casting process

A variety of interacting complex phenomena takes place during the casting of metallic components. Here molten metal is poured into a mould cavity where it flows, cools, solidifies and then deforms in its solid state. As the metal cools, thermal gradients will promote thermal convection which will redistribute the heat around the component (usually from feeders or risers) towards the solidification front and mushy zone. Also, as the evolving solid regions of the cast component deform they will form gap at the cast-mould interface. This gap may change the rate of solidification in certain parts the casting, hence affecting the manner in which the cast component solidifies. Interaction between a cast component and its surrounding mould will also govern stress magnitudes in both the cast and mould -these may lead to defects such as cracks. This paper presents a multiphysics modelling approach to this complex process. Emphasis will be placed on the interacting phenomena taking place during the process and the modelling strategy used. Comparisons with plant data are also be given.

Investigating wetting characteristics on microstructured surfaces for superhydrophobicity and metal microcasting

The engineering of liquid behavior on surfaces is important for infrastructure, transportation, manufacturing, and sensing. Surfaces can be rendered superhydrophobic by microstructuring, and superhydrophobic devices could lead to practical corrosion inhibition, self-cleaning, fluid flow control, and surface drag reduction. To more fully understand how liquid interacts with microstructured surfaces, this dissertation introduces a direct method for determining droplet solid-liquid-vapor interfacial geometry on microstructured surfaces. The technique performs metrology on molten metal droplets deposited onto microstructured surfaces and then frozen. Unlike other techniques, this visualization technique can be used on large areas of curved and opaque microstructured surfaces to determine contact line. This dissertation also presents measurements and models for how curvature and flexing of microstructured polymers affects hydrophobicity. Increasing curvature of microstructured surfaces leads to decreased slide angle for liquid droplets suspended on the surface asperities. For a surface with regularly spaced asperities, as curvature becomes more positive, droplets suspended on the tops of asperities are suspended on fewer asperities. Curvature affects superhydrophobicity because microscopic curvature changes solid-liquid interaction, pitch is altered, and curvature changes the shape of the three phase contact line. This dissertation presents a model of droplet interactions with curved microstructured surfaces that can be used to design microstructure geometries that maintain the suspension of a droplet when curved surfaces are covered with microstructured polymers. Controlling droplet dynamics could improve microfluidic devices and the shedding of liquids from expensive equipment, preventing corrosion and detrimental performance. This dissertation demonstrates redirection of dynamic droplet spray with anisotropic microstructures. Superhydrophobic microstructured surfaces can be economically fabricated using metal embossing masters, so this dissertation describes casting-based microfabrication of metal microstructures and nanostructures. Low melting temperature metal was cast into flexible silicone molds which were themselves cast from microfabricated silicon templates. The flexibility of the silicone mold permits casting of curved surfaces, which this dissertation demonstrates by fabricating a cylindrical metal roller with microstructures. The metal microstructures can be in turn used as a reusable molding tool. This dissertation also describes an industrial investment casting process to produce aluminum molds having integrated microstructures. Unlike conventional micromolding tools, the aluminum mold was large and had complex curved surfaces. The aluminum was cast into curved microstructured ceramic molds which were themselves cast from curved microstructured rubber. Many structures were successfully cast into the aluminum with excellent replication fidelity, including circular, square, and triangular holes. This dissertation demonstrates molding of large, curved surfaces having surface microstructures using the aluminum mold. This work contributes a more full understanding of the phenomenon of superhydrophobicity and techniques for the economic fabrication of superhydrophobic microstructures.

Studies of turbulent flows in continuous casting of steel with and without magnetic field

This thesis develops and tests various transient and steady-state computational models such as direct numerical simulation (DNS), large eddy simulation (LES), filtered unsteady Reynolds-averaged Navier-Stokes (URANS) and steady Reynolds-averaged Navier-Stokes (RANS) with and without magnetic field to investigate turbulent flows in canonical as well as in the nozzle and mold geometries of the continuous casting process. The direct numerical simulations are first performed in channel, square and 2:1 aspect rectangular ducts to investigate the effect of magnetic field on turbulent flows. The rectangular duct is a more practical geometry for continuous casting nozzle and mold and has the option of applying magnetic field either perpendicular to broader side or shorter side. This work forms the part of a graphic processing unit (GPU) based CFD code (CU-FLOW) development for magnetohydrodynamic (MHD) turbulent flows. The DNS results revealed interesting effects of the magnetic field and its orientation on primary, secondary flows (instantaneous and mean), Reynolds stresses, turbulent kinetic energy (TKE) budgets, momentum budgets and frictional losses, besides providing DNS database for two-wall bounded square and rectangular duct MHD turbulent flows. Further, the low- and high-Reynolds number RANS models (k-ε and Reynolds stress models) are developed and tested with DNS databases for channel and square duct flows with and without magnetic field. The MHD sink terms in k- and ε-equations are implemented as proposed by Kenjereš and Hanjalić using a user defined function (UDF) in FLUENT. This work revealed varying accuracies of different RANS models at different levels. This work is useful for industry to understand the accuracies of these models, including continuous casting. After realizing the accuracy and computational cost of RANS models, the steady-state k-ε model is then combined with the particle image velocimetry (PIV) and impeller probe velocity measurements in a 1/3rd scale water model to study the flow quality coming out of the well- and mountain-bottom nozzles and the effect of stopper-rod misalignment on fluid flow. The mountain-bottom nozzle was found more prone to the longtime asymmetries and higher surface velocities. The left misalignment of stopper gave higher surface velocity on the right leading to significantly large number of vortices forming behind the nozzle on the left. Later, the transient and steady-state models such as LES, filtered URANS and steady RANS models are combined with ultrasonic Doppler velocimetry (UDV) measurements in a GaInSn model of typical continuous casting process. LES-CU-LOW is the fastest and the most accurate model owing to much finer mesh and a smaller timestep. This work provided a good understanding on the performance of these models. The behavior of instantaneous flows, Reynolds stresses and proper orthogonal decomposition (POD) analysis quantified the nozzle bottom swirl and its importance on the turbulent flow in the mold. Afterwards, the aforementioned work in GaInSn model is extended with electromagnetic braking (EMBr) to help optimize a ruler-type brake and its location for the continuous casting process. The magnetic field suppressed turbulence and promoted vortical structures with their axis aligned with the magnetic field suggesting tendency towards 2-d turbulence. The stronger magnetic field at the nozzle well and around the jet region created large scale and lower frequency flow behavior by suppressing nozzle bottom swirl and its front-back alternation. Based on this work, it is advised to avoid stronger magnetic field around jet and nozzle bottom to get more stable and less defect prone flow.

Chaotic mixing of a competitive-consecutive reaction

The evolution of a competitive-consecutive chemical reaction is computed numerically in a two-dimensional chaotic fluid flow with initially segregated reactants. Results from numerical simulations are used to evaluate a variety of reduced models commonly adopted to model the full advection-reaction-diffusion problem. Particular emphasis is placed upon fast reactions, where the yield varies most significantly with Peclet number (the ratio of diffusive to advective time scales). When effects of the fluid mechanical mixing are strongest, we find that the yield of the reaction is underestimated by a one-dimensional lamellar model that ignores the effects of fluid mixing, but overestimated by two other lamellar models that include fluid mixing.

The study of the sedimentation mechanism in channels under the impact of tidal currents

The purpose of this research is to study sedimentation mechanism by mathematical modeling in access channels which are affected by tidal currents. The most important factor for recognizing sedimentation process in every water environment is the flow pattern of that environment. It is noteworthy that the flow pattern is affected by the geometry and the shape of the environment as well as the type of existing affects in area. The area under the study in this thesis is located in Bushehr Gulf and the access channels (inner and outer). The study utilizes the hydrodynamic modeling with unstructured triangular and non-overlapping grids, using the finite volume, From method analysis in two scale sizes: large scale (200 m to 7.5km) and small scale (50m to 7.5km) in two different time durations of 15 days and 3.5 days to obtain the flow patterns. The 2D governing equations used in the model are the Depth-Averaged Shallow Water Equations. Turbulence Modeling is required to calculate the Eddy Viscosity Coefficient using the Smagorinsky Model with coefficient of 0.3. In addition to the flow modeling in two different scales and the use of the data of 3.5 day tidal current modeling have been considered to study the effects of the sediments equilibrium in the area and the channels. This model is capable of covering the area which is being settled and eroded and to identify the effects of tidal current of these processes. The required data of the above mentioned models such as current and sediments data have been obtained by the measurements in Bushehr Gulf and the access channels which was one of the PSO's (Port and Shipping Organization) project-titled, "The Sedimentation Modeling in Bushehr Port" in 1379. Hydrographic data have been obtained from Admiralty maps (2003) and Cartography Organization (1378, 1379). The results of the modeling includes: cross shore currents in northern and north western coasts of Bushehr Gulf during the neap tide and also the same current in northern and north eastern coasts of the Gulf during the spring tide. These currents wash and carry fine particles (silt, clay, and mud) from the coastal bed of which are generally made of mud and clay with some silts. In this regard, the role of sediments in the islands of this area and the islands made of depot of dredged sediments should not be ignored. The result of using 3.5 day modeling is that the cross channels currents leads to settlement places in inner and outer channels in tidal period. In neap tide the current enters the channel from upside bend of the two channels and outer channel. Then it crosses the channel oblique in some places of the outer channel. Also the oblique currents or even almost perpendicular current from up slope of inner channel between No. 15 and No. 18 buoys interact between the parallel currents in the channel and made secondary oblique currents which exit as a down-slope current in the channel and causes deposit of sediments as well as settling the suspended sediments carried by these currents. In addition in outer channel the speed of parallel currents in the bend of the channel which is naturally deeper increases. Therefore, it leads to erosion and suspension of sediments in this area. The speed of suspended sediments carried by this current which is parallel to the channel axis decreases when they pass through the shallower part of the channel where it is in the buoys No.7 and 8 to 5 and 6 are located. Therefore, the suspended sediment settles and because of this process these places will be even shallower. Furthermore, the passing of oblique upstream leads to settlement of the sediments in the up-slope and has an additional effect on the process of decreasing the depth of these locations. On the contrary, in the down-slope channel, as the results of sediments and current modeling indicates the speed of current increases and the currents make the particles of down-slope channel suspended and be carried away. Thus, in a vast area of downstream of both channels, the sediments have settled. At the end of the neap tide, the process along with circulations in this area produces eddies which causes sedimentation in the area. During spring some parts of this active location for sedimentation will enter both channels in a reverse process. The above mentioned processes and the places of sedimentation and erosion in inner and outer channels are validated by the sediments equilibrium modeling. This model will be able to estimate the suspended, bed load and the boundary layer thickness in each point of both channels and in the modeled area.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.

Influência da lubrirrefrigeração na qualidade superficial de metais retificados

Na fabricação de componentes mecânicos precisos, que necessitam de alta resistência mecânica e ao desgaste, utiliza-se o processo de retificação, para conferir o acabamento final desejado e, também, para eliminar as deformações ocorridas durante a têmpera do aço. No entanto, as condições de retificação devem ser adequadas, para que não sejam introduzidas falhas na peça. Novos conceitos de lubrificação e refrigeração, para o processo de retificação, estão sendo pesquisados, de forma a diminuir os custos e os danos ambientais causados pelos fluidos de corte. Nesse trabalho, é analisada a influência das técnicas de mínima quantidade de lubrificante (MQL), refrigeração otimizada e refrigeração convencional, com diferentes vazões e velocidade de aplicação do fluido de corte, na qualidade das peças produzidas com aço ABNT 4340 endurecido, no processo de retificação cilíndrica externa de mergulho com a utilização de rebolos de CBN. O Aço ABNT 4340 apresenta várias aplicações industriais, sendo considerado de uso aeronáutico devido, sua alta resistência mecânica sem aumentar o peso dos componentes que o utilizam. A análise da qualidade das peças foi realizada com a verficação das rugosidades e com a análise de microscopias eletrônicas de varredura. Verificou-se, ainda, a força tangencial de corte. em relação às diferentes formas de aplicação do fluido de corte, notou-se o melhor desempenho da aplicação otimizada, para maiores velocidades, mostrando a eficiência do bocal utilizado. O processo otimizado e o processo MQL foram capazes de manter a integridade superficial das peças produzidas. Exceção somente para a condição MQL com vazão de fluido de corte de 40ml/h, que produziu trincas e queima superficial. Rebolos com baixa concentração de CBN, conseqüentemente mais baratos, proporcionaram bons resultados, quando associados com técnicas mais eficientes de aplicação de fluido de corte apresentando desgaste reduzido.

A systems approach to cardiovascular health and disease with a focus on aortic wave dynamics

Cardiovascular diseases (CVDs) have reached an epidemic proportion in the US and worldwide with serious consequences in terms of human suffering and economic impact. More than one third of American adults are suffering from CVDs. The total direct and indirect costs of CVDs are more than $500 billion per year. Therefore, there is an urgent need to develop noninvasive diagnostics methods, to design minimally invasive assist devices, and to develop economical and easy-to-use monitoring systems for cardiovascular diseases. In order to achieve these goals, it is necessary to gain a better understanding of the subsystems that constitute the cardiovascular system. The aorta is one of these subsystems whose role in cardiovascular functioning has been underestimated. Traditionally, the aorta and its branches have been viewed as resistive conduits connected to an active pump (left ventricle of the heart). However, this perception fails to explain many observed physiological results. My goal in this thesis is to demonstrate the subtle but important role of the aorta as a system, with focus on the wave dynamics in the aorta.

The operation of a healthy heart is based on an optimized balance between its pumping characteristics and the hemodynamics of the aorta and vascular branches. The delicate balance between the aorta and heart can be impaired due to aging, smoking, or disease. The heart generates pulsatile flow that produces pressure and flow waves as it enters into the compliant aorta. These aortic waves propagate and reflect from reflection sites (bifurcations and tapering). They can act constructively and assist the blood circulation. However, they may act destructively, promoting diseases or initiating sudden cardiac death. These waves also carry information about the diseases of the heart, vascular disease, and coupling of heart and aorta. In order to elucidate the role of the aorta as a dynamic system, the interplay between the dominant wave dynamic parameters is investigated in this study. These parameters are heart rate, aortic compliance (wave speed), and locations of reflection sites. Both computational and experimental approaches have been used in this research. In some cases, the results are further explained using theoretical models.

The main findings of this study are as follows: (i) developing a physiologically realistic outflow boundary condition for blood flow modeling in a compliant vasculature; (ii) demonstrating that pulse pressure as a single index cannot predict the true level of pulsatile workload on the left ventricle; (iii) proving that there is an optimum heart rate in which the pulsatile workload of the heart is minimized and that the optimum heart rate shifts to a higher value as aortic rigidity increases; (iv) introducing a simple bio-inspired device for correction and optimization of aortic wave reflection that reduces the workload on the heart; (v) deriving a non-dimensional number that can predict the optimum wave dynamic state in a mammalian cardiovascular system; (vi) demonstrating that waves can create a pumping effect in the aorta; (vii) introducing a system parameter and a new medical index, Intrinsic Frequency, that can be used for noninvasive diagnosis of heart and vascular diseases; and (viii) proposing a new medical hypothesis for sudden cardiac death in young athletes.

Adaptive coatings based on polyaniline for direct 2D observation of diffusion processes in microfluidic systems

[EN] Therefore the understanding and proper evaluation of the flow and mixing behaviour at microscale becomes a very important issue. In this study, the diffusion behaviour of two reacting solutions of HCI and NaOH were directly observed in a glass/polydimethylsiloxane microfluidic device using adaptive coatings based on the conductive polymer polyaniline that are covalently attached to the microchannel walls. The two liquid streams were combined at the junction of a Y-shaped microchannel, and allowed to diffuse into each other and react. The results showed excellent correlation between optical observation of the diffusion process and the numerical results. A numerical model which is based on finite volume method (FVM) discretisation of steady Navier-Stokes (fluid flow) equations and mass transport equations without reactions was used to calculate the flow variables at discrete points in the finite volume mesh element. The high correlation between theory and practical data indicates the potential of such coatings to monitor diffusion processes and mixing behaviour inside microfluidic channels in a dye free environment.

Inverse Spectral Methods in Acoustic Normal Mode Velocimetry of High Reynolds Number Spherical Couette Flows

Experimental geophysical fluid dynamics often examines regimes of fluid flow infeasible for computer simulations. Velocimetry of zonal flows present in these regimes brings many challenges when the fluid is opaque and vigorously rotating; spherical Couette flows with molten metals are one such example. The fine structure of the acoustic spectrum can be related to the fluid’s velocity field, and inverse spectral methods can be used to predict and, with sufficient acoustic data, mathematically reconstruct the velocity field. The methods are to some extent inherited from helioseismology. This work develops a Finite Element Method suitable to matching the geometries of experimental setups, as well as modelling the acoustics based on that geometry and zonal flows therein. As an application, this work uses the 60-cm setup Dynamo 3.5 at the University of Maryland Nonlinear Dynamics Laboratory. Additionally, results obtained using a small acoustic data set from recent experiments in air are provided.