A discrete simulation of 2-D fluid flow on TERASYS
Mullins, P.G.; Krolak, P.D.
1995-12-01
A discrete simulation of two-dimensional (2-D) fluid flow, on a recently designed novel architecture called TERASYS is presented. The simulation uses a cellular automaton approach, implemented in a new language called data-parallel bit C (dbC). A performance comparison between our implementation on TERASYS and an implementation on the Connection Machine is discussed. We comment briefly on the suitability of the TERASYS system for modeling fluid flow using cellular automata.
An Integrative Model of Excitation Driven Fluid Flow in a 2D Uterine Channel
NASA Astrophysics Data System (ADS)
Maggio, Charles; Fauci, Lisa; Chrispell, John
2009-11-01
We present a model of intra-uterine fluid flow in a sagittal cross-section of the uterus by inducing peristalsis in a 2D channel. This is an integrative multiscale computational model that takes as input fluid viscosity, passive tissue properties of the uterine channel and a prescribed wave of membrane depolarization. This voltage pulse is coupled to a model of calcium dynamics inside a uterine smooth muscle cell, which in turn drives a kinetic model of myosin phosphorylation governing contractile muscle forces. Using the immersed boundary method, these muscle forces are communicated to a fluid domain to simulate the contractions which occur in a human uterus. An analysis of the effects of model parameters on the flow properties and emergent geometry of the peristaltic channel will be presented.
A case study of fluid flow in fractured rock mass based on 2-D DFN modeling
NASA Astrophysics Data System (ADS)
Han, Jisu; Noh, Young-Hwan; Um, Jeong-Gi; Choi, Yosoon
2014-05-01
A two dimensional steady-state fluid flow through fractured rock mass of an abandoned copper mine in Korea is addressed based on discrete fracture network modeling. An injection well and three observation wells were installed at the field site to monitor the variations of total heads induced by injection of fresh water. A series of packer tests were performed to estimate the rock mass permeability. First, the two dimensional stochastic fracture network model was built and validated for a granitic rock mass using the geometrical and statistical data obtained from surface exposures and borehole logs. This validated fracture network model was combined with the fracture data observed on boreholes to generate a stochastic-deterministic fracture network system. Estimated apertures for each of the fracture sets using permeability data obtained from borehole packer tests were discussed next. Finally, a systematic procedure for fluid flow modeling in fractured rock mass in two dimensional domain was presented to estimate the conductance, flow quantity and nodal head in 2-D conceptual linear pipe channel network. The results obtained in this study clearly show that fracture geometry parameters (orientation, density and size) play an important role in the hydraulic behavior of fractured rock masses.
NASA Astrophysics Data System (ADS)
Shaari, M. F.; Abu Bakar, H.; Nordin, N.; Saw, S. K.; Samad, Z.
2013-12-01
Contractile body is an alternative mechanism instead of rotating blade propeller to generate water jet for locomotion. The oscillating motion of the actuator at different frequencies varies the pressure and volume of the pressure chamber in time to draw in and jet out the water at a certain mass flow rate. The aim of this research was to analyze the influence of the actuating frequency of the fluid flow in the pressure chamber of the thruster during this inflation-deflation process. A 70mm × 70mm × 18mm (L × W × T) 2D water jet thruster was fabricated for this purpose. The contractile function was driven using two lateral pneumatic actuators where the fluid flow analysis was focused on the X-Y plane vector. Observation was carried out using a video camera and Matlab image measurement technique to determine the volume of the flowing mass. The result demonstrated that the greater actuating frequency decreases the fluid flow rate and the Reynolds number. This observation shows that the higher frequency would give a higher mass flow rate during water jet generation.
Two-Phase Acto-Cytosolic Fluid Flow in a Moving Keratocyte: A 2D Continuum Model.
Nikmaneshi, M R; Firoozabadi, B; Saidi, M S
2015-09-01
The F-actin network and cytosol in the lamellipodia of crawling cells flow in a centripetal pattern and spout-like form, respectively. We have numerically studied this two-phase flow in the realistic geometry of a moving keratocyte. Cytosol has been treated as a low viscosity Newtonian fluid flowing through the high viscosity porous medium of F-actin network. Other involved phenomena including myosin activity, adhesion friction, and interphase interaction are also discussed to provide an overall view of this problem. Adopting a two-phase coupled model by myosin concentration, we have found new accurate perspectives of acto-cytosolic flow and pressure fields, myosin distribution, as well as the distribution of effective forces across the lamellipodia of a keratocyte with stationary shape. The order of magnitude method is also used to determine the contribution of forces in the internal dynamics of lamellipodia. PMID:26403420
NASA Astrophysics Data System (ADS)
Kim, Ho Jun; Lee, Hae June
2016-06-01
The wide applicability of capacitively coupled plasma (CCP) deposition has increased the interest in developing comprehensive numerical models, but CCP imposes a tremendous computational cost when conducting a transient analysis in a three-dimensional (3D) model which reflects the real geometry of reactors. In particular, the detailed flow features of reactive gases induced by 3D geometric effects need to be considered for the precise calculation of radical distribution of reactive species. Thus, an alternative inclusive method for the numerical simulation of CCP deposition is proposed to simulate a two-dimensional (2D) CCP model based on the 3D gas flow results by simulating flow, temperature, and species fields in a 3D space at first without calculating the plasma chemistry. A numerical study of a cylindrical showerhead-electrode CCP reactor was conducted for particular cases of SiH4/NH3/N2/He gas mixture to deposit a hydrogenated silicon nitride (SiN x H y ) film. The proposed methodology produces numerical results for a 300 mm wafer deposition reactor which agree very well with the deposition rate profile measured experimentally along the wafer radius.
Scaglione, S; Wendt, D; Miggino, S; Papadimitropoulos, A; Fato, M; Quarto, R; Martin, I
2008-08-01
In this study, we investigated the effect of the long-term (10 days) application of a defined and uniform level of fluid flow (uniform shear stress of 1.2 x 10(-3) N/m(2)) on human bone marrow stromal cells (BMSC) cultured on different substrates (i.e., uncoated glass or calcium phosphate coated glass, Osteologictrade mark) in a 2D parallel plate model. Both exposure to flow and culture on Osteologic significantly reduced the number of cell doublings. BMSC cultured under flow were more intensely stained for collagen type I and by von Kossa for mineralized matrix. BMSC exposed to flow displayed an increased osteogenic commitment (i.e., higher mRNA expression of cbfa-1 and osterix), although phenotype changes in response to flow (i.e., mRNA expression of osteopontin, osteocalcin and bone sialoprotein) were dependent on the substrate used. These findings highlight the importance of the combination of physical forces and culture substrate to determine the functional state of differentiating osteoblastic cells. The results obtained using a simple and controlled 2D model system may help to interpret the long-term effects of BMSC culture under perfusion within 3D porous scaffolds, where multiple experimental variables cannot be easily studied independently, and shear stresses cannot be precisely computed. PMID:17969030
An IPOT meshless method using DC PSE approximation for fluid flow equations in 2D and 3D geometries
NASA Astrophysics Data System (ADS)
Bourantas, G. C.; Loukopoulos, V. C.; Skouras, E. D.; Burganos, V. N.; Nikiforidis, G. C.
2016-06-01
Navier-Stokes (N-S) equations, in their primitive variable (u-v-p) formulation, are numerically solved using the Implicit Potential (IPOT) numerical scheme in the context of strong form Meshless Point Collocation (MPC) method. The unknown field functions are computed using the Discretization Correction Particle Strength Exchange (DC PSE) approximation method. The latter makes use of discrete moment conditions to derive the operator kernels, which leads to low condition number for the moment matrix compared to other meshless interpolation methods and increased stability for the numerical solution. The proposed meshless scheme is applied on 2D and 3D spatial domains, using uniform or irregular set of nodes to represent the domain. The numerical results obtained are compared against those obtained using well-established methods.
ENERGY LANDSCAPE OF 2D FLUID FORMS
Y. JIANG; ET AL
2000-04-01
The equilibrium states of 2D non-coarsening fluid foams, which consist of bubbles with fixed areas, correspond to local minima of the total perimeter. (1) The authors find an approximate value of the global minimum, and determine directly from an image how far a foam is from its ground state. (2) For (small) area disorder, small bubbles tend to sort inwards and large bubbles outwards. (3) Topological charges of the same sign repel while charges of opposite sign attract. (4) They discuss boundary conditions and the uniqueness of the pattern for fixed topology.
Rheological Properties of Quasi-2D Fluids in Microgravity
NASA Technical Reports Server (NTRS)
Stannarius, Ralf; Trittel, Torsten; Eremin, Alexey; Harth, Kirsten; Clark, Noel; Maclennan, Joseph; Glaser, Matthew; Park, Cheol; Hall, Nancy; Tin, Padetha
2015-01-01
In recent years, research on complex fluids and fluids in restricted geometries has attracted much attention in the scientific community. This can be attributed not only to the development of novel materials based on complex fluids but also to a variety of important physical phenomena which have barely been explored. One example is the behavior of membranes and thin fluid films, which can be described by two-dimensional (2D) rheology behavior that is quite different from 3D fluids. In this study, we have investigated the rheological properties of freely suspended films of a thermotropic liquid crystal in microgravity experiments. This model system mimics isotropic and anisotropic quasi 2D fluids [46]. We use inkjet printing technology to dispense small droplets (inclusions) onto the film surface. The motion of these inclusions provides information on the rheological properties of the films and allows the study of a variety of flow instabilities. Flat films have been investigated on a sub-orbital rocket flight and curved films (bubbles) have been studied in the ISS project OASIS. Microgravity is essential when the films are curved in order to avoid sedimentation. The experiments yield the mobility of the droplets in the films as well as the mutual mobility of pairs of particles. Experimental results will be presented for 2D-isotropic (smectic-A) and 2D-nematic (smectic-C) phases.
Lappala, E.G.; Healy, R.W.; Weeks, E.P.
1987-01-01
This report documents FORTRAN computer code for solving problems involving variably saturated single-phase flow in porous media. The flow equation is written with total hydraulic potential as the dependent variable, which allows straightforward treatment of both saturated and unsaturated conditions. The spatial derivatives in the flow equation are approximated by central differences, and time derivatives are approximated either by a fully implicit backward or by a centered-difference scheme. Nonlinear conductance and storage terms may be linearized using either an explicit method or an implicit Newton-Raphson method. Relative hydraulic conductivity is evaluated at cell boundaries by using either full upstream weighting, the arithmetic mean, or the geometric mean of values from adjacent cells. Nonlinear boundary conditions treated by the code include infiltration, evaporation, and seepage faces. Extraction by plant roots that is caused by atmospheric demand is included as a nonlinear sink term. These nonlinear boundary and sink terms are linearized implicitly. The code has been verified for several one-dimensional linear problems for which analytical solutions exist and against two nonlinear problems that have been simulated with other numerical models. A complete listing of data-entry requirements and data entry and results for three example problems are provided. (USGS)
NASA Astrophysics Data System (ADS)
Baudon, Catherine; Gillet, Hervé; Cremer, Michel
2013-04-01
High-quality bathymetric, 2D seismic and Chirp data located in the southern parts of the Bay of Biscay, France, collected by the University of Bordeaux 1 (Cruises ITSAS 2, 2001; PROSECAN 3, 2006 and SARGASS, 2010) have recently been compiled. The survey area widely covers the Capbreton Canyon, which lies on the boundary between two major structural zones: the Aquitanian passive margin to the North, and the Basque-Cantabrian margin to the South which corresponds to the offshore Pyrenean front. The dataset revealed a large number of key seafloor features potentially associated with focused fluid-flow processes and subsurface sediment-remobilization. Focused fluid migration through sub-seabed sediments is a common phenomenon on continental margins worldwide and has widespread implications from both industrial and fundamental perspectives, from seafloor marine environmental issues to petroleum exploration and hazard assessments. Our study analyses the relationships between seafloor features, deeper structures and fluid migration through the Plio-Quaternary sedimentary pile. The geometrical characteristics, mechanisms of formation and kinematics of four main groups of seabed features have been investigated. (i) A 150km2 field of pockmarks can be observed on the Basque margin. These features are cone-shaped circular or elliptical depressions that are either randomly distributed as small pockmarks (diameter < 20m) or aligned in trains of large pockmarks (ranging from 200 to 600m in diameter) along shallow troughs leading downstream to the Capbreton Canyon. Seismic data show that most pockmarks reach the seabed through vertically staked V-shaped features but some are buried and show evidence of lateral migration through time. (ii) A second field of widely-spaced groups of pockmarks pierce the upper slope of the Aquitanian margin. These depressions are typically a few hundred meters in diameter and seem to be preferentially located in the troughs or on the stoss sides of
A 2D electrohydrodynamic model for electrorotation of fluid drops.
Feng, James Q
2002-02-01
A theoretical analysis of spontaneous electrorotation of deformable fluid drops in a DC electric field is presented with a 2D electrohydrodynamic model. The fluids in the system are assumed to be leaky dielectric and Newtonian. If the rotating flow is dominant over the cellular convection type of electrohydrodynamic flow, closed-form solutions for drops of small deformations can be obtained. Because the governing equations are in general nonlinear even when drop deformations are ignored, the general solution for even undeformed drop takes a form of infinite series and can only be evaluated by numerical means. Both closed-form solutions for special cases and numerical solutions for more general cases are obtained here to describe steady-state field variables and first-order drop deformations. In a DC electric field of strength beyond the threshold value, spontaneous electrorotation of a drop is shown to occur when charge relaxation in the surrounding fluid is faster than the fluid inside the drop. With increasing the strength of the applied electric field from the threshold for onset of electrorotation, the axis of drop contraction deviates from from that of the applied electric field in the direction of the rotating flow with an angle increasing with the field strength. PMID:16290391
Collective excitations in 2D hard-disc fluid.
Huerta, Adrian; Bryk, Taras; Trokhymchuk, Andrij
2015-07-01
Collective dynamics of a two-dimensional (2D) hard-disc fluid was studied by molecular dynamics simulations in the range of packing fractions that covers states up to the freezing. Some striking features concerning collective excitations in this system were observed. In particular, the short-wavelength shear waves while being absent at low packing fractions were observed in the range of high packing fractions, just before the freezing transition in a 2D hard-disc fluid. In contrast, the so-called "positive sound dispersion" typically observed in dense Lennard-Jones-like fluids, was not detected for the 2D hard-disc fluid. The ratio of specific heats in the 2D hard-disc fluid shows a monotonic increase with density approaching the freezing, resembling in this way the similar behavior in the vicinity of the Widom line in the case of supercritical fluids. PMID:25595625
A new inversion method for (T2, D) 2D NMR logging and fluid typing
NASA Astrophysics Data System (ADS)
Tan, Maojin; Zou, Youlong; Zhou, Cancan
2013-02-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology has some significant limitations in fluid typing. However, not only can two-dimensional nuclear magnetic resonance (2D NMR) provide some accurate porosity parameters, but it can also identify fluids more accurately than 1D NMR. In this paper, based on the relaxation mechanism of (T2, D) 2D NMR in a gradient magnetic field, a hybrid inversion method that combines least-squares-based QR decomposition (LSQR) and truncated singular value decomposition (TSVD) is examined in the 2D NMR inversion of various fluid models. The forward modeling and inversion tests are performed in detail with different acquisition parameters, such as magnetic field gradients (G) and echo spacing (TE) groups. The simulated results are discussed and described in detail, the influence of the above-mentioned observation parameters on the inversion accuracy is investigated and analyzed, and the observation parameters in multi-TE activation are optimized. Furthermore, the hybrid inversion can be applied to quantitatively determine the fluid saturation. To study the effects of noise level on the hybrid method and inversion results, the numerical simulation experiments are performed using different signal-to-noise-ratios (SNRs), and the effect of different SNRs on fluid typing using three fluid models are discussed and analyzed in detail.
Flow past 2-D Hemispherical Rigid Canopies
NASA Astrophysics Data System (ADS)
Carnasciali, Maria-Isabel
2013-11-01
The flow past a 2-dimensional rigid hemispherical shape is investigated using PIV. Flow field measurements and images were generated with the use of a Thermoflow® apparatus. Results of this study are compared to prior work (APS DFD 2012 Session E9.00003) which employed CFD to investigate the flow in the near wake of hemispherical parachutes. The various sized gaps/open areas were positioned at distinct locations. The work presented here is part of a larger research project to investigate flow fields in deceleration devices and parachutes. Understanding the pitch-stability of parachutes is essential for accurate design and implementation of these deceleration devices but they present a difficult system to analyze. The flexibility of the parachute fabric results in large variations in the parachute geometry leading to complex fluid-structure interactions. Such flow, combined with flow through gaps and open areas, has been postulated to shed alternating vortices causing pitching/oscillations of the canopy. The results presented provide some insight into which geometric features affect vortex shedding and may enable the redesign of the baseline parachute to minimize instabilities.
Topological phase transition in 2D porous media flows
NASA Astrophysics Data System (ADS)
Waisbord, Nicolas; Stoop, Norbert; Kantsler, Vasily; Guasto, Jeffrey S.; Dunkel, Jorn; Guasto Team; Dunkel Team; Kantsler Team
2015-11-01
Since the establishment of Darcy's law, analysis of porous-media flows has focused primarily on linking macroscopic transport properties, such as mean flow rate and dispersion, to the pore statistics of the material matrix. Despite intense efforts to understand the fluid velocity statistics from the porous-media structure, a qualitative and quantitative connection remains elusive. Here, we combine precisely controlled experiments with theory to quantify how geometric disorder in the matrix affects the flow statistics and transport in a quasi-2D microfluidic channel. Experimentally measured velocity fields for a range of different microstructure configurations are found to be in excellent agreement with large-scale numerical simulations. By successively increasing the matrix disorder, we study the transition from periodic flow structures to transport networks consisting of extended high-velocity channels. Morse-Smale complex analysis of the flow patterns reveals a topological phase transition that is linked to a qualitative change in the physical transport properties. This work demonstrates that topological flow analysis provides a mathematically well-defined, broadly applicable framework for understanding and quantifying fluid transport in complex geometries.
Advecting Procedural Textures for 2D Flow Animation
NASA Technical Reports Server (NTRS)
Kao, David; Pang, Alex; Moran, Pat (Technical Monitor)
2001-01-01
This paper proposes the use of specially generated 3D procedural textures for visualizing steady state 2D flow fields. We use the flow field to advect and animate the texture over time. However, using standard texture advection techniques and arbitrary textures will introduce some undesirable effects such as: (a) expanding texture from a critical source point, (b) streaking pattern from the boundary of the flowfield, (c) crowding of advected textures near an attracting spiral or sink, and (d) absent or lack of textures in some regions of the flow. This paper proposes a number of strategies to solve these problems. We demonstrate how the technique works using both synthetic data and computational fluid dynamics data.
NASA Astrophysics Data System (ADS)
Juez, C.; Caviedes-Voullième, D.; Murillo, J.; García-Navarro, P.
2014-12-01
Dense granular flows are present in geophysics and in several industrial processes, which has lead to an increasing interest for the knowledge and understanding of the physics which govern their propagation. For this reason, a wide range of laboratory experiments on gravity-driven flows have been carried out during the last two decades. The present work is focused on geomorphological processes and, following previous work, a series of laboratory studies which constitute a further step in mimicking natural phenomena are described and simulated. Three situations are considered with some common properties: a two-dimensional configuration, variable slope of the topography and the presence of obstacles. The setup and measurement technique employed during the development of these experiments are deeply explained in the companion work. The first experiment is based on a single obstacle, the second one is performed against multiple obstacles and the third one studies the influence of a dike on which overtopping occurs. Due to the impact of the flow against the obstacles, fast moving shocks appear, and a variety of secondary waves emerge. In order to delve into the physics of these types of phenomena, a shock-capturing numerical scheme is used to simulate the cases. The suitability of the mathematical models employed in this work has been previously validated. Comparisons between computed and experimental data are presented for the three cases. The computed results show that the numerical tool is able to predict faithfully the overall behavior of this type of complex dense granular flow.
NASA Astrophysics Data System (ADS)
Erwee, M. W.; Reynolds, Q. G.; Zietsman, J. H.
2016-03-01
Furnace tap-holes vary in design depending on the type of furnace and process involved, but they share one common trait: The tap-hole must be opened and closed periodically. In general, tap-holes are plugged with refractory clay after tapping, thereby stopping the flow of molten material. Once a furnace is ready to be tapped, drilling and/or lancing with oxygen are typically used to remove tap-hole clay from the tap-hole. Lancing with oxygen is an energy-intensive, mostly manual process, which affects the performance and longevity of the tap-hole refractory material as well as the processes inside the furnace. Computational modeling offers an opportunity to gain insight into the possible effects of oxygen lancing on various aspects of furnace operation.
NASA Astrophysics Data System (ADS)
Erwee, M. W.; Reynolds, Q. G.; Zietsman, J. H.
2016-06-01
Furnace tap-holes vary in design depending on the type of furnace and process involved, but they share one common trait: The tap-hole must be opened and closed periodically. In general, tap-holes are plugged with refractory clay after tapping, thereby stopping the flow of molten material. Once a furnace is ready to be tapped, drilling and/or lancing with oxygen are typically used to remove tap-hole clay from the tap-hole. Lancing with oxygen is an energy-intensive, mostly manual process, which affects the performance and longevity of the tap-hole refractory material as well as the processes inside the furnace. Computational modeling offers an opportunity to gain insight into the possible effects of oxygen lancing on various aspects of furnace operation.
Mean flow and anisotropic cascades in decaying 2D turbulence
NASA Astrophysics Data System (ADS)
Liu, Chien-Chia; Cerbus, Rory; Gioia, Gustavo; Chakraborty, Pinaki
2015-11-01
Many large-scale atmospheric and oceanic flows are decaying 2D turbulent flows embedded in a non-uniform mean flow. Despite its importance for large-scale weather systems, the affect of non-uniform mean flows on decaying 2D turbulence remains unknown. In the absence of mean flow it is well known that decaying 2D turbulent flows exhibit the enstrophy cascade. More generally, for any 2D turbulent flow, all computational, experimental and field data amassed to date indicate that the spectrum of longitudinal and transverse velocity fluctuations correspond to the same cascade, signifying isotropy of cascades. Here we report experiments on decaying 2D turbulence in soap films with a non-uniform mean flow. We find that the flow transitions from the usual isotropic enstrophy cascade to a series of unusual and, to our knowledge, never before observed or predicted, anisotropic cascades where the longitudinal and transverse spectra are mutually independent. We discuss implications of our results for decaying geophysical turbulence.
Billeter, Thomas R.; Philipp, Lee D.; Schemmel, Richard R.
1976-01-01
A microwave fluid flow meter is described utilizing two spaced microwave sensors positioned along a fluid flow path. Each sensor includes a microwave cavity having a frequency of resonance dependent upon the static pressure of the fluid at the sensor locations. The resonant response of each cavity with respect to a variation in pressure of the monitored fluid is represented by a corresponding electrical output which can be calibrated into a direct pressure reading. The pressure drop between sensor locations is then correlated as a measure of fluid velocity. In the preferred embodiment the individual sensor cavities are strategically positioned outside the path of fluid flow and are designed to resonate in two distinct frequency modes yielding a measure of temperature as well as pressure. The temperature response can then be used in correcting for pressure responses of the microwave cavity encountered due to temperature fluctuations.
Geophysical fluid flow experiment
NASA Technical Reports Server (NTRS)
Broome, B. G.; Fichtl, G.; Fowlis, W.
1979-01-01
The essential fluid flow processes associated with the solar and Jovian atmospheres will be examined in a laboratory experiment scheduled for performance on Spacelab Missions One and Three. The experimental instrumentation required to generate and to record convective fluid flow is described. Details of the optical system configuration, the lens design, and the optical coatings are described. Measurement of thermal gradient fields by schlieren techniques and measurement of fluid flow velocity fields by photochromic dye tracers is achieved with a common optical system which utilizes photographic film for data recording. Generation of the photochromic dye tracers is described, and data annotation of experimental parameters on the film record is discussed.
Validation of a 2-D semi-coupled numerical model for fluid-structure-seabed interaction
NASA Astrophysics Data System (ADS)
Ye, Jianhong; Jeng, Dongsheng; Wang, Ren; Zhu, Changqi
2013-10-01
A 2-D semi-coupled model PORO-WSSI 2D (also be referred as FSSI-CAS 2D) for the Fluid-Structure-Seabed Interaction (FSSI) has been developed by employing RANS equations for wave motion in fluid domain, VARANS equations for porous flow in porous structures; and taking the dynamic Biot's equations (known as "u - p" approximation) for soil as the governing equations. The finite difference two-step projection method and the forward time difference method are adopted to solve the RANS, VARANS equations; and the finite element method is adopted to solve the "u - p" approximation. A data exchange port is developed to couple the RANS, VARANS equations and the dynamic Biot's equations together. The analytical solution proposed by Hsu and Jeng (1994) and some experiments conducted in wave flume or geotechnical centrifuge in which various waves involved are used to validate the developed semi-coupled numerical model. The sandy bed involved in these experiments is poro-elastic or poro-elastoplastic. The inclusion of the interaction between fluid, marine structures and poro-elastoplastic seabed foundation is a special point and highlight in this paper, which is essentially different with other previous coupled models The excellent agreement between the numerical results and the experiment data indicates that the developed coupled model is highly reliablefor the FSSI problem.
2D vortex interaction in a non-uniform flow
NASA Astrophysics Data System (ADS)
Perrot, Xavier; Carton, X.
2010-03-01
In a two-dimensional incompressible fluid, we study the interaction of two like-signed Rankine vortices embedded in a steady shear/strain flow. The numerical results of vortex evolutions are compared with the analytical results for point vortices. We show the existence of vortex equilibria, and of merger for initial distances larger than those without external flow. The evolutions depend on the initial orientation of the vortices in the external flow.
Hierarchical streamline bundles for visualizing 2D flow fields.
Shene, Ching-Kuang; Wang, Chaoli; Yu, Hongfeng; Chen, Jacqueline H.
2010-08-01
We present hierarchical streamline bundles, a new approach to simplifying and visualizing 2D flow fields. Our method first densely seeds a flow field and produces a large number of streamlines that capture important flow features such as critical points. Then, we group spatially neighboring and geometrically similar streamlines to construct a hierarchy from which we extract streamline bundles at different levels of detail. Streamline bundles highlight multiscale flow features and patterns through a clustered yet non-cluttered display. This selective visualization strategy effectively accentuates visual foci and therefore is able to convey the desired insight into the flow fields. The hierarchical streamline bundles we have introduced offer a new way to characterize and visualize the flow structure and patterns in multiscale fashion. Streamline bundles highlight critical points clearly and concisely. Exploring the hierarchy allows a complete visualization of important flow features. Thanks to selective streamline display and flexible LOD refinement, our multiresolution technique is scalable and is promising for viewing large and complex flow fields. In the future, we would like to seek a cost-effective way to generate streamlines without enforcing the dense seeding condition. We will also extend this approach to handle real-world 3D complex flow fields.
Magnetically stimulated fluid flow patterns
Martin, Jim; Solis, Kyle
2014-08-06
Sandia National Laboratories' Jim Martin and Kyle Solis explain research on the effects of magnetic fields on fluid flows and how they stimulate vigorous flows. Fluid flow is a necessary phenomenon in everything from reactors to cooling engines in cars.
Magnetically stimulated fluid flow patterns
Martin, Jim; Solis, Kyle
2014-03-06
Sandia National Laboratories' Jim Martin and Kyle Solis explain research on the effects of magnetic fields on fluid flows and how they stimulate vigorous flows. Fluid flow is a necessary phenomenon in everything from reactors to cooling engines in cars.
General Transient Fluid Flow Algorithm
Energy Science and Technology Software Center (ESTSC)
1992-03-12
SALE2D calculates two-dimensional fluid flows at all speeds, from the incompressible limit to highly supersonic. An implicit treatment of the pressure calculation similar to that in the Implicit Continuous-fluid Eulerian (ICE) technique provides this flow speed flexibility. In addition, the computing mesh may move with the fluid in a typical Lagrangian fashion, be held fixed in an Eulerian manner, or move in some arbitrarily specified way to provide a continuous rezoning capability. This latitude resultsmore » from use of an Arbitrary Lagrangian-Eulerian (ALE) treatment of the mesh. The partial differential equations solved are the Navier-Stokes equations and the mass and internal energy equations. The fluid pressure is determined from an equation of state and supplemented with an artificial viscous pressure for the computation of shock waves. The computing mesh consists of a two-dimensional network of quadrilateral cells for either cylindrical or Cartesian coordinates, and a variety of user-selectable boundary conditions are provided in the program.« less
2D Mixed Convection Thermal Incompressible Viscous Flows
NASA Astrophysics Data System (ADS)
Bermudez, Blanca; Nicolas, Alfredo
2005-11-01
Mixed convection thermal incomprressible viscous fluid flows in rectangular cavities are presented. These kind of flows may be governed by the time-dependent Boussinesq approximation in terms of the stream function-vorticity variables formulation. The results are obtained with a simple numerical scheme based mainly on a fixed point iterative process applied to the non-linear system of elliptic equations that is obtained after a second order time discretization. Numerical experiments are reported for the problem of a cavity with fluid boundary motion on the top. Some results correspond to validation examples and others, to the best of our knowledge, correspond to new results. To show that the new results are correct, a mesh size and time independence studies are carried out, and the acceptable errors are measured point-wise. For the optimal mesh size and time step the final times when the steady state is reached, as solution from the unsteady problem, are reported; it should be seen that they are larger than the ones for natural convection which, physically speaking, show the agreement that mixed convection flows are more active than those of natural convection due to the fluid boundary motion on the top of the cavity. The flow parameters are: the Reynolds number, the Grashof number and the aspect ratio.
CAS2D- NONROTATING BLADE-TO-BLADE, STEADY, POTENTIAL TRANSONIC CASCADE FLOW ANALYSIS CODE
NASA Technical Reports Server (NTRS)
Dulikravich, D. S.
1994-01-01
An exact, full-potential-equation model for the steady, irrotational, homoentropic, and homoenergetic flow of a compressible, inviscid fluid through a two-dimensional planar cascade together with its appropriate boundary conditions has been derived. The CAS2D computer program numerically solves an artificially time-dependent form of the actual full-potential-equation, providing a nonrotating blade-to-blade, steady, potential transonic cascade flow analysis code. Comparisons of results with test data and theoretical solutions indicate very good agreement. In CAS2D, the governing equation is discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field is discretized by providing a boundary-fitted, nonuniform computational mesh. This mesh is generated by using a sequence of conformal mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the full-potential equation is solved iteratively by using successive line over relaxation. Possible isentropic shocks are captured by the explicit addition of an artificial viscosity in a conservative form. In addition, a four-level, consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two-dimensional cascade flows. The results from CAS2D are not directly applicable to three-dimensional, potential, rotating flows through a cascade of blades because CAS2D does not consider the effects of the Coriolis force that would be present in the three-dimensional case. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 200K of 8 bit bytes. The CAS2D program was developed in 1980.
2-D traveling-wave patterns in binary fluid convection
Surko, C.M.; Porta, A.L.
1996-12-31
An overview is presented of recent experiments designed to study two-dimensional traveling-wave convection in binary fluid convection in a large aspect ratio container. Disordered patterns are observed when convection is initiated. As time proceeds, they evolve to more ordered patterns, consisting of several domains of traveling-waves separated by well-defined domain boundaries. The detailed character of the patterns depends sensitively on the Rayleigh number. Numerical techniques are described which were developed to provide a quantitative characterization of the traveling-wave patterns. Applications of complex demodulation techniques are also described, which make a detailed study of the structure and dynamics of the domain boundaries possible.
Improvement of a 2D numerical model of lava flows
NASA Astrophysics Data System (ADS)
Ishimine, Y.
2013-12-01
I propose an improved procedure that reduces an improper dependence of lava flow directions on the orientation of Digital Elevation Model (DEM) in two-dimensional simulations based on Ishihara et al. (in Lava Flows and Domes, Fink, JH eds., 1990). The numerical model for lava flow simulations proposed by Ishihara et al. (1990) is based on two-dimensional shallow water model combined with a constitutive equation for a Bingham fluid. It is simple but useful because it properly reproduces distributions of actual lava flows. Thus, it has been regarded as one of pioneer work of numerical simulations of lava flows and it is still now widely used in practical hazard prediction map for civil defense officials in Japan. However, the model include an improper dependence of lava flow directions on the orientation of DEM because the model separately assigns the condition for the lava flow to stop due to yield stress for each of two orthogonal axes of rectangular calculating grid based on DEM. This procedure brings a diamond-shaped distribution as shown in Fig. 1 when calculating a lava flow supplied from a point source on a virtual flat plane although the distribution should be circle-shaped. To improve the drawback, I proposed a modified procedure that uses the absolute value of yield stress derived from both components of two orthogonal directions of the slope steepness to assign the condition for lava flows to stop. This brings a better result as shown in Fig. 2. Fig. 1. (a) Contour plots calculated with the original model of Ishihara et al. (1990). (b) Contour plots calculated with a proposed model.
Gradient-Driven Vortex Motion in Nonneutral Plasmas and Ideal 2D Fluids
NASA Astrophysics Data System (ADS)
Schecter, David A.
2000-10-01
Two-dimensional (2D) turbulent flows can relax to metastable patterns without dissipation of kinetic energy. This ``rapid'' relaxation has been observed in computer simulations of ideal 2D fluids, and more recently in experiments with pure electron plasmas, which can obey similar dynamics. The late stage of relaxation often involves small vortices moving in a larger ``background'' shear-flow.(X.P. Huang et al., Phys. Rev. Lett. 74), 4424 (1995). In time, positive vortices (rotating counter-clockwise) move to peaks in background vorticity, whereas negative vortices (rotating clockwise) move to minima.(C.G. Rossby, J. Mar. Res. 7), 175 (1948); C.H. Liu and L. Ting, Comp. & Fluids 15, 77 (1987). In general, the rate of this migration increases with the magnitude of the background vorticity gradient, whereas it decreases as the background shear intensifies.\\vspace12pt Positive and negative vortices can also be classified as either prograde or retrograde, depending on whether they rotate with or against the local background shear. Surprisingly, a retrograde vortex moves up or down a background vorticity gradient orders of magnitude faster than a prograde vortex of equal strength.(D.A. Schecter and D.H.E. Dubin, Phys. Rev. Lett. 83), 2191 (1999). An accurate expression for the velocity of a weak retrograde vortex is obtained from an analytic calculation, in which the response of the background flow to the vortex is linearized. However, this linear theory fails for prograde vortices of any strength. Interestingly, the velocity of a prograde vortex can be obtained from a simple estimate, which accounts for the nonlinear ``trapping'' of background fluid around the vortex. The analytic expressions for the velocities of both prograde and retrograde vortices are in good quantitative agreement with vortex-in-cell simulations, and with electron plasma experiments, when the background shear is below a critical level. When the ratio of background shear to background vorticity
Experimental and Computational Study of Multiphase Flow Hydrodynamics in 2D Trickle Bed Reactors
NASA Astrophysics Data System (ADS)
Nadeem, H.; Ben Salem, I.; Kurnia, J. C.; Rabbani, S.; Shamim, T.; Sassi, M.
2014-12-01
Trickle bed reactors are largely used in the refining processes. Co-current heavy oil and hydrogen gas flow downward on catalytic particle bed. Fine particles in the heavy oil and/or soot formed by the exothermic catalytic reactions deposit on the bed and clog the flow channels. This work is funded by the refining company of Abu Dhabi and aims at mitigating pressure buildup due to fine deposition in the TBR. In this work, we focus on meso-scale experimental and computational investigations of the interplay between flow regimes and the various parameters that affect them. A 2D experimental apparatus has been built to investigate the flow regimes with an average pore diameter close to the values encountered in trickle beds. A parametric study is done for the development of flow regimes and the transition between them when the geometry and arrangement of the particles within the porous medium are varied. Liquid and gas flow velocities have also been varied to capture the different flow regimes. Real time images of the multiphase flow are captured using a high speed camera, which were then used to characterize the transition between the different flow regimes. A diffused light source was used behind the 2D Trickle Bed Reactor to enhance visualizations. Experimental data shows very good agreement with the published literature. The computational study focuses on the hydrodynamics of multiphase flow and to identify the flow regime developed inside TBRs using the ANSYS Fluent Software package. Multiphase flow inside TBRs is investigated using the "discrete particle" approach together with Volume of Fluid (VoF) multiphase flow modeling. The effect of the bed particle diameter, spacing, and arrangement are presented that may be used to provide guidelines for designing trickle bed reactors.
Lagrangian statistics and flow topology in forced 2-D turbulence
Kadoch, B.; Del-Castillo-Negrete, Diego B; Bos, W.J.T.; Schneider, Kai
2011-01-01
A study of the relationship between Lagrangian statistics and flow topology in fluid turbulence is presented. The topology is characterized using the Weiss criterion, which provides a conceptually simple tool to partition the flow into topologically different regions: elliptic (vortex dominated), hyperbolic (deformation dominated), and intermediate (turbulent background). The flow corresponds to forced two-dimensional Navier-Stokes turbulence in doubly periodic and circular bounded domains, the latter with no-slip boundary conditions. In the double periodic domain, the probability density function (pdf) of the Weiss field exhibits a negative skewness consistent with the fact that in periodic domains the flow is dominated by coherent vortex structures. On the other hand, in the circular domain, the elliptic and hyperbolic regions seem to be statistically similar. We follow a Lagrangian approach and obtain the statistics by tracking large ensembles of passively advected tracers. The pdfs of residence time in the topologically different regions are computed introducing the Lagrangian Weiss field, i.e., the Weiss field computed along the particles' trajectories. In elliptic and hyperbolic regions, the pdfs of the residence time have self-similar algebraic decaying tails. In contrast, in the intermediate regions the pdf has exponential decaying tails. The conditional pdfs (with respect to the flow topology) of the Lagrangian velocity exhibit Gaussian-like behavior in the periodic and in the bounded domains. In contrast to the freely decaying turbulence case, the conditional pdfs of the Lagrangian acceleration in forced turbulence show a comparable level of intermittency in both the periodic and the bounded domains. The conditional pdfs of the Lagrangian curvature are characterized, in all cases, by self-similar power-law behavior with a decay exponent of order - 2.
In situ fluid typing and quantification with 1D and 2D NMR logging.
Sun, Boqin
2007-05-01
In situ nuclear magnetic resonance (NMR) fluid typing has recently gained momentum due to data acquisition and inversion algorithm enhancement of NMR logging tools. T(2) distributions derived from NMR logging contain information on bulk fluids and pore size distributions. However, the accuracy of fluid typing is greatly overshadowed by the overlap between T(2) peaks arising from different fluids with similar apparent T(2) relaxation times. Nevertheless, the shapes of T(2) distributions from different fluid components are often different and can be predetermined. Inversion with predetermined T(2) distributions allows us to perform fluid component decomposition to yield individual fluid volume ratios. Another effective method for in situ fluid typing is two-dimensional (2D) NMR logging, which results in proton population distribution as a function of T(2) relaxation time and fluid diffusion coefficient (or T(1) relaxation time). Since diffusion coefficients (or T(1) relaxation time) for different fluid components can be very different, it is relatively easy to separate oil (especially heavy oil) from water signal in a 2D NMR map and to perform accurate fluid typing. Combining NMR logging with resistivity and/or neutron/density logs provides a third method for in situ fluid typing. We shall describe these techniques with field examples. PMID:17466778
Turbulent flow over a surface-mounted 2-D block in thermally-stratified boundary layers
NASA Astrophysics Data System (ADS)
Zhang, W.; Markfort, C. D.; Porte-Agel, F.
2013-12-01
Turbulent boundary-layer flows over complex topography have been of great interest in the atmospheric sciences and wind engineering communities. The geometry of the topography, surface characteristics and atmospheric thermal stability play important roles in determining momentum and scalar flux distribution. Studies of turbulent flow over simplified topography, such as 2-D or 3-D blocks and 2-D or 3-D sinusoidal hills, conducted under neutrally stratified boundary-layer conditions have provided insightful information of fluid dynamics. However, atmospheric thermal stability has rarely been incorporated into laboratory simulations, in particular, wind-tunnel experiments. Extension of such studies in thermally-stratified wind tunnels will fill this gap and advance our understanding of the underlying physics of flow over complex topography. Additionally, experimental data are useful for the development of new parameterizations for surface fluxes and validation of numerical models such as Large-Eddy Simulation (LES). A series of experiments involving neutral and thermally-stratified boundary-layer flows over a surface-mounted 2-D block, conducted at the Saint Anthony Falls Laboratory boundary-layer wind tunnel, will be presented. The 2-D block, with a width to height ratio of 2:1, occupied the lowest 25% of the turbulent boundary layer. Thermal stratification of the boundary layer was achieved by independently controlling the temperature of both the airflow, the test section floor and block surfaces. Measurements using high-resolution PIV, x-wire/cold-wire anemometry, thermal-couples and surface heat flux sensors were made to identify and quantify the turbulent flow properties, including the size of the recirculation zone, coherent vortex structures and the subsequent boundary layer recovery. Emphasis will be put on addressing thermal stability effects on momentum and scalar flux distribution.
Numerical simulation of ( T 2, T 1) 2D NMR and fluid responses
NASA Astrophysics Data System (ADS)
Tan, Mao-Jin; Zou, You-Long; Zhang, Jin-Yan; Zhao, Xin
2012-12-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time ( T 1) and transverse relaxation time ( T 2) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times ( T W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings ( T E s) and T W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.
Eulerian and Lagrangian methods for vortex tracking in 2D and 3D flows
NASA Astrophysics Data System (ADS)
Huang, Yangzi; Green, Melissa
2014-11-01
Coherent structures are a key component of unsteady flows in shear layers. Improvement of experimental techniques has led to larger amounts of data and requires of automated procedures for vortex tracking. Many vortex criteria are Eulerian, and identify the structures by an instantaneous local swirling motion in the field, which are indicated by closed or spiral streamlines or pathlines in a reference frame. Alternatively, a Lagrangian Coherent Structures (LCS) analysis is a Lagrangian method based on the quantities calculated along fluid particle trajectories. In the current work, vortex detection is demonstrated on data from the simulation of two cases: a 2D flow with a flat plate undergoing a 45 ° pitch-up maneuver and a 3D wall-bounded turbulence channel flow. Vortices are visualized and tracked by their centers and boundaries using Γ1, the Q criterion, and LCS saddle points. In the cases of 2D flow, saddle points trace showed a rapid acceleration of the structure which indicates the shedding from the plate. For channel flow, saddle points trace shows that average structure convection speed exhibits a similar trend as a function of wall-normal distance as the mean velocity profile, and leads to statistical quantities of vortex dynamics. Dr. Jeff Eldredge and his research group at UCLA are gratefully acknowledged for sharing the database of simulation for the current research. This work was supported by the Air Force Office of Scientific Research under AFOSR Award No. FA9550-14-1-0210.
2-D Three Fluid Simulation of Upstreaming Ions Above Auroral Precipitation
NASA Astrophysics Data System (ADS)
Danielides, M. A.; Lummerzheim, D.; Otto, A.; Stevens, R. J.
2006-12-01
The ionosphere is a rich reservoir of charged particles from which a variable fraction is transported to the magnetosphere. An important transport phenomena is the formation of upward ion flow above auroral structure. A primary region of the outflow is not known, but contributions come from polar cap, dayside cusp/cleft region, auroral oval, or even from mid-latitudes. In the past global magnetospheric models and fluid codes were used to simulate large scale ion outflow above, e.g., the polar-cap aurora. However, satellites orbiting at low- altitudes have repeatingly detected localized ion outflow above the auroral oval. Ionosphere-magnetosphere coupling simulations gave first insides into the small-scale dynamics of aurora. The aim of this study is the investigation of coupled plasma and neutral dynamics in smaller scale aurora to explain the generation, structure, and dynamics of vertical ion upstream. We consider auroral electron precipitation at ionospheric heights in a 2-D three fluid ionospheric-magnetospheric coupling code (Otto and Zhu, 2003). Specially we examine the effects of the electron precipitation, heat conduction and heating in field- aligned current through coulomb collisions or turbulence causing: i) electron heating, ii) electron pressure gradients, and iii) upstreaming of ions through a resulting ambipolar electric field. Our first case studies are performed for different boundary conditions and for different auroral electron precipitation parameters (variation in characteristic auroral energy, auroral energy flux and horizontal scale). The results shall clarify how auroral precipitation can drive ions upwards. Finally we discuss the effect of ion drag and the interaction of the upstreaming ions with a stable neutral constituent. Otto, O. and H. Zhu, Fluid plasma simulation of coupled systems: Ionosphere and magnetosphere, Space Plasma Simulation. Edited by J. Buechner, C. Dum, and M. Scholer., Lecture Notes in Physics, vol. 615, p.193
McKay, M.D.; Sweeney, C.E.; Spangler, B.S. Jr.
1993-11-30
A flow meter and temperature measuring device are described comprising a tube with a body centered therein for restricting flow and a sleeve at the upper end of the tube to carry several channels formed longitudinally in the sleeve to the appropriate axial location where they penetrate the tube to allow pressure measurements and temperature measurements with thermocouples. The high pressure measurement is made using a channel penetrating the tube away from the body and the low pressure measurement is made at a location at the widest part of the body. An end plug seals the end of the device and holes at its upper end allow fluid to pass from the interior of the tube into a plenum. The channels are made by cutting grooves in the sleeve, the grooves widened at the surface of the sleeve and then a strip of sleeve material is welded to the grooves closing the channels. Preferably the sleeve is packed with powdered graphite before cutting the grooves and welding the strips. 7 figures.
McKay, Mark D.; Sweeney, Chad E.; Spangler, Jr., B. Samuel
1993-01-01
A flow meter and temperature measuring device comprising a tube with a body centered therein for restricting flow and a sleeve at the upper end of the tube to carry several channels formed longitudinally in the sleeve to the appropriate axial location where they penetrate the tube to allow pressure measurements and temperature measurements with thermocouples. The high pressure measurement is made using a channel penetrating the tube away from the body and the low pressure measurement is made at a location at the widest part of the body. An end plug seals the end of the device and holes at its upper end allow fluid to pass from the interior of the tube into a plenum. The channels are made by cutting grooves in the sleeve, the grooves widened at the surface of the sleeve and then a strip of sleeve material is welded to the grooves closing the channels. Preferably the sleeve is packed with powdered graphite before cutting the grooves and welding the strips.
Acoustic concentration of particles in fluid flow
Ward, Michael D.; Kaduchak, Gregory
2010-11-23
An apparatus for acoustic concentration of particles in a fluid flow includes a substantially acoustically transparent membrane and a vibration generator that define a fluid flow path therebetween. The fluid flow path is in fluid communication with a fluid source and a fluid outlet and the vibration generator is disposed adjacent the fluid flow path and is capable of producing an acoustic field in the fluid flow path. The acoustic field produces at least one pressure minima in the fluid flow path at a predetermined location within the fluid flow path and forces predetermined particles in the fluid flow path to the at least one pressure minima.
Foam flows in 2D porous media: intermittency and bubble fragmentation
NASA Astrophysics Data System (ADS)
Meheust, Y.; Géraud, B.; Jones, S. A.; Dollet, B.; Cantat, I.
2014-12-01
Flowing foams are used in many engineering and technical applications. A well-known application is oil recovery. Another one is the remediation of polluted soils: the foam is injected into the ground in order to mobilize chemical species present in the medium. Apart from potential interesting physico-chemical and biochemical properties, foams have peculiar flow properties that applications might benefit of. In particular, viscous dissipation arises mostly from the contact zones between the soap films and the walls, which results in peculiar friction laws allowing the foam to invade narrow pores more efficiently than Newtonian fluids would. In most experimental studies no local information of the foam structure is possible, and only global quantities such as the effective viscosity can be measured. Using two-dimensional transparent flow cells, we have previously shown that foam structural [1] and elastic [2] effects significantly impact the flow of foams in porous media. We now present an investigation of foam flow through a two-dimensional (2D) porous medium consisting of circular obstacles positioned randomly in a Hele-Shaw cell (see figure). The foam structure is recorded in time by a video camera and subsequently analyzed by image processing, which provides us with the velocity field and spatial distribution of bubble sizes. The flow exhibits a rich phenomenology, including flow irreversibility, preferential flow paths, local flow intermittency/non-stationarity (despite the imposed permanent global flow rate). Moreover, the medium impacts the nature of the flowing fluid by selecting the bubble size through bubble fragmentation. We investigate how preferential flow paths and intermittency depend on the imposed global flow rate and foam quality (the water content), and show that the spatial distribution of bubble sizes is to some extent correlated with the velocity field. We furthermore measure the evolution, along the flow direction, of the probability density
Integrated Coupling of Surface and Subsurface Flow with HYDRUS-2D
NASA Astrophysics Data System (ADS)
Hartmann, Anne; Šimůnek, Jirka; Wöhling, Thomas; Schütze, Niels
2016-04-01
of various model parameters on simulated overland flow (while considering or neglecting the effects of subsurface flow) was carried out to verify the applicability of the model to different problems. The model produced reasonable results in describing the diffusion wave approximation and its interactions with subsurface flow processes. The model could handle coupled surface-subsurface processes for conditions involving runoff generated by infiltration excess, saturation excess, or run-on, as well as a combination of these runoff generating processes. Several standard features of the HYDRUS 2D model, such as root water uptake and evaporation from the soil surface, as well as evaporation from runoff layer, can still be considered by the new model. The code required relatively small time steps when overland flow was active, resulting in long simulation times, and sometimes produced poor mass balance. The model nevertheless showed potential to be a useful tool for addressing various issues related to irrigation research and to natural generation of overland flow at the hillslope scale. Maxwell, R., Putti, M., Meyerhoff, S., Delf, J., Ferguson, I., Ivanov, V., Kim, J., Kolditz, O., Kollet, S., Kumar, M., Lopez, S., Niu, J., Paniconi, C., Park, Y.-J., Phanikumar, M., Shen, C., Sudicky, E., and Sulis, M. (2014). Surface-subsurface model intercomparison: A first set of benchmark results to diagnose integrated hydrology and feedbacks. Water Resourc. Res., 50:1531-1549. Šimůnek, J., van Genuchten, M. T., and Šejna, M. (2011). The HYDRUS Software Package for Simulating Two- and Three-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media. Technical Manual, Version 2.0, PC Progress, Prague, Czech Republic. Takizawa, K., Bazilevs Y., Tezduyar, T. E., Long, C.C., Marsden, A. L. and Schjodt.K., Patient-Specific Cardiovascular Fluid Mechanics Analysis with the ST and ALE-VMS Method in Idelsohn, S. R. (2014). Numerical Simulations of Coupled Problems
Geophysical Fluid Flow Cell Simulation
NASA Technical Reports Server (NTRS)
1998-01-01
Computer simulation of atmospheric flow corresponds well to imges taken during the second Geophysical Fluid Flow Cell (BFFC) mission. The top shows a view from the pole, while the bottom shows a view from the equator. Red corresponds to hot fluid rising while blue shows cold fluid falling. This simulation was developed by Anil Deane of the University of Maryland, College Park and Paul Fischer of Argorne National Laboratory. Credit: NASA/Goddard Space Flight Center
Template Matching Using a Fluid Flow Model
NASA Astrophysics Data System (ADS)
Newman, William Curtis
Template matching is successfully used in machine recognition of isolated spoken words. In these systems a word is broken into frames (20 millisecond time slices) and the spectral characteristics of each frame are found. Thus, each word is represented as a 2-dimensional (2-D) function of spectral characteristic and frame number. An unknown word is recognized by matching its 2-D representation to previously stored example words, or templates, also in this 2-D form. A new model for this matching step will be introduced. The 2-D representations of the template and unknown are used to determine the shape of a volume of viscous fluid. This volume is broken up into many small elements. The unknown is changed into the template by allowing flows between the element boundaries. Finally the match between the template and unknown is determined by calculating a weighted squared sum of the flow values. The model also allows the relative flow resistance between the element boundaries to be changed. This is useful for characterizing the important features of a given template. The flow resistances are changed according to the gradient of a simple performance function. This performance function is evaluated using a set of training samples provided by the user. The model is applied to isolated word and single character recognition tasks. Results indicate the applications where this model works best.
Simulation of abrasive flow machining process for 2D and 3D mixture models
NASA Astrophysics Data System (ADS)
Dash, Rupalika; Maity, Kalipada
2015-12-01
Improvement of surface finish and material removal has been quite a challenge in a finishing operation such as abrasive flow machining (AFM). Factors that affect the surface finish and material removal are media viscosity, extrusion pressure, piston velocity, and particle size in abrasive flow machining process. Performing experiments for all the parameters and accurately obtaining an optimized parameter in a short time are difficult to accomplish because the operation requires a precise finish. Computational fluid dynamics (CFD) simulation was employed to accurately determine optimum parameters. In the current work, a 2D model was designed, and the flow analysis, force calculation, and material removal prediction were performed and compared with the available experimental data. Another 3D model for a swaging die finishing using AFM was simulated at different viscosities of the media to study the effects on the controlling parameters. A CFD simulation was performed by using commercially available ANSYS FLUENT. Two phases were considered for the flow analysis, and multiphase mixture model was taken into account. The fluid was considered to be a
CAS2D: FORTRAN program for nonrotating blade-to-blade, steady, potential transonic cascade flows
NASA Technical Reports Server (NTRS)
Dulikravich, D. S.
1980-01-01
An exact, full-potential-equation (FPE) model for the steady, irrotational, homentropic and homoenergetic flow of a compressible, homocompositional, inviscid fluid through two dimensional planar cascades of airfoils was derived, together with its appropriate boundary conditions. A computer program, CAS2D, was developed that numerically solves an artificially time-dependent form of the actual FPE. The governing equation was discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field was discretized by providing a boundary-fitted, nonuniform computational mesh. The mesh was generated by using a sequence of conforming mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the FPE was solved iteratively by using successive line overrelaxation. The possible isentropic shocks were correctly captured by adding explicitly an artificial viscosity in a conservative form. In addition, a three-level consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two dimensional cascade flows.
Programming fluid flow with microstructures
NASA Astrophysics Data System (ADS)
Amini, Hamed; Masaeli, Mahdokht; di Carlo, Dino
2011-11-01
Flow control and fluid interface manipulation in microfluidic platforms are of great importance in a variety of applications. Current approaches to manipulate fluids generally rely on complex designs, difficult-to-fabricate 3D platforms or use of active methods. Here we show that in the presence of simple cylindrical obstacles (i.e. pillars) in a microchannel, at moderate to high flow rates, streamlines tend to turn and stretch in a manner that, unlike intuition for Stokes flow, does not precisely reverse after passing the pillar. The asymmetric flow behavior up- and down-stream of the pillar due to fluid inertia manifests itself as a total deformation of the topology of streamlines that effectively creates a net secondary flow which resembles the recirculating Dean flow in curving channels. Confocal images were taken to investigate the secondary flow for a variety of microstructure settings. We also developed a numerical technique to map the fluid motion in the channel which is utilized to characterize the secondary flow as well as to engineer the fluid patterns within the channel. This passive method creates the possibility of exceptional control of the 3D structure of the fluid within a microfluidic platform which can significantly advance applications requiring fluid interface control (e.g. optofluidics), ultrafast mixing and solution control around cells.
Relaminarization of fluid flows
NASA Technical Reports Server (NTRS)
Narasimha, R.; Sreenivasan, K. R.
1979-01-01
The mechanisms of the relaminarization of turbulent flows are investigated with a view to establishing any general principles that might govern them. Three basic archetypes of reverting flows are considered: the dissipative type, the absorptive type, and the Richardson type exemplified by a turbulent boundary layer subjected to severe acceleration. A number of other different reverting flows are then considered in the light of the analysis of these archetypes, including radial Poiseuille flow, convex boundary layers, flows reverting by rotation, injection, and suction, as well as heated horizontal and vertical gas flows. Magnetohydrodynamic duct flows are also examined. Applications of flow reversion for turbulence control are discussed.
Fluid Flow Phenomena during Welding
Zhang, Wei
2011-01-01
MOLTEN WELD POOLS are dynamic. Liquid in the weld pool in acted on by several strong forces, which can result in high-velocity fluid motion. Fluid flow velocities exceeding 1 m/s (3.3 ft/s) have been observed in gas tungsten arc (GTA) welds under ordinary welding conditions, and higher velocities have been measured in submerged arc welds. Fluid flow is important because it affects weld shape and is related to the formation of a variety of weld defects. Moving liquid transports heat and often dominates heat transport in the weld pool. Because heat transport by mass flow depends on the direction and speed of fluid motion, weld pool shape can differ dramatically from that predicted by conductive heat flow. Temperature gradients are also altered by fluid flow, which can affect weld microstructure. A number of defects in GTA welds have been attributed to fluid flow or changes in fluid flow, including lack of penetration, top bead roughness, humped beads, finger penetration, and undercutting. Instabilities in the liquid film around the keyhole in electron beam and laser welds are responsible for the uneven penetration (spiking) characteristic of these types of welds.
A 2D pendulum submitted to an incoming flow: drag acting like gravity and new instabilities
NASA Astrophysics Data System (ADS)
Fani, Andrea; Gallaire, Francois
2013-11-01
Flow induced oscillations of slender bodies facing an incoming flow are relevant in a large number of engineering applications, such as the design of tubular structures of offshore platforms, heat exchangers and energy harvesting. Numerical simulations and experiments available in literature often consider a circular cylinder in an uniform flow which can move only transversally with respect to the flow direction. In a recent work Semin et al. (JFM, 2011) studied a tethered 2D cylinder strongly confined between two parallel plane walls. It is shown that confinement alters significantly the flow dynamics, with a new instability, denoted confinement induced vibration (CIV), which occur at a Reynolds number much lower than the vortex induced vibration (VIV) critical one. In the present work we characterize the instability scenario of a confined tethered cylinder by means of a global stability analysis of the fluid-structure problem. In strongly confined channels, a periodic unstable mode, related to CIV vibrations, is observed, while for moderated confinement a new steady diverging instability is founded.
Insertable fluid flow passage bridgepiece and method
Jones, Daniel O.
2000-01-01
A fluid flow passage bridgepiece for insertion into an open-face fluid flow channel of a fluid flow plate is provided. The bridgepiece provides a sealed passage from a columnar fluid flow manifold to the flow channel, thereby preventing undesirable leakage into and out of the columnar fluid flow manifold. When deployed in the various fluid flow plates that are used in a Proton Exchange Membrane (PEM) fuel cell, bridgepieces of this invention prevent mixing of reactant gases, leakage of coolant or humidification water, and occlusion of the fluid flow channel by gasket material. The invention also provides a fluid flow plate assembly including an insertable bridgepiece, a fluid flow plate adapted for use with an insertable bridgepiece, and a method of manufacturing a fluid flow plate with an insertable fluid flow passage bridgepiece.
COYOTE: A computer program for 2-D reactive flow simulations
Cloutman, L.D.
1990-04-01
We describe the numerical algorithm used in the COYOTE two- dimensional, transient, Eulerian hydrodynamics program for reactive flows. The program has a variety of options that provide capabilities for a wide range of applications, and it is designed to be robust and relatively easy to use while maintaining adequate accuracy and efficiency to solve realistic problems. It is based on the ICE method, and it includes a general species and chemical reaction network for simulating reactive flows. It also includes swirl, turbulence transport models, and a nonuniform mesh capability. We describe several applications of the program. 33 refs., 4 figs.
2D VARIABLY SATURATED FLOWS: PHYSICAL SCALING AND BAYESIAN ESTIMATION
A novel dimensionless formulation for water flow in two-dimensional variably saturated media is presented. It shows that scaling physical systems requires conservation of the ratio between capillary forces and gravity forces. A direct result of this finding is that for two phys...
Chang, F.H.; Santee, G.E. Jr.; Mortensen, G.A.; Brockett, G.F.; Gross, M.B.; Silling, S.A.; Belytschko, T.
1981-03-01
This report, the second in a series of reports for RP-1065, describes the second step in the stepwise approach for developing the three-dimensional, nonlinear, fluid/structure interaction methodology to assess the hydroloads on a large PWR during the subcooled portions of a hypothetical LOCA. The second step in the methodology considers enhancements and special modifications to the 2D STEALTH-HYDRO computer program and the 2D WHAMSE computer program. The 2D STEALTH-HYDRO enhancements consist of a fluid-fluid coupling control-volume model and an orifice control-volume model. The enhancements to 2D WHAMSE include elimination of the implicit integration routines, material models, and structural elements not required for the hydroloads application. In addition the logic for coupling the 2D STEALTH-HYDRO computer program to the 2D WHAMSE computer program is discussed.
Fluid Physics of Foam Evolution and Flow
NASA Technical Reports Server (NTRS)
Aref, H.; Thoroddsen, S. T.; Sullivan, J. M.
2003-01-01
The grant supported theoretical, numerical and experimental work focused on the elucidation of the fluid physics of foam structure, evolution and flow. The experimental work concentrated on these subject areas: (a) Measurements of the speed of reconnections within a foam; (b) statistics of bubble rearrangements; and (c) three-dimensional reconstruction of the foam structure. On the numerical simulation and theory side our efforts concentrated on the subjects: (a) simulation techniques for 2D and 3D foams; (b) phase transition in a compressible foam; and (c) TCP structures.
Roton Excitations and the Fluid-Solid Phase Transition in Superfluid 2D Yukawa Bosons
NASA Astrophysics Data System (ADS)
Molinelli, S.; Galli, D. E.; Reatto, L.; Motta, M.
2016-05-01
We compute several ground-state properties and the dynamical structure factor of a zero-temperature system of Bosons interacting with the 2D screened Coulomb (2D-SC) potential. We resort to the exact shadow path integral ground state (SPIGS) quantum Monte Carlo method to compute the imaginary-time correlation function of the model, and to the genetic algorithm via falsification of theories (GIFT) to retrieve the dynamical structure factor. We provide a detailed comparison of ground-state properties and collective excitations of 2D-SC and ^4 He atoms. The roton energy of the 2D-SC system is an increasing function of density, and not a decreasing one as in ^4 He. This result is in contrast with the view that the roton is the soft mode of the fluid-solid transition. We uncover a remarkable quasi-universality of backflow and of other properties when expressed in terms of the amount of short-range order as quantified by the height of the first peak of the static structure factor.
Rotational fluid flow experiment
NASA Technical Reports Server (NTRS)
1991-01-01
This project which began in 1986 as part of the Worcester Polytechnic Institute (WPI) Advanced Space Design Program focuses on the design and implementation of an electromechanical system for studying vortex behavior in a microgravity environment. Most of the existing equipment was revised and redesigned by this project team, as necessary. Emphasis was placed on documentation and integration of the electrical and mechanical subsystems. Project results include reconfiguration and thorough testing of all hardware subsystems, implementation of an infrared gas entrainment detector, new signal processing circuitry for the ultrasonic fluid circulation device, improved prototype interface circuits, and software for overall control of experiment operation.
NASA Technical Reports Server (NTRS)
Kapoor, Kamlesh; Anderson, Bernhard H.; Shaw, Robert J.
1994-01-01
A two-dimensional computational code, PRLUS2D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for two-dimensional shock-wave/turbulent-boundary-layer interactions. The problem of compression corners at supersonic speeds was solved using the RPLUS2D code. To validate the RPLUS2D code for hypersonic speeds, it was applied to a realistic hypersonic inlet geometry. Both the Baldwin-Lomax and the Chien two-equation turbulence models were used. Computational results showed that the RPLUS2D code compared very well with experimentally obtained data for supersonic compression corner flows, except in the case of large separated flows resulting from the interactions between the shock wave and turbulent boundary layer. The computational results compared well with the experiment results in a hypersonic NASA P8 inlet case, with the Chien two-equation turbulence model performing better than the Baldwin-Lomax model.
NASA Astrophysics Data System (ADS)
Stone, James M.; Norman, Michael L.
1992-06-01
A detailed description of ZEUS-2D, a numerical code for the simulation of fluid dynamical flows including a self-consistent treatment of the effects of magnetic fields and radiation transfer is presented. Attention is given to the hydrodynamic (HD) algorithms which form the foundation for the more complex MHD and radiation HD algorithms. The effect of self-gravity on the flow dynamics is accounted for by an iterative solution of the sparse-banded matrix resulting from discretizing the Poisson equation in multidimensions. The results of an extensive series of HD test problems are presented. A detailed description of the MHD algorithms in ZEUS-2D is presented. A new method of computing the electromotive force is developed using the method of characteristics (MOC). It is demonstrated through the results of an extensive series of MHD test problems that the resulting hybrid MOC-constrained transport method provides for the accurate evolution of all modes of MHD wave families.
2D numerical modelling of fluid percolation in the subduction zone
NASA Astrophysics Data System (ADS)
Dymkova, D.; Gerya, T.; Podladchikov, Y.
2012-04-01
Subducting slab dehydration and resulting aqueous fluid percolation triggers partial melting in the mantle wedge and is accompanied with the further melt percolation through the porous space to the region above the slab. This problem is a complex coupled chemical, thermal and mechanical process responsible for the magmatic arcs formation and change of the mantle wedge properties. We have created a two-dimensional model of a two-phase flow in a porous media solving a coupled Darcy-Stokes system of equations for two incompressible media for the case of nonlinear visco-plastic rheology of solid matrix. Our system of equation is expanded for the high-porosity limits and stabilized for the case of high porosity contrasts. We use a finite-difference method with fully staggered grid in a combination with marker-in-cell technique for advection of fluid and solid phase. We performed a comparison with a benchmark of a thermal convection in a porous media in a bottom-heated box to verify the interdependency of Rayleigh and Nusselt numbers with earlier obtained ones (Cherkaoui & Wilcock, 1999). We have demonstrated the stability and robustness of the algorithm in case of strongly non-linear visco-plastic rheology of solid including cases with localization of both deformation and porous flow along spontaneously forming shear bands. We have checked our model for the forming of localized porous channels under a simple shear stress (Katz et al, 2006). We have developed a setup of a self-initiating due to gravitational instability subduction. With our coupled fluid-solid flow we have achieved a self-consistent water downward suction by a slab bending predicted by the other models with a simplified fluid kinematical motion implementation (Faccenda et al, 2009). With this setup we have obtained a self-consistent upper crust weakening by a porous fluid pressure which was theoretically assumed in the previously existing subduction models (Gerya & Meilick, 2011; Faccenda et al, 2009
A multispeed Discrete Boltzmann Model for transcritical 2D shallow water flows
NASA Astrophysics Data System (ADS)
La Rocca, Michele; Montessori, Andrea; Prestininzi, Pietro; Succi, Sauro
2015-03-01
In this work a Discrete Boltzmann Model for the solution of transcritical 2D shallow water flows is presented and validated. In order to provide the model with transcritical capabilities, a particular multispeed velocity set has been employed for the discretization of the Boltzmann equation. It is shown that this particular set naturally yields a simple and closed procedure to determine higher order equilibrium distribution functions needed to simulate transcritical flow. The model is validated through several classical benchmarks and is proven to correctly and accurately simulate both 1D and 2D transitions between the two flow regimes.
Lee, Shern-Long; Yuan, Zhongyi; Chen, Long; Mali, Kunal S; Müllen, Klaus; De Feyter, Steven
2014-05-28
Controlling crystal polymorphism constitutes a formidable challenge in contemporary chemistry. Two-dimensional (2D) crystals often provide model systems to decipher the complications in 3D crystals. In this contribution, we explore a unique way of governing 2D polymorphism at the organic liquid-solid interface. We demonstrate that a directional solvent flow could be used to stabilize crystalline monolayers of a metastable polymorph. Furthermore, flow fields active within the applied flow generate millimeter-sized domains of either polymorph in a controlled and reproducible fashion. PMID:24867142
2D-simulation of wet steam flow in a steam turbine with spontaneous condensation
NASA Astrophysics Data System (ADS)
Sun, Lan-Xin; Zheng, Qun; Liu, Shun-Long
2007-06-01
Removal of condensates from wet steam flow in the last stages of steam turbines significantly promotes stage efficiency and prevents erosion of rotors. In this paper, homogeneous spontaneous condensation in transonic steam flow in the 2-D rotor-tip section of a stage turbine is investigated. Calculated results agree with experimental data reasonably well. On the basis of the above work, a 2-D numerical simulation of wet steam flow in adjacent root sections of a complex steam turbine stage was carried out. Computational results were analyzed and provide insights into effective removal of humidity.
Fluctuating Pressure Data from 2-D Nozzle Cold Flow Tests (Dual Bell)
NASA Technical Reports Server (NTRS)
Nesman, Tomas E.
2001-01-01
Rocket engines nozzle performance changes as a vehicle climbs through the atmosphere. An altitude compensating nozzle, ACN, is intended to improve on a fixed geometry bell nozzle that performs at optimum at only one trajectory point. In addition to nozzle performance, nozzle transient loads are an important consideration. Any nozzle experiences large transient toads when shocks pass through the nozzle at start and shutdown. Additional transient toads will occur at transitional flow conditions. The objectives of cold flow nozzle testing at MSFC are CFD benchmark / calibration and Unsteady flow / sideloads. Initial testing performed with 2-D inserts to 14" transonic wind tunnel. Recent review of 2-D data in preparation for nozzle test facility 3-D testing. This presentation shows fluctuating pressure data and some observations from 2-D dual-bell nozzle cold flow tests.
Turbulence modeling for subsonic separated flows over 2-D airfoils and 3-D wings
NASA Astrophysics Data System (ADS)
Rosen, Aaron M.
Accurate predictions of turbulent boundary layers and flow separation through computational fluid dynamics (CFD) are becoming more and more essential for the prediction of loads in the design of aerodynamic flight components. Standard eddy viscosity models used in many commercial codes today do not capture the nonequilibrium effects seen in a separated flow and thus do not generally make accurate separation predictions. Part of the reason for this is that under nonequilibrium conditions such as a strong adverse pressure gradient, the history effects of the flow play an important role in the growth and decay of turbulence. More recent turbulence models such as Olsen and Coakley's Lag model and Lillard's lagRST model seek to simulate these effects by lagging the turbulent variables when nonequilibrium effects become important. The purpose of the current research is to assess how these nonequilibrium turbulence models capture the separated regions on various 2-D airfoils and 3-D wings. Nonequilibrium models including the Lag model and the lagRST model are evaluated in comparison with three baseline models (Spalart-Allmaras, Wilcox's k-omega, and Menter's SST) using a modified version of the OVERFLOW code. Tuning the model coefficients of the Lag and lagRST models is also explored. Results show that the various lagRST formulations display an improvement in velocity profile predictions over the standard RANS models, but have trouble capturing the edge of the boundary layer. Experimental separation location measurements were not available, but several trends are noted which may be useful to tuning the model coefficients in the future.
NASA Astrophysics Data System (ADS)
Zou, B.; Li, D. F.; Hu, H. J.; Zhang, H. W.; Lou, L. H.; Chen, M.; Lv, Z. Y.
Based on the verified two dimensional(2D) finite element model for river flow simulation, the effect of estuary training levees on the water flow and sediment movement in the Yellow River estuary is analyzed. For disclosing the effect of setting the two training levees on the flow and sediment motion, the calculation and analysis for the two projects, (one is no levees, the other is setting up two no levees) are given. The results show that when setting up two training levees, water flow is bound by levees and the water flows become more concentrated. As a result, velocity increases in the main channel, sediment carrying capacity of water flow increases correspondingly.
Fluid flow electrophoresis in space
NASA Technical Reports Server (NTRS)
Griffin, R. N.
1975-01-01
Four areas relating to free-flow electrophoresis in space were investigated. The first was the degree of improvement over earthbound operations that might be expected. The second area of investigation covered the problems in developing a flowing buffer electrophoresis apparatus. The third area of investigation was the problem of testing on the ground equipment designed for use in space. The fourth area of investigation was the improvement to be expected in space for purification of biologicals. The results of some ground-based experiments are described. Other studies included cooling requirements in space, fluid sealing techniques, and measurement of voltage drop across membranes.
Ferroelectric Fluid Flow Control Valve
NASA Technical Reports Server (NTRS)
Jalink, Antony, Jr. (Inventor); Hellbaum, Richard F. (Inventor); Rohrbach, Wayne W. (Inventor)
1999-01-01
An active valve is controlled and driven by external electrical actuation of a ferroelectric actuator to provide for improved passage of the fluid during certain time periods and to provide positive closure of the valve during other time periods. The valve provides improved passage in the direction of flow and positive closure in the direction against the flow. The actuator is a dome shaped internally prestressed ferroelectric actuator having a curvature, said dome shaped actuator having a rim and an apex. and a dome height measured from a plane through said rim said apex that varies with an electric voltage applied between an inside and an outside surface of said dome shaped actuator.
Cerebrospinal fluid flow in adults.
Bradley, William G; Haughton, Victor; Mardal, Kent-Andre
2016-01-01
This chapter uses magnetic resonance imaging phase-contrast cerebrospinal fluid (CSF) flow measurements to predict which clinical normal-pressure hydrocephalus (NPH) patients will respond to shunting as well as which patients with Chiari I are likely to develop symptoms of syringomyelia. Symptomatic NPH patients with CSF flow (measured as the aqueductal CSF stroke volume) which is shown to be hyperdynamic (defined as twice normal) are quite likely to respond to ventriculoperitoneal shunting. The hyperdynamic CSF flow results from normal systolic brain expansion compressing the enlarged ventricles. When atrophy occurs, there is less brain expansion, decreased aqueductal CSF flow, and less likelihood of responding to shunting. It appears that NPH is a "two-hit" disease, starting as benign external hydrocephalus in infancy, followed by deep white-matter ischemia in late adulthood, which causes increased resistance to CSF outflow through the extracellular space of the brain. Using computational flow dynamics (CFD), CSF flow can be modeled at the foramen magnum and in the upper cervical spine. As in the case of NPH, hyperdynamic CSF flow appears to cause the signs and symptoms in Chiari I and can provide an additional indication for surgical decompression. CFD can also predict CSF pressures over the cardiac cycle. It has been hypothesized that elevated pressure pulses may be a significant etiologic factor in some cases of syringomyelia. PMID:27432684
An F2D analysis of the Flow Instability Test (FIT) experiment
Suo-Anttila, A.
1993-10-01
The F2D code was used to analyze the Flow-Instability-Test (FIT) experiments at Brookhaven National Laboratories. A one-dimensional analysis of the experiment indicated that at the higher temperature levels the element should be unstable. The experimental data corroborated this theory. The two-dimensional simulation behaved in a manner that was very similar to the experimentally measured behavior. In conclusion, the FIT experimental analysis yields partial code validation of F2D, and it also validates the methodology that is used in analyzing thermal flow stability.
Extension and application of the Preissmann slot model to 2D transient mixed flows
NASA Astrophysics Data System (ADS)
Maranzoni, Andrea; Dazzi, Susanna; Aureli, Francesca; Mignosa, Paolo
2015-08-01
This paper presents an extension of the Preissmann slot concept for the modeling of highly transient two-dimensional (2D) mixed flows. The classic conservative formulation of the 2D shallow water equations for free surface flows is adapted by assuming that two fictitious vertical slots, aligned along the two Cartesian plane directions and normally intersecting, are added on the ceiling of each integration element. Accordingly, transitions between free surface and pressurized flow can be handled in a natural and straightforward way by using the same set of governing equations. The opportunity of coupling free surface and pressurized flows is actually useful not only in one-dimensional (1D) problems concerning sewer systems but also for modeling 2D flooding phenomena in which the pressurization of bridges, culverts, or other crossing hydraulic structures can be expected. Numerical simulations are performed by using a shock-capturing MUSCL-Hancock finite volume scheme combined with the FORCE (First-Order Centred) solver for the evaluation of the numerical fluxes. The validation of the mathematical model is accomplished on the basis of both exact solutions of 1D discontinuous initial value problems and reference radial solutions of idealized test cases with cylindrical symmetry. Furthermore, the capability of the model to deal with practical field-scale applications is assessed by simulating the transit of a bore under an arch bridge. Numerical results show that the proposed model is suitable for the prediction of highly transient 2D mixed flows.
Combination 3D TOP with 2D PC MRA Techique for cerebral blood flow volume measurement.
Guo, G; Wu, R H; Zhang, Y P; Guan, J T; Guo, Y L; Cheng, Y; terBrugge, K; Mikulis, D J
2006-01-01
To demonstrate the discrepancy of cerebral blood flow volume (BFV) estimation with 2D phase-contrast (2D PC) MRA guided with 3D time-of-flight (3D TOF) MR localization by using an "internal" standard. 20 groups of the common (CCA), internal (ICA), and external (ECA) carotid arteries in 10 healthy subjects were examined with 2D PC MRA guided by 3D TOF MR angiograms. The sum BFV of the internal and external carotid arteries was then compared with the ipsilateral common carotid artery flow. An accurate technique would demonstrate no difference. The difference was therefore a measure of accuracy of the method. 3D TOF MRA localization is presented to allow the determination of a slice orientation to improve the accuracy of 2D PC MRA in estimate the BFV. By using the combined protocols, there was better correlation in BFV estimate between the sum of ICA+ECA with the ipsilateral CCA (R2=0.729, P=0.000). The inconsistency (mean +/- SD) was found to be 6.95 +/- 5.95% for estimate the BFV in ICA+ECA and ipsilateral CCA. The main inconsistency was contributed to the ECA and its branches. Guided with 3D TOF MRA localization, 2D PC MRA is more accurate in the determination of blood flow volume in the carotid arteries. PMID:17946401
Transient Wellbore Fluid Flow Model
Energy Science and Technology Software Center (ESTSC)
1982-04-06
WELBORE is a code to solve transient, one-dimensional two-phase or single-phase non-isothermal fluid flow in a wellbore. The primary thermodynamic variables used in solving the equations are the pressure and specific energy. An equation of state subroutine provides the density, quality, and temperature. The heat loss out of the wellbore is calculated by solving a radial diffusion equation for the temperature changes outside the bore. The calculation is done at each node point in themore » wellbore.« less
Use of finite volume radiation for predicting the Knudsen minimum in 2D channel flow
Malhotra, Chetan P.; Mahajan, Roop L.
2014-12-09
In an earlier paper we employed an analogy between surface-to-surface radiation and free-molecular flow to model Knudsen flow through tubes and onto planes. In the current paper we extend the analogy between thermal radiation and molecular flow to model the flow of a gas in a 2D channel across all regimes of rarefaction. To accomplish this, we break down the problem of gaseous flow into three sub-problems (self-diffusion, mass-motion and generation of pressure gradient) and use the finite volume method for modeling radiation through participating media to model the transport in each sub-problem as a radiation problem. We first model molecular self-diffusion in the stationary gas by modeling the transport of the molecular number density through the gas starting from the analytical asymptote for free-molecular flow to the kinetic theory limit of gaseous self-diffusion. We then model the transport of momentum through the gas at unit pressure gradient to predict Poiseuille flow and slip flow in the 2D gas. Lastly, we predict the generation of pressure gradient within the gas due to molecular collisions by modeling the transport of the forces generated due to collisions per unit volume of gas. We then proceed to combine the three radiation problems to predict flow of the gas over the entire Knudsen number regime from free-molecular to transition to continuum flow and successfully capture the Knudsen minimum at Kn ∼ 1.
Evaluation of 2D shallow-water model for spillway flow with a complex geometry
Technology Transfer Automated Retrieval System (TEKTRAN)
Although the two-dimensional (2D) shallow water model is formulated based on several assumptions such as hydrostatic pressure distribution and vertical velocity is negligible, as a simple alternative to the complex 3D model, it has been used to compute water flows in which these assumptions may be ...
Parallelized CCHE2D flow model with CUDA Fortran on Graphics Process Units
Technology Transfer Automated Retrieval System (TEKTRAN)
This paper presents the CCHE2D implicit flow model parallelized using CUDA Fortran programming technique on Graphics Processing Units (GPUs). A parallelized implicit Alternating Direction Implicit (ADI) solver using Parallel Cyclic Reduction (PCR) algorithm on GPU is developed and tested. This solve...
Application of a generalized minimal residual method to 2D unsteady flows
NASA Technical Reports Server (NTRS)
Hixon, Ray; Sankar, L. N.
1992-01-01
A generalized minimum residual scheme (GMRES), previously developed for solving nonlinear and linear systems of equations, has been applied to the numerical solution of 2D unsteady compressible flows. It is found that the use of GMRES significantly increases the time step that may be used, compared to noniterative implicit schemes. The feasibility of reducing the memory requirements of the GMRES scheme using a multigrid strategy has also been explored. Several sample steady and unsteady viscous flow applications are presented.
2D fluid simulations of acoustic waves in pulsed ICP discharges: Comparison with experiments
NASA Astrophysics Data System (ADS)
Despiau-Pujo, Emilie; Cunge, Gilles; Sadeghi, Nader; Braithwaite, N. St. J.
2012-10-01
Neutral depletion, which is mostly caused by gas heating under typical material processing conditions, is an important phenomenon in high-density plasmas. In low pressure pulsed discharges, experiments show that additional depletion due to electron pressure (Pe) may have a non-negligible influence on radical transport [1]. To evaluate this effect, comparisons between 2D fluid simulations and measurements of gas convection in Ar/Cl2 pulsed ICP plasmas are reported. In the afterglow, Pe drops rapidly by electron cooling which generates a neutral pressure gradient between the plasma bulk and the reactor walls. This in turn forces the cold surrounding gas to move rapidly towards the center, thus launching an acoustic wave in the reactor. Time-resolved measurements of atoms drift velocity and gas temperature by LIF and LAS in the early afterglow are consistent with gas drifting at acoustic wave velocity followed by rapid gas cooling. Similar results are predicted by the model. The ion flux at the reactor walls is also shown to oscillate in phase with the acoustic wave due to ion-neutral friction forces. Finally, during plasma ignition, experiments show opposite phenomena when Pe rises.[4pt] [1] Cunge et al, APL 96, 131501 (2010)
Fluid flows around nanoelectromechanical resonators
NASA Astrophysics Data System (ADS)
Svitelskiy, O.; Sauer, V.; Liu, N.; Vick, D.; Cheng, K. M.; Freeman, M. R.; Hiebert, W. K.
2012-02-01
To explore properties of fluids on a nanosize scale, we fabricated by a standard top down technique a series of nanoelectromechanical resonators (cantilevers and bridges) with widths w and thicknesses t from 100 to 500 nm; lengths l from 0.5 to 12 micron; and resonant frequencies f from 10 to 400 MHz. For the sake of purity of the experiment, the undercut in the widest (w=500 nm) devices was eliminated using the focused ion beam. To model the fluidic environment the devices were placed in the atmosphere of compressed gases (He, N2, CO2, Ar, H2) at pressures from vacuum up to 20 MPa, and in liquid CO2; their properties were studied by the real time stroboscopic optical interferometry. Thus, we fully explored the Newtonian and non-Newtonian flow damping models. Observing free molecular flow extending above atmospheric pressure, we find the fluid relaxation time model to be the best approximation throughout, but not beyond, the non-Newtonian regime, and both, vibrating spheres model and the model based on Knudsen number, to be valid in the viscous limit.
Simulation of the flow and mass transfer for KDP crystals undergoing 2D translation during growth
NASA Astrophysics Data System (ADS)
Zhou, Chuan; Li, Mingwei; Hu, Zhitao; Yin, Huawei; Wang, Bangguo; Cui, Qidong
2016-09-01
In this study, a novel motion mode for crystals during growth, i.e., 2D translation, is proposed. Numerical simulations of flow and mass transfer are conducted for the growth of large-scale potassium dihydrogen phosphate (KDP) crystals subjected to the new motion mode. Surface supersaturation and shear stress are obtained as functions of the translational velocity, distance, size, orientation of crystals. The dependence of these two parameters on the flow fields around the crystals is also discussed. The thicknesses of the solute boundary layer varied with translational velocity are described. The characteristics of solution flow and surface supersaturation distribution are summarized, where it suggests that the morphological stability of a crystal surface can be enhanced if the proposed 2D translation is applied to crystal growth.
Approximate analytic solutions to 3D unconfined groundwater flow within regional 2D models
NASA Astrophysics Data System (ADS)
Luther, K.; Haitjema, H. M.
2000-04-01
We present methods for finding approximate analytic solutions to three-dimensional (3D) unconfined steady state groundwater flow near partially penetrating and horizontal wells, and for combining those solutions with regional two-dimensional (2D) models. The 3D solutions use distributed singularities (analytic elements) to enforce boundary conditions on the phreatic surface and seepage faces at vertical wells, and to maintain fixed-head boundary conditions, obtained from the 2D model, at the perimeter of the 3D model. The approximate 3D solutions are analytic (continuous and differentiable) everywhere, including on the phreatic surface itself. While continuity of flow is satisfied exactly in the infinite 3D flow domain, water balance errors can occur across the phreatic surface.
Investigation of flow rate in a quasi-2D hopper with two symmetric outlets
NASA Astrophysics Data System (ADS)
Zhang, Xuezhi; Zhang, Sheng; Yang, Guanghui; Lin, Ping; Tian, Yuan; Wan, Jiang-Feng; Yang, Lei
2016-03-01
In this study the flow rate in a quasi-2D hopper with two symmetric outlets was investigated numerically. When two outlets merge into one, the flow rate is obviously bigger than that with a large enough distance between them. To explain the influence of the interaction between these two outlets, it was found that the velocity field above one outlet will be affected by the neighboring outlet by reducing the stagnant zone. With increasing distance, this interaction becomes weaker. Finally the formula of flow rate considering various outlet sizes and distances is fitted to the simulation data.
Fluid Flow Experiment for Undergraduate Laboratory.
ERIC Educational Resources Information Center
Vilimpochapornkul, Viroj; Obot, Nsima T.
1986-01-01
The undergraduate fluid mechanics laboratory at Clarkson University consists of three experiments: mixing; drag measurements; and fluid flow and pressure drop measurements. The latter experiment is described, considering equipment needed, procedures used, and typical results obtained. (JN)
Fluid flow, mineral reactions, and metasomatism
Ferry, J.M.; Dipple, G.M. )
1991-03-01
A general model that relates fluid flow along a temnperature gradient to chemical reaction in rocks can be used to quantitatively interpret petrologic and geochemical data on metasomatism from ancient flow systems in terms of flow direction and time-integrated fluid flux. The model is applied to regional metamorphism, quartz veins, and a metasomatized ductile fault zone.
Anomalous diffusion of an ellipsoid in quasi-2D active fluids
NASA Astrophysics Data System (ADS)
Peng, Yi; Yang, Ou; Tang, Chao; Cheng, Xiang
Enhanced diffusion of a tracer particle is a unique feature in active fluids. Here, we studied the diffusion of an ellipsoid in a free-standing film of E. coli. Particle diffusion is linearly enhanced at low bacterial concentrations, whereas a non-linear enhancement is observed at high bacterial concentrations due to the giant fluctuation. More importantly, we uncover an anomalous coupling between the translational and rotational degrees of freedom that is strictly prohibited in the classical Brownian diffusion. Combining experiments with theoretical modeling, we show that such an anomaly arises from the stretching flow induced by the force dipole of swimming bacteria. Our work illustrates a novel universal feature of active matter and transforms the understanding of fundamental transport processes in microbiological systems. ACS Petroleum Research Fund #54168-DNI9, NSF Faculty Early Career Development Program, DMR-1452180.
NASA Astrophysics Data System (ADS)
Kolb, Evelyne; Algarra, Nicolas; Vandembroucq, Damien; Lazarus, Arnaud
2015-11-01
We propose a new fluid/structure interaction in the unusual case of a dense granular medium flowing against an elastic fibre acting as a flexible intruder. We experimentally studied the deflection of a mylar flexible beam clamped at one side, the other free side facing a 2D granular flow in a horizontal cell moving at a constant velocity. We investigated the reconfiguration of the fibre as a function of the fibre's rigidity and of the granular packing fraction close but below the jamming in 2D. Imposing the fibre geometry like its length or thickness sets the critical buckling force the fibre is able to resist if it was not supported by lateral grains, while increasing the granular packing fraction might laterally consolidate the fibre and prevent it from buckling. But on the other side, the approach to jamming transition by increasing the granular packing fraction will be characterized by a dramatically increasing size of the cluster of connected grains forming a solid block acting against the fibre, which might promote the fibre's deflection. Thus, we investigated the granular flow fields, the fibre's deflexion as well as the forces experienced by the fibre and compared them with theoretical predictions from elastica for different loadings along the fibre. PMMH, CNRS UMR 7636, UPMC, ESPCI-ParisTech, 10 rue Vauquelin, 75231 Paris Cedex 05, France.
Investigation of the effect of wall friction on the flow rate in 2D and 3D Granular Flow
NASA Astrophysics Data System (ADS)
Carballo-Ramirez, Brenda; Pleau, Mollie; Easwar, Nalini; Birwa, Sumit; Shah, Neil; Tewari, Shubha
We have measured the mass flow rate of spherical steel spheres under gravity in vertical, straight-walled 2 and 3-dimensional hoppers, where the flow velocity is controlled by the opening size. Our measurements focus on the role of friction and its placement along the walls of the hopper. In the 2D case, an increase in the coefficient of static friction from μ = 0.2 to 0.6 is seen to decrease the flow rate significantly. We have changed the placement of frictional boundaries/regions from the front and back walls of the 2D hopper to the side walls and floor to investigate the relative importance of the different regions in determining the flow rate. Fits to the Beverloo equation show significant departure from the expected exponent of 1.5 in the case of 2D flow. In contrast, 3D flow rates do not show much dependence on wall friction and its placement. We compare the experimental data to numerical simulations of gravity driven hopper granular flow with varying frictional walls constructed using LAMMPS*. *http://lammps.sandia.gov Supported by NSF MRSEC DMR 0820506.
Moment Invariants for 2D Flow Fields via Normalization in Detail.
Bujack, Roxana; Hotz, Ingrid; Scheuermann, Gerik; Hitzer, Eckhard
2015-08-01
The analysis of 2D flow data is often guided by the search for characteristic structures with semantic meaning. One way to approach this question is to identify structures of interest by a human observer, with the goal of finding similar structures in the same or other datasets. The major challenges related to this task are to specify the notion of similarity and define respective pattern descriptors. While the descriptors should be invariant to certain transformations, such as rotation and scaling, they should provide a similarity measure with respect to other transformations, such as deformations. In this paper, we propose to use moment invariants as pattern descriptors for flow fields. Moment invariants are one of the most popular techniques for the description of objects in the field of image recognition. They have recently also been applied to identify 2D vector patterns limited to the directional properties of flow fields. Moreover, we discuss which transformations should be considered for the application to flow analysis. In contrast to previous work, we follow the intuitive approach of moment normalization, which results in a complete and independent set of translation, rotation, and scaling invariant flow field descriptors. They also allow to distinguish flow features with different velocity profiles. We apply the moment invariants in a pattern recognition algorithm to a real world dataset and show that the theoretical results can be extended to discrete functions in a robust way. PMID:26357255
Flow Quantification from 2D Phase Contrast MRI in Renal Arteries Using Clustering
NASA Astrophysics Data System (ADS)
Zöllner, Frank G.; Monnsen, Jan Ankar; Lundervold, Arvid; Rørvik, Jarle
We present an approach based on clustering to segment renal arteries from 2D PC Cine MR images to measure blood velocity and flow. Such information are important in grading renal artery stenosis and support the decision on surgical interventions like percutan transluminal angioplasty. Results show that the renal arteries could be extracted automatically and the corresponding velocity profiles could be calculated. Furthermore, the clustering could detect possible phase wrap effects automatically as well as differences in the blood flow patterns within the vessel.
Thermal stability effects on the separated flow over a steep 2-D hill
NASA Astrophysics Data System (ADS)
Zhang, W.; Porte-Agel, F.
2012-12-01
Transport of momentum and scalars in turbulent boundary-layer flows over complex topography has been of great interest in the atmospheric sciences and wind engineering communities. Applications include but are not limited to weather forecasting, air pollution dispersion, aviation safety control, and wind energy project planning. Linear models have been well accepted to predict boundary-layer flows over topography with gentle slope. However, once the slope of the topography is sufficientlyo steep that flow separation occurs, linear models are not applicable. Modeling the turbulent transport of momentum and scalars in such flows has to be achieved through non-linear models, such as Reynolds-averaged Navier-Stokes solvers and large-eddy simulations (LES). Dynamics of the separated boundary-layer flows over steep topography is affected by the shape and size of the topography, surface characteristics (e.g., roughness and temperature) and atmospheric thermal stability. Most wind-tunnel experiments of boundary-layer flows over idealized topography (e.g. 2-D or 3-D hills, axisymmetric bumps) do not take thermal stability effects into account due to difficulty of physical simulation. We conducted comprehensive experimental investigation of stably- and unstably- stratified boundary layers over a steep 2-D hill in the thermally-controlled boundary-layer wind tunnel at the Saint Anthony Falls Laboratory. The 2-D model hill has a steepest slope of 0.73 and its shape follows a cosine square function: h=Hcos^2 (πx/L) for -L/2 ≤ x ≤ L/2 , where the maximum height H is 7 cm and the total width L is 15 cm. High-resolution Particle Image Velocimetry (PIV) provides dynamic information of the separated shear layer, the recirculation zone and flow reattachment. Turbulent momentum and scalar (heat) fluxes were characterized up to the top of the thermal boundary layer using a triple-wire (cross-wire and cold-wire) anemometer. Results indicate that promoted and suppressed turbulence
Development of models for the two-dimensional, two-fluid code for sodium boiling NATOF-2D. [LMFBR
Zielinski, R.G.; Kazimi, M.S.
1981-09-01
Several features were incorporated into NATOF-2D, a two-dimensional, two fluid code developed at MIT for the purpose of analysis of sodium boiling transients under LMFBR conditions. They include improved interfacial mass, momentum and energy exchange rate models, and a cell-to-cell radial heat conduction mechanism which was calibrated by simulation of Westinghouse Blanket Heat Transfer Test Program Runs 544 and 545. Finally, a direct method of pressure field solution was implemented into a direct method of pressure field solution was implemented into NATOF-2D, replacing the iterative technique previously available, and resulted in substantially reduced computational costs.
Fluid dynamics and vibration of tube banks in fluid flow
Zukauskas, A.; Ulinskas, R.; Katinas, V.
1988-01-01
This work presents results derived in fluid dynamics, hydraulic drag and flow-induced vibrations within transverse and yawed tube banks. The studies encompass banks of smooth, rough and finned tubes at Reynolds numbers from 1 to 2x10/sup 6/. Highlighted in the text are fluid dynamic parameters of tube banks measured at inter-tube spaces and tube surfaces.
Value for controlling flow of cryogenic fluid
Knapp, Philip A.
1996-01-01
A valve is provided for accurately controlling the flow of cryogenic fluids such as liquid nitrogen. The valve comprises a combination of disc and needle valves affixed to a valve stem in such a manner that the disc and needle are free to rotate about the stem, but are constrained in lateral and vertical movements. This arrangement provides accurate and precise fluid flow control and positive fluid isolation.
Instability of plasma waves during relaxation of 2D turbulent flows
NASA Astrophysics Data System (ADS)
Kabantsev, A. A.; Drsicoll, C. F.
2015-11-01
We observe strong excitation of novel low-frequency z-dependent plasma waves (mθ = 0 ,kz = 1) , occurring during the nominally 2D relaxation of turbulent initial conditions (10 -100 interacting vortices) in strongly magnetized electron columns. This initial relaxation often results in ``2D vortex crystal'' states. Here we describe experiments showing the concomitant growth of ill-understood low-frequency plasma waves, probably due to ``leakage'' of 2D turbulent potential energy into z-dependent fluctuations. With plasma injection, the lowest regular Trivelpiece- Gould mode (mθ = 0 ,kz = 1) is observed at fTG (t) ~ 2 . 8 MHz and exponential decay time τTG ~ 1 msec. Also, we observe rapid exponential growth of a novel low-frequency mode with fLF (t) ~ 0 . 3 MHz, nominally also with mθ = 0 ,kz = 1 . In a few milliseconds (several tens of rotation times at B = 10kG), the LF-mode becomes highly nonlinear, developing up to a dozen temporal harmonics. When a LF-harmonic resonates with the decaying TG-mode, LF-mode energy is transferred into the TG-mode, and both modes remain at moderate amplitudes until the 2D turbulent relaxation abates (hundreds of rotation times). The ill-understood fLF is independent of B, even though the growth and duration times follow scale as B1 from the 2D flows. Supported by National Science Foundation Grant PHY-1414570, Department of Energy Grants DE-SC0008693.
Fluid Flow Control with Transformation Media
NASA Astrophysics Data System (ADS)
Urzhumov, Yaroslav A.; Smith, David R.
2011-08-01
We introduce a new concept for the manipulation of fluid flow around three-dimensional bodies. Inspired by transformation optics, the concept is based on a mathematical idea of coordinate transformations and physically implemented with anisotropic porous media permeable to the flow of fluids. In two situations—for an impermeable object placed either in a free-flowing fluid or in a fluid-filled porous medium—we show that the object can be coated with an inhomogeneous, anisotropic permeable medium, such as to preserve the flow that would have existed in the absence of the object. The proposed fluid flow cloak eliminates downstream wake and compensates viscous drag, hinting at the possibility of novel propulsion techniques.
TRENT2D WG: a smart web infrastructure for debris-flow modelling and hazard assessment
NASA Astrophysics Data System (ADS)
Zorzi, Nadia; Rosatti, Giorgio; Zugliani, Daniel; Rizzi, Alessandro; Piffer, Stefano
2016-04-01
Mountain regions are naturally exposed to geomorphic flows, which involve large amounts of sediments and induce significant morphological modifications. The physical complexity of this class of phenomena represents a challenging issue for modelling, leading to elaborate theoretical frameworks and sophisticated numerical techniques. In general, geomorphic-flows models proved to be valid tools in hazard assessment and management. However, model complexity seems to represent one of the main obstacles to the diffusion of advanced modelling tools between practitioners and stakeholders, although the UE Flood Directive (2007/60/EC) requires risk management and assessment to be based on "best practices and best available technologies". Furthermore, several cutting-edge models are not particularly user-friendly and multiple stand-alone software are needed to pre- and post-process modelling data. For all these reasons, users often resort to quicker and rougher approaches, leading possibly to unreliable results. Therefore, some effort seems to be necessary to overcome these drawbacks, with the purpose of supporting and encouraging a widespread diffusion of the most reliable, although sophisticated, modelling tools. With this aim, this work presents TRENT2D WG, a new smart modelling solution for the state-of-the-art model TRENT2D (Armanini et al., 2009, Rosatti and Begnudelli, 2013), which simulates debris flows and hyperconcentrated flows adopting a two-phase description over a mobile bed. TRENT2D WG is a web infrastructure joining advantages offered by the software-delivering model SaaS (Software as a Service) and by WebGIS technology and hosting a complete and user-friendly working environment for modelling. In order to develop TRENT2D WG, the model TRENT2D was converted into a service and exposed on a cloud server, transferring computational burdens from the user hardware to a high-performing server and reducing computational time. Then, the system was equipped with an
A fast and accurate method to predict 2D and 3D aerodynamic boundary layer flows
NASA Astrophysics Data System (ADS)
Bijleveld, H. A.; Veldman, A. E. P.
2014-12-01
A quasi-simultaneous interaction method is applied to predict 2D and 3D aerodynamic flows. This method is suitable for offshore wind turbine design software as it is a very accurate and computationally reasonably cheap method. This study shows the results for a NACA 0012 airfoil. The two applied solvers converge to the experimental values when the grid is refined. We also show that in separation the eigenvalues remain positive thus avoiding the Goldstein singularity at separation. In 3D we show a flow over a dent in which separation occurs. A rotating flat plat is used to show the applicability of the method for rotating flows. The shown capabilities of the method indicate that the quasi-simultaneous interaction method is suitable for design methods for offshore wind turbine blades.
Field Evaluation of a Novel 2D Preferential Flow Snowpack Hydrology Model
NASA Astrophysics Data System (ADS)
Leroux, N.; Pomeroy, J. W.; Kinar, N. J.
2015-12-01
Accurate estimation of snowmelt flux is of primary importance for runoff hydrograph prediction, which is used for water management and flood forecasting. Lateral flows and preferential flow pathways in porous media flow have proven critical for improving soil and groundwater flow models, but though many physically-based layered snowmelt models have been developed, only 1D matrix flow is accounted for in these models. Therefore, there is a need for snowmelt models that include these processes so as to examine the potential to improve snowmelt hydrological modelling. A 2D model is proposed that enables an improved understanding of energy and water flows within deep heterogeneous snowpacks, including those on slopes. A dual pathway theory is presented that simulates the formation of preferential flow paths, vertical and lateral water flows through the snow matrix and flow fingers, internal energy fluxes, melt, wet snow metamorphism, and internal refreezing. The dual pathway model utilizes an explicit finite volume method to solve for the energy and water flux equations over a non-orthogonal grid. It was run and evaluated using in-situ data collected from snowpit - accessed gravimetric, thermometric, photographic, and dielectric observations and novel non-invasive acoustic observations of layering, temperature, flowpath geometry, density and wetness at the Fortress Mountain Snow Laboratory, Alberta, Canada. The melt of a natural snowpack was artificially generated after detailed observation of snowpack initial conditions such as snow layer properties, temperature, and liquid water content. Snowpack ablation and liquid water content distribution over time were then measured and used for model parameterization and validation. Energy available at the snow surface and soil slope angle were set as mondel inputs. Model verification was based on snowpack property evolution. The heterogeneous flow model can be an important tool to help understand snowmelt flow processes, how
VISCOPLASTIC FLUID MODEL FOR DEBRIS FLOW ROUTING.
Chen, Cheng-lung
1986-01-01
This paper describes how a generalized viscoplastic fluid model, which was developed based on non-Newtonian fluid mechanics, can be successfully applied to routing a debris flow down a channel. The one-dimensional dynamic equations developed for unsteady clear-water flow can be used for debris flow routing if the flow parameters, such as the momentum (or energy) correction factor and the resistance coefficient, can be accurately evaluated. The writer's generalized viscoplastic fluid model can be used to express such flow parameters in terms of the rheological parameters for debris flow in wide channels. A preliminary analysis of the theoretical solutions reveals the importance of the flow behavior index and the so-called modified Froude number for uniformly progressive flow in snout profile modeling.
Frictional flow characteristics of microconvective flow for variable fluid properties
NASA Astrophysics Data System (ADS)
Kumar, Rajan; Mahulikar, Shripad P.
2015-12-01
The present work investigates the frictional flow characteristics of water flowing through a circular microchannel with variable fluid properties. The computational analysis reveals the importance of physical mechanisms due to variations in thermophysical fluid properties such as viscosity μ(T), thermal conductivity k(T) and density ρ(T) and also their contribution in the characteristics of frictional flow. Various combinations of thermophysical fluid properties have been used to find their effects on fluid friction. It is observed that the fluid friction attains the maximum value in the vicinity of the inlet and diminishes along the flow. The main reasons are attributed to this, (1) near the inlet, there is a flow undevelopment (the reverse process of flow development) due to μ(T) variation. (2) The viscosity of the water decreases with increasing temperature, which reduces fluid friction along the flow. It is noted that the skin friction coefficient (cf) reduces with increasing fluid mean velocity for a same value of constant wall heat flux ({q}{{w}}\\prime\\prime ). In the vicinity of the inlet, the deviation of Poiseuille number (Po) from 64 (constant properties solution) is also investigated in this paper. Additionally, the relationship between Reynolds number (Re) and cf, Po and Re have been proposed for different combinations of thermophysical fluid properties. This investigation also shows that the effect of fluid property variations on pressure drop is highly significant for microconvective water flow.
NASA Astrophysics Data System (ADS)
Sharapov, V. N.; Cherepanov, A. N.; Popov, V. N.; Bykova, V. G.
2012-11-01
A model describing two-dimensional (2D) dynamics of heat transfer in the fluid systems with a localized sink of a magmatic fluid into local fractured zones above the roof of crystallizing crustal intrusions is suggested. Numerical modeling of the migration of the phase boundaries in 2D intrusive chambers under retrograde boiling of magma with relatively high initial water content in the melt shows that, depending on the character of heat dissipation from a magmatic fluid into the host rock, two types of fluid magmatic systems can arise. (1) At high heat losses, the zoning of fluidogenic ore formation is determined by the changes in temperature of the rocks within the contact aureole of the intrusive bodies. These temperature variations are controlled by the migration of the phase boundaries in the cooling melt towards the center of the magmatic bodies from their contacts. (2) In the case of a localized sink of the magmatic fluid in different parts of the top of the intrusive chambers, a specific characteristic scenario of cooling of the magmatic bodies is probably implemented. In 2D systems with a heat transfer coefficient α k < 5 × 104 W/m2 K, an area with quasi-stationary phase boundaries develops close to the region of fluid drainage through the fractured zone in the intrusion. Therefore, as the phase boundaries contract to the sink zone of a fluid, specific thermal tubes arise, whose characteristics depend on the width of the fluid-conductive zone and the heat losses into the side rocks. (3) The time required for the intrusion to solidify varies depending on the particular position of the fluid conductor above the top of the magmatic body.
SAGE 2D and 3D Simulations of the Explosive Venting of Supercritical Fluids Through Porous Media
NASA Astrophysics Data System (ADS)
Weaver, R.; Gisler, G.; Svensen, H.; Mazzini, A.
2008-12-01
Magmatic intrusive events in large igneous provinces heat sedimentary country rock leading to the eventual release of volatiles. This has been proposed as a contributor to climate change and other environmental impacts. By means of numerical simulations, we examine ways in which these volatiles can be released explosively from depth. Gases and fluids cooked out of country rock by metamorphic heating may be confined for a time by impermeable clays or other barriers, developing high pressures and supercritical fluids. If confinement is suddenly breached (by an earthquake for example) in such a way that the fluid has access to porous sediments, a violent eruption of a non-magmatic mixture of fluid and sediment may result. Surface manifestations of these events could be hydrothermal vent complexes, kimberlite pipes, pockmarks, or mud volcanoes. These are widespread on Earth, especially in large igneous provinces, as in the Karoo Basin of South Africa, the North Sea off the Norwegian margin, and the Siberian Traps. We have performed 2D and 3D simulations with the Sage hydrocode (from Los Alamos and Science Applications International) of supercritical venting in a variety of geometries and configurations. The simulations show several different patterns of propagation and fracturing in porous or otherwise weakened overburden, dependent on depth, source conditions (fluid availability, temperature, and pressure), and manner of confinement breach. Results will be given for a variety of 2D and 3D simulations of these events exploring the release of volatiles into the atmosphere.
Measurement of Diffusion in Flowing Complex Fluids
Leonard, Edward F.; Aucoin, Christian P.; Nanne, Edgar E.
2006-01-01
A microfluidic device for the measurement of solute diffusion as well as particle diffusion and migration in flowing complex fluids is described. The device is particularly suited to obtaining diffusivities in such fluids, which require a desired flow state to be maintained during measurement. A method based on the Loschmidt diffusion theory and short times of exposure is presented to allow calculation of diffusivities from concentration differences in the flow streams leaving the cell. PMID:18560469
Two Dimensional Fluid Flow Models Offshore Southwestern Taiwan
NASA Astrophysics Data System (ADS)
Chen, L. W.; Wu, S. K.; Chi, W. C.; Liu, C. S.; Shyu, C. T.; Wang, Y. S.
2012-04-01
Fluid migration rates are important parameters for understanding the structural characteristics and evolution of the crustal tectonics and hydrocarbon exploration. However, they are difficult to measure on the seafloor. Dense distribution of bottom-simulating reflectors (BSRs) as the index of fluid existence to shed light on our study of the fluid migration. In this study, We acquired 2D fluid flow patterns in two potential gas hydrate prospect sites offshore southwestern Taiwan, and respectively modeled across Yung-An and Formosa ridge in N-S and E-W direction southwestern Taiwan. Temperature field in the shallow crust is used as a tracer to examine the fluid flow patterns. We use thermal information directly measured by thermal probes and topography data to develop the theoretical 2D temperature field using a thermal conduction model, which was derived from a finite element method. The discrepancy between the observed temperature data and the conductive model is attributed to advection heat transfer due to fluid migration. For Yung-An Ridge, we found the BSR-based temperatures are about 2oC higher than the conduction model in the following zones: (1) near a fault zone, (2) on the eastern flank where there are strong seismic reflectors in a pseudo 3D seismic dataset, (3) a seismic chimney zone. We interpret that there is possible active dewatering inside the accretionary prism to allow fluid to migrate upward here. For Formosa Ridge in the passive margin, the BSR-based temperatures are about 2oC colder than the theoretical model, especially on the flanks. We interpret that cold seawater is moving into the ridge from the flanks, cooling the ridge, and then some of the fluid is expelled at the ridge top. The shallow temperature fields are strongly affected by 2D or even 3D bathymetry effects. But we can still gain much information regarding fluid flow patterns through modeling. In the near future, we will extend such study into 3D. Keywords: fluid migration
Lattice Boltzmann simulations of 2D laminar flows past two tandem cylinders
NASA Astrophysics Data System (ADS)
Mussa, Alberto; Asinari, Pietro; Luo, Li-Shi
2009-03-01
We apply the lattice Boltzmann equation (LBE) with multiple-relaxation-time (MRT) collision model to simulate laminar flows in two-dimensions (2D). In order to simulate flows in an unbounded domain with the LBE method, we need to address two issues: stretched non-uniform mesh and inflow and outflow boundary conditions. We use the interpolated grid stretching method to address the need of non-uniform mesh. We demonstrate that various inflow and outflow boundary conditions can be easily and consistently realized with the MRT-LBE. The MRT-LBE with non-uniform stretched grids is first validated with a number of test cases: the Poiseuille flow, the flow past a cylinder asymmetrically placed in a channel, and the flow past a cylinder in an unbounded domain. We use the LBE method to simulate the flow past two tandem cylinders in an unbounded domain with Re = 100. Our results agree well with existing ones. Through this work we demonstrate the effectiveness of the MRT-LBE method with grid stretching.
Assessment of the Impacts of Compensation Flow Changes Upon Instream Habitat Using 2D Modelling
NASA Astrophysics Data System (ADS)
Mould, D. C.; Lane, S. N.; Christmas, M.
2004-05-01
Many millstone-grit rivers in northern England are impounded. In such cases the water company in the area has to release compensation flows from the reservoirs, traditionally to meet industrial needs: these flows are rarely set with ecology in mind; and have commonly involved constant flow. Dam overtopping may create spates, but spawning in many fish species is prompted by a spate flow in the early autumn when dams are rarely full enough to overtop. Such flows are important for fine sediment flushing and controlling the wetted useable area for spawning. Classical physical habitat modelling for instream habitat has been largely reliant upon 1D approaches, such as the Instream Flow Incremental Methodology (IFIM). Here we use a 2D finite element model (FESWMS), to simulate changes in instream habitat with variations in the compensation flow regimes. The spatial resolution of 2D models can be adapted to the scale of fish habitats so providing better representation of the reach-scale flow processes (such as slack water in the margins, wetting and drying) than the 1D case. The model is applied to the Rivers Rivelin and Loxley in Sheffield, Northern England. At the confluence of the two rivers, the compensation flow level is set at 30.6 Thousand Cubic Metres per Day (TCMD). Due to historical reasons, the compensation is not divided equally, as the Loxley receives 28 TCMD whilst the Rivelin receives only 2.6 TCMD. The model is used to simulate a transfer of 6 TCMD from the Loxley to the Rivelin. After validation, model predictions are combined with available habitat requirement data (e.g. velocity and depth needs) to develop an index of change in habitat suitability in terms of first order variables (e.g. velocity, depth and wetted useable area). This suggests that the change in compensation may significantly improve instream ecology in relation to macroinvertebrates, brown trout (Salmo trutta) and bullhead (Cottus gobio) in the Rivelin without causing detrimental impacts
NASA Astrophysics Data System (ADS)
Stone, James M.; Norman, Michael L.
1992-06-01
In this, the second of a series of three papers, we continue a detailed description of ZEUS-2D, a numerical code for the simulation of fluid dynamical flows in astrophysics including a self-consistent treatment of the effects of magnetic fields and radiation transfer. In this paper, we give a detailed description of the magnetohydrodynamical (MHD) algorithms in ZEUS-2D. The recently developed constrained transport (CT) algorithm is implemented for the numerical evolution of the components of the magnetic field for MHD simulations. This formalism guarantees the numerically evolved field components will satisfy the divergence-free constraint at all times. We find, however, that the method used to compute the electromotive forces must be chosen carefully to propagate accurately all modes of MHD wave families (in particular shear Alfvén waves). A new method of computing the electromotive force is developed using the method of characteristics (MOC). It is demonstrated through the results of an extensive series of MHD test problems that the resulting hybrid MOC-CT method provides for the accurate evolution of all modes of MHD wave families.
Visualizing vector field topology in fluid flows
NASA Technical Reports Server (NTRS)
Helman, James L.; Hesselink, Lambertus
1991-01-01
Methods of automating the analysis and display of vector field topology in general and flow topology in particular are discussed. Two-dimensional vector field topology is reviewed as the basis for the examination of topology in three-dimensional separated flows. The use of tangent surfaces and clipping in visualizing vector field topology in fluid flows is addressed.
Method and Apparatus for Measuring Fluid Flow
NASA Technical Reports Server (NTRS)
Arndt, G. Dickey (Inventor); Nguyen, Thanh X. (Inventor); Carl, James R. (Inventor)
1997-01-01
Method and apparatus for making measurements on fluids related to their complex permeability are disclosed. A microwave probe is provided for exposure to the fluids. The probe can be non-intrusive or can also be positioned at the location where measurements are to be made. The impedance of the probe is determined. in part. by the complex dielectric constant of the fluids at the probe. A radio frequency signal is transmitted to the probe and the reflected signal is phase and amplitude detected at a rapid rate for the purpose of identifying the fluids. Multiple probes may be selectively positioned to monitor the behavior of the fluids including their flow rate. Fluids may be identified as between two or more different fluids as well as multiple phases of the same fluid based on differences between their complex permittivities.
Untangling tracer trajectories and clarifying coherence in 2D flows using braid theory
NASA Astrophysics Data System (ADS)
Filippi, Margaux; Atis, Séverine; Thiffeault, Jean-Luc; Budišić, Marko; Allshouse, Michael; Peacock, Thomas
2014-11-01
Interpreting ocean surface transport is crucial to many areas of oceanography, ranging from marine ecology to pollution management. To better understand surface mixing, we investigate a braid theory method to detect transport barriers bounding coherent structures in two-dimensional fluid flows. Whereas most existing techniques rely on an extensive spatiotemporal knowledge of the flow field, we seek to identify these structures from sparse data sets involving trajectories of a few tracer particles or floats. We present the results of model and laboratory experimental studies to test the robustness and applicability of the braid theory method, and discuss the potential applicability to oceanic data sets.
Magnetic field applied to thermochemical non-equilibrium reentry flows in 2D - five species
NASA Astrophysics Data System (ADS)
Sávio de Góes Maciel, Edisson
2015-07-01
In this work, a study involving magnetic field actuation over reentry flows in thermochemical non-equilibrium is performed. The Euler and Navier-Stokes equations are studied. The proposed numerical algorithm is centred and second-order accurate. The hypersonic flow around a blunt body is simulated. Three time integration methods are tested. The reactive simulations involve Earth atmosphere of five species. The work of Gaitonde is the reference to couple the fluid dynamics and Maxwell equations of electromagnetism. The results have indicated that the Maciel scheme, using the Mavriplis dissipation model, yields the best prediction of the stagnation pressure.
Wavelet diagnostics of the flow control of unsteady separation on a 2D Wind Turbine Airfoil
NASA Astrophysics Data System (ADS)
Bai, Zhe; Lewalle, Jacques; Wang, Guannan; Glauser, Mark
2013-11-01
We investigated the aerodynamic characteristics of a 2D wind turbine airfoil. Unsteadiness was associated with the wake of a cylinder upstream of the airfoil. The experiments were conducted in both the baseline case, and with active closed-loop control on the suction surface of the airfoil. The data consisted of surface pressure time series. Continuous wavelet analysis gave the phase, band-pass filtered signals and envelope of harmonics of the fundamental shedding frequency. Coherence of pairs of signals was also used to map the flow characteristics. For the baseline and controlled case, we will report on the relation between phase of the leading edge fluctuations, unsteady flow separation and lift and drag coefficients. Our goal is to develop a more effective controller. The experiment was funded by DoE through University of Minnesota Wind Energy Consortium. Thanks for the support from the MAE department of Syracuse University.
Quasi-simultaneous interaction method for solving 2D boundary layer flows over plates and airfoils
NASA Astrophysics Data System (ADS)
Bijleveld, H. A.; Veldman, A. E. P.
2012-11-01
This paper studies unsteady 2D boundary layer flows over dented plates and a NACA 0012 airfoil. An inviscid flow is assumed to exist outside the boundary layer and is solved iteratively with the boundary layer flow together with the interaction method until a matching solution is achieved. Hereto a quasi-simultaneous interaction method is applied, in which the integral boundary layer equations are solved together with an interaction-law equation. The interaction-law equation is an approximation of the external flow and based on thin-airfoil theory. It is an algebraic relation between the velocity and displacement thickness. The interaction-law equation ensures that the eigenvalues of the system of equations do not have a sign change and that no singularities occur. Three numerical schemes are used to solve the boundary layer flow with the interaction method. These are: a standard scheme, a splitting method and a characteristics solver. All schemes use a finite difference discretization. The three schemes yield comparable results for the simulations carried out. The standard scheme is deviating most from the splitting and characteristics solvers. The results show that the eigenvalues remain positive, even in separation. As expected, the addition of the interaction-law equation prevents a sign change of the eigenvalues. The quasi-simultaneous interaction scheme is applicable to the three numerical schemes tested.
Self-Organization in 2D Traffic Flow Model with Jam-Avoiding Drive
NASA Astrophysics Data System (ADS)
Nagatani, Takashi
1995-04-01
A stochastic cellular automaton (CA) model is presented to investigate the traffic jam by self-organization in the two-dimensional (2D) traffic flow. The CA model is the extended version of the 2D asymmetric exclusion model to take into account jam-avoiding drive. Each site contains either a car moving to the up, a car moving to the right, or is empty. A up car can shift right with probability p ja if it is blocked ahead by other cars. It is shown that the three phases (the low-density phase, the intermediate-density phase and the high-density phase) appear in the traffic flow. The intermediate-density phase is characterized by the right moving of up cars. The jamming transition to the high-density jamming phase occurs with higher density of cars than that without jam-avoiding drive. The jamming transition point p 2c increases with the shifting probability p ja. In the deterministic limit of p ja=1, it is found that a new jamming transition occurs from the low-density synchronized-shifting phase to the high-density moving phase with increasing density of cars. In the synchronized-shifting phase, all up cars do not move to the up but shift to the right by synchronizing with the move of right cars. We show that the jam-avoiding drive has an important effect on the dynamical jamming transition.
NASA Astrophysics Data System (ADS)
LeVeque, R. J.; Motley, M. R.
2015-12-01
A series of tsunami wave basin experiments of flow through a scale model of Seaside, Oregon have been used as validation data for a 2015 benchmarking workshop hosted by the National Tsunami Mitigation Program, which focused on better understanding the ability of tsunami models to predict flow velocities and inundation depths following a coastal inundation event. As researchers begin to assess the safety of coastal infrastructures, proper assessment of tsunami-induced forces on coastal structures is critical. Hydrodynamic forces on these structures are fundamentally proportional to the local momentum flux of the fluid, and experimental data included momentum flux measurements at many instrumented gauge locations. The GeoClaw tsunami model, which solves the two-dimensional shallow water equations, was compared against other codes during the benchmarking workshop, and more recently a three-dimensional computational fluid dynamics model using the open-source OpenFOAM software has been developed and results from this model are being compared with both the experimental data and the 2D GeoClaw results. In addition, the 3D model allows for computation of fluid forces on the faces of structures, permitting an investigation of the common use of momentum flux as a proxy for these forces. This work aims to assess the potential to apply these momentum flux predictions locally within the model to determine tsunami-induced forces on critical structures. Difficulties in working with these data sets and cross-model comparisons will be discussed. Ultimately, application of the more computationally efficient GeoClaw model, informed by the 3D OpenFOAM models, to predict forces on structures at the community scale can be expected to improve the safety and resilience of coastal communities.
NASA Astrophysics Data System (ADS)
Pan, Li-Hua; Hou, Peng-Fei; Chen, Jia-Yun
2016-08-01
The 2D steady-state solutions regarding the expressions of stress and strain for fluid-saturated, orthotropic, poroelastic plane are derived in this paper. For this object, the general solutions of the corresponding governing equation are first obtained and expressed in harmonic functions. Based on these compact general solutions, the suitable harmonic functions with undetermined constants for line fluid source in the interior of infinite poroelastic body and a line fluid source on the surface of semi-infinite poroelastic body are presented, respectively. The fundamental solutions can be obtained by substituting these functions into the general solution, and the undetermined constants can be obtained by the continuous conditions, equilibrium conditions and boundary conditions.
Viscoelastic fluid flow in inhomogeneous porous media
Siginer, D.A.; Bakhtiyarov, S.I.
1996-09-01
The flow of inelastic and viscoelastic fluids in two porous media of different permeabilities and same priority arranged in series has been investigated both theoretically and experimentally. The fluids are an oil field spacer fluid and aqueous solutions of polyacrylamide. The porous medium is represented by a cylindrical tube randomly packed with glass spheres. Expressions for the friction factor and the resistance coefficient as a function of the Reynolds number have been developed both for shear thinning and viscoelastic fluids based on the linear fluidity and eight constant Oldroyd models, respectively. The authors show that the energy loss is higher if the viscoelastic fluid flows first through the porous medium with the smaller permeability rather than through the section of the cylinder with the larger permeability. This effect is not observed for Newtonian and shear thinning fluids flowing through the same configuration. Energy requirements for the same volume flow rate are much higher than a Newtonian fluid of the same zero shear viscosity as the polymeric solution. Energy loss increases with increasing Reynolds number at a fixed concentration. At a fixed Reynolds number, the loss is a strong function of the concentration and increases with increasing concentration. The behavior of all fluids is predicted qualitatively except the difference in energy requirements.
Method and device for measuring fluid flow
Atherton, Richard; Marinkovich, Phillip S.; Spadaro, Peter R.; Stout, J. Wilson
1976-11-23
This invention is a fluid flow measuring device for determining the coolant flow at the entrance to a specific nuclear reactor fuel region. The device comprises a plurality of venturis having the upstream inlet and throat pressure of each respectively manifolded together to provide one static pressure signal for each region monitored. The device provides accurate flow measurement with low pressure losses and uniform entrance and discharge flow distribution.
Debris Flow Hazard Map Simulation using FLO-2D For Selected Areas in the Philippines
NASA Astrophysics Data System (ADS)
Khallil Ferrer, Peter; Llanes, Francesca; dela Resma, Marvee; Realino, Victoriano, II; Obrique, Julius; Ortiz, Iris Jill; Aquino, Dakila; Narod Eco, Rodrigo; Mahar Francisco Lagmay, Alfredo
2014-05-01
On December 4, 2012, Super Typhoon Bopha wreaked havoc in the southern region of Mindanao, leaving 1,067 people dead and causing USD 800 million worth of damage. Classified as a Category 5 typhoon by the Joint Typhoon Warning Center (JTWC), Bopha brought intense rainfall and strong winds that triggered landslides and debris flows, particularly in Barangay (village) Andap, New Bataan municipality, in the southern Philippine province of Compostela Valley. The debris flow destroyed school buildings and covered courts and an evacuation center. Compostela Valley also suffered the most casualties of any province: 612 out of a total of 1,067. In light of the disaster in Compostela, measures were immediately devised to improve available geohazard maps to raise public awareness about landslides and debris flows. A debris flow is a very rapid to extremely rapid flow of saturated non-plastic debris in a steep channel. They are generated when heavy rainfall saturates sediments, causing them to flow down river channels within an alluvial fan situated at the base of the slope of a mountain drainage network. Many rural communities in the Philippines, such as Barangay Andap, are situated at the apex of alluvial fans and in the path of potential debris flows. In this study, we conducted simulations of debris flows to assess the risks in inhabited areas throughout the Philippines and validated the results in the field, focusing on the provinces of Pangasinan and Aurora as primary examples. Watersheds that drain in an alluvial fan using a 10-m resolution Synthetic Aperture Radar (SAR)-derived Digital Elevation Model (DEM) was first delineated, and then a 1 in 100-year rain return rainfall scenario for the watershed was used to simulate debris flows using FLO-2D, a flood-routing software. The resulting simulations were used to generate debris flow hazard maps which are consistent with danger zones in alluvial fans delineated previously from satellite imagery and available DEMs. The
Vlasov Fluid stability of a 2-D plasma with a linear magnetic field null
Kim, J.S.
1984-01-01
Vlasov Fluid stability of a 2-dimensional plasma near an O type magnetic null is investigated. Specifically, an elongated Z-pinch is considered, and applied to Field Reversed Configurations at Los Alamos National Laboratory by making a cylindrical approximation of the compact torus. The orbits near an elliptical O type null are found to be very complicated; the orbits are large and some are stochastic. The kinetic corrections to magnetohydrodynamics (MHD) are investigated by evaluating the expectation values of the growth rates of a Vlasov Fluid dispersion functional by using a set of trial functions based on ideal MHD. The dispersion functional involves fluid parts and orbit dependent parts. The latter involves phase integral of two time correlations. The phase integral is replaced by the time integral both for the regular and for the stochastic orbits. Two trial functions are used; one has a large displacement near the null and the other away from the null.
Engineering fluid flow using sequenced microstructures
NASA Astrophysics Data System (ADS)
Amini, Hamed; Sollier, Elodie; Masaeli, Mahdokht; Xie, Yu; Ganapathysubramanian, Baskar; Stone, Howard A.; di Carlo, Dino
2013-05-01
Controlling the shape of fluid streams is important across scales: from industrial processing to control of biomolecular interactions. Previous approaches to control fluid streams have focused mainly on creating chaotic flows to enhance mixing. Here we develop an approach to apply order using sequences of fluid transformations rather than enhancing chaos. We investigate the inertial flow deformations around a library of single cylindrical pillars within a microfluidic channel and assemble these net fluid transformations to engineer fluid streams. As these transformations provide a deterministic mapping of fluid elements from upstream to downstream of a pillar, we can sequentially arrange pillars to apply the associated nested maps and, therefore, create complex fluid structures without additional numerical simulation. To show the range of capabilities, we present sequences that sculpt the cross-sectional shape of a stream into complex geometries, move and split a fluid stream, perform solution exchange and achieve particle separation. A general strategy to engineer fluid streams into a broad class of defined configurations in which the complexity of the nonlinear equations of fluid motion are abstracted from the user is a first step to programming streams of any desired shape, which would be useful for biological, chemical and materials automation.
NASA Astrophysics Data System (ADS)
Magri, F.; Inbar, N.; Raggad, M.; Möller, S.; Siebert, C.; Möller, P.; Kuehn, M.
2014-12-01
Lake Kinneret (Lake Tiberias or Sea of Galilee) is the most important freshwater reservoir in the Northern Jordan Valley. Simulations that couple fluid flow, heat and mass transport are built to understand the mechanisms responsible for the salinization of this important resource. Here the effects of permeability distribution on 2D and 3D convective patterns are compared. 2D simulations indicate that thermal brine in Haon and some springs in the Yamourk Gorge (YG) are the result of mixed convection, i.e. the interaction between the regional flow from the bordering heights and thermally-driven flow (Magri et al., 2014). Calibration of the calculated temperature profiles suggests that the faults in Haon and the YG provides paths for ascending hot waters, whereas the fault in the Golan recirculates water between 1 and 2 km depths. At higher depths, faults induce 2D layered convection in the surrounding units. The 2D assumption for a faulted basin can oversimplify the system, and the conclusions might not be fully correct. The 3D results also point to mixed convection as the main mechanism for the thermal anomalies. However, in 3D the convective structures are more complex allowing for longer flow paths and residence times. In the fault planes, hydrothermal convection develops in a finger regime enhancing inflow and outflow of heat in the system. Hot springs can form locally at the surface along the fault trace. By contrast, the layered cells extending from the faults into the surrounding sediments are preserved and are similar to those simulated in 2D. The results are consistent with the theory from Zhao et al. (2003), which predicts that 2D and 3D patterns have the same probability to develop given the permeability and temperature ranges encountered in geothermal fields. The 3D approach has to be preferred to the 2D in order to capture all patterns of convective flow, particularly in the case of planar high permeability regions such as faults. Magri, F., et al., 2014
Fluid/structure interactions. Internal flows
NASA Astrophysics Data System (ADS)
Weaver, D. S.
1991-05-01
Flow-induced vibrations are found wherever structures are exposed to high velocity fluid flows. Internal flows are usually characterized by the close proximity of solid boundaries. There are surfaces against which separated flows may reattach, or from which pressure disturbances may be reflected resulting in acoustic resonance. When the fluid is a liquid, the close proximity of solid boundaries to a vibrating component can produce very high added mass effects. This paper presents three different experimental studies of flow-induced vibration problems associated with internal flows. The emphasis was on experimental techniques developed for understanding excitation mechanisms. In difficult flow-induced vibration problems, a useful experimental technique is flow visualization using a large scale model and strobe light triggered by the phenomenon being observed. This should be supported by point measurements of velocity and frequency spectra. When the flow excitation is associated with acoustic resonance, the sound can be fed back to enhance or eliminate the instability. This is potentially a very useful tool for studying and controlling fluid-structure interaction problems. Some flow-induced vibration problems involve a number of different excitation mechanisms and care must be taken to ensure that the mechanisms are properly identified. Artificially imposing structural vibrations or acoustic fields may induce flow structures not naturally present in the system.
Apparatus for measuring fluid flow
Smith, J.E.; Thomas, D.G.
Flow measuring apparatus includes a support loop having strain gages mounted thereon and a drag means which is attached to one end of the support loop and which bends the sides of the support loop and induces strains in the strain gages when a flow stream impacts thereon.
Apparatus for measuring fluid flow
Smith, Jack E.; Thomas, David G.
1984-01-01
Flow measuring apparatus includes a support loop having strain gages mounted thereon and a drag means which is attached to one end of the support loop and which bends the sides of the support loop and induces strains in the strain gages when a flow stream impacts thereon.
Directed flow fluid rinse trough
Kempka, Steven N.; Walters, Robert N.
1996-01-01
Novel rinse troughs accomplish thorough uniform rinsing. The tanks are suitable for one or more essentially planar items having substantially the same shape. The troughs ensure that each surface is rinsed uniformly. The new troughs also require less rinse fluid to accomplish a thorough rinse than prior art troughs.
Directed flow fluid rinse trough
Kempka, S.N.; Walters, R.N.
1996-07-02
Novel rinse troughs accomplish thorough uniform rinsing. The tanks are suitable for one or more essentially planar items having substantially the same shape. The troughs ensure that each surface is rinsed uniformly. The new troughs also require less rinse fluid to accomplish a thorough rinse than prior art troughs. 9 figs.
Fluid Flow Within Fractured Porous Media
Crandall, D.M.; Ahmadi, G.; Smith, D.H.; Bromhal, G.S.
2006-10-01
Fractures provide preferential flow paths to subterranean fluid flows. In reservoir scale modeling of geologic flows fractures must be approximated by fairly simple formulations. Often this is accomplished by assuming fractures are parallel plates subjected to an applied pressure gradient. This is known as the cubic law. An induced fracture in Berea sandstone has been digitized to perform numerical flow simulations. A commercially available computational fluid dynamics software package has been used to solve the flow through this model. Single phase flows have been compared to experimental works in the literature to evaluate the accuracy with which this model can be applied. Common methods of fracture geometry classification are also calculated and compared to experimentally obtained values. Flow through regions of the fracture where the upper and lower fracture walls meet (zero aperture) are shown to induce a strong channeling effect on the flow. This model is expanded to include a domain of surrounding porous media through which the flow can travel. The inclusion of a realistic permeability in this media shows that the regions of small and zero apertures contribute to the greatest pressure losses over the fracture length and flow through the porous media is most prevalent in these regions. The flow through the fracture is shown to be the largest contributor to the net flow through the media. From this work, a novel flow relationship is proposed for flow through fractured media.
NASA Astrophysics Data System (ADS)
Tucciarelli, T.
2012-12-01
A new methodology for the solution of irrotational 2D flow problems in domains with strongly unstructured meshes is presented. A fractional time step procedure is applied to the original governing equations, solving consecutively a convective prediction system and a diffusive corrective system. The non linear components of the problem are concentrated in the prediction step, while the correction step leads to the solution of a linear system, of the order of the number of computational cells. A MArching in Space and Time (MAST) approach is applied for the solution of the convective prediction step. The major advantages of the model, as well as its ability to maintain the solution monotonicity even in strongly irregular meshes, are briefly described. The algorithm is applied to the solution of diffusive shallow water equations in a simple domain.
An Experimental Study of Flow Separation over 2D Transverse Grooves
NASA Astrophysics Data System (ADS)
Jones, Emily; Lang, Amy; Afroz, Farhana; Wheelus, Jennifer; Smith, Drew
2011-11-01
A shark's scales help to reduce drag over its body by controlling boundary layer separation over its skin. It is theorized that the scales bristle when encountering a reversing flow, thereby trapping vortices between the scales, creating a partial slip condition over the surface and inducing turbulence augmentation in the boundary layer. In an attempt to replicate and study these effects, a spinning cylinder was used in a water tunnel to induce separation over a flat plate with 2 millimeter square 2D transverse grooves. The results were compared to separation occurring over a flat plate without grooves using DPIV. The angular speed of the cylinder was varied. The observed delays in separation, changes in separation bubble shedding frequency and other effects upon the boundary layer are discussed.
Flow-induced protein crystallization: Macroscopic effects on 2D crystals
NASA Astrophysics Data System (ADS)
Young, James; Posada, David; Hirsa, Amir; Lopez, Juan
2012-11-01
Proteins must first be crystallized before their molecular structure can be studied in detail. However, crystallizing protein is a challenging task which is often met with limited success. Although 2-D protein crystals at the air/water interface are usually obtained under quiescent conditions, it was recently shown that crystallization can be enhanced by a shearing flow. Here we examine the relationship between Reynolds number and the crystal growth process using the deep-channel surface viscometer geometry. It consists of an annular region bounded by stationary inner and outer cylinders and driven by a constant rotation of the floor. The interfacial velocity measurements are compared to Navier-Stokes computations with the Boussinesq-Scriven surface model. The interfacial film is lifted onto a solid substrate, and the protein crystals are observed via optical and atomic force microscopy. For a particular protein surface concentration, a Reynolds number threshold has been identified for flow-induced crystallization. This flow geometry also allows for the determination of the surface shear viscosity, which provides a quantitative measure of the mesoscale interactions associated with protein crystallization.
Micro PIV measurements of turbulent flow over 2D structured roughness
NASA Astrophysics Data System (ADS)
Hartenberger, Joel; Perlin, Marc
2015-11-01
We investigate the turbulent boundary layer over surfaces with 2D spanwise square and triangular protrusions having nominal heights of 100 - 300 microns for Reynolds numbers ranging from Reτ ~ 1500 through Reτ ~ 4500 using a high speed, high magnification imaging system. Micro PIV analysis gives finely resolved velocity fields of the flow (on the order of 10 microns between vectors) enabling a detailed look at the inner region as well as the flow in the immediate vicinity of the roughness elements. Additionally, planar PIV with lower resolution is performed to capture the remainder of the boundary layer to the freestream flow. Varying the streamwise distance between individual roughness elements from one to ten times the nominal heights allows investigation of k-type and d-type roughness in both the transitionally rough and fully rough regimes. Preliminary results show a shift in the mean velocity profile similar to the results of previous studies. Turbulent statistics will be presented also. The authors would like to acknowledge the support of NAVSEA which funded this project through the Naval Engineering Education Center (NEEC).
Analysis of High-Speed Rotating Flow in 2D Polar (r - θ)Coordinate
NASA Astrophysics Data System (ADS)
Pradhan, S.
2016-03-01
The generalized analytical model for the radial boundary layer in a high-speed rotating cylinder is formulated for studying the gas flow field due to insertion of mass, momentum and energy into the rotating cylinder in the polar (r - θ) plane. The analytical solution includes the sixth order differential equation for the radial boundary layer at the cylindrical curved surface in terms of master potential (χ) , which is derived from the equations of motion in a polar (r - θ) plane. The linearization approximation (Wood & Morton, J. Fluid Mech-1980; Pradhan & Kumaran, J. Fluid Mech-2011; Kumaran & Pradhan, J. Fluid Mech-2014) is used, where the equations of motion are truncated at linear order in the velocity and pressure disturbances to the base flow, which is a solid-body rotation. Additional assumptions in the analytical model include constant temperature in the base state (isothermal condition), and high Reynolds number, but there is no limitation on the stratification parameter. In this limit, the gas flow is restricted to a boundary layer of thickness (Re (1 / 3) R) at the wall of the cylinder. Here, the stratification parameter A = √ ((mΩ 2R2) / (2kB T)) . This parameter Ais the ratio of the peripheral speed, ΩR , to the most probable molecular speed, √(2 k_B T/m), the Reynolds number Re = (ρ _w ΩR2 / μ) , where m is the molecular mass, Ω and R are the rotational speed and radius of the cylinder, k_B is the Boltzmann constant, T is the gas temperature, ρ_w is the gas density at wall, and μ is the gas viscosity. The analytical solutions are then compared with direct simulation Monte Carlo (DSMC) simulations.
Pattern formation in flowing electrorheological fluids.
von Pfeil, Karl; Graham, Michael D; Klingenberg, Daniel J; Morris, Jeffrey F
2002-05-01
A two-fluid continuum model is developed to describe mass transport in electro- and magnetorheological suspensions. The particle flux is related to the field-induced stresses. Solutions of the resulting mass balance show column formation in the absence of flow, and stripe formation when a suspension is subjected simultaneously to an applied electric field and shear flow. PMID:12005727
Instrument continuously measures density of flowing fluids
NASA Technical Reports Server (NTRS)
Jacobs, R. B.; Macinko, J.; Miller, C. E.
1967-01-01
Electromechanical densitometer continuously measures the densities of either single-phase or two-phase flowing cryogenic fluids. Measurement is made on actual flow. The instrument operates on the principle that the mass of any vibrating system is a primary factor in determining the dynamic characteristics of the system.
Electromagnetic probe technique for fluid flow measurements
NASA Technical Reports Server (NTRS)
Arndt, G. D.; Carl, J. R.
1994-01-01
The probes described herein, in various configurations, permit the measurement of the volume fraction of two or more fluids flowing through a pipe. Each probe measures the instantaneous relative dielectric constant of the fluid in immediate proximity. As long as separation of the relative dielectric constant of the fluid is possible, several or even many fluids can be measured in the same flow stream. By using multiple probes, the velocity of each fluid can generally be determined as well as the distribution of each constituent in the pipe. The values are determined by statistical computation. There are many potential applications for probes of this type in industry and government. Possible NASA applications include measurements of helium/hydrazine flow during rocket tests at White Sands, liquid/gas flow in hydrogen or oxygen lines in Orbiter engines, and liquid/gaseous Freon flow in zero gravity tests with the KS135 aircraft at JSC. Much interest has been shown recently by the oil industry. In this a good method is needed to measure the fractions of oil, water, and natural gas flowing in a pipeline and the velocity of each. This particular problem involves an extension of what has been developed to date and our plans to solve this problem will be discussed herein.
Fundamental Processes of Atomization in Fluid-Fluid Flows
NASA Technical Reports Server (NTRS)
McCready, M. J.; Chang, H.-C.; Leighton, D. T.
2001-01-01
This report outlines the major results of the grant "Fundamental Processes of Atomization in Fluid-Fluid Flows." These include: 1) the demonstration that atomization in liquid/liquid shear flow is driven by a viscous shear instability that triggers the formation of a long thin sheet; 2) discovery of a new mode of interfacial instability for oscillatory two-layer systems whereby a mode that originates within the less viscous liquid phase causes interfacial deformation as the oscillation proceeds; 3) the demonstration that rivulet formation from gravity front occurs because the local front shape specified by gravity and surface tension changes from a nose to a wedge geometry, thus triggering a large increase in viscous resistance; and 4) extension of the studies on nonlinear wave evolution on falling films and in stratified flow, particularly the evolution towards large-amplitude solitary waves that tend to generate drops.
Manifest: A computer program for 2-D flow modeling in Stirling machines
NASA Technical Reports Server (NTRS)
Gedeon, David
1989-01-01
A computer program named Manifest is discussed. Manifest is a program one might want to use to model the fluid dynamics in the manifolds commonly found between the heat exchangers and regenerators of Stirling machines; but not just in the manifolds - in the regenerators as well. And in all sorts of other places too, such as: in heaters or coolers, or perhaps even in cylinder spaces. There are probably nonStirling uses for Manifest also. In broad strokes, Manifest will: (1) model oscillating internal compressible laminar fluid flow in a wide range of two-dimensional regions, either filled with porous materials or empty; (2) present a graphics-based user-friendly interface, allowing easy selection and modification of region shape and boundary condition specification; (3) run on a personal computer, or optionally (in the case of its number-crunching module) on a supercomputer; and (4) allow interactive examination of the solution output so the user can view vector plots of flow velocity, contour plots of pressure and temperature at various locations and tabulate energy-related integrals of interest.
Comparing a 2D fluid model of the DC planar magnetron cathode to experiments
Garcia, M.
1996-05-01
Planar magnetron cathodes have arching magnetic field lines which concentrate plasma density near the electrode surface. This enhances the ion bombardment of the surface and the yield of sputtered atoms. Magnetron cathodes are used in the Plasma Electrode Pockels Cell (PEPC) devices of the Laser Program because they provide for significantly higher conduction than do glow discharges. An essential feature of magnetron cathodes is that the vector product of the perpendicular electric field, E[sub y], with the parallel component of the magnetic field, B[sub x], forms a closed track with a circulating current along the cathode surface. An analytical, 2D, two component, quasi-neutral, continuum model yields formulas for the plasma density, the total and component current densities, the electric field, and the positive electrical potential, between the cathode surface and a distant, uniform plasma. For a specific gas, the free parameters are electron temperature, gas number density, and total current. The model is applied to the interpretation of experimental data from the PEPC device, as well as a small vacuum facility for testing magnetron cathodes. Finally, the model has been applied to generate cross sectional views of a PEPC magnetron cathode track.
Thermochemical Nonequilibrium 2D Modeling of Nitrogen Inductively Coupled Plasma Flow
NASA Astrophysics Data System (ADS)
Yu, Minghao; Yusuke, Takahashi; Hisashi, Kihara; Ken-ichi, Abe; Kazuhiko, Yamada; Takashi, Abe; Satoshi, Miyatani
2015-09-01
Two-dimensional (2D) numerical simulations of thermochemical nonequilibrium inductively coupled plasma (ICP) flows inside a 10-kW inductively coupled plasma wind tunnel (ICPWT) were carried out with nitrogen as the working gas. Compressible axisymmetric Navier-Stokes (N-S) equations coupled with magnetic vector potential equations were solved. A four-temperature model including an improved electron-vibration relaxation time was used to model the internal energy exchange between electron and heavy particles. The third-order accuracy electron transport properties (3rd AETP) were applied to the simulations. A hybrid chemical kinetic model was adopted to model the chemical nonequilibrium process. The flow characteristics such as thermal nonequilibrium, inductive discharge, effects of Lorentz force were made clear through the present study. It was clarified that the thermal nonequilibrium model played an important role in properly predicting the temperature field. The prediction accuracy can be improved by applying the 3rd AETP to the simulation for this ICPWT. supported by Grant-in-Aid for Scientific Research (No. 23560954), sponsored by the Japan Society for the Promotion of Science
Hydrodynamics of embedded planets' first atmospheres - I. A centrifugal growth barrier for 2D flows
NASA Astrophysics Data System (ADS)
Ormel, Chris W.; Kuiper, Rolf; Shi, Ji-Ming
2015-01-01
In the core accretion paradigm of planet formation, gas giants only form a massive atmosphere after their progenitors exceeded a threshold mass: the critical core mass. Most (exo)planets, being smaller and rock/ice-dominated, never crossed this line. Nevertheless, they were massive enough to attract substantial amounts of gas from the disc, while their atmospheres remained in pressure-equilibrium with the disc. Our goal is to characterize the hydrodynamical properties of the atmospheres of such embedded planets and the implications for their (long-term) evolution. In this paper - the first in series - we start to investigate the properties of an isothermal and inviscid flow past a small, embedded planet by conducting local, 2D hydrodynamical simulations. Using the PLUTO code, we confirm that the flow is steady and bound. This steady outcome is most apparent for the log-polar grid (with the grid spacing proportional to the distance from the planet). For low-mass planets, Cartesian grids are somewhat less efficient as they have difficulty to follow the circular, large speeds in the deep atmosphere. Relating the amount of rotation to the gas fraction of the atmosphere, we find that more massive atmospheres rotate faster - a finding consistent with Kelvin's circulation theorem. Rotation therefore limits the amount of gas that planets can acquire from the nebula. Dependent on the Toomre-Q parameter of the circumstellar disc, the planet's atmosphere will reach Keplerian rotation before self-gravity starts to become important.
Accurate Angle Estimator for High-Frame-Rate 2-D Vector Flow Imaging.
Villagomez Hoyos, Carlos Armando; Stuart, Matthias Bo; Hansen, Kristoffer Lindskov; Nielsen, Michael Bachmann; Jensen, Jorgen Arendt
2016-06-01
This paper presents a novel approach for estimating 2-D flow angles using a high-frame-rate ultrasound method. The angle estimator features high accuracy and low standard deviation (SD) over the full 360° range. The method is validated on Field II simulations and phantom measurements using the experimental ultrasound scanner SARUS and a flow rig before being tested in vivo. An 8-MHz linear array transducer is used with defocused beam emissions. In the simulations of a spinning disk phantom, a 360° uniform behavior on the angle estimation is observed with a median angle bias of 1.01° and a median angle SD of 1.8°. Similar results are obtained on a straight vessel for both simulations and measurements, where the obtained angle biases are below 1.5° with SDs around 1°. Estimated velocity magnitudes are also kept under 10% bias and 5% relative SD in both simulations and measurements. An in vivo measurement is performed on a carotid bifurcation of a healthy individual. A 3-s acquisition during three heart cycles is captured. A consistent and repetitive vortex is observed in the carotid bulb during systoles. PMID:27093598
Fluid flow nozzle energy harvesters
NASA Astrophysics Data System (ADS)
Sherrit, Stewart; Lee, Hyeong Jae; Walkemeyer, Phillip; Winn, Tyler; Tosi, Luis Phillipe; Colonius, Tim
2015-04-01
Power generation schemes that could be used downhole in an oil well to produce about 1 Watt average power with long-life (decades) are actively being developed. A variety of proposed energy harvesting schemes could be used to extract energy from this environment but each of these has their own limitations that limit their practical use. Since vibrating piezoelectric structures are solid state and can be driven below their fatigue limit, harvesters based on these structures are capable of operating for very long lifetimes (decades); thereby, possibly overcoming a principle limitation of existing technology based on rotating turbo-machinery. An initial survey [1] identified that spline nozzle configurations can be used to excite a vibrating piezoelectric structure in such a way as to convert the abundant flow energy into useful amounts of electrical power. This paper presents current flow energy harvesting designs and experimental results of specific spline nozzle/ bimorph design configurations which have generated suitable power per nozzle at or above well production analogous flow rates. Theoretical models for non-dimensional analysis and constitutive electromechanical model are also presented in this paper to optimize the flow harvesting system.
Heat Flow Partitioning Between Continents and Oceans - from 2D to 3D
NASA Astrophysics Data System (ADS)
Moresi, L. N.; Cooper, C. M.; Lenardic, A.
2010-12-01
Scalings derived from thermal network theory explain how the presence of continents can influence the Earth’s overall heat loss. Intuitively, it may seem that increasing the proportion of a planet’s surface area covered by continents would decrease the efficiency of heat transfer given that continents do not participate in convective overturn. However, this ignores the potential feedback between the insulating effect of continents and the temperature-dependent viscosity of the mantle (Lenardic et al, 2005, Cooper et al, 2007). When this feedback is considered, a clear regime exists in which the partial stagnation and insulation of the surface by buoyant continental crust can lead to an increase in heat flow compared to the uninsulated case. The numerical results used to verify the scalings have mostly been conducted in two dimensions in order to cover a very wide range of Rayleigh number, fraction of continental coverage, and continental thickness. However as more recent results show that the configuration of the crust also plays a role in determining the heat flow partitioning and global heat flow (See Lenardic et al, “Continents, Super-Continents, Mantle Thermal Mixing, and Mantle Thermal Isolation” in this session), we have begun to repeat this exhaustive and exhausting 2D study in 3D. Cooper, C.M., A. Lenardic, and L.-N. Moresi "Effects of continental insulation and the partioning of heat producing elements on the Earth's heat loss." Geophys. Res. Lett., 33 ,10.1029, 2006. Lenardic, A., L.-N. Moresi, A.M. Jellinek, and M. Manga "Continental insulation, mantle cooling, and the surface area of oceans and continents." Earth Planet. Sci. Lett., 234 ,317-333, 2005.
MEANS FOR VISUALIZING FLUID FLOW PATTERNS
Lynch, F.E.; Palmer, L.D.; Poppendick, H.F.; Winn, G.M.
1961-05-16
An apparatus is given for determining both the absolute and relative velocities of a phosphorescent fluid flowing through a transparent conduit. The apparatus includes a source for exciting a narrow trsnsverse band of the fluid to phosphorescence, detecting means such as a camera located downstream from the exciting source to record the shape of the phosphorescent band as it passes, and a timer to measure the time elapsed between operation of the exciting source and operation of the camera.
Two-fluid equilibrium with flow: FLOW2
Guazzotto, L.; Betti, R.
2015-09-15
The effects of finite macroscopic velocities on axisymmetric ideal equilibria are examined using the two-fluid (ions and electrons) model. A new equilibrium solver, the code FLOW2, is introduced for the two-fluid model and used to investigate the importance of various flow patterns on the equilibrium of tight aspect ratio (NSTX) and regular tokamak (DIII-D) configurations. Several improvements to the understanding and calculation of two-fluid equilibria are presented, including an analytical and numerical proof of the single-fluid and static limits of the two-fluid model, a discussion of boundary conditions, a user-friendly free-function formulation, and the explicit evaluation of velocity components normal to magnetic surfaces.
Electro-osmotic flow in bicomponent fluids
NASA Astrophysics Data System (ADS)
Bazarenko, Andrei; Sega, Marcello
The electroosmotic flow (EOF) is a widely used technique that uses the action of external electric fields on solvated ions to move fluids around in microfluidics devices. For homogeneous fluids, the characteristics of the flow can be well approximated by simple analytical models, but in multicomponent systems such as oil-in-water droplets one has to rely to numerical simulations. The purpose of this study is to investigate physical properties of the EOF in a bicomponent fluid by solving the coupled equations of motions of explicit ions in interaction with a continuous model of the flow. To do so we couple the hydrodynamics equations as solved by a Shan-Chen Lattice-Boltzmann method to the molecular dynamics of the ions. The presence of explicit ions allows us to go beyond the simple Poisson-Boltzmann approximations, and investigate a variety of EOF regimes. ETN-COLLDENSE (H2020-MCSA-ITN-2014, Grant No. 642774).
Tomographic reconstruction of stratified fluid flow.
Winters, K B; Rouseff, D
1993-01-01
A method for imaging a moving fluid is proposed and evaluated by numerical simulation. A cross-section of a three-dimensional fluid is probed by high-frequency acoustic waves from several different directions. Assuming straight-ray geometric acoustics, the time of flight depends on both the scaler sound speed and the vector fluid velocity. By appropriately combining travel times, projections of both the sound speed and the velocity are isolated. The sound speed is reconstructed using the standard filtered backprojection algorithm. Though complete inversion of velocity is not possible, sufficient information is available to recover the component of fluid vorticity transverse to the plane of insonification. A new filtered backprojection algorithm for vorticity is developed and implemented. To demonstrate the inversion procedure, a 3-D stratified fluid is simulated and travel time data are calculated by path integration. These data are then inverted to recover both the scaler sound speed and the vorticity of the evolving flow. PMID:18263153
Gabbour, Maya; Schnell, Susanne; Jarvis, Kelly; Robinson, Joshua D.; Markl, Michael
2015-01-01
Background Doppler echocardiography (echo) is the reference standard for blood flow velocity analysis, and two-dimensional (2-D) phase-contrast magnetic resonance imaging (MRI) is considered the reference standard for quantitative blood flow assessment. However, both clinical standard-of-care techniques are limited by 2-D acquisitions and single-direction velocity encoding and may make them inadequate to assess the complex three-dimensional hemodynamics seen in congenital heart disease. Four-dimensional flow MRI (4-D flow) enables qualitative and quantitative analysis of complex blood flow in the heart and great arteries. Objectives The objectives of this study are to compare 4-D flow with 2-D phase-contrast MRI for quantification of aortic and pulmonary flow and to evaluate the advantage of 4-D flow-based volumetric flow analysis compared to 2-D phase-contrast MRI and echo for peak velocity assessment in children and young adults. Materials and methods Two-dimensional phase-contrast MRI of the aortic root, main pulmonary artery (MPA), and right and left pulmonary arteries (RPA, LPA) and 4-D flow with volumetric coverage of the aorta and pulmonary arteries were performed in 50 patients (mean age: 13.1±6.4 years). Four-dimensional flow analyses included calculation of net flow and regurgitant fraction with 4-D flow analysis planes similarly positioned to 2-D planes. In addition, 4-D flow volumetric assessment of aortic root/ascending aorta and MPA peak velocities was performed and compared to 2-D phase-contrast MRI and echo. Results Excellent correlation and agreement were found between 2-D phase-contrast MRI and 4-D flow for net flow (r=0.97, P<0.001) and excellent correlation with good agreement was found for regurgitant fraction (r= 0.88, P<0.001) in all vessels. Two-dimensional phase-contrast MRI significantly underestimated aortic (P= 0.032) and MPA (P<0.001) peak velocities compared to echo, while volumetric 4-D flow analysis resulted in higher (aortic: P=0
Full 2D observation of water surface elevation from SWOT under different flow conditions
NASA Astrophysics Data System (ADS)
Domeneghetti, Alessio; Schumann, Guy; Rui, Wei; Durand, Michael; Pavelsky, Tamlin
2016-04-01
The upcoming Surface Water and Ocean Topography (SWOT) satellite mission is a joint project of NASA, Centre National d'Etudes Spatiales (CNES, France), the Canadian Space Agency, and the Space Agency of the UK that will provide a first global, high-resolution observation of ocean and terrestrial water surface heights. Characterized by an observation swath of 120 km and an orbit repeat interval of about 21 days, SWOT will provide unprecedented bi-dimensional observations of rivers wider than 50-100 m. Despite many research activities that have investigated potential uses of remotely sensed data from SWOT, potentials and limitations of the spatial observations provided by the satellite mission for flood modeling still remain poorly understood and investigated. In this study we present a first analysis of the spatial observation of water surface elevation that is expected from SWOT for a 140 km reach of the middle-lower portion of the Po River, in Northern Italy. The river stretch is characterized by a main channel varying from 200-500 m in width and a floodplain that can be as wide as 5 km and that is delimited by a system of major embankments. The reconstruction of the hydraulic behavior of the Po River is performed by means of a quasi-2d model built with detailed topographic and bathymetric information (LiDAR, 2 m resolution), while the simulation of the spatial observation sensed by SWOT is performed with a SWOT simulator that mimics the satellite sensor characteristics. Referring to water surface elevations associated with different flow conditions (maximum, minimum and average flow reproduced by means of the quasi-2d numerical model) this work provides a first characterization of the spatial observations provided by SWOT and highlights the strengths and limitations of the expected products. By referring to a real river reach the analysis provides a credible example of the type of spatial observations that will be available after launch of SWOT and offers a first
Method and Apparatus for Measuring Fluid Flow
NASA Technical Reports Server (NTRS)
Arndt, G. Dickey (Inventor); Nguyen, Than X. (Inventor); Carl, James R. (Inventor)
1995-01-01
The invention is a method and apparatus for monitoring the presence, concentration, and the movement of fluids. It is based on utilizing electromagnetic measurements of the complex permittivity of the fluids for detecting and monitoring the fluid. More particularly the apparatus uses one or more microwave probes which are placed at the locations where the measurements are to be made. A radio frequency signal is transmitted to the probe and the reflected signal is phase and amplitude detected at a rapid rate for the purpose of identifying the fluids, based on their dielectric constant at the probe. The apparatus can be used for multiple purposes including measures of flow rates, turbulence, dispersion, fluid identification, and changes in flow conditions of multiple fluids or multiple states of a single fluid in a flowline or a holding container. The apparatus includes a probe consisting of two electrical conductors separated by an insulator. A radio frequency signal is communicated to the probe and is reflected back from the portion of the probe exposed to the fluid. The radio frequency signal also provides a reference signal. An oscillator generates a second signal which combined with each of the reference signal and the reflected signal to produce signals of lower frequencies to facilitate filtering and amplifying those signals. The two signals are then mixed in a detector to produce an output signal that is representative of the phase and amplitude change caused by the reflection of the signal at the probe exposed to the fluid. The detector may be a dual phase detector that provides two such output signals that are in phase quadrature. A phase shifter may be provided for selectively changing the phase of the reference signal to improve the sensitivity of at least one of the output signals for more accurate readings and/or for calibration purposes. The two outputs that are in quadrature with respect to each other may be simultaneously monitored to account for
NASA Astrophysics Data System (ADS)
Choquard, Ph.; Vuffray, M.
2014-10-01
The coupling between dilatation and vorticity, two coexisting and fundamental processes in fluid dynamics (Wu et al., 2006, pp. 3, 6) is investigated here, in the simplest cases of inviscid 2D isotropic Burgers and pressureless Euler-Coriolis fluids respectively modeled by single vortices confined in compressible, local, inertial and global, rotating, environments. The field equations are established, inductively, starting from the equations of the characteristics solved with an initial Helmholtz decomposition of the velocity fields namely a vorticity free and a divergence free part (Wu et al., 2006, Sects. 2.3.2, 2.3.3) and, deductively, by means of a canonical Hamiltonian Clebsch like formalism (Clebsch, 1857, 1859), implying two pairs of conjugate variables. Two vector valued fields are constants of the motion: the velocity field in the Burgers case and the momentum field per unit mass in the Euler-Coriolis one. Taking advantage of this property, a class of solutions for the mass densities of the fluids is given by the Jacobian of their sum with respect to the actual coordinates. Implementation of the isotropy hypothesis entails a radial dependence of the velocity potentials and of the stream functions associated to the compressible and to the rotational part of the fluids and results in the cancellation of the dilatation-rotational cross terms in the Jacobian. A simple expression is obtained for all the radially symmetric Jacobians occurring in the theory. Representative examples of regular and singular solutions are shown and the competition between dilatation and vorticity is illustrated. Inspired by thermodynamical, mean field theoretical analogies, a genuine variational formula is proposed which yields unique measure solutions for the radially symmetric fluid densities investigated. We stress that this variational formula, unlike the Hopf-Lax formula, enables us to treat systems which are both compressible and rotational. Moreover in the one
Fluid Flow in An Evaporating Droplet
NASA Technical Reports Server (NTRS)
Hu, H.; Larson, R.
1999-01-01
Droplet evaporation is a common phenomenon in everyday life. For example, when a droplet of coffee or salt solution is dropped onto a surface and the droplet dries out, a ring of coffee or salt particles is left on the surface. This phenomenon exists not only in everyday life, but also in many practical industrial processes and scientific research and could also be used to assist in DNA sequence analysis, if the flow field in the droplet produced by the evaporation could be understood and predicted in detail. In order to measure the fluid flow in a droplet, small particles can be suspended into the fluid as tracers. From the ratio of gravitational force to Brownian force a(exp 4)(delta rho)(g)/k(sub B)T, we find that particle's tendency to settle is proportional to a(exp 4) (a is particle radius). So, to keep the particles from settling, the droplet size should be chosen to be in a range 0.1 -1.0 microns in experiments. For such small particles, the Brownian force will affect the motion of the particle preventing accurate measurement of the flow field. This problem could be overcome by using larger particles as tracers to measure fluid flow under microgravity since the gravitational acceleration g is then very small. For larger particles, Brownian force would hardly affect the motion of the particles. Therefore, accurate flow field could be determined from experiments in microgravity. In this paper, we will investigate the fluid flow in an evaporating droplet under normal gravity, and compare experiments to theories. Then, we will present our ideas about the experimental measurement of fluid flow in an evaporating droplet under microgravity.
Messinian Salinity Crisis and basin fluid flow
NASA Astrophysics Data System (ADS)
Bertoni, Claudia; Cartwight, Joe
2014-05-01
Syn- and post-depositional movement of fluids through sediments is one of the least understood aspects in the evolution of a basin. The conventional hydrostratigraphic view on marine sedimentary basins assumes that compactional and meteoric groundwater fluid circulation drives fluid movement and defines its timing. However, in the past few years, several examples of instantaneous and catastrophic release of fluids have been observed even through low-permeability sediments. A particularly complex case-study involves the presence of giant salt bodies in the depocentres of marine basins. Evaporites dramatically change the hydrostratigraphy and fluid-dynamics of the basin, and influence the P/T regimes, e.g. through changes in the geothermal gradient and in the compaction of underlying sediments. Our paper reviews the impact of the Messinian Salinity Crisis (MSC) and evaporites on fluid flow in the Mediterranean sub-basins. The analysis of geological and geophysical sub-surface data provides examples from this basin, and the comparison with analogues in other well-known evaporitic provinces. During the MSC, massive sea-level changes occurred in a relatively limited time interval, and affected the balance of fluid dynamics, e.g. with sudden release or unusual trapping of fluids. Fluid expulsion events are here analysed and classified in relation to the long and short-term effects of the MSC. Our main aim is to build a framework for the correct identification of the fluid flow-related events, and their genetic mechanisms. On basin margins, where evaporites are thin or absent, the sea-level changes associated with the MSC force a rapid basinward shift of the mixing zone of meteoric/gravity flow and saline/compactional flow, 100s-km away from its pre-MSC position. This phenomenon changes the geometry of converging flows, creates hydraulic traps for fluids, and triggers specific diagenetic reactions in pre-MSC deep marine sediments. In basin-centre settings, unloading and
Fluid flow meter for measuring the rate of fluid flow in a conduit
NASA Technical Reports Server (NTRS)
White, P. R. (Inventor)
1986-01-01
A tube fluid flow rate meter consists of a reservoir divided by flexible diaphragm into two separate isolated compartments. The incoming and outgoing tubes open into the compartments. The orifice is sized to allow maximum tube fluid flow. Opposing compression springs are secured within the two compartments on opposite sides of the orifice to maintain orifice position when the tube fluid pressure is zero. A tapered element is centered in, and extends through the orifice into the compartment, leaving an annular opening between the element and the perimeter of the oriface. The size varies as the diaphragm flexes with changes in the tube fluid pressure to change the fluid flow through the opening. The light source directs light upon the element which in turn scatters the light through the opening into the compartment. The light detector in the compartment senses the scattered light to generate a signal indicating the amount of fluid.
Time-Dependent 2D Modeling of Magnetron Plasma Torch in Turbulent Flow
NASA Astrophysics Data System (ADS)
Li, Lincun; Xia, Weidong
2008-06-01
A theoretical model is presented to describe the electromagnetic, heat transfer and fluid flow phenomena within a magnetron plasma torch and in the resultant plume, by using a commercial computational fluid dynamics (CFD) code FLUENT. Specific calculations are presented for a pure argon system (i.e., an argon plasma discharging into an argon environment), operated in a turbulent mode. An important finding of this work is that the external axial magnetic field (AMF) may have a significant effect on the behavior of arc plasma and thus affects the resulting plume. The AMF impels the plasma to retract axially and expand radially. As a result, the plasma intensity distribution on the cross section of torch seems to be more uniform. Numerical results also show that with AMF, the highest plasma temperature decreases and the anode arc root moves upstream significantly, while the current density distribution at the anode is more concentrated with a higher peak value. In addition, the use of AMF then induces a strong backflow at the torch spout and its magnitude increases with the AMF strength but decreases with the inlet gas velocity.
Method and apparatus for controlling fluid flow
Miller, J.R.
1980-06-27
A method and apparatus for precisely controlling the rate (and hence amount) of fluid flow are given. The controlled flow rate is finely adjustable, can be extremely small (on the order of microliter-atmospheres per second), can be adjusted to zero (flow stopped), and is stable to better than 1% with time. The dead volume of the valve can be made arbitrarily small, in fact essentially zero. The valve employs no wearing mechanical parts (including springs, stems, or seals). The valve is finely adjustable, has a flow rate dynamic range of many decades, can be made compatible with any fluid, and is suitable for incorporation into an open or closed loop servo-control system.
Fluid flow along faults in carbonate rocks
NASA Astrophysics Data System (ADS)
Romano, Valentina; Battaglia, Maurizio; Bigi, Sabina
2015-04-01
The study of fluid flow in fractured rocks plays a key role in reservoir management, including CO2 sequestration and waste isolation. We present a mathematical model of fluid flow in a fault zone, based on field data acquired in Majella Mountain, in the Central Apennines (Italy). The Majella is a thrust related, asymmetric, box shaped anticline. The mountain carbonate outcrops are part of a lower Cretaceous-Miocene succession, covered by a siliciclastic sequence of lower Pliocene age. We study a fault zone located in the Bolognano Formation (Oligo-Miocene age) and exposed in the Roman Valley Quarry near the town of Lettomanoppello, in the northern sector of the Majella Mountain. This is one of the best places in the Apennines to investigate a fault zone and has been the subject of numerous field studies. Faults are mechanical and permeability heterogeneities in the upper crust, so they strongly influence fluid flow. The distribution of the main components (core, damage zone) can lead a fault zone to act as a conduit, a barrier or a combined conduit-barrier system. We integrated existing and our own structural surveys of the area to better identify the major fault features (e.g., kind of fractures, statistical properties, geometry and pertrophysical characteristics). Our analytical model describe the Bolognano Formation using a dual porosity/dual permeability model: global flow occurs through the fracture network only, while rock matrix contain the majority of fluid storage and provide fluid drainage to the fractures. Pressure behavior is analyzed by examining the pressure drawdown curves, the derivative plots and the effects of the characteristic parameters. The analytical model has been calibrated against published data on fluid flow and pressure distribution in the Bolognano Formation.
Titan2D Based Pyroclastic Flows Hazard Maps for Santa Ana Volcano, El Salvador
NASA Astrophysics Data System (ADS)
Bajo, J. V.; Martinez-Hackert, B.; Escobar, C. D.; Gutierrez, R. E.
2009-05-01
Santa Ana Volcano is located in the Apaneca Volcanic Field located to the west of El Salvador, Central America. It is one the six active volcanoes monitor by the Servicios Nacionales de Estudios Territoriales (SNET) in El Salvador, out of twenty that are considered active in this small country by Smithsonian definition. The Santa Ana Volcano is surrounded by rural communities in its proximal areas and in its close distal areas by the second largest city of the country. On October 1st 2005, after a few months of increased fumarolic and seismic activity, it erupted generating a 10 km high steam and ash plume, reportedly seen by some aircraft and estimated using photography by SNET members. Ash was deposited to the west, north-west part of the country, following typical wind pattern for the region, as well as small pyroclastic flows and major lahars in its eastern part. Coffee plantations were lost, as was some crop of coffee in the following season. However, to the west the ash fertilized the land and resulted in an enhanced harvest of coffee beans. Only 2 people were killed from the Blast, thanks to the auto evacuation of proximal communities. Whilst the last eruption had a relatively low human life toll, a stronger eruption spells havoc almost certainly for the region. At this moment no exhaustive study and understanding exists of the pyroclastic flows generated by the Santa Ana Volcano nor a map for this particular hazard. This study proposes the use of Titan2D for those two purposes, using a DEM generated by the SNET using topographic maps as well as DEMs generated using Advanced Spaceborne Thermal Emission and Reflection Radiometer Images (ASTER).
NASA Astrophysics Data System (ADS)
Fairman, Stephen J.; Johnson, Joseph A.; Walkiewicz, Thomas A.
2003-09-01
Graphical analysis of experimental data that exhibit exponential behavior is typically postponed at many institutions until students are able to understand the theory underlying the concept of radioactive decay or of RC time constants in ac circuits. In 1960 Smithson and Pinkston described a laboratory exercise that uses the flow of water from a vertical column through a long horizontal capillary tube as a source of data that models radioactive decay. Many institutions have used this experiment simply as an early introduction to exponential behavior without reference to radioactive decay or ac circuits. Greenslade2 recently described a modification of this experiment to demonstrate the concept of secular equilibrium in radioactive decay. This paper presents results of similar experiments, but visual measurements are replaced in this work by data obtained with modern sensors interfaced to a computer. Experiments are described from simple exponential decay to an analogue of the complex interactions of three nuclides in a radioactive-series decay chain.
NASA Astrophysics Data System (ADS)
Zavolzhenskii, M. V.
1982-09-01
Boussinesq equations are used in studying the spectral problem of the stability loss in the equilibrium state of a rotating layer of viscous fluid subjected to temperature inversion. It is shown that this loss can take the form of eddy flows localized around the axis of rotation. It is noted that flows of this type have properties similar to those of waterspouts, tornados, and other vortices.
Automated analysis for fluid flow topology
NASA Technical Reports Server (NTRS)
Helman, James; Hesselink, Lambertus
1989-01-01
A new approach for visualizing vector data sets was developed by reducing the original vector field to a set of critical points and their connections, and was applied to fluid flow data sets. The critical point representation allows for considerable reduction in the data complexity. The representations are displayed as surfaces which are much simpler than the original data set, yet retain all the pertinent flow topology information. It is suggested that topological representations may be useful for database comparison.
Maximal mixing by incompressible fluid flows
NASA Astrophysics Data System (ADS)
Seis, Christian
2013-12-01
We consider a model for mixing binary viscous fluids under an incompressible flow. We prove the impossibility of perfect mixing in finite time for flows with finite viscous dissipation. As measures of mixedness we consider a Monge-Kantorovich-Rubinstein transportation distance and, more classically, the H-1 norm. We derive rigorous a priori lower bounds on these mixing norms which show that mixing cannot proceed faster than exponentially in time. The rate of the exponential decay is uniform in the initial data.
Advanced designs for fluid flow visualization
NASA Technical Reports Server (NTRS)
1978-01-01
Research was carried out on existing and new designs for minimally intrusive measurement of flow fields in the Geophysical Fluid Flow Cell and the proposed Atmospheric General Circulation Experiment. The following topics are discussed: (1) identification and removal of foreign particles, (2) search for higher dielectric photochromic solutions, (3) selection of uv light source, (4) analysis of refractive techniques and (5) examination of fresnel lens applicability.
Geophysical Fluid Flow Cell (GFFC) Simulation
NASA Technical Reports Server (NTRS)
1999-01-01
These simulations of atmospheric flow use the same experimental parameters but started with slightly different initial conditions in the model. The simulations were part of data analysis for the Geophysical Fluid Flow Cell (GFFC), a planet in a test tube apparatus flown on Spacelab to mimic the atmospheres on gas giant planets and stars. (Credit: Dr. Tim Miller of Global Hydrology and Climate Center at the Marshall Space Flight Center)
Classification of the stratified fluid flows regimes around a square cylinder
NASA Astrophysics Data System (ADS)
Gushchin, V. A.; Matyushin, P. V.
2015-10-01
The 2D density stratified (in vertical direction) viscous fluid flows around a square cylinder with diameter d (moving in horizontal direction with the velocity U) have been simulated on the basis of the Navier-Stokes equations in the Boussinesq approximation. For solving of the Navier-Stokes equations the Splitting on physical factors Method for Incompressible Fluid flows (SMIF) with hybrid explicit finite difference scheme (second-order accuracy in space, minimum scheme viscosity and dispersion, monotonous) has been used. The numerical method SMIF has been successfully applied for solving of the different problems: 2D and 3D separated homogeneous and stratified fluid flows around a sphere and a circular cylinder; the flows with free surface including regimes with broken surface wave; the air, heat and mass transfer in the clean rooms. At the present paper the original refined classification of 2D stratified viscous fluid flow regimes around a square cylinder at Re ≤ 200 has been obtained and the interesting fluid flows with two hanging vortices in the wake and with two wavy hanging sheets of density (connected with two hanging vortices) have been investigated in details at Fr = 0.1, Re = 50, where Re = U.d/ν is the Reynolds number, Fr = U/(N.d) is the internal Froude number, ν is the kinematical viscosity coefficient, N is the buoyancy frequency.
Finite scale equations for compressible fluid flow
Margolin, Len G
2008-01-01
Finite-scale equations (FSE) describe the evolution of finite volumes of fluid over time. We discuss the FSE for a one-dimensional compressible fluid, whose every point is governed by the Navier-Stokes equations. The FSE contain new momentum and internal energy transport terms. These are similar to terms added in numerical simulation for high-speed flows (e.g. artificial viscosity) and for turbulent flows (e.g. subgrid scale models). These similarities suggest that the FSE may provide new insight as a basis for computational fluid dynamics. Our analysis of the FS continuity equation leads to a physical interpretation of the new transport terms, and indicates the need to carefully distinguish between volume-averaged and mass-averaged velocities in numerical simulation. We make preliminary connections to the other recent work reformulating Navier-Stokes equations.
The Geophysical Fluid Flow Cell Experiment
NASA Technical Reports Server (NTRS)
Hart, J. E.; Ohlsen, D.; Kittleman, S.; Borhani, N.; Leslie, F.; Miller, T.
1999-01-01
The Geophysical Fluid Flow Cell (GFFC) experiment performed visualizations of thermal convection in a rotating differentially heated spherical shell of fluid. In these experiments dielectric polarization forces are used to generate a radially directed buoyancy force. This enables the laboratory simulation of a number of geophysically and astrophysically important situations in which sphericity and rotation both impose strong constraints on global scale fluid motions. During USML-2 a large set of experiments with spherically symmetric heating were carried out. These enabled the determination of critical points for the transition to various forms of nonaxisymmetric convection and, for highly turbulent flows, the transition latitudes separating the different modes of motion. This paper presents a first analysis of these experiments as well as data on the general performance of the instrument during the USML-2 flight.
Cryogenic fluid flow instabilities in heat exchangers
NASA Technical Reports Server (NTRS)
Fleming, R. B.; Staub, F. W.
1969-01-01
Analytical and experimental investigation determines the nature of oscillations and instabilities that occur in the flow of two-phase cryogenic fluids at both subcritical and supercritical pressures in heat exchangers. Test results with varying system parameters suggest certain design approaches with regard to heat exchanger geometry.
Volcanic termor: Nonlinear excitation by fluid flow
NASA Astrophysics Data System (ADS)
Julian, Bruce R.
1994-06-01
A nonlinear process analogous to the excitation mechanism of musical wind instruments and human vocal cords can explain many characteristics of volcanic tremor, including (1) periodic and 'chaotic' oscillations, with peaked and irregular spectra respectively, (2) rapid pulsations in eruptions occurring at the same frequency as tremor, (3) systematic changes in tremor amplitude as channel geometry evolves during an eruption, (4) the period doubling reported for Hawaiian deep tremor, and (5) the fact that the onset of termor can be either gradual or abrupt. Volcanic 'long-period' earthquakes can be explained as oscillations excited by transient disturbances produced by nearby earthquakes, fluid heterogeneity, or changes in channel geometry, when the magma flow rate is too low to excite continuous tremor. A simple lumped-parameter tremor model involving the flow of an incompressible viscous fluid through a channel with movable elastic walls leads to a third-order system of nonlinear ordinary differential equations. For different driving fluid pressures, numerical solutions exhibit steady flow, simple limit-cycle oscillations, a cascade of period-doubling subharmonic bifurcations, and chaotic oscillations controlled by a strange attractor of Rossler type. In this model, tremor occurs most easily at local constrictions, and fluid discharge is lower than would occur in unstable steady flow.
14 CFR 23.1095 - Carburetor deicing fluid flow rate.
Code of Federal Regulations, 2014 CFR
2014-01-01
... 14 Aeronautics and Space 1 2014-01-01 2014-01-01 false Carburetor deicing fluid flow rate. 23.1095... Induction System § 23.1095 Carburetor deicing fluid flow rate. (a) If a carburetor deicing fluid system is used, it must be able to simultaneously supply each engine with a rate of fluid flow, expressed...
14 CFR 23.1095 - Carburetor deicing fluid flow rate.
Code of Federal Regulations, 2013 CFR
2013-01-01
... 14 Aeronautics and Space 1 2013-01-01 2013-01-01 false Carburetor deicing fluid flow rate. 23.1095... Induction System § 23.1095 Carburetor deicing fluid flow rate. (a) If a carburetor deicing fluid system is used, it must be able to simultaneously supply each engine with a rate of fluid flow, expressed...
14 CFR 23.1095 - Carburetor deicing fluid flow rate.
Code of Federal Regulations, 2012 CFR
2012-01-01
... 14 Aeronautics and Space 1 2012-01-01 2012-01-01 false Carburetor deicing fluid flow rate. 23.1095... Induction System § 23.1095 Carburetor deicing fluid flow rate. (a) If a carburetor deicing fluid system is used, it must be able to simultaneously supply each engine with a rate of fluid flow, expressed...
14 CFR 23.1095 - Carburetor deicing fluid flow rate.
Code of Federal Regulations, 2011 CFR
2011-01-01
... 14 Aeronautics and Space 1 2011-01-01 2011-01-01 false Carburetor deicing fluid flow rate. 23.1095... Induction System § 23.1095 Carburetor deicing fluid flow rate. (a) If a carburetor deicing fluid system is used, it must be able to simultaneously supply each engine with a rate of fluid flow, expressed...
14 CFR 23.1095 - Carburetor deicing fluid flow rate.
Code of Federal Regulations, 2010 CFR
2010-01-01
... 14 Aeronautics and Space 1 2010-01-01 2010-01-01 false Carburetor deicing fluid flow rate. 23.1095... Induction System § 23.1095 Carburetor deicing fluid flow rate. (a) If a carburetor deicing fluid system is used, it must be able to simultaneously supply each engine with a rate of fluid flow, expressed...
A qualitative view of cryogenic fluid injection into high speed flows
NASA Technical Reports Server (NTRS)
Hendricks, R. C.; Schlumberger, J.; Proctor, M.
1991-01-01
The injection of supercritical pressure, subcritical temperature fluids, into a 2-D, ambient, static temperature and static pressure supersonic tunnel and free jet supersonic nitrogen flow field was observed. Observed patterns with fluid air were the same as those observed for fluid nitrogen injected into the tunnel at 90 deg to the supersonic flow. The nominal injection pressure was of 6.9 MPa and tunnel Mach number was 2.7. When injected directly into and opposing the tunnel exhaust flow, the observed patterns with fluid air were similar to those observed for fluid nitrogen but appeared more diffusive. Cryogenic injection creates a high density region within the bow shock wake but the standoff distance remains unchanged from the gaseous value. However, as the temperature reaches a critical value, the shock faded and advanced into the supersonic stream. For both fluids, nitrogen and air, the phenomena was completely reversible.
Enhanced fluid flow through nanoscale carbon pipes.
Whitby, Max; Cagnon, Laurent; Thanou, Maya; Quirke, Nick
2008-09-01
Recent experimental and theoretical studies demonstrate that pressure driven flow of fluids through nanoscale ( d < 10 nm) carbon pores occurs 4 to 5 orders of magnitude faster than predicted by extrapolation from conventional theory. Here, we report experimental results for flow of water, ethanol, and decane through carbon nanopipes with larger inner diameters (43 +/- 3 nm) than previously investigated. We find enhanced transport up to 45 times theoretical predictions. In contrast to previous work, in our systems, decane flows faster than water. These nanopipes were composed of amorphous carbon deposited from ethylene vapor in alumina templates using a single step fabrication process. PMID:18680352
Onset of turbulence from the receptivity stage of fluid flows.
Sengupta, T K; Bhaumik, S
2011-10-01
The traditional viewpoint of fluid flow considers the transition to turbulence to occur by the secondary and nonlinear instability of wave packets, which have been created experimentally by localized harmonic excitation. The boundary layer has been shown theoretically to support spatiotemporal growing wave fronts by Sengupta, Rao, and Venkatasubbaiah [Phys. Rev. Lett. 96, 224504 (2006)] by a linear mechanism, which is shown here to grow continuously, causing the transition to turbulence. Here, we track spatiotemporal wave fronts to a nonlinear turbulent state by solving the full 2D Navier-Stokes equation, without any limiting assumptions. Thus, this is the only demonstration of deterministic disturbances evolving from a receptivity stage to the full turbulent flow. This is despite the prevalent competing conjectures of the event being three-dimensional and/or stochastic in nature. PMID:22107294
Granular Material Flows with Interstitial Fluid Effects
NASA Technical Reports Server (NTRS)
Hunt, Melany L.; Brennen, Christopher E.
2004-01-01
The research focused on experimental measurements of the rheological properties of liquid-solid and granular flows. In these flows, the viscous effects of the interstitial fluid, the inertia of the fluid and particles, and the collisional interactions of the particles may all contribute to the flow mechanics. These multiphase flows include industrial problems such as coal slurry pipelines, hydraulic fracturing processes, fluidized beds, mining and milling operation, abrasive water jet machining, and polishing and surface erosion technologies. In addition, there are a wide range of geophysical flows such as debris flows, landslides and sediment transport. In extraterrestrial applications, the study of transport of particulate materials is fundamental to the mining and processing of lunar and Martian soils and the transport of atmospheric dust (National Research Council 2000). The recent images from Mars Global Surveyor spacecraft dramatically depict the complex sand and dust flows on Mars, including dune formation and dust avalanches on the slip-face of dune surfaces. These Aeolian features involve a complex interaction of the prevailing winds and deposition or erosion of the sediment layer; these features make a good test bed for the verification of global circulation models of the Martian atmosphere.
Behaviour of a rimmed elliptical inclusion in 2D slow incompressible viscous flow
NASA Astrophysics Data System (ADS)
Mancktelow, N. S.
2012-04-01
The shape preferred orientation of natural populations of inclusions (or "porphyroclasts") is often inconsistent with predictions from established analytical theory for inclusions with coherent boundaries (e.g., Pennacchioni et al. 2001). A totally incoherent or slipping interface can explain observed stable back-rotated (or antithetic) orientations but not the observed cut-off axial ratio, below which inclusions still rotate. However, this behaviour is reproduced by a rimmed inclusion with a rim viscosity that is not infinitely weak but still weaker than the matrix (e.g., Schmid and Podladchikov 2005; Johnson et al. 2009). In this study, finite-element numerical modelling (FEM) is employed to investigate this system in 2D over a very wide parameter space, from a viscosity ratio (relative to the matrix) of the inclusion from 106 to 1, the rim from 10-6 to 1, the axial ratio from 1.00025 to 20, and the rim thickness from 5% to 20%. Theoretical consideration of a concentric elliptical inclusion and ellipse reduces the number of scalar values to be determined to fully characterize the system to two: one for the rate of stretch of the inclusion and one for the rate of rotation. From these two values, the rotation and stretching rate can be calculated for any orientation and 2D background flow field. For effectively rigid particles, the cut-off axial ratio between rotation and stabilization is determined by the remaining two parameters, namely the rim viscosity and the thickness, with low rim viscosity or thick rims promoting stabilization. The shape fabric of a population of particles in a high strain shear zone, presented as a typical Rf/φ plot, can be forward modelled using an initial value Ordinary Differential Equation (ODE) approach. Because the rim does not remain elliptical to high strain, this method cannot accurately model the behaviour of individual inclusions. However, a statistical approach, allowing variation in rim viscosity, which is also a proxy for
Approaches to Modeling Coupled Flow and Reaction in a 2-D Cementation Experiment
Steefel, Carl; Cochepin, B.; Trotignon, L.; Bildstein, O.; Steefel, C.; Lagneau, V.; van der Lee, J.
2008-04-01
Porosity evolution at reactive interfaces is a key process that governs the evolution and performances of many engineered systems that have important applications in earth and environmental sciences. This is the case, for example, at the interface between cement structures and clays in deep geological nuclear waste disposals. Although in a different transport regime, similar questions arise for permeable reactive barriers used for biogeochemical remediation in surface environments. The COMEDIE project aims at investigating the coupling between transport, hydrodynamics and chemistry when significant variations of porosity occur. The present work focuses on a numerical benchmark used as a design exercise for the future COMEDIE-2D experiment. The use of reactive transport simulation tools like Hytec and Crunch provides predictions of the physico-chemical evolutions that are expected during the future experiments in laboratory. Focus is given in this paper on the evolution during the simulated experiment of precipitate, permeability and porosity fields. A first case is considered in which the porosity is constant. Results obtained with Crunch and Hytec are in relatively good agreement. Differences are attributable to the models of reactive surface area taken into account for dissolution/precipitation processes. Crunch and Hytec simulations taking into account porosity variations are then presented and compared. Results given by the two codes are in qualitative agreement, with differences attributable in part to the models of reactive surface area for dissolution/precipitation processes. As a consequence, the localization of secondary precipitates predicted by Crunch leads to lower local porosities than for predictions obtained by Hytec and thus to a stronger coupling between flow and chemistry. This benchmark highlights the importance of the surface area model employed to describe systems in which strong porosity variations occur as a result of dissolution
Using heteroclinic orbits to quantify topological entropy in fluid flows
NASA Astrophysics Data System (ADS)
Sattari, Sulimon; Chen, Qianting; Mitchell, Kevin A.
2016-03-01
Topological approaches to mixing are important tools to understand chaotic fluid flows, ranging from oceanic transport to the design of micro-mixers. Typically, topological entropy, the exponential growth rate of material lines, is used to quantify topological mixing. Computing topological entropy from the direct stretching rate is computationally expensive and sheds little light on the source of the mixing. Earlier approaches emphasized that topological entropy could be viewed as generated by the braiding of virtual, or "ghost," rods stirring the fluid in a periodic manner. Here, we demonstrate that topological entropy can also be viewed as generated by the braiding of ghost rods following heteroclinic orbits instead. We use the machinery of homotopic lobe dynamics, which extracts symbolic dynamics from finite-length pieces of stable and unstable manifolds attached to fixed points of the fluid flow. As an example, we focus on the topological entropy of a bounded, chaotic, two-dimensional, double-vortex cavity flow. Over a certain parameter range, the topological entropy is primarily due to the braiding of a period-three orbit. However, this orbit does not explain the topological entropy for parameter values where it does not exist, nor does it explain the excess of topological entropy for the entire range of its existence. We show that braiding by heteroclinic orbits provides an accurate computation of topological entropy when the period-three orbit does not exist, and that it provides an explanation for some of the excess topological entropy when the period-three orbit does exist. Furthermore, the computation of symbolic dynamics using heteroclinic orbits has been automated and can be used to compute topological entropy for a general 2D fluid flow.
Using heteroclinic orbits to quantify topological entropy in fluid flows.
Sattari, Sulimon; Chen, Qianting; Mitchell, Kevin A
2016-03-01
Topological approaches to mixing are important tools to understand chaotic fluid flows, ranging from oceanic transport to the design of micro-mixers. Typically, topological entropy, the exponential growth rate of material lines, is used to quantify topological mixing. Computing topological entropy from the direct stretching rate is computationally expensive and sheds little light on the source of the mixing. Earlier approaches emphasized that topological entropy could be viewed as generated by the braiding of virtual, or "ghost," rods stirring the fluid in a periodic manner. Here, we demonstrate that topological entropy can also be viewed as generated by the braiding of ghost rods following heteroclinic orbits instead. We use the machinery of homotopic lobe dynamics, which extracts symbolic dynamics from finite-length pieces of stable and unstable manifolds attached to fixed points of the fluid flow. As an example, we focus on the topological entropy of a bounded, chaotic, two-dimensional, double-vortex cavity flow. Over a certain parameter range, the topological entropy is primarily due to the braiding of a period-three orbit. However, this orbit does not explain the topological entropy for parameter values where it does not exist, nor does it explain the excess of topological entropy for the entire range of its existence. We show that braiding by heteroclinic orbits provides an accurate computation of topological entropy when the period-three orbit does not exist, and that it provides an explanation for some of the excess topological entropy when the period-three orbit does exist. Furthermore, the computation of symbolic dynamics using heteroclinic orbits has been automated and can be used to compute topological entropy for a general 2D fluid flow. PMID:27036190
Flow behaviour of extremely bidisperse magnetizable fluids
NASA Astrophysics Data System (ADS)
Susan-Resiga, Daniela; Bica, Doina; Vékás, L.
2010-10-01
In this paper we investigated the rheological and magnetorheological behaviours of an extremely bidisperse (nano-micro) magnetizable fluid (sample D1) for comparison of a commercial magnetorheological fluid (MRF-140CG; LORD Co. (USA)) with the same magnetic solid volume fraction, using the Physica MCR-300 rheometer with a 20 mm diameter plate-plate magnetorheological cell (MRD180). D1 sample is a suspension of micrometer range Fe particles in a transformer oil based magnetic fluid as carrier. For both types of samples, the experimental data for zero and non-zero magnetic field conditions were fitted to equations derived from the Newtonian and Cross type flow equations, as well as the Herschel-Bulkley model. The main advantage of both rheological equations for the quantitative description of the magnetic field behaviour of samples is that they can be used in regular CFD codes to compute the flow properties of the magnetorheological fluid and of the bidisperse magnetizable fluid for practical applications.
Ultrasonic fluid flow measurement method and apparatus
Kronberg, J.W.
1993-10-12
An apparatus for measuring the flow of a fluid in a pipe using ultrasonic waves. The apparatus comprises an ultrasonic generator, a lens for focusing the sound energy produced by the generator, and means for directing the focused energy into the side of the pipe through an opening and in a direction close to parallel to the long axis of the pipe. A cone carries the sound energy to the lens from the generator. Depending on the choice of materials, there may be a quarter-wave, acoustic impedance matching section between the generator and the cone to reduce the reflections of energy at the cone boundary. The lens material has an acoustic impedance similar to that of the cone material but a different sonic velocity so that the lens can converge the sound waves in the fluid. A transition section between the lens and the fluid helps to couple the energy to the fluid and assures it is directed as close to parallel to the fluid flow direction as possible. 3 figures.
Ultrasonic fluid flow measurement method and apparatus
Kronberg, James W.
1993-01-01
An apparatus for measuring the flow of a fluid in a pipe using ultrasonic waves. The apparatus comprises an ultrasonic generator, a lens for focusing the sound energy produced by the generator, and means for directing the focused energy into the side of the pipe through an opening and in a direction close to parallel to the long axis of the pipe. A cone carries the sound energy to the lens from the generator. Depending on the choice of materials, there may be a quarter-wave, acoustic impedance matching section between the generator and the cone to reduce the reflections of energy at the cone boundary. The lens material has an acoustic impedance similar to that of the cone material but a different sonic velocity so that the lens can converge the sound waves in the fluid. A transition section between the lens and the fluid helps to couple the energy to the fluid and assures it is directed as close to parallel to the fluid flow direction as possible.
2-Phase Fluid Flow & Heat Transport
Energy Science and Technology Software Center (ESTSC)
1993-03-13
GEOTHER is a three-dimensional, geothermal reservoir simulation code. The model describes heat transport and flow of a single component, two-phase fluid in porous media. It is based on the continuity equations for steam and water, which are reduced to two nonlinear partial differential equations in which the dependent variables are fluid pressure and enthalpy. GEOTHER can be used to simulate the fluid-thermal interaction in rock that can be approximated by a porous media representation. Itmore » can simulate heat transport and the flow of compressed water, two-phase mixtures, and superheated steam in porous media over a temperature range of 10 to 300 degrees C. In addition, it can treat the conversion from single to two-phase flow, and vice versa. It can be used for evaluation of a near repository spatial scale and a time scale of a few years to thousands of years. The model can be used to investigate temperature and fluid pressure changes in response to thermal loading by waste materials.« less
NASA Astrophysics Data System (ADS)
Lopes Filho, Milton C.; Nussenzveig Lopes, Helena J.; Titi, Edriss S.; Zang, Aibin
2015-06-01
The second-grade fluid equations are a model for viscoelastic fluids, with two parameters: α > 0, corresponding to the elastic response, and , corresponding to viscosity. Formally setting these parameters to 0 reduces the equations to the incompressible Euler equations of ideal fluid flow. In this article we study the limits of solutions of the second-grade fluid system, in a smooth, bounded, two-dimensional domain with no-slip boundary conditions. This class of problems interpolates between the Euler- α model (), for which the authors recently proved convergence to the solution of the incompressible Euler equations, and the Navier-Stokes case ( α = 0), for which the vanishing viscosity limit is an important open problem. We prove three results. First, we establish convergence of the solutions of the second-grade model to those of the Euler equations provided , as α → 0, extending the main result in (Lopes Filho et al., Physica D 292(293):51-61, 2015). Second, we prove equivalence between convergence (of the second-grade fluid equations to the Euler equations) and vanishing of the energy dissipation in a suitably thin region near the boundary, in the asymptotic regime , as α → 0. This amounts to a convergence criterion similar to the well-known Kato criterion for the vanishing viscosity limit of the Navier-Stokes equations to the Euler equations. Finally, we obtain an extension of Kato's classical criterion to the second-grade fluid model, valid if , as . The proof of all these results relies on energy estimates and boundary correctors, following the original idea by Kato.
Monotectic composite growth with fluid flow
NASA Astrophysics Data System (ADS)
Stöcker, C.; Ratke, L.
It is a well-known fluid-mechanical phenomenon that thermocapillary forces induce surface convection on a fluid-fluid interface. This so-called Marangoni convection depends on the variation of the surface energy along the interface. In our present work we focus our attention on the evolution of a fibrous monotectic microstructure with liquid L 2 fibers. We will show, that the Marangoni convection has a strong influence on the transport of solute in front of the solidification front, despite the flow induced by density differences. The resulting flow field affects the constitutional undercooling and therefore the mean undercooling of a monotectic solidification front. In a previous paper we discussed qualitatively the influence of fluid flow on the microstructure evolution of composite monotetic growth (C. Stöcker, L. Ratke, J. Crystal Growth 203 (1999) 582). We introduced an analytical model that takes the density differences of the phases and the surface convection on the L 1-L 2 surface into consideration. With this extended Jackson and Hunt theory for composite monotectic growth we derived a characteristic equation for the inter-rod distance depending on solidification velocity and temperature gradient. In this paper we develop a more accurate model. We solve numerically the diffusion equation coupled with the Navier-Stokes equation in the L 1 phase to find the minimal undercooling for a given velocity and temperature gradient. We derive a Jackson and Hunt diagram and show that the fluid flow leads to a strong dependence of the inter-rod distance on the temperature gradient opposite to eutectic solidification.
Modeling Tools Predict Flow in Fluid Dynamics
NASA Technical Reports Server (NTRS)
2010-01-01
"Because rocket engines operate under extreme temperature and pressure, they present a unique challenge to designers who must test and simulate the technology. To this end, CRAFT Tech Inc., of Pipersville, Pennsylvania, won Small Business Innovation Research (SBIR) contracts from Marshall Space Flight Center to develop software to simulate cryogenic fluid flows and related phenomena. CRAFT Tech enhanced its CRUNCH CFD (computational fluid dynamics) software to simulate phenomena in various liquid propulsion components and systems. Today, both government and industry clients in the aerospace, utilities, and petrochemical industries use the software for analyzing existing systems as well as designing new ones."
Slip mechanisms in complex fluid flows.
Hatzikiriakos, Savvas G
2015-10-28
The classical no-slip boundary condition of fluid mechanics is not always a valid assumption for the flow of several classes of complex fluids including polymer melts, their blends, polymer solutions, microgels, glasses, suspensions and pastes. In fact, it appears that slip effect in these systems is the rule and not the exemption. The occurrence of slip complicates the analysis of rheological data, although it provides new opportunities to understand their behavior in restricted environments delineating additional molecular mechanisms i.e. entropic restrictions due to limitations in the number of molecular conformations. This article discusses these complexities and provides future research opportunities. PMID:26345121
SIMULATIONS OF 2D AND 3D THERMOCAPILLARY FLOWS BY A LEAST-SQUARES FINITE ELEMENT METHOD. (R825200)
Numerical results for time-dependent 2D and 3D thermocapillary flows are presented in this work. The numerical algorithm is based on the Crank-Nicolson scheme for time integration, Newton's method for linearization, and a least-squares finite element method, together with a matri...
2-D transmitral flows simulation by means of the immersed boundary method on unstructured grids
NASA Astrophysics Data System (ADS)
Denaro, F. M.; Sarghini, F.
2002-04-01
Interaction between computational fluid dynamics and clinical researches recently allowed a deeper understanding of the physiology of complex phenomena involving cardio-vascular mechanisms. The aim of this paper is to develop a simplified numerical model based on the Immersed Boundary Method and to perform numerical simulations in order to study the cardiac diastolic phase during which the left ventricle is filled with blood flowing from the atrium throughout the mitral valve. As one of the diagnostic problems to be faced by clinicians is the lack of a univocal definition of the diastolic performance from the velocity measurements obtained by Eco-Doppler techniques, numerical simulations are supposed to provide an insight both into the physics of the diastole and into the interpretation of experimental data. An innovative application of the Immersed Boundary Method on unstructured grids is presented, fulfilling accuracy requirements related to the development of a thin boundary layer along the moving immersed boundary. It appears that this coupling between unstructured meshes and the Immersed Boundary Method is a promising technique when a wide range of spatial scales is involved together with a moving boundary. Numerical simulations are performed in a range of physiological parameters and a qualitative comparison with experimental data is presented, in order to demonstrate that, despite the simplified model, the main physiological characteristics of the diastole are well represented. Copyright
Gao, Hang; Bijnens, Nathalie; Coisne, Damien; Lugiez, Mathieu; Rutten, Marcel; D'hooge, Jan
2015-01-01
Despite the availability of multiple ultrasound approaches to left ventricular (LV) flow characterization in two dimensions, this technique remains in its childhood and further developments seem warranted. This article describes a new methodology for tracking the 2-D LV flow field based on ultrasound data. Hereto, a standard speckle tracking algorithm was modified by using a dynamic kernel embedding Navier-Stokes-based regularization in an iterative manner. The performance of the proposed approach was first quantified in synthetic ultrasound data based on a computational fluid dynamics model of LV flow. Next, an experimental flow phantom setup mimicking the normal human heart was used for experimental validation by employing simultaneous optical particle image velocimetry as a standard reference technique. Finally, the applicability of the approach was tested in a clinical setting. On the basis of the simulated data, pointwise evaluation of the estimated velocity vectors correlated well (mean r = 0.84) with the computational fluid dynamics measurement. During the filling period of the left ventricle, the properties of the main vortex obtained from the proposed method were also measured, and their correlations with the reference measurement were also calculated (radius, r = 0.96; circulation, r = 0.85; weighted center, r = 0.81). In vitro results at 60 bpm during one cardiac cycle confirmed that the algorithm properly measures typical characteristics of the vortex (radius, r = 0.60; circulation, r = 0.81; weighted center, r = 0.92). Preliminary qualitative results on clinical data revealed physiologic flow fields. PMID:25438850
Two-fluid Hydrodynamic Model for Fluid-Flow Simulation in Fluid-Solids Systems
Energy Science and Technology Software Center (ESTSC)
1994-06-20
FLUFIX is a two-dimensional , transient, Eulerian, and finite-difference program, based on a two-fluid hydrodynamic model, for fluid flow simulation in fluid-solids systems. The software is written in a modular form using the Implicit Multi-Field (IMF) numerical technique. Quantities computed are the spatial distribution of solids loading, gas and solids velocities, pressure, and temperatures. Predicted are bubble formation, bed frequencies, and solids recirculation. Applications include bubbling and circulating atmospheric and pressurized fluidized bed reactors, combustors,more » gasifiers, and FCC (Fluid Catalytic Cracker) reactors.« less
Analysis of Fluid Flow over a Surface
NASA Technical Reports Server (NTRS)
McCloud, Peter L. (Inventor)
2013-01-01
A method, apparatus, and computer program product for modeling heat radiated by a structure. The flow of a fluid over a surface of a model of the structure is simulated. The surface has a plurality of surface elements. Heat radiated by the plurality of surface elements in response to the fluid flowing over the surface of the model of the structure is identified. An effect of heat radiated by at least a portion of the plurality of surface elements on each other is identified. A model of the heat radiated by the structure is created using the heat radiated by the plurality of surface elements and the effect of the heat radiated by at least a portion of the plurality of surface elements on each other.
Visualization of working fluid flow in gravity assisted heat pipe
NASA Astrophysics Data System (ADS)
Nemec, Patrik; Malcho, Milan
2015-05-01
Heat pipe is device working with phase changes of working fluid inside hermetically closed pipe at specific pressure. The phase changes of working fluid from fluid to vapor and vice versa help heat pipe to transport high heat flux. The article deal about construction and processes casing in heat pipe during operation. Experiment visualization of working fluid flow is performed with glass heat pipe filed with ethanol. The visualization of working fluid flow explains the phenomena as working fluid boiling, nucleation of bubbles, vapor flow, vapor condensation on the wall, vapor and condensate flow interaction, flow down condensate film thickness on the wall, occurred during the heat pipe operation.
Bernoulli theorem generalized to rheologically complex viscous fluid flow
NASA Astrophysics Data System (ADS)
Brutyan, M. A.; Krapivskii, P. L.
1992-08-01
The Bernoulli theorem is generalized to two-dimensional and axisymmetric micropolar incompressible fluid flows. It is shown that the approach developed is also applicable to magnetohydrodynamic flows of a viscous Newtonian fluid.
Numerical Simulations of High-Frequency Respiratory Flows in 2D and 3D Lung Bifurcation Models
NASA Astrophysics Data System (ADS)
Chen, Zixi; Parameswaran, Shamini; Hu, Yingying; He, Zhaoming; Raj, Rishi; Parameswaran, Siva
2014-07-01
To better understand the human pulmonary system and optimize the high-frequency oscillatory ventilation (HFOV) design, numerical simulations were conducted under normal breathing frequency and HFOV condition using a CFD code Ansys Fluent and its user-defined C programs. 2D and 3D double bifurcating lung models were created, and the geometry corresponds to fifth to seventh generations of airways with the dimensions based on the Weibel's pulmonary model. Computations were carried out for different Reynolds numbers (Re = 400 and 1000) and Womersley numbers (α = 4 and 16) to study the air flow fields, gas transportation, and wall shear stresses in the lung airways. Flow structure was compared with experimental results. Both 2D and 3D numerical models successfully reproduced many results observed in the experiment. The oxygen concentration distribution in the lung model was investigated to analyze the influence of flow oscillation on gas transport inside the lung model.
NASA Astrophysics Data System (ADS)
Stahr, Donald W.; Law, Richard D.
2014-11-01
We model the development of shape preferred orientation (SPO) of a large population of two- and three-dimensional (2D and 3D) rigid clasts suspended in a linear viscous matrix deformed by superposed steady and continuously non-steady plane strain flows to investigate the sensitivity of clasts to changing boundary conditions during a single or superposed deformation events. Resultant clast SPOs are compared to one developed by an identical initial population that experienced a steady flow history of constant kinematic vorticity and reached an identical finite strain state, allowing examination of SPO sensitivity to deformation path. Rotation paths of individual triaxial inclusions are complex, even for steady plane strain flow histories. It has been suggested that the 3D nature of the system renders predictions based on 2D models inadequate for applied clast-based kinematic vorticity gauges. We demonstrate that for a large population of clasts, simplification to a 2D model does provide a good approximation to the SPO predicted by full 3D analysis for steady and non-steady plane strain deformation paths. Predictions of shape fabric development from 2D models are not only qualitatively similar to the more complex 3D analysis, but they display the same limitations of techniques based on clast SPO commonly used as a quantitative kinematic vorticity gauge. Our model results from steady, superposed, and non-steady flow histories with a significant pure shearing component at a wide range of finite strain resemble predictions for an identical initial population that experienced a single steady simple shearing deformation. We conclude that individual 2D and 3D clasts respond instantaneously to changes in boundary conditions, however, in aggregate, the SPO of a population of rigid inclusions does not reflect the late-stage kinematics of deformation, nor is it an indicator of the unique 'mean' kinematic vorticity experienced by a deformed rock volume.
Visualization periodic flows in a continuously stratified fluid.
NASA Astrophysics Data System (ADS)
Bardakov, R.; Vasiliev, A.
2012-04-01
To visualize the flow pattern of viscous continuously stratified fluid both experimental and computational methods were developed. Computational procedures were based on exact solutions of set of the fundamental equations. Solutions of the problems of flows producing by periodically oscillating disk (linear and torsion oscillations) were visualized with a high resolutions to distinguish small-scale the singular components on the background of strong internal waves. Numerical algorithm of visualization allows to represent both the scalar and vector fields, such as velocity, density, pressure, vorticity, stream function. The size of the source, buoyancy and oscillation frequency, kinematic viscosity of the medium effects were traced in 2D an 3D posing problems. Precision schlieren instrument was used to visualize the flow pattern produced by linear and torsion oscillations of strip and disk in a continuously stratified fluid. Uniform stratification was created by the continuous displacement method. The buoyancy period ranged from 7.5 to 14 s. In the experiments disks with diameters from 9 to 30 cm and a thickness of 1 mm to 10 mm were used. Different schlieren methods that are conventional vertical slit - Foucault knife, vertical slit - filament (Maksoutov's method) and horizontal slit - horizontal grating (natural "rainbow" schlieren method) help to produce supplementing flow patterns. Both internal wave beams and fine flow components were visualized in vicinity and far from the source. Intensity of high gradient envelopes increased proportionally the amplitude of the source. In domains of envelopes convergence isolated small scale vortices and extended mushroom like jets were formed. Experiments have shown that in the case of torsion oscillations pattern of currents is more complicated than in case of forced linear oscillations. Comparison with known theoretical model shows that nonlinear interactions between the regular and singular flow components must be taken
Piezoelectric energy harvesting in internal fluid flow.
Lee, Hyeong Jae; Sherrit, Stewart; Tosi, Luis Phillipe; Walkemeyer, Phillip; Colonius, Tim
2015-01-01
We consider piezoelectric flow energy harvesting in an internal flow environment with the ultimate goal powering systems such as sensors in deep oil well applications. Fluid motion is coupled to structural vibration via a cantilever beam placed in a converging-diverging flow channel. Two designs were considered for the electromechanical coupling: first; the cantilever itself is a piezoelectric bimorph; second; the cantilever is mounted on a pair of flextensional actuators. We experimentally investigated varying the geometry of the flow passage and the flow rate. Experimental results revealed that the power generated from both designs was similar; producing as much as 20 mW at a flow rate of 20 L/min. The bimorph designs were prone to failure at the extremes of flow rates tested. Finite element analysis (FEA) showed fatigue failure was imminent due to stress concentrations near the bimorph's clamped region; and that robustness could be improved with a stepped-joint mounting design. A similar FEA model showed the flextensional-based harvester had a resonant frequency of around 375 Hz and an electromechanical coupling of 0.23 between the cantilever and flextensional actuators in a vacuum. These values; along with the power levels demonstrated; are significant steps toward building a system design that can eventually deliver power in the Watts range to devices down within a well. PMID:26473879
Piezoelectric Energy Harvesting in Internal Fluid Flow
Lee, Hyeong Jae; Sherrit, Stewart; Tosi, Luis Phillipe; Walkemeyer, Phillip; Colonius, Tim
2015-01-01
We consider piezoelectric flow energy harvesting in an internal flow environment with the ultimate goal powering systems such as sensors in deep oil well applications. Fluid motion is coupled to structural vibration via a cantilever beam placed in a converging-diverging flow channel. Two designs were considered for the electromechanical coupling: first; the cantilever itself is a piezoelectric bimorph; second; the cantilever is mounted on a pair of flextensional actuators. We experimentally investigated varying the geometry of the flow passage and the flow rate. Experimental results revealed that the power generated from both designs was similar; producing as much as 20 mW at a flow rate of 20 L/min. The bimorph designs were prone to failure at the extremes of flow rates tested. Finite element analysis (FEA) showed fatigue failure was imminent due to stress concentrations near the bimorph’s clamped region; and that robustness could be improved with a stepped-joint mounting design. A similar FEA model showed the flextensional-based harvester had a resonant frequency of around 375 Hz and an electromechanical coupling of 0.23 between the cantilever and flextensional actuators in a vacuum. These values; along with the power levels demonstrated; are significant steps toward building a system design that can eventually deliver power in the Watts range to devices down within a well. PMID:26473879
NASA Astrophysics Data System (ADS)
Kawamura, E.; Lichtenberg, A. J.; Lieberman, M. A.; Marakhtanov, A. M.
2016-06-01
A fast 2D axisymmetric fluid-analytical multifrequency capacitively coupled plasma (CCP) reactor code is used to study center high nonuniformity in a low pressure electronegative chlorine discharge. In the code, a time-independent Helmholtz wave equation is used to solve for the capacitive fields in the linearized frequency domain. This eliminates the time dependence from the electromagnetic (EM) solve, greatly speeding up the simulations at the cost of neglecting higher harmonics. However, since the code allows up to three driving frequencies, we can add the two most important harmonics to the CCP simulations as the second and third input frequencies. The amplitude and phase of these harmonics are estimated by using a recently developed 1D radial nonlinear transmission line (TL) model of a highly asymmetric cylindrical discharge (Lieberman et al 2015 Plasma Sources Sci. Technol. 24 055011). We find that at higher applied frequencies, the higher harmonics contribute significantly to the center high nonuniformity due to their shorter plasma wavelengths.
Physical aspects of computing the flow of a viscous fluid
NASA Technical Reports Server (NTRS)
Mehta, U. B.
1984-01-01
One of the main themes in fluid dynamics at present and in the future is going to be computational fluid dynamics with the primary focus on the determination of drag, flow separation, vortex flows, and unsteady flows. A computation of the flow of a viscous fluid requires an understanding and consideration of the physical aspects of the flow. This is done by identifying the flow regimes and the scales of fluid motion, and the sources of vorticity. Discussions of flow regimes deal with conditions of incompressibility, transitional and turbulent flows, Navier-Stokes and non-Navier-Stokes regimes, shock waves, and strain fields. Discussions of the scales of fluid motion consider transitional and turbulent flows, thin- and slender-shear layers, triple- and four-deck regions, viscous-inviscid interactions, shock waves, strain rates, and temporal scales. In addition, the significance and generation of vorticity are discussed. These physical aspects mainly guide computations of the flow of a viscous fluid.
NASA Astrophysics Data System (ADS)
Shen, D.; Wang, J.; Cheng, X.; Rui, Y.; Ye, S.
2015-02-01
The rapid progress of Light Detection And Ranging (LiDAR) technology has made acquirement and application of high-resolution digital elevation model (DEM) data increasingly popular, especially with regards to the study of floodplain flow modeling. High-resolution DEM data include many redundant interpolation points, needs a high amount of calculation, and does not match the size of computational mesh. These disadvantages are a common problem for floodplain flow modeling studies. Two-dimensional (2-D) hydraulic modeling, a popular method of analyzing floodplain flow, offers high precision of elevation parameterization for computational mesh while ignoring much micro-topographic information of the DEM data itself. We offer a flood simulation method that integrates 2-D hydraulic model results and high-resolution DEM data, enabling the calculation of flood water levels in DEM grid cells through local inverse distance weighted interpolation. To get rid of the false inundation areas during interpolation, it employs the run-length encoding method to mark the inundated DEM grid cells and determine the real inundation areas through the run-length boundary tracing technique, which solves the complicated problem of the connectivity between DEM grid cells. We constructed a 2-D hydraulic model for the Gongshuangcha polder, a flood storage area of Dongting Lake, using our integrated method to simulate the floodplain flow. The results demonstrate that this method can solve DEM associated problems efficiently and simulate flooding processes with greater accuracy than DEM only simulations.
A central difference type approximation of convection in fluid flow
NASA Astrophysics Data System (ADS)
Tzanos, Constantine P.
In fluid flow and heat transfer problems, to eliminate instability and oscillations resulting from the numerical approximation of the convection term, it is common practice to use first order upwind differencing. For example, first order upwind differencing is used in codes like COMMIX and PHOENICS. However, in problems characterized by high Peclet numbers, this approach is highly inaccurate unless an excessive number of nodes is used. To eliminate this difficulty, a method of grid adaptation using central differencing in one-dimensional problems was presented. In the effort to extend this approach to 2-D and 3-D problems a new method was developed for multidimensional problems that is superior to first order upwind differencing even for a few non-adaptive grid points. This method is presented in this work.
Fluid flow in solidifying monotectic alloys
NASA Technical Reports Server (NTRS)
Ecker, A.; Frazier, D. O.; Alexander, J. Iwan D.
1989-01-01
Use of a two-wavelength holographic technique results in a simultaneous determination of temperature and composition profiles during directional solidification in a system with a miscibility gap. The relationships among fluid flow, phase separation, and mass transport during the solidification of the monotectic alloy are discussed. The primary sources of fluid motion in this system are buoyancy and thermocapillary forces. These forces act together when phase separation results in the formation of droplets (this occurs at the solid-liquid interface and in the bulk melt). In the absence of phase separation, buoyancy results from density gradients related to temperature and compositional gradients in the single-phase bulk melt. The effects of buoyancy are especially evident in association with water- or ethanol-rich volumes created at the solid-liquid growth interface.
Testing the Markov hypothesis in fluid flows
NASA Astrophysics Data System (ADS)
Meyer, Daniel W.; Saggini, Frédéric
2016-05-01
Stochastic Markov processes are used very frequently to model, for example, processes in turbulence and subsurface flow and transport. Based on the weak Chapman-Kolmogorov equation and the strong Markov condition, we present methods to test the Markov hypothesis that is at the heart of these models. We demonstrate the capabilities of our methodology by testing the Markov hypothesis for fluid and inertial particles in turbulence, and fluid particles in the heterogeneous subsurface. In the context of subsurface macrodispersion, we find that depending on the heterogeneity level, Markov models work well above a certain scale of interest for media with different log-conductivity correlation structures. Moreover, we find surprising similarities in the velocity dynamics of the different media considered.
NASA Technical Reports Server (NTRS)
Thompson David S.; Soni, Bharat K.
2001-01-01
An integrated geometry/grid/simulation software package, ICEG2D, is being developed to automate computational fluid dynamics (CFD) simulations for single- and multi-element airfoils with ice accretions. The current version, ICEG213 (v2.0), was designed to automatically perform four primary functions: (1) generate a grid-ready surface definition based on the geometrical characteristics of the iced airfoil surface, (2) generate high-quality structured and generalized grids starting from a defined surface definition, (3) generate the input and restart files needed to run the structured grid CFD solver NPARC or the generalized grid CFD solver HYBFL2D, and (4) using the flow solutions, generate solution-adaptive grids. ICEG2D (v2.0) can be operated in either a batch mode using a script file or in an interactive mode by entering directives from a command line within a Unix shell. This report summarizes activities completed in the first two years of a three-year research and development program to address automation issues related to CFD simulations for airfoils with ice accretions. As well as describing the technology employed in the software, this document serves as a users manual providing installation and operating instructions. An evaluation of the software is also presented.
NASA Astrophysics Data System (ADS)
Li, Y.; Kazemifar, F.; Blois, G.; Christensen, K. T.
2015-12-01
Multiphase flow of water and supercritical carbon dioxide (CO2) in porous media is central to geological sequestration of CO2 into saline aquifers. However, our fundamental understanding of the coupled flow dynamics of CO2 and water in complex geologic media still remains limited, especially at the pore scale. Recently, studies have been carried out in 2D homogeneous models with the micro-PIV technique, yielding very interesting observations of pore-scale flow transport. The primary aim of this work is to leverage this experimental protocol to quantify the pore-scale flow of water and liquid/supercritical CO2 in 2D heterogeneous porous micromodels under reservoir-relevant conditions. The goal is to capture the dynamics of this multi-phase flow in a porous matrix that mimics the heterogeneity of natural rock. Fluorescent microscopy and the micro-PIV technique are employed to simultaneously measure the spatially-resolved instantaneous velocity field in the water and quantify the instantaneous spatial configuration of both phases. The results for heterogeneous micromodels will be presented and compared with those for homogeneous micromodels, yielding valuable insight into flow processes at the pore scale in natural rock.
Adaptive clutter filter in 2-D color flow imaging based on in vivo I/Q signal.
Zhou, Xiaoming; Zhang, Congyao; Liu, Dong C
2014-01-01
Color flow imaging has been well applied in clinical diagnosis. For the high quality color flow images, clutter filter is important to separate the Doppler signals from blood and tissue. Traditional clutter filters, such as finite impulse response, infinite impulse response and regression filters, were applied, which are based on the hypothesis that the clutter signal is stationary or tissue moves slowly. However, in realistic clinic color flow imaging, the signals are non-stationary signals because of accelerated moving tissue. For most related papers, simulated RF signals are widely used without in vivo I/Q signal. Hence, in this paper, adaptive polynomial regression filter, which is down mixing with instantaneous clutter frequency, was proposed based on in vivo carotid I/Q signal in realistic color flow imaging. To get the best performance, the optimal polynomial order of polynomial regression filter and the optimal polynomial order for estimation of instantaneous clutter frequency respectively were confirmed. Finally, compared with the mean blood velocity and quality of 2-D color flow image, the experiment results show that adaptive polynomial regression filter, which is down mixing with instantaneous clutter frequency, can significantly enhance the mean blood velocity and get high quality 2-D color flow image. PMID:24211911
Fluid flow in solidifying monotectic alloys
NASA Astrophysics Data System (ADS)
Ecker, A.; Frazier, D. O.; Alexander, J. Iwan D.
1989-11-01
Use of a two-wavelength holographic technique results in a simultaneous determination of temperature and composition profiles during directional solidification in a system with a miscibility gap. A shadowgraph technique is employed for flow visualization. By these methods, flow regimes are identified and related to particular melt compositions. We discuss the relationships among fluid flow, phase separation, and mass transport during the solidification of the monotectic alloy. The primary sources of fluid motion in this system are buoyancy and thermocapillary forces. These forces act together when phase separation results in the formation of droplets (this occurs at the solid-liquid interface and in the bulk melt). While buoyancy forces arise due to density differences between the droplet and the host phase, thermocapillary forces (associated with temperature gradients in the droplet surface) may predominate. In the absence of phase separation, buoyancy results from density gradients related to temperature and compositional gradients in the single-phase bulk melt. The effects of buoyancy are especially evident in association with water- or ethanol-rich volumes created at the solid-liquid growth interface.
Fluid flow through packings of rotating obstacles
NASA Astrophysics Data System (ADS)
Oliveira, Rafael S.; Andrade, José S.; Andrade, Roberto F. S.
2015-03-01
We investigate through numerical simulation the nonstationary flow of a Newtonian fluid through a two-dimensional channel filled with an array of circular obstacles of distinct sizes. The disks may rotate around their respective centers, modeling a nonstationary, inhomogeneous porous medium. Obstacle sizes and positions are defined by the geometry of an Apollonian packing (AP). To allow for fluid flow, the radii of the disks are uniformly reduced by a factor 0.6 ≤s ≤0.8 for assemblies corresponding to the four first AP generations. The investigation is targeted to elucidate the main features of the rotating regime as compared to the fixed disk condition. It comprises the evaluation of the region of validity of Darcy's law as well as the study of the nonlinear hydraulic resistance as a function of the channel Reynolds number, the reduction factor s , and the AP generation. Depending on a combination of these factors, the resistance of rotating disks may be larger or smaller than that of the corresponding static case. We also analyze the flow redistribution in the interdisk channels as a result of the rotation pattern and characterize the angular velocity of the disks. Here, the striking feature is the emergence of a stable oscillatory behavior of the angular velocity for almost all disks that are inserted into the assemblies after the second generation.
Fluid flow dynamics in MAS systems
NASA Astrophysics Data System (ADS)
Wilhelm, Dirk; Purea, Armin; Engelke, Frank
2015-08-01
The turbine system and the radial bearing of a high performance magic angle spinning (MAS) probe with 1.3 mm-rotor diameter has been analyzed for spinning rates up to 67 kHz. We focused mainly on the fluid flow properties of the MAS system. Therefore, computational fluid dynamics (CFD) simulations and fluid measurements of the turbine and the radial bearings have been performed. CFD simulation and measurement results of the 1.3 mm-MAS rotor system show relatively low efficiency (about 25%) compared to standard turbo machines outside the realm of MAS. However, in particular, MAS turbines are mainly optimized for speed and stability instead of efficiency. We have compared MAS systems for rotor diameter of 1.3-7 mm converted to dimensionless values with classical turbomachinery systems showing that the operation parameters (rotor diameter, inlet mass flow, spinning rate) are in the favorable range. This dimensionless analysis also supports radial turbines for low speed MAS probes and diagonal turbines for high speed MAS probes. Consequently, a change from Pelton type MAS turbines to diagonal turbines might be worth considering for high speed applications. CFD simulations of the radial bearings have been compared with basic theoretical values proposing considerably smaller frictional loss values. The discrepancies might be due to the simple linear flow profile employed for the theoretical model. Frictional losses generated inside the radial bearings result in undesired heat-up of the rotor. The rotor surface temperature distribution computed by CFD simulations show a large temperature gradient over the rotor.
Fluid flow dynamics in MAS systems.
Wilhelm, Dirk; Purea, Armin; Engelke, Frank
2015-08-01
The turbine system and the radial bearing of a high performance magic angle spinning (MAS) probe with 1.3mm-rotor diameter has been analyzed for spinning rates up to 67kHz. We focused mainly on the fluid flow properties of the MAS system. Therefore, computational fluid dynamics (CFD) simulations and fluid measurements of the turbine and the radial bearings have been performed. CFD simulation and measurement results of the 1.3mm-MAS rotor system show relatively low efficiency (about 25%) compared to standard turbo machines outside the realm of MAS. However, in particular, MAS turbines are mainly optimized for speed and stability instead of efficiency. We have compared MAS systems for rotor diameter of 1.3-7mm converted to dimensionless values with classical turbomachinery systems showing that the operation parameters (rotor diameter, inlet mass flow, spinning rate) are in the favorable range. This dimensionless analysis also supports radial turbines for low speed MAS probes and diagonal turbines for high speed MAS probes. Consequently, a change from Pelton type MAS turbines to diagonal turbines might be worth considering for high speed applications. CFD simulations of the radial bearings have been compared with basic theoretical values proposing considerably smaller frictional loss values. The discrepancies might be due to the simple linear flow profile employed for the theoretical model. Frictional losses generated inside the radial bearings result in undesired heat-up of the rotor. The rotor surface temperature distribution computed by CFD simulations show a large temperature gradient over the rotor. PMID:26073599
A 2-D oscillating flow analysis in Stirling engine heat exchangers
NASA Technical Reports Server (NTRS)
Ahn, Kyung H.; Ibrahim, Mounir B.
1991-01-01
A two-dimensional oscillating flow analysis was conducted, simulating the gas flow inside Stirling heat exchangers. Both laminar and turbulent oscillating pipe flow were investigated numerically for Re(max) = 1920 (Va = 80), 10800 (Va = 272), 19300 (Va = 272), and 60800 (Va = 126). The results are compared with experimental results of previous investigators. Also, predictions of the flow regime on present oscillating flow conditions were checked by comparing velocity amplitudes and phase differences with those from laminar theory and quasi-steady profile. A high Reynolds number k-epsilon turbulence model was used for turbulent oscillating pipe flow. Finally, performance evaluation of the K-epsilon model was made to explore the applicability of quasi-steady turbulent models to unsteady oscillating flow analysis.
A 2-D oscillating flow analysis in Stirling engine heat exchangers
NASA Technical Reports Server (NTRS)
Ahn, Kyung H.; Ibrahim, Mounir B.
1991-01-01
A two dimensional oscillating flow analysis was conducted, simulating the gas flow inside Stirling heat exchangers. Both laminar and turbulent oscillating pipe flow were investigated numerically for Re(max) = 1920 (Va = 80), 10800 (Va = 272), 19300 (Va = 272), and 60800 (Va = 126). The results are compared with experimental results of previous investigators. Also, predictions of the flow regime on present oscillating flow conditions were checked by comparing velocity amplitudes and phase differences with those from laminar theory and quasi-steady profile. A high Reynolds number k-epsilon turbulence model was used for turbulent oscillating pipe flow. Finally, performance evaluation of the K-epsilon model was made to explore the applicability of quasi-steady turbulent models to unsteady oscillating flow analysis.
Flow Diode and Method for Controlling Fluid Flow Origin of the Invention
NASA Technical Reports Server (NTRS)
Dyson, Rodger W (Inventor)
2015-01-01
A flow diode configured to permit fluid flow in a first direction while preventing fluid flow in a second direction opposite the first direction is disclosed. The flow diode prevents fluid flow without use of mechanical closures or moving parts. The flow diode utilizes a bypass flowline whereby all fluid flow in the second direction moves into the bypass flowline having a plurality of tortuous portions providing high fluidic resistance. The portions decrease in diameter such that debris in the fluid is trapped. As fluid only travels in one direction through the portions, the debris remains trapped in the portions.
Flow rate measurement in aggressive conductive fluids
NASA Astrophysics Data System (ADS)
Dubovikova, Nataliia; Kolesnikov, Yuri; Karcher, Christian
2014-03-01
Two non-contact experimental methods of flow rate measurements for aggressive conductive liquids are described. The techniques are based on electromagnetic forces and Faraday's law: Lorentz force is induced inside moving conductive liquid under influence of variable magnetic field of permanent magnets. They are mounted along a liquid metal channel or (in case of the second method) inserted into rotated metal wheels. The force acts in the opposite of fluids' velocity direction and hence it is possible to measure reaction force of it that takes place according to Newton's law on magnetic field source - permanent magnets. And by knowing the force, which linearly depends on velocity, one can calculate mean flow rate of liquid. In addition experimental "dry" calibration and its results are described for one of the measurements' techniques.
Geophysical Fluid Flow Cell (GFFC) Cross Section
NASA Technical Reports Server (NTRS)
1995-01-01
This drawing shows a cross-section view of the test cell at the heart of the Geophysical Fluid Flow Cell (GFFC) that flew on two Spacelab missions. The middle and lower drawings depict the volume of the silicone oil layer that served as the atmosphere as the steel ball rotated and an electrostatic field pulled the oil inward to mimic gravity's effects during the experiments. The GFFC thus produced flow patterns that simulated conditions inside the atmospheres of Jupiter and the Sun and other stars. The principal investigator was John Hart of the University of Colorado at Boulder. It was managed by NASA's Marshall Space Flight Center (MSFC). An Acrobat PDF copy of this drawing is available at http://microgravity.nasa.gov/gallery. (Credit: NASA/Marshall Space Flight Center)
Bubble Formation in Yield Stress Fluids Using Flow-Focusing and T -Junction Devices
NASA Astrophysics Data System (ADS)
Laborie, Benoit; Rouyer, Florence; Angelescu, Dan E.; Lorenceau, Elise
2015-05-01
We study the production of bubbles inside yield stress fluids (YSFs) in axisymmetric T -junction and flow-focusing devices. Taking advantage of yield stress over capillary stress, we exhibit a robust break-up mechanism reminiscent of the geometrical operating regime in 2D flow-focusing devices for Newtonian fluids. We report that when the gas is pressure driven, the dynamics is unsteady due to hydrodynamic feedback and YSF deposition on the walls of the channels. However, the present study also identifies pathways for potential steady-state production of bubbly YSFs at large scale.
Microscale imaging of cilia-driven fluid flow
Huang, Brendan K.; Choma, Michael A.
2015-01-01
Cilia-driven fluid flow is important for multiple processes in the body, including respiratory mucus clearance, gamete transport in the oviduct, right-left patterning in the embryonic node, and cerebrospinal fluid circulation. Multiple imaging techniques have been applied towards quantifying ciliary flow. Here we review common velocimetry methods of quantifying fluid flow. We then discuss four important optical modalities, including light microscopy, epifluorescence, confocal microscopy, and optical coherence tomography, that have been used to investigate cilia-driven flow. PMID:25417211
Unified slip boundary condition for fluid flows.
Thalakkottor, Joseph John; Mohseni, Kamran
2016-08-01
Determining the correct matching boundary condition is fundamental to our understanding of several everyday problems. Despite over a century of scientific work, existing velocity boundary conditions are unable to consistently explain and capture the complete physics associated with certain common but complex problems, such as moving contact lines and corner flows. The widely used Maxwell and Navier slip boundary conditions make an implicit assumption that velocity varies only in the wall normal direction. This makes their boundary condition inapplicable in the vicinity of contact lines and corner points, where velocity gradient exists both in the wall normal and wall tangential directions. In this paper, by identifying this implicit assumption we are able to extend Maxwell's slip model. Here, we present a generalized velocity boundary condition that shows that slip velocity is a function of not only the shear rate but also the linear strain rate. In addition, we present a universal relation for slip length, which shows that, for a general flow, slip length is a function of the principal strain rate. The universal relation for slip length along with the generalized velocity boundary condition provides a unified slip boundary condition to model a wide range of steady Newtonian fluid flows. We validate the unified slip boundary for simple Newtonian liquids by using molecular dynamics simulations and studying both the moving contact line and corner flow problems. PMID:27627398
VFLOW2D - A Vorte-Based Code for Computing Flow Over Elastically Supported Tubes and Tube Arrays
WOLFE,WALTER P.; STRICKLAND,JAMES H.; HOMICZ,GREGORY F.; GOSSLER,ALBERT A.
2000-10-11
A numerical flow model is developed to simulate two-dimensional fluid flow past immersed, elastically supported tube arrays. This work is motivated by the objective of predicting forces and motion associated with both deep-water drilling and production risers in the oil industry. This work has other engineering applications including simulation of flow past tubular heat exchangers or submarine-towed sensor arrays and the flow about parachute ribbons. In the present work, a vortex method is used for solving the unsteady flow field. This method demonstrates inherent advantages over more conventional grid-based computational fluid dynamics. The vortex method is non-iterative, does not require artificial viscosity for stability, displays minimal numerical diffusion, can easily treat moving boundaries, and allows a greatly reduced computational domain since vorticity occupies only a small fraction of the fluid volume. A gridless approach is used in the flow sufficiently distant from surfaces. A Lagrangian remap scheme is used near surfaces to calculate diffusion and convection of vorticity. A fast multipole technique is utilized for efficient calculation of velocity from the vorticity field. The ability of the method to correctly predict lift and drag forces on simple stationary geometries over a broad range of Reynolds numbers is presented.
Fluid Flow and Solidification Under Combined Action of Magnetic Fields and Microgravity
NASA Technical Reports Server (NTRS)
Li, B. Q.; Shu, Y.; Li, K.; deGroh, H. C.
2002-01-01
Mathematical models, both 2-D and 3-D, are developed to represent g-jitter induced fluid flows and their effects on solidification under combined action of magnetic fields and microgravity. The numerical model development is based on the finite element solution of governing equations describing the transient g-jitter driven fluid flows, heat transfer and solutal transport during crystal growth with and without an applied magnetic field in space vehicles. To validate the model predictions, a ground-based g-jitter simulator is developed using the oscillating wall temperatures where timely oscillating fluid flows are measured using a laser PIV system. The measurements are compared well with numerical results obtained from the numerical models. Results show that a combined action derived from magnetic damping and microgravity can be an effective means to control the melt flow and solutal transport in space single crystal growth systems.
A numerical method for computing unsteady 2-D boundary layer flows
NASA Technical Reports Server (NTRS)
Krainer, Andreas
1988-01-01
A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.
2D and 3D potential flows with rotational source terms in turbomachines
NASA Astrophysics Data System (ADS)
Alkalai, K.; Leboeuf, F.
A computational method capable of treating two- and three-dimensional potential flows is developed which includes blade effects and viscosity in source terms determined over the entire flowfield considered. Details of the mathematical and numerical formulations are given, and grid generation and density calculation are discussed. Preliminary results obtained with the codes developed here are then presented, and the possibility of applying the method to more complex flows is examined.
Study of Unsteady Flow Actuation Produced by Surface Plasma Actuator on 2-D Airfoil
NASA Astrophysics Data System (ADS)
Phan, Minh Khang; Shin, Jichul
2014-10-01
Effect of flow actuation driven by low current continuous or pulsed DC surface glow discharge plasma actuator is studied. Schlieren image of induced flow on flat plate taken at a high repetition rate reveals that the actuation is mostly initiated near the cathode. Assuming that the actuation is mostly achieved by ions in the cathode sheath region, numerical model for the source of flow actuation is obtained by analytical estimation of ion pressure force created in DC plasma sheath near the cathode and added in momentum equation as a body force term. Modeled plasma flow actuator is simulated with NACA0012 airfoil oscillating over a certain range of angle of attack (AoA) at specific reduced frequencies of airfoil. By changing actuation authority according to the change in AoA, stabilization of unsteady flow field is improved and hence steady aerodynamic performance can be maintained. Computational result shows that plasma actuation is only effective in modifying aerodynamic characteristics of separated flow. It turns out that plasma pulse frequency should be tuned for optimal performance depending on phase angle and rotating speed. The actuation authority can be parameterized by a ratio between plasma pulse frequency and reduced frequency.
Space Coffee Cup: Capillary Flow Driven Fluids in Space
Interested in learning more about how fluids react in Space? In this video, Professor Mark Weislogel, and Dr. Marshall Porterfield will discuss the Space Coffee Cup and Capillary Flow Driven Fluids...
Thermal and Fluid Flow Brazing Simulations
HOSKING, FLOYD MICHAEL; GIANOULAKIS,STEVEN E.; GIVLER,RICHARD C.; SCHUNK,P. RANDALL
1999-12-15
The thermal response of fixtured parts in a batch-type brazing furnace can require numerous, time-consuming development runs before an acceptable furnace schedule or joint design is established. Powerful computational simulation tools are being developed to minimize the required number of verification experiments, improve furnace throughput, and increase product yields. Typical furnace simulations are based on thermal, fluid flow, and structural codes that incorporate the fundamental physics of the brazing process. The use of massively parallel computing to predict furnace and joint-level responses is presented. Measured and computed data are compared. Temperature values are within 1-270 of the expected peak brazing temperature for different loading conditions. Sensitivity studies reveal that the thermal response is more sensitive to the thermal boundary conditions of the heating enclosure than variability y in the materials data. Braze flow simulations predict fillet geometry and free surface joint defects. Dynamic wetting conditions, interfacial reactions, and solidification structure add a high degree of uncertainty to the flow results.
An experimental study of flow separation over a flat plate with 2D transverse grooves
NASA Astrophysics Data System (ADS)
Jones, Emily Michelle
Nature has long been an inspiration for research in engineering. In particular, the biological surfaces of aquatic swimmers have been studied for their potential as drag reducing surfaces. The hydrodynamic benefit of riblets, or grooves embedded parallel to the flow, which appear on many aquatic biological surfaces, have been well documented and implemented in practical engineering applications. However the skin of dolphins is embedded with grooves that run perpendicular to the flow of water over their bodies. It is theorized that the transverse grooves present on dolphin skin trap vortices between them, creating a partial slip condition over the surface and inducing turbulence augmentation in the boundary layer, thus controlling boundary layer separation over the dolphin's skin. Similarly, sharks are covered with scales that are flexible at the base and capable of bristling, forming grooves running transverse to the flow. It is theorized that the scales bristle when encountering a reversing flow, thereby trapping vortices between the scales and, similarly, delaying boundary layer separation. In an attempt to test this hypothesis and study these affects, a spinning cylinder was used in a water tunnel to induce separation over a flat plate with 2 mm, rectangular transverse grooves and sinusoidal grooves of similar scaling. The results were compared to tripped, turbulent boundary layer separation occurring over a flat plate without grooves using time-resolved particle image velocimetry. The strength of the adverse pressure gradient was varied, and the observed delay in flow separation and other affects upon the boundary layer are discussed.
2D Flow patterning in Hele-Shaw configurations using Non-Uniform Electroosmotic Slip
NASA Astrophysics Data System (ADS)
Boyko, Evgeniy; Rubin, Shimon; Gat, Amir; Bercovici, Moran
2015-11-01
We present an analytical study, validated by numerical simulations, of electroosmotic flow in a Hele-Shaw configuration with non-uniform zeta potential distribution. Applying the lubrication approximation and assuming thin electric double layer, we derive a pair of uncoupled Poisson equations for the pressure and the stream function, and show that the inhomogeneous parts in these equations are governed by gradients in zeta potential parallel and perpendicular to the applied electric field, respectively. We obtain a solution for the case of a disk with uniform zeta potential and show that the flow field created is an exact dipole, even in the immediate vicinity of the disk. We then illustrate the ability to generate complex flow fields using superposition of such disks. Furthermore, we study the inverse problem in which we define the desired flow pattern and solve for the zeta potential distribution required in order to establish it. We demonstrate that such inverse problem solutions can be used to create directional flows confined within narrow regions, without physical walls. We show that these solutions can be assembled to create complex microfluidic networks, composed of intersecting channels and turns, which are basic building blocks in microfluidic devices.
NASA Astrophysics Data System (ADS)
Tong, Mingming; Browne, David J.
2012-01-01
Smoothed particle hydrodynamics is employed, for the first time, to develop a numerical model for the melting and fluid flow during laser welding process. In this meshlessLagrangian method the gas-melt two phase flow, heat transfer, surface tension, and melting of solid parent material are considered. This model was used to study the evolution of temperature field and fluid flow in the case study of laser spot welding in 2D. The simulation results show a strong influence of the melting process on the flow of liquid metal and a clear influence of the Marangoni flow on the heat transfer is also found.
NASA Astrophysics Data System (ADS)
Zhang, Wei; Markfort, Corey; Porté-Agel, Fernando
2014-05-01
Turbulent boundary-layer flows over complex topography have been extensively studied in the atmospheric sciences and wind engineering communities. The upwind turbulence level, the atmospheric thermal stability and the shape of the topography as well as surface characteristics play important roles in turbulent transport of momentum and scalar fluxes. However, to the best of our knowledge, atmospheric thermal stability has rarely been taken into account in laboratory simulations, particularly in wind-tunnel experiments. Extension of such studies in thermally-stratified wind tunnels will substantially advance our understanding of thermal stability effects on the physics of flow over complex topography. Additionally, high-resolution experimental data can be used for development of new parameterization of surface fluxes and validation of numerical models such as Large-Eddy Simulation (LES). A series of experiments of neutral and thermally-stratified boundary-layer flows over a wall-mounted 2-D block were conducted at the Saint Anthony Falls Laboratory boundary-layer wind tunnel. The 2-D block, with a width to height ratio of 2:1, occupied the lowest 25% of the turbulent boundary layer. Stable and convective boundary layers were simulated by independently controlling the temperature of air flow, the test section floor, and the wall-mounted block surfaces. Measurements using high-resolution Particle Image Velocimetry (PIV), x-wire/cold-wire anemometry, thermal-couples and surface heat flux sensors were made to quantify the turbulent properties and surface fluxes in distinct macroscopic flow regions, including the separation/recirculation zones, evolving shear layer and the asymptotic far wake. Emphasis will be put on addressing thermal stability effects on the spatial distribution of turbulent kinetic energy (TKE) and turbulent fluxes of momentum and scalar from the near to far wake region. Terms of the TKE budget equation are also inferred from measurements and
NASA Astrophysics Data System (ADS)
Zhang, Wei; Markfort, Corey; Porté-Agel, Fernando
2014-11-01
Turbulent flows over complex surface topography have been of great interest in the atmospheric science and wind engineering communities. The geometry of the topography, surface roughness and temperature characteristics as well as the atmospheric thermal stability play important roles in determining momentum and scalar flux distribution. Studies of turbulent flow over simplified topography models, under neutrally stratified boundary-layer conditions, have provided insights into fluid dynamics. However, atmospheric thermal stability has rarely been considered in laboratory experiments, e.g., wind-tunnel experiments. Series of wind-tunnel experiments of thermally-stratified boundary-layer flow over a surface-mounted 2-D block, in a well-controlled boundary-layer wind tunnel, will be presented. Measurements using high-resolution PIV, x-wire/cold-wire anemometry and surface heat flux sensors were conducted to quantify the turbulent flow properties, including the size of the recirculation zone, coherent vortex structures and the subsequent boundary layer recovery. Results will be shown to address thermal stability effects on momentum and scalar flux distribution in the wake, as well as dominant mechanism of turbulent kinetic energy generation and consumption. The authors gratefully acknowledge funding from the Swiss National Foundation (Grant 200021-132122), the National Science Foundation (Grant ATM-0854766) and NASA (Grant NNG06GE256).
On stability and turbulence of fluid flows
NASA Technical Reports Server (NTRS)
Heisenberg, Werner
1951-01-01
This investigation is divided into two parts, the treatment of the stability problem of fluid flows on the one hand, and that of the turbulent motion on the other. The first part summarizes all previous investigations under a unified point of view, that is, sets up as generally as possible the conditions under which a profile possesses unstable or stable characteristics, and indicates the methods for solution of the stability equation for any arbitrary velocity profile and for calculation of the critical Reynolds number for unstable profiles. In the second part, under certain greatly idealizing assumptions, differential equations for the turbulent motions are derived and from them qualitative information about several properties of the turbulent velocity distribution is obtained.
NASA Astrophysics Data System (ADS)
Shen, D.; Wang, J.; Cheng, X.; Rui, Y.; Ye, S.
2015-08-01
The rapid progress of lidar technology has made the acquirement and application of high-resolution digital elevation model (DEM) data increasingly popular, especially in regards to the study of floodplain flow. However, high-resolution DEM data pose several disadvantages for floodplain modeling studies; e.g., the data sets contain many redundant interpolation points, large numbers of calculations are required to work with data, and the data do not match the size of the computational mesh. Two-dimensional (2-D) hydraulic modeling, which is a popular method for analyzing floodplain flow, offers highly precise elevation parameterization for computational mesh while ignoring much of the micro-topographic information of the DEM data itself. We offer a flood simulation method that integrates 2-D hydraulic model results and high-resolution DEM data, thus enabling the calculation of flood water levels in DEM grid cells through local inverse distance-weighted interpolation. To get rid of the false inundation areas during interpolation, it employs the run-length encoding method to mark the inundated DEM grid cells and determine the real inundation areas through the run-length boundary tracing technique, which solves the complicated problem of connectivity between DEM grid cells. We constructed a 2-D hydraulic model for the Gongshuangcha detention basin, which is a flood storage area of Dongting Lake in China, by using our integrated method to simulate the floodplain flow. The results demonstrate that this method can solve DEM associated problems efficiently and simulate flooding processes with greater accuracy than simulations only with DEM.
MEAN FLOW AND TURBULENCE MEASUREMENTS AROUND A 2-D ARRAY OF BUILDINGS IN A WIND TUNNEL
In order to predict the dispersion of harmful materials released in or near an urban environment, it is important to first understand the complex flow patterns which result from the interaction of the wind with buildings and, more commonly, clusters of buildings. Recent advanc...
Turbulence Measurements on a 2D NACA 0036 with Synthetic Jet Flow Control
NASA Technical Reports Server (NTRS)
Wilson, J. S.
2006-01-01
An active flow control experiment was conducted on a 2-ft chord NACA 0036 airfoil in a 3-ft by 4-ft Wind Tunnel at Re = 1 x 10(exp 6). The model was equipped with synthetic jet actuators at x/c = 0.30 and 0.65 that provided 120 Hz periodic excitation at a C(sub mu) 0.86% through 0.06-in wide slots. Three different slot con gurations were tested, including a baseline with no slots. Surface pressure data was collected to compare to previous tests and to combine with turbulence data to aid future CFD modeling efforts. Turbulence data, measured by hot-wire, was compared with and without flow control. Pressure data corroborates previous test data and provides more points for CFD validation. Hot-wire results showed ow control reduced the separated wake size and brought the high Reynolds stress shear layer closer to the airfoil surface. The position of this layer to the surface was altered more significantly than the magnitude of the peak stresses. Flow control was shown to increase turbulent energy in the attached boundary layer downstream of the slot but to have little effect upstream. These results provide further justification to continue assessing the potential of active flow control to reduce drag of helicopter airframe components.
Saffer, D.M.; Bekins, B.A.
1998-01-01
Down-hole geochemical anomalies encountered in active accretionary systems can be used to constrain the timing, rates, and localization of fluid flow. Here we combine a coupled flow and solute transport model with a kinetic model for smectite dehydration to better understand and quantify fluid flow in the Nankai accretionary complex offshore of Japan. Compaction of sediments and clay dehydration provide fluid sources which drive the model flow system. We explicitly include the consolidation rate of underthrust sediments in our calculations to evaluate the impact that variations in this unknown quantity have on pressure and chloride distribution. Sensitivity analysis of steady state pressure solutions constrains bulk and flow conduit permeabilities. Steady state simulations with 30% smectite in the incoming sedimentary sequence result in minimum chloride concentrations at site 808 of 550 mM, but measured chlorinity is as low as 447 mM. We simulate the transient effects of hydrofracture or a strain event by assuming an instantaneous permeability increase of 3-4 orders of magnitude along a flow conduit (in this case the de??collement), using steady state results as initial conditions. Transient results with an increase in de??collement permeability from 10-16 m2 to 10-13 m2 and 20% smectite reproduce the observed chloride profile at site 808 after 80-160 kyr. Modeled chloride concentrations are highly sensitive to the consolidation rate of underthrust sediments, such that rapid compaction of underthrust material leads to increased freshening. Pressures within the de??collement during transient simulations rise rapidly to a significant fraction of lithostatic and remain high for at least 160 kyr, providing a mechanism for maintaining high permeability. Flow rates at the deformation front for transient simulations are in good agreement with direct measurements, but steady state flow rates are 2-3 orders of magnitude smaller than observed. Fluid budget calculations
NASA Astrophysics Data System (ADS)
Ague, J. J.
2004-12-01
Fluids are generally expected to be driven upward in the deep parts of orogens, but permeability heterogeneity and anisotropy must also be considered to properly interpret fluid infiltration and kinetic reaction histories preserved in the rock record. This paper focuses on new 2-D models of Darcian fluid flow incorporating permeability contrasts between rock units, the permeability tensor, and reactive fluid sources (e.g., dehydration). Factor of ten contrasts between the minimum and maximum permeability values in anisotropic rocks can strongly divert flow, but contrasts of as little as a factor of two still influence flow behavior. The first example considers fluid flow in subduction zone mélange, Syros, Greece. Geochemical evidence suggests that the interiors of meta-mafic blocks of oceanic crust in the mélange underwent limited fluid-rock reaction, despite extensive dehydration and decarbonation of the subduction complex. Modeling shows that if the blocks have lower permeability than the surrounding serpentine-rich matrix, then flow is diverted around the blocks resulting in little infiltration except at block margins, consistent with field relations. In this way, the subducted oceanic crust could preserve little evidence of fluid infiltration, even though considerable flow occurred through the mélange. The largest fluid fluxes are concentrated in matrix where blocks are in close proximity, and this effect increases as the anisotropy of the matrix increases. The lack of fluid infiltration into blocks could account for the observed limited metamorphism and strong kinetic overstepping of reactions that in some cases allowed preservation of ocean-floor mineral assemblages even at blueschist-eclogite facies conditions. The second example examines fluid flow through a folded sequence in which the direction of maximum permeability is parallel to the folded layering, and is based on field relations of Barrovian metamorphic sequences in CT, USA, and Scotland. As the
Fluid flow dynamics under location uncertainty
NASA Astrophysics Data System (ADS)
Mémin, Etienne
2014-03-01
We present a derivation of a stochastic model of Navier Stokes equations that relies on a decomposition of the velocity fields into a differentiable drift component and a time uncorrelated uncertainty random term. This type of decomposition is reminiscent in spirit to the classical Reynolds decomposition. However, the random velocity fluctuations considered here are not differentiable with respect to time, and they must be handled through stochastic calculus. The dynamics associated with the differentiable drift component is derived from a stochastic version of the Reynolds transport theorem. It includes in its general form an uncertainty dependent "subgrid" bulk formula that cannot be immediately related to the usual Boussinesq eddy viscosity assumption constructed from thermal molecular agitation analogy. This formulation, emerging from uncertainties on the fluid parcels location, explains with another viewpoint some subgrid eddy diffusion models currently used in computational fluid dynamics or in geophysical sciences and paves the way for new large-scales flow modelling. We finally describe an applications of our formalism to the derivation of stochastic versions of the Shallow water equations or to the definition of reduced order dynamical systems.
A coupled model of fluid flow in jointed rock
Swenson, Daniel; Martineau, Rick; James, Mark; Brown, Don
1991-01-01
We present a fully coupled model of fluid flow in jointed rock, where the fluid flow depends on the joint openings and the joint openings depend on the fluid pressure. The joints and rock blocks are modeled discretely using the finite element method. Solutions for the fluid and rock are obtained and iteration is performed until both solutions converge. Example applications include an examination of the effects of back-pressure on flow in a geothermal reservoir and transient fluid injection into a reservoir.
A knowledge-based approach to automated flow-field zoning for computational fluid dynamics
NASA Technical Reports Server (NTRS)
Vogel, Alison Andrews
1989-01-01
An automated three-dimensional zonal grid generation capability for computational fluid dynamics is shown through the development of a demonstration computer program capable of automatically zoning the flow field of representative two-dimensional (2-D) aerodynamic configurations. The applicability of a knowledge-based programming approach to the domain of flow-field zoning is examined. Several aspects of flow-field zoning make the application of knowledge-based techniques challenging: the need for perceptual information, the role of individual bias in the design and evaluation of zonings, and the fact that the zoning process is modeled as a constructive, design-type task (for which there are relatively few examples of successful knowledge-based systems in any domain). Engineering solutions to the problems arising from these aspects are developed, and a demonstration system is implemented which can design, generate, and output flow-field zonings for representative 2-D aerodynamic configurations.