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
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.
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.
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."
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.
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
Microfluidic flow switching design using volume of fluid model.
Chein, Reiyu; Tsai, S H
2004-03-01
In this study, a volume of fluid (VOF) model was employed for microfluidic switch design. The VOF model validity in predicting the interface between fluid streams with different viscosities co-flowing in a microchannel was first verified by experimental observation. It was then extended to microfluidic flow switch design. Two specific flow switches, one with a guided fluid to one of five desired outlet ports, and another with a guided fluid flows into one, two, or three outlet ports equally distributed along the outlet channel of a Y-shaped channel. The flow switching was achieved by controlling the flow rate ratios between tested and buffer fluids. The numerical results showed that the VOF model could successfully predict the flow switching phenomena in these flow switches. The numerical results also showed that the flow rate ratio required for flow switching depends on the viscosity ratio between the tested and buffer fluids. The numerical simulation was verified by experimental study and the agreement was good. PMID:15307449
Modelling fluid flow in a reciprocating compressor
NASA Astrophysics Data System (ADS)
Tuhovcak, Jan; Hejčík, Jiří; Jícha, Miroslav
2015-05-01
Efficiency of reciprocating compressor is strongly dependent on the valves characteristics, which affects the flow through the suction and discharge line. Understanding the phenomenon inside the compressor is necessary step in development process. Commercial CFD tools offer wide capabilities to simulate the flow inside the reciprocating compressor, however they are too complicated in terms of computational time and mesh creation. Several parameters describing compressor could be therefore examined without the CFD analysis, such is valve characteristic, flow through the cycle and heat transfer. The aim of this paper is to show a numerical tool for reciprocating compressor based on the energy balance through the cycle, which provides valve characteristics, flow through the cycle and heat losses from the cylinder. Spring-damping-mass model was used for the valve description. Boundary conditions were extracted from the performance test of 4-cylinder semihermetic compressor and numerical tool validation was performed with indicated p-V diagram comparison.
A two-fluid model for avalanche and debris flows.
Pitman, E Bruce; Le, Long
2005-07-15
Geophysical mass flows--debris flows, avalanches, landslides--can contain O(10(6)-10(10)) m(3) or more of material, often a mixture of soil and rocks with a significant quantity of interstitial fluid. These flows can be tens of meters in depth and hundreds of meters in length. The range of scales and the rheology of this mixture presents significant modelling and computational challenges. This paper describes a depth-averaged 'thin layer' model of geophysical mass flows containing a mixture of solid material and fluid. The model is derived from a 'two-phase' or 'two-fluid' system of equations commonly used in engineering research. Phenomenological modelling and depth averaging combine to yield a tractable set of equations, a hyperbolic system that describes the motion of the two constituent phases. If the fluid inertia is small, a reduced model system that is easier to solve may be derived. PMID:16011934
Surface tension driven flow in glass melts and model fluids
NASA Technical Reports Server (NTRS)
Mcneil, T. J.; Cole, R.; Subramanian, R. S.
1982-01-01
Surface tension driven flow has been investigated analytically and experimentally using an apparatus where a free column of molten glass or model fluids was supported at its top and bottom faces by solid surfaces. The glass used in the experiments was sodium diborate, and the model fluids were silicone oils. In both the model fluid and glass melt experiments, conclusive evidence was obtained to prove that the observed flow was driven primarily by surface tension forces. The experimental observations are in qualitative agreement with predictions from the theoretical model.
Pulmonary fluid flow challenges for experimental and mathematical modeling.
Levy, Rachel; Hill, David B; Forest, M Gregory; Grotberg, James B
2014-12-01
Modeling the flow of fluid in the lungs, even under baseline healthy conditions, presents many challenges. The complex rheology of the fluids, interaction between fluids and structures, and complicated multi-scale geometry all add to the complexity of the problem. We provide a brief overview of approaches used to model three aspects of pulmonary fluid and flow: the surfactant layer in the deep branches of the lung, the mucus layer in the upper airway branches, and closure/reopening of the airway. We discuss models of each aspect, the potential to capture biological and therapeutic information, and open questions worthy of further investigation. We hope to promote multi-disciplinary collaboration by providing insights into mathematical descriptions of fluid-mechanics in the lung and the kinds of predictions these models can make. PMID:25096289
Pulmonary Fluid Flow Challenges for Experimental and Mathematical Modeling
Levy, Rachel; Hill, David B.; Forest, M. Gregory; Grotberg, James B.
2014-01-01
Modeling the flow of fluid in the lungs, even under baseline healthy conditions, presents many challenges. The complex rheology of the fluids, interaction between fluids and structures, and complicated multi-scale geometry all add to the complexity of the problem. We provide a brief overview of approaches used to model three aspects of pulmonary fluid and flow: the surfactant layer in the deep branches of the lung, the mucus layer in the upper airway branches, and closure/reopening of the airway. We discuss models of each aspect, the potential to capture biological and therapeutic information, and open questions worthy of further investigation. We hope to promote multi-disciplinary collaboration by providing insights into mathematical descriptions of fluid-mechanics in the lung and the kinds of predictions these models can make. PMID:25096289
Benchmarking computational fluid dynamics models for lava flow simulation
NASA Astrophysics Data System (ADS)
Dietterich, Hannah; Lev, Einat; Chen, Jiangzhi
2016-04-01
Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, and COMSOL. Using the new benchmark scenarios defined in Cordonnier et al. (Geol Soc SP, 2015) as a guide, we model viscous, cooling, and solidifying flows over horizontal and sloping surfaces, topographic obstacles, and digital elevation models of natural topography. We compare model results to analytical theory, analogue and molten basalt experiments, and measurements from natural lava flows. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We can apply these models to reconstruct past lava flows in Hawai'i and Saudi Arabia using parameters assembled from morphology, textural analysis, and eruption observations as natural test cases. Our study highlights the strengths and weaknesses of each code, including accuracy and computational costs, and provides insights regarding code selection.
On the coupling between fluid flow and mesh motion in the modelling of fluid structure interaction
NASA Astrophysics Data System (ADS)
Dettmer, Wulf G.; Perić, Djordje
2008-12-01
Partitioned Newton type solution strategies for the strongly coupled system of equations arising in the computational modelling of fluid solid interaction require the evaluation of various coupling terms. An essential part of all ALE type solution strategies is the fluid mesh motion. In this paper, we investigate the effect of the terms which couple the fluid flow with the fluid mesh motion on the convergence behaviour of the overall solution procedure. We show that the computational efficiency of the simulation of many fluid solid interaction processes, including fluid flow through flexible pipes, can be increased significantly if some of these coupling terms are calculated exactly.
A preliminary study to Assess Model Uncertainties in Fluid Flows
Marc Oliver Delchini; Jean C. Ragusa
2009-09-01
The goal of this study is to assess the impact of various flow models for a simplified primary coolant loop of a light water nuclear reactor. The various fluid flow models are based on the Euler equations with an additional friction term, gravity term, momentum source, and energy source. The geometric model is purposefully chosen simple and consists of a one-dimensional (1D) loop system in order to focus the study on the validity of various fluid flow approximations. The 1D loop system is represented by a rectangle; the fluid is heated up along one of the vertical legs and cooled down along the opposite leg. A pressurizer and a pump are included in the horizontal legs. The amount of energy transferred and removed from the system is equal in absolute value along the two vertical legs. The various fluid flow approximations are compressible vs. incompressible, and complete momentum equation vs. Darcy’s approximation. The ultimate goal is to compute the fluid flow models’ uncertainties and, if possible, to generate validity ranges for these models when applied to reactor analysis. We also limit this study to single phase flows with low-Mach numbers. As a result, sound waves carry a very small amount of energy in this particular case. A standard finite volume method is used for the spatial discretization of the system.
Hydromechanical Modeling of Fluid Flow in the Lower Crust
NASA Astrophysics Data System (ADS)
Connolly, J.
2011-12-01
The lower crust lies within an ambiguous rheological regime between the brittle upper crust and ductile sub-lithospheric mantle. This ambiguity has allowed two schools of thought to develop concerning the nature of fluid flow in the lower crust. The classical school holds that lower crustal rocks are inviscid and that any fluid generated by metamorphic devolatilization is squeezed out of rocks as rapidly as it is produced. According to this school, permeability is a dynamic property and fluid flow is upward. In contrast, the modern school uses concepts from upper crustal hydrology that presume implicitly, if not explicitly, that rocks are rigid or, at most, brittle. For the modern school, the details of crustal permeability determine fluid flow and as these details are poorly known almost anything is possible. Reality, to the extent that it is reflected by inference from field studies, offers some support to both schools. In particular, evidence of significant lateral and channelized fluid flow are consistent with flow in rigid media, while evidence for short (104 - 105 y) grain-scale fluid-rock interaction during much longer metamorphic events, suggests that reaction-generated grain-scale permeability is sealed rapidly by compaction; a phenomenon that is also essential to prevent extensive retrograde metamorphism. These observations provide a compelling argument for recognizing in conceptual models of lower crustal fluid flow that rocks are neither inviscid nor rigid, but compact by viscous mechanisms on a finite time-scale. This presentation will review the principle consequences of, and obstacles to, incorporating compaction in such models. The role of viscous compaction in the lower crust is extraordinarily uncertain, but ignoring this uncertainty in models of lower crustal fluid flow does not make the models any more certain. Models inevitably invoke an initial steady state hydraulic regime. This initial steady state is critical to model outcomes because it
A numerical model for dynamic crustal-scale fluid flow
NASA Astrophysics Data System (ADS)
Sachau, Till; Bons, Paul; Gomez-Rivas, Enrique; Koehn, Daniel
2015-04-01
Fluid flow in the crust is often envisaged and modeled as continuous, yet minimal flow, which occurs over large geological times. This is a suitable approximation for flow as long as it is solely controlled by the matrix permeability of rocks, which in turn is controlled by viscous compaction of the pore space. However, strong evidence (hydrothermal veins and ore deposits) exists that a significant part of fluid flow in the crust occurs strongly localized in both space and time, controlled by the opening and sealing of hydrofractures. We developed, tested and applied a novel computer code, which considers this dynamic behavior and couples it with steady, Darcian flow controlled by the matrix permeability. In this dual-porosity model, fractures open depending on the fluid pressure relative to the solid pressure. Fractures form when matrix permeability is insufficient to accommodate fluid flow resulting from compaction, decompression (Staude et al. 2009) or metamorphic dehydration reactions (Weisheit et al. 2013). Open fractures can close when the contained fluid either seeps into the matrix or escapes by fracture propagation: mobile hydrofractures (Bons, 2001). In the model, closing and sealing of fractures is controlled by a time-dependent viscous law, which is based on the effective stress and on either Newtonian or non-Newtonian viscosity. Our simulations indicate that the bulk of crustal fluid flow in the middle to lower upper crust is intermittent, highly self-organized, and occurs as mobile hydrofractures. This is due to the low matrix porosity and permeability, combined with a low matrix viscosity and, hence, fast sealing of fractures. Stable fracture networks, generated by fluid overpressure, are restricted to the uppermost crust. Semi-stable fracture networks can develop in an intermediate zone, if a critical overpressure is reached. Flow rates in mobile hydrofractures exceed those in the matrix porosity and fracture networks by orders of magnitude
Particle hopping vs. fluid-dynamical models for traffic flow
Nagel, K.
1995-12-31
Although particle hopping models have been introduced into traffic science in the 19509, their systematic use has only started recently. Two reasons for this are, that they are advantageous on modem computers, and that recent theoretical developments allow analytical understanding of their properties and therefore more confidence for their use. In principle, particle hopping models fit between microscopic models for driving and fluiddynamical models for traffic flow. In this sense, they also help closing the conceptual gap between these two. This paper shows connections between particle hopping models and traffic flow theory. It shows that the hydrodynamical limits of certain particle hopping models correspond to the Lighthill-Whitham theory for traffic flow, and that only slightly more complex particle hopping models produce already the correct traffic jam dynamics, consistent with recent fluid-dynamical models for traffic flow. By doing so, this paper establishes that, on the macroscopic level, particle hopping models are at least as good as fluid-dynamical models. Yet, particle hopping models have at least two advantages over fluid-dynamical models: they straightforwardly allow microscopic simulations, and they include stochasticity.
Modeling of fluid and heat flow in fractured geothermal reservoirs
Pruess, K.
1988-08-01
In most geothermal reservoirs large-scale permeability is dominated by fractures, while most of the heat and fluid reserves are stored in the rock matrix. Early-time fluid production comes mostly from the readily accessible fracture volume, while reservoir behavior at later time depends upon the ease with which fluid and heat can be transferred from the rock matrix to the fractures. Methods for modeling flow in fractured porous media must be able to deal with this matrix-fracture exchange, the so-called interporosity flow. This paper reviews recent work at Lawrence Berkeley Laboratory on numerical modeling of nonisothermal multiphase flow in fractured porous media. We also give a brief summary of simulation applications to problems in geothermal production and reinjection. 29 refs., 1 fig.
Numerical modeling of fluid flow with rafts: An application to lava flows
NASA Astrophysics Data System (ADS)
Tsepelev, Igor; Ismail-Zadeh, Alik; Melnik, Oleg; Korotkii, Alexander
2016-07-01
Although volcanic lava flows do not significantly affect the life of people, its hazard is not negligible as hot lava kills vegetation, destroys infrastructure, and may trigger a flood due to melting of snow/ice. The lava flow hazard can be reduced if the flow patterns are known, and the complexity of the flow with debris is analyzed to assist in disaster risk mitigation. In this paper we develop three-dimensional numerical models of a gravitational flow of multi-phase fluid with rafts (mimicking rigid lava-crust fragments) on a horizontal and topographic surfaces to explore the dynamics and the interaction of lava flows. We have obtained various flow patterns and spatial distribution of rafts depending on conditions at the surface of fluid spreading, obstacles on the way of a fluid flow, raft landing scenarios, and the size of rafts. Furthermore, we analyze two numerical models related to specific lava flows: (i) a model of fluid flow with rafts inside an inclined channel, and (ii) a model of fluid flow from a single vent on an artificial topography, when the fluid density, its viscosity, and the effusion rate vary with time. Although the studied models do not account for lava solidification, crust formation, and its rupture, the results of the modeling may be used for understanding of flows with breccias before a significant lava cooling.
Beyond poiseuille: preservation fluid flow in an experimental model.
Singh, Saurabh; Randle, Lucy V; Callaghan, Paul T; Watson, Christopher J E; Callaghan, Chris J
2013-01-01
Poiseuille's equation describes the relationship between fluid viscosity, pressure, tubing diameter, and flow, yet it is not known if cold organ perfusion systems follow this equation. We investigated these relationships in an ex vivo model and aimed to offer some rationale for equipment selection. Increasing the cannula size from 14 to 20 Fr increased flow rate by a mean (SD) of 13 (12)%. Marshall's hyperosmolar citrate was three times less viscous than UW solution, but flows were only 45% faster. Doubling the bag pressure led to a mean (SD) flow rate increase of only 19 (13)%, not twice the rate. When external pressure devices were used, 100 mmHg of continuous pressure increased flow by a mean (SD) of 43 (17)% when compared to the same pressure applied initially only. Poiseuille's equation was not followed; this is most likely due to "slipping" of preservation fluid within the plastic tubing. Cannula size made little difference over the ranges examined; flows are primarily determined by bag pressure and fluid viscosity. External infusor devices require continuous pressurisation to deliver high flow. Future studies examining the impact of perfusion variables on graft outcomes should include detailed equipment descriptions. PMID:24062943
Beyond Poiseuille: Preservation Fluid Flow in an Experimental Model
Singh, Saurabh; Randle, Lucy V.; Callaghan, Paul T.; Watson, Christopher J. E.; Callaghan, Chris J.
2013-01-01
Poiseuille's equation describes the relationship between fluid viscosity, pressure, tubing diameter, and flow, yet it is not known if cold organ perfusion systems follow this equation. We investigated these relationships in an ex vivo model and aimed to offer some rationale for equipment selection. Increasing the cannula size from 14 to 20 Fr increased flow rate by a mean (SD) of 13 (12)%. Marshall's hyperosmolar citrate was three times less viscous than UW solution, but flows were only 45% faster. Doubling the bag pressure led to a mean (SD) flow rate increase of only 19 (13)%, not twice the rate. When external pressure devices were used, 100 mmHg of continuous pressure increased flow by a mean (SD) of 43 (17)% when compared to the same pressure applied initially only. Poiseuille's equation was not followed; this is most likely due to “slipping” of preservation fluid within the plastic tubing. Cannula size made little difference over the ranges examined; flows are primarily determined by bag pressure and fluid viscosity. External infusor devices require continuous pressurisation to deliver high flow. Future studies examining the impact of perfusion variables on graft outcomes should include detailed equipment descriptions. PMID:24062943
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
ANFIS modeling for prediction of particle motions in fluid flows
NASA Astrophysics Data System (ADS)
Safdari, Arman; Kim, Kyung Chun
2015-11-01
Accurate dynamic analysis of parcel of solid particles driven in fluid flow system is of interest for many natural and industrial applications such as sedimentation process, study of cloud particles in atmosphere, etc. In this paper, numerical modeling of solid particles in incompressible flow using Eulerian-Lagrangian approach is carried out to investigate the dynamic behavior of particles in different flow conditions; channel and cavity flow. Although modern computers have been well developed, the high computational time and costs for this kind of problems are still demanded. The Lattice Boltzmann Method (LBM) is used to simulate fluid flows and combined with the Lagrangian approach to predict the motion of particles in the range of masses. Some particles are selected, and subjected to Adaptive-network-based fuzzy inference system (ANFIS) to predict the trajectory of moving solid particles. Using a hybrid learning procedure from computational particle movement, the ANFIS can construct an input-output mapping based on fuzzy if-then rules and stipulated computational fluid dynamics prediction pairs. The obtained results from ANFIS algorithm is validated and compared with the set of benchmark data provided based on point-like approach coupled with the LBM method.
Particle-fluid two-phase flow modeling
Mortensen, G.A.; Trapp, J.A. |
1992-09-01
This paper describes a numerical scheme and computer program, DISCON, for the calculation of two-phase flows that does not require the use of flow regime maps. This model is intermediate between-thermal instantaneous and the averaged two-fluid model. It solves the Eulerian continuity, momentum, and energy equations for each liquid control volume, and the Lagrangian mass, momentum, energy, and position equations for each bubble. The bubbles are modeled individually using a large representative number of bubbles thus avoiding the numerical diffusion associated with Eulerian models. DISCON has been used to calculate the bubbling of air through a column of water and the subcooled boiling of water in a flow channel. The results of these calculations are presented.
Particle-fluid two-phase flow modeling
Mortensen, G.A. ); Trapp, J.A. Idaho National Engineering Lab., Idaho Falls, ID )
1992-01-01
This paper describes a numerical scheme and computer program, DISCON, for the calculation of two-phase flows that does not require the use of flow regime maps. This model is intermediate between-thermal instantaneous and the averaged two-fluid model. It solves the Eulerian continuity, momentum, and energy equations for each liquid control volume, and the Lagrangian mass, momentum, energy, and position equations for each bubble. The bubbles are modeled individually using a large representative number of bubbles thus avoiding the numerical diffusion associated with Eulerian models. DISCON has been used to calculate the bubbling of air through a column of water and the subcooled boiling of water in a flow channel. The results of these calculations are presented.
Particle-fluid two-phase flow modeling
NASA Astrophysics Data System (ADS)
Mortensen, G. A.; Trapp, J. A.
This paper describes a numerical scheme and computer program, DISCON, for the calculation of two-phase flows that does not require the use of flow regime maps. This model is intermediate between-thermal instantaneous and the averaged two-fluid model. It solves the Eulerian continuity, momentum, and energy equations for each liquid control volume, and the Lagrangian mass, momentum, energy, and position equations for each bubble. The bubbles are modeled individually using a large representative number of bubbles, thus, avoiding the numerical diffusion associated with Eulerian models. DISCON has been used to calculate the bubbling of air through a column of water and the subcooled boiling of water in a flow channel. The results of these calculations are presented.
GEOSIM: A numerical model for geophysical fluid flow simulation
NASA Technical Reports Server (NTRS)
Butler, Karen A.; Miller, Timothy L.; Lu, Huei-Iin
1991-01-01
A numerical model which simulates geophysical fluid flow in a wide range of problems is described in detail, and comparisons of some of the model's results are made with previous experimental and numerical studies. The model is based upon the Boussinesq Navier-Stokes equations in spherical coordinates, which can be reduced to a cylindrical system when latitudinal walls are used near the pole and the ratio of latitudinal length to the radius of the sphere is small. The equations are approximated by finite differences in the meridional plane and spectral decomposition in the azimuthal direction. The user can specify a variety of boundary and initial conditions, and there are five different spectral truncation options. The results of five validation cases are presented: (1) the transition between axisymmetric flow and baroclinic wave flow in the side heated annulus; (2) the steady baroclinic wave of the side heated annulus; (3) the wave amplitude vacillation of the side heated annulus; (4) transition to baroclinic wave flow in a bottom heated annulus; and (5) the Spacelab Geophysical Fluid Flow Cell (spherical) experiment.
Modelling couplings between reaction, fluid flow and deformation: Kinetics
NASA Astrophysics Data System (ADS)
Malvoisin, Benjamin; Podladchikov, Yury Y.; Connolly, James A. D.
2016-04-01
Mineral assemblages out of equilibrium are commonly found in metamorphic rocks testifying of the critical role of kinetics for metamorphic reactions. As experimentally determined reaction rates in fluid-saturated systems generally indicate complete reaction in less than several years, i.e. several orders of magnitude faster than field-based estimates, metamorphic reaction kinetics are generally thought to be controlled by transport rather than by processes at the mineral surface. However, some geological processes like earthquakes or slow-slip events have shorter characteristic timescales, and transport processes can be intimately related to mineral surface processes. Therefore, it is important to take into account the kinetics of mineral surface processes for modelling fluid/rock interactions. Here, a model coupling reaction, fluid flow and deformation was improved by introducing a delay in the achievement of equilibrium. The classical formalism for dissolution/precipitation reactions was used to consider the influence of the distance from equilibrium and of temperature on the reaction rate, and a dependence on porosity was introduced to model evolution of reacting surface area during reaction. The fitting of experimental data for three reactions typically occurring in metamorphic systems (serpentine dehydration, muscovite dehydration and calcite decarbonation) indicates a systematic faster kinetics close from equilibrium on the dehydration side than on the hydration side. This effect is amplified through the porosity term in the reaction rate since porosity is formed during dehydration. Numerical modelling indicates that this difference in reaction rate close from equilibrium plays a key role in microtextures formation. The developed model can be used in a wide variety of geological systems where couplings between reaction, deformation and fluid flow have to be considered.
Two-fluid model for two-phase flow
NASA Astrophysics Data System (ADS)
Ishii, M.
1987-06-01
The two-fluid model formulation is discussed in detail. The emphasis of the paper is on the three-dimensional formulation and the closure issues. The origin of the interfacial and turbulent transfer terms in the averaged formulation is explained and their original mathematical forms are examined. The interfacial transfer of mass, momentum, and energy is proportional to the interfacial area and driving force. This is not a postulate but a result of the careful examination of the mathematical form of the exact interfacial terms. These two effects are considered separately. Since all the interfacial transfer terms involve the interfacial area concentration, the accurate modeling of the local interfacial area concentration is the first step to be taken for a development of a reliable two-fluid model closure relations. The interfacial momentum interaction has been studied in terms of the standard-drag, lift, virtual mass, and Basset forces. Available analytical and semi-empirical correlations and closure relations are reviewed and existing shortcomings are pointed out. The other major area of importance is the modeling of turbulent transfer in two-phase flow. The two-phase flow turbulence problem is coupled with the phase separation problem even in a steady-state fully developed flow. Thus the two-phase turbulence cannot be understood without understanding the interfacial drag and lift forces accurately. There are some indications that the mixing length type model may not be sufficient to describe the three-dimensional turbulent and flow structures. Although it is a very difficult challenge, the two-phase flow turbulence should be investigated both experimentally and analytically with long time-scale research.
A Generalized Fluid System Simulation Program to Model Flow Distribution in Fluid Networks
NASA Technical Reports Server (NTRS)
Majumdar, Alok; Bailey, John W.; Schallhorn, Paul; Steadman, Todd
1998-01-01
This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.
Complex fluid flow modeling with SPH on GPU
NASA Astrophysics Data System (ADS)
Bilotta, Giuseppe; Hérault, Alexis; Del Negro, Ciro; Russo, Giovanni; Vicari, Annamaria
2010-05-01
We describe an implementation of the Smoothed Particle Hydrodynamics (SPH) method for the simulation of complex fluid flows. The algorithm is entirely executed on Graphic Processing Units (GPUs) using the Compute Unified Device Architecture (CUDA) developed by NVIDIA and fully exploiting their computational power. An increase of one to two orders of magnitude in simulation speed over equivalent CPU code is achieved. A complete modeling of the flow of a complex fluid such as lava is challenging from the modelistic, numerical and computational points of view. The natural topography irregularities, the dynamic free boundaries and phenomena such as solidification, presence of floating solid bodies or other obstacles and their eventual fragmentation make the problem difficult to solve using traditional numerical methods (finite volumes, finite elements): the need to refine the discretization grid in correspondence of high gradients, when possible, is computationally expensive and with an often inadequate control of the error; for real-world applications, moreover, the information needed by the grid refinement may not be available (e.g. because the Digital Elevation Models are too coarse); boundary tracking is also problematic with Eulerian discretizations, more so with complex fluids due to the presence of internal boundaries given by fluid inhomogeneity and presence of solidification fronts. An alternative approach is offered by mesh-free particle methods, that solve most of the problems connected to the dynamics of complex fluids in a natural way. Particle methods discretize the fluid using nodes which are not forced on a given topological structure: boundary treatment is therefore implicit and automatic; the movement freedom of the particles also permits the treatment of deformations without incurring in any significant penalty; finally, the accuracy is easily controlled by the insertion of new particles where needed. Our team has developed a new model based on the
NASA Astrophysics Data System (ADS)
Huang, Yu-Ning
In this work, we derive necessary and sufficient conditions for turbulent secondary flows of a Newtonian fluid and necessary and sufficient conditions for laminar steady secondary flows of a Non-Newtonian fluid in a straight tube. It is found that there is a striking similarity between them. This similarity motivates the assumption used in developing a generalized non-linear K- epsilon model. Based on an analogy that exists between the constitutive relations for turbulent mean flows of a Newtonian fluid and that for laminar flows of a Non -Newtonian fluid, and making use of the constitutive framework of extended thermodynamics, we develop a generalized non -linear K-epsilon model with the same relaxation time as that which appears in the turbulence model proposed by Yakhot, Orszag, Thangam, Gatski and Speziale in 1992. We show that the non-linear K-epsilon model developed by Speziale in 1987 is unable to predict the relaxation phenomena of the Reynolds stresses because of involving no K and dotepsilon , and a coefficient of which leads to a negative relaxation time for the Reynolds stresses. To correct this deficiency, we resort to making use of the relaxation time in the model of Yakhot et al.. The approximate form of our generalized non-linear K-epsilon model, which can predict the relaxation phenomena of the Reynolds stresses and is frame indifferent, is an extension of the standard K-epsilon model and the non-linear K-epsilon model of Speziale.
Draft: Modeling Two-Phase Flow in Porous Media Including Fluid-Fluid Interfacial Area
Crandall, Dustin; Niessner, Jennifer; Hassanizadeh, S Majid
2008-01-01
We present a new numerical model for macro-scale twophase flow in porous media which is based on a physically consistent theory of multi-phase flow.The standard approach for modeling the flow of two fluid phases in a porous medium consists of a continuity equation for each phase, an extended form of Darcy’s law as well as constitutive relationships for relative permeability and capillary pressure. This approach is known to have a number of important shortcomings and, in particular, it does not account for the presence and role of fluid - fluid interfaces. An alternative is to use an extended model which is founded on thermodynamic principles and is physically consistent. In addition to the standard equations, the model uses a balance equation for specific interfacial area. The constitutive relationship for capillary pressure involves not only saturation, but also specific interfacial area. We show how parameters can be obtained for the alternative model using experimental data from a new kind of flow cell and present results of a numerical modeling study
Advances in modelling of biomimetic fluid flow at different scales
2011-01-01
The biomimetic flow at different scales has been discussed at length. The need of looking into the biological surfaces and morphologies and both geometrical and physical similarities to imitate the technological products and processes has been emphasized. The complex fluid flow and heat transfer problems, the fluid-interface and the physics involved at multiscale and macro-, meso-, micro- and nano-scales have been discussed. The flow and heat transfer simulation is done by various CFD solvers including Navier-Stokes and energy equations, lattice Boltzmann method and molecular dynamics method. Combined continuum-molecular dynamics method is also reviewed. PMID:21711847
Wu, Binxin
2010-12-01
In this paper, 12 turbulence models for single-phase non-newtonian fluid flow in a pipe are evaluated by comparing the frictional pressure drops obtained from computational fluid dynamics (CFD) with those from three friction factor correlations. The turbulence models studied are (1) three high-Reynolds-number k-ε models, (2) six low-Reynolds-number k-ε models, (3) two k-ω models, and (4) the Reynolds stress model. The simulation results indicate that the Chang-Hsieh-Chen version of the low-Reynolds-number k-ε model performs better than the other models in predicting the frictional pressure drops while the standard k-ω model has an acceptable accuracy and a low computing cost. In the model applications, CFD simulation of mixing in a full-scale anaerobic digester with pumped circulation is performed to propose an improvement in the effective mixing standards recommended by the U.S. EPA based on the effect of rheology on the flow fields. Characterization of the velocity gradient is conducted to quantify the growth or breakage of an assumed floc size. Placement of two discharge nozzles in the digester is analyzed to show that spacing two nozzles 180° apart with each one discharging at an angle of 45° off the wall is the most efficient. Moreover, the similarity rules of geometry and mixing energy are checked for scaling up the digester. PMID:21047058
Fluid flow induced calcium response in osteoblasts: mathematical modeling.
Su, J H; Xu, F; Lu, X L; Lu, T J
2011-07-28
Fluid flow in the bone lacuno-canalicular network can induce dynamic fluctuation of intracellular calcium concentration ([Ca(2+)](i)) in osteoblasts, which plays an important role in bone remodeling. There has been limited progress in the mathematical modeling of this process probably due to its complexity, which is controlled by various factors such as Ca(2+) channels and extracellular messengers. In this study we developed a mathematical model to describe [Ca(2+)](i) response induced by fluid shear stress (SS) by integrating the major factors involved and analyzed the effects of different experimental setups (e.g. [Ca(2+)](i) baseline, pretreatment with ATP). In this model we considered the ATP release process and the activities of multiple ion channels and purinergic receptors. The model was further verified quantitatively by comparing the simulation results with experimental data reported in literature. The results showed that: (i) extracellular ATP concentration has more significant effect on [Ca(2+)](i) baseline (73% increase in [Ca(2+)](i) with extracellular ATP concentration varying between 0 and 10 μM), as compared to that induced by SS (25% variation in [Ca(2+)](i) with SS varying from 0 to 3.5 Pa); (ii) Pretreatment with ATP-medium results in different [Ca(2+)](i) response as compared to the control group (ATP-free medium) under SS; (iii) Relative [Ca(2+)](i) fluctuation over baseline is more reliable to show the [Ca(2+)](i) response process than the absolute [Ca(2+)](i) response peak. The developed model may improve the experimental design and facilitate our understanding of the mechanotransduction process in osteoblasts. PMID:21665208
Modelling Fault Zone Evolution: Implications for fluid flow.
NASA Astrophysics Data System (ADS)
Moir, H.; Lunn, R. J.; Shipton, Z. K.
2009-04-01
Flow simulation models are of major interest to many industries including hydrocarbon, nuclear waste, sequestering of carbon dioxide and mining. One of the major uncertainties in these models is in predicting the permeability of faults, principally in the detailed structure of the fault zone. Studying the detailed structure of a fault zone is difficult because of the inaccessible nature of sub-surface faults and also because of their highly complex nature; fault zones show a high degree of spatial and temporal heterogeneity i.e. the properties of the fault change as you move along the fault, they also change with time. It is well understood that faults influence fluid flow characteristics. They may act as a conduit or a barrier or even as both by blocking flow across the fault while promoting flow along it. Controls on fault hydraulic properties include cementation, stress field orientation, fault zone components and fault zone geometry. Within brittle rocks, such as granite, fracture networks are limited but provide the dominant pathway for flow within this rock type. Research at the EU's Soultz-sous-Forệt Hot Dry Rock test site [Evans et al., 2005] showed that 95% of flow into the borehole was associated with a single fault zone at 3490m depth, and that 10 open fractures account for the majority of flow within the zone. These data underline the critical role of faults in deep flow systems and the importance of achieving a predictive understanding of fault hydraulic properties. To improve estimates of fault zone permeability, it is important to understand the underlying hydro-mechanical processes of fault zone formation. In this research, we explore the spatial and temporal evolution of fault zones in brittle rock through development and application of a 2D hydro-mechanical finite element model, MOPEDZ. The authors have previously presented numerical simulations of the development of fault linkage structures from two or three pre-existing joints, the results of
Mathematical modeling of slope flows with entrainment as flows of non-Newtonian fluids
NASA Astrophysics Data System (ADS)
Zayko, Julia; Eglit, Margarita
2015-04-01
Non-Newtonian fluids in which the shear stresses are nonlinear functions of the shear strain rates are used to model slope flows such as snow avalanches, mudflows, debris flows. The entrainment of bottom material is included into the model basing on the assumption that in entraining flows the bed friction is equal to the shear stress of the bottom material (Issler et al, 2011). Unsteady motion down long homogeneous slopes with constant inclines is studied numerically for different flow rheologies and different slope angles. Variation of the velocity profile, increase of the flow depth and velocity due to entrainment as well as the value of the entrainment rate is calculated. Asymptotic formulae for the entrainment rate are derived for unsteady flows of different rheological properties. REFERENCES Chowdhury M., Testik F., 2011. Laboratory testing of mathematical models for high-concentration fluid mud turbidity currents. Ocean Engineering 38, 256-270. Eglit, M.E., Demidov, K.S., 2005. Mathematical modeling of snow entrainment in avalanche motion. Cold Reg. Sci. Technol. 43 (1-2), 10-23. Eglit M. E., Yakubenko A. E., 2012, Mathematical Modeling of slope flows entraining bottom material. Eglit M. E., Yakubenko A. E., 2014, Numerical modeling of slope flows entraining bottom material. Cold Reg. Sci. Technol. 108, 139-148. Issler D, M. Pastor Peréz. 2011. Interplay of entrainment and rheology in snow avalanches; a numerical study. Annals of Glaciology, 52(58), pp.143-147 Kern M. A., Tiefenbacher F., McElwaine J., N., 2004. The rheology of snow in large chute flows. Cold Regions Science and Technology, 39, 181 -192. Naaim, M., Faug, T., Naaim-Bouvet, F., 2003. Dry granular flow modelling including erosion and deposition. Surv. Geophys. 24, 569-585. Naaim, M., Naaim-Bouvet, F., Faug, T., Bouchet, A., 2004. Dense snow avalanche modeling: flow, erosion, deposition and obstacle effects. Cold Reg. Sci. Technol. 39, 193-204. Rougier, J & Kern, M 2010, 'Predicting snow
Spatial and Temporal Low-Dimensional Models for Fluid Flow
NASA Technical Reports Server (NTRS)
Kalb, Virginia
2008-01-01
A document discusses work that obtains a low-dimensional model that captures both temporal and spatial flow by constructing spatial and temporal four-mode models for two classic flow problems. The models are based on the proper orthogonal decomposition at two reference Reynolds numbers. Model predictions are made at an intermediate Reynolds number and compared with direct numerical simulation results at the new Reynolds number.
Winters, W.S.
1984-01-01
An overview of the computer code TOPAZ (Transient-One-Dimensional Pipe Flow Analyzer) is presented. TOPAZ models the flow of compressible and incompressible fluids through complex and arbitrary arrangements of pipes, valves, flow branches and vessels. Heat transfer to and from the fluid containment structures (i.e. vessel and pipe walls) can also be modeled. This document includes discussions of the fluid flow equations and containment heat conduction equations. The modeling philosophy, numerical integration technique, code architecture, and methods for generating the computational mesh are also discussed.
Modeling of unsteady-state flows of viscoelastic plastic fluids
Shulman, Z.P.; Dornyak, O.R.; Khusid, B.M.; Ryklina, I.L.; Zal'tsgendler, E.A.
1989-04-01
Unsteady-state flows of media that possess a complex of rheological properties such as elasticity, viscosity and plasticity are studied. The fluid is assumed to exhibit elastic properties at stresses below the fluidity limit. A Trikomi-type boundary-value problem for describing the unsteady-state forced flows in these media is formulated. A difference scheme for the unobstructed calculation is constructed, and the conditions of its efficiency are investigated. The numerical results obtained illustrate the essential effect of elastic properties in the region where the stress is below the fluidity limit; in particular, those cases are studied wherein the period of the elastic shear wave is commensurable with the characteristic hydrodynamic time of the process.
Chang, F.C.; Hull, J.R.; Wang, Y.H.; Blazek, K.E.
1996-02-01
A computer model was developed to predict eddy currents and fluid flows in molten steel. The model was verified by comparing predictions with experimental results of liquid-metal containment and fluid flow in electromagnetic (EM) edge dams (EMDs) designed at Inland Steel for twin-roll casting. The model can optimize the EMD design so it is suitable for application, and minimize expensive, time-consuming full-scale testing. Numerical simulation was performed by coupling a three-dimensional (3-D) finite-element EM code (ELEKTRA) and a 3-D finite-difference fluids code (CaPS-EM) to solve heat transfer, fluid flow, and turbulence transport in a casting process that involves EM fields. ELEKTRA is able to predict the eddy- current distribution and the electromagnetic forces in complex geometries. CaPS-EM is capable of modeling fluid flows with free surfaces. Results of the numerical simulation compared well with measurements obtained from a static test.
Shaded computer graphic techniques for visualizing and interpreting analytic fluid flow models
NASA Technical Reports Server (NTRS)
Parke, F. I.
1981-01-01
Mathematical models which predict the behavior of fluid flow in different experiments are simulated using digital computers. The simulations predict values of parameters of the fluid flow (pressure, temperature and velocity vector) at many points in the fluid. Visualization of the spatial variation in the value of these parameters is important to comprehend and check the data generated, to identify the regions of interest in the flow, and for effectively communicating information about the flow to others. The state of the art imaging techniques developed in the field of three dimensional shaded computer graphics is applied to visualization of fluid flow. Use of an imaging technique known as 'SCAN' for visualizing fluid flow, is studied and the results are presented.
Predicting multidimensional annular flow with a locally based two-fluid model
Antal, S.P.; Edwards, D.P.; Strayer, T.D.
1998-06-01
The purpose of this work was to: develop a methodology to predict annular flows using a multidimensional four-field, two-fluid Computational Fluid Dynamics (CFD) computer code; develop closure models which use the CFD predicted local velocities, phasic volume fractions, etc...; implement a numerical method which allows the discretized equations to have the same characteristics as the differential form; and compare predicted results to local flow field data taken in a R-134a working fluid test section.
Regional Fluid Flow and Basin Modeling in Northern Alaska
Kelley, Karen D., (Edited By)
2007-01-01
INTRODUCTION The foothills of the Brooks Range contain an enormous accumulation of zinc (Zn) in the form of zinc sulfide and barium (Ba) in the form of barite in Carboniferous shale, chert, and mudstone. Most of the resources and reserves of Zn occur in the Red Dog deposit and others in the Red Dog district; these resources and reserves surpass those of most deposits worldwide in terms of size and grade. In addition to zinc and lead sulfides (which contain silver, Ag) and barite, correlative strata host phosphate deposits. Furthermore, prolific hydrocarbon source rocks of Carboniferous and Triassic to Early Jurassic age generated considerable amounts of petroleum that may have contributed to the world-class petroleum resources of the North Slope. Deposits of Zn-Pb-Ag or barite as large as those in the Brooks Range are very rare on a global basis and, accordingly, multiple coincident favorable factors must be invoked to explain their origins. To improve our understanding of these factors and to contribute to more effective assessments of resources in sedimentary basins of northern Alaska and throughout the world, the Mineral Resources Program and the Energy Resources Program of the U.S. Geological Survey (USGS) initiated a project that was aimed at understanding the petroleum maturation and mineralization history of parts of the Brooks Range that were previously poorly characterized. The project, titled ?Regional Fluid Flow and Basin Modeling in Northern Alaska,? was undertaken in collaboration with industry, academia, and other government agencies. This Circular contains papers that describe the results of the recently completed project. The studies that are highlighted in these papers have led to a better understanding of the following: *The complex sedimentary facies relationships and depositional settings and the geochemistry of the sedimentary rocks that host the deposits (sections 2 and 3). *The factors responsible for formation of the barite and zinc deposits
Multi-fluid pipe flow model for analysis of wellbore dynamics
NASA Astrophysics Data System (ADS)
Krasnopolsky, B.; Starostin, A.; Osiptsov, A.
2012-09-01
We present the results of developing a multi-fluid model for a transient multiphase flow in a pipe with application to oil and gas well cleanup after drilling or hydraulic fracturing. The model is an extension of the earlier work on multi-fluid models for near-horizontal pipes to cover any inclination angle and include additional phases. The 1D multi-fluid model governs a transient isothermal flow of three compressible fluids in a horizontal, deviated, or a vertical pipe. The multi-fluid model research code is based on well-known SIMPLE schemes, and features implementation of either a geometry or a mass conservation based algorithms. The multi-fluid model is able to predict the flow regimes with a single interface between gas and liquid mixture: stratified flow, flow with rolling waves, and flow with slugs. The code is validated against an analytical solutions and flow pattern maps. Agreement with theory is found in predicting the transition between the regimes. The final application of the model will be in transient wellbore dynamics with flow regime transitions.
E-1 Dynamic Fluid-Flow Model Update: EASY/ROCETS Enhancement and Model Development Support
NASA Technical Reports Server (NTRS)
Follett, Randolph F.; Taylor, Robert P.
1998-01-01
This report documents the research conducted to update computer models for dynamic fluid flow simulation of the E-1 test stand subsystems at te NASA John C. Stennis Space Center.Work also involved significant upgrades to the capabilities of EASY/ROCKETS library through the inclusion of the NIST-12 thermodynamic property database and development of new control system modules.
Fluid Structure Modelling of Blood Flow in Vessels.
Moatamedi, M; Souli, M; Al-Bahkali, E
2014-12-01
This paper describes the capabilities of fluid structure interaction based multi-physics numerical modelling in solving problems related to vascular biomechanics. In this research work, the onset of a pressure pulse was simulated at the entrance of a three dimensional straight segment of the blood vessel like circular tube and the resulting dynamic response in the form of a propagating pulse wave through the wall was analysed and compared. Good agreement was found between the numerical results and the theoretical description of an idealized artery. Work has also been done on implementing the material constitutive models specific for vascular applications. PMID:26336693
Field and modelling studies of immiscible fluid flow above a contaminated water-table aquifer
Herkelrath, W.N.; Essaid, H.I.; Hess, K.M.
1991-01-01
A method was developed for measuring the spatial distribution of immiscible liquid contaminants in the subsurface. Fluid saturation distributions measured at a crude-oil spill site were used to test a numerical multiphase flow model.
A fully coupled porous flow and geomechanics model for fluid driven cracks: a peridynamics approach
NASA Astrophysics Data System (ADS)
Ouchi, Hisanao; Katiyar, Amit; York, Jason; Foster, John T.; Sharma, Mukul M.
2015-03-01
A state-based non-local peridynamic formulation is presented for simulating fluid driven fractures in an arbitrary heterogeneous poroelastic medium. A recently developed peridynamic formulation of porous flow has been coupled with the existing peridynamic formulation of solid and fracture mechanics resulting in a peridynamic model that for the first time simulates poroelasticity and fluid-driven fracture propagation. This coupling is achieved by modeling the role of pore pressure on the deformation of porous media and vice versa through porosity variation with medium deformation, pore pressure and total mean stress. The poroelastic model is verified by simulating the one-dimensional consolidation of fluid saturated rock. An additional porous flow equation with material permeability dependent on fracture width is solved to simulate fluid flow in the fractured region. Finally, single fluid-driven fracture propagation with a two-dimensional plane strain assumption is simulated and verified against the corresponding classical analytical solution.
Modelling of fluid-structure interaction with multiphase viscous flows using an immersed-body method
NASA Astrophysics Data System (ADS)
Yang, P.; Xiang, J.; Fang, F.; Pavlidis, D.; Latham, J.-P.; Pain, C. C.
2016-09-01
An immersed-body method is developed here to model fluid-structure interaction for multiphase viscous flows. It does this by coupling a finite element multiphase fluid model and a combined finite-discrete element solid model. A coupling term containing the fluid stresses is introduced within a thin shell mesh surrounding the solid surface. The thin shell mesh acts as a numerical delta function in order to help apply the solid-fluid boundary conditions. When used with an advanced interface capturing method, the immersed-body method has the capability to solve problems with fluid-solid interfaces in the presence of multiphase fluid-fluid interfaces. Importantly, the solid-fluid coupling terms are treated implicitly to enable larger time steps to be used. This two-way coupling method has been validated by three numerical test cases: a free falling cylinder in a fluid at rest, elastic membrane and a collapsing column of water moving an initially stationary solid square. A fourth simulation example is of a water-air interface with a floating solid square being moved around by complex hydrodynamic flows including wave breaking. The results show that the immersed-body method is an effective approach for two-way solid-fluid coupling in multiphase viscous flows.
Hydraulic modeling of unsteady debris-flow surges with solid-fluid interactions
Iverson, Richard M.
1997-01-01
Interactions of solid and fluid constituents produce the unique style of motion that typifies debris flows. To simulate this motion, a new hydraulic model represents debris flows as deforming masses of granular solids variably liquefied by viscous pore fluid. The momentum equation of the model describes how internal and boundary forces change as coarse-grained surge heads dominated by grain-contact friction grade into muddy debris-flow bodies more strongly influenced by fluid viscosity and pressure. Scaling analysis reveals that pore-pressure variations can cause flow resistance in surge heads to surpass that in debris-flow bodies by orders of magnitude. Numerical solutions of the coupled momentum and continuity equations provide good predictions of unsteady, nonuniform motion of experimental debris flows from initiation through deposition.
An accurate and comprehensive model of thin fluid flows with inertia on curved substrates
NASA Astrophysics Data System (ADS)
Roberts, A. J.; Li, Zhenquan
2006-04-01
Consider the three-dimensional flow of a viscous Newtonian fluid upon a curved two-dimensional substrate when the fluid film is thin, as occurs in many draining, coating and biological flows. We derive a comprehensive model of the dynamics of the film, the model being expressed in terms of the film thickness eta and the average lateral velocity bar{bm u}. Centre manifold theory assures us that the model accurately and systematically includes the effects of the curvature of substrate, gravitational body force, fluid inertia and dissipation. The model resolves wavelike phenomena in the dynamics of viscous fluid flows over arbitrarily curved substrates such as cylinders, tubes and spheres. We briefly illustrate its use in simulating drop formation on cylindrical fibres, wave transitions, three-dimensional instabilities, Faraday waves, viscous hydraulic jumps, flow vortices in a compound channel and flow down and up a step. These models are the most complete models for thin-film flow of a Newtonian fluid; many other thin-film models can be obtained by different restrictions and truncations of the model derived here.
Modeling deformation-induced fluid flow in cortical bone's canalicular-lacunar system.
Gururaja, S; Kim, H J; Swan, C C; Brand, R A; Lakes, R S
2005-01-01
To explore the potential role that load-induced fluid flow plays as a mechano-transduction mechanism in bone adaptation, a lacunar-canalicular scale bone poroelasticity model is developed and implemented. The model uses micromechanics to homogenize the pericanalicular bone matrix, a system of straight circular cylinders in the bone matrix through which bone fluids can flow, as a locally anisotropic poroelastic medium. In this work, a simplified two-dimensional model of a periodic array of lacunae and their surrounding systems of canaliculi is used to quantify local fluid flow characteristics in the vicinity of a single lacuna. When the cortical bone model is loaded, microscale stress, and strain concentrations occur in the vicinity of individual lacunae and give rise to microscale spatial variations in the pore fluid pressure field. Furthermore, loading of the bone matrix containing canaliculi generates fluid pressures in the contained fluids. Consequently, loading of cortical bone induces fluid flow in the canaliculi and exchange of fluid between canaliculi and lacunae. For realistic bone morphology parameters, and a range of loading frequencies, fluid pressures and fluid-solid drag forces in the canalicular bone are computed and the associated energy dissipation in the models compared to that measured in physical in vitro experiments on human cortical bone. The proposed model indicates that deformation-induced fluid pressures in the lacunar-canalicular system have relaxation times on the order of milliseconds as opposed to the much shorter times (hundredths of milliseconds) associated with deformation-induced pressures in the Haversian system. PMID:15709702
Direct pore-level modeling of incompressible fluid flow in porous media
Ovaysi, Saeed; Piri, Mohammad
2010-09-20
We present a dynamic particle-based model for direct pore-level modeling of incompressible viscous fluid flow in disordered porous media. The model is capable of simulating flow directly in three-dimensional high-resolution micro-CT images of rock samples. It is based on moving particle semi-implicit (MPS) method. We modify this technique in order to improve its stability for flow in porous media problems. Using the micro-CT image of a rock sample, the entire medium, i.e., solid and fluid, is discretized into particles. The incompressible Navier-Stokes equations are then solved for each particle using the MPS summations. The model handles highly irregular fluid-solid boundaries effectively. An algorithm to split and merge fluid particles is also introduced. To handle the computational load, we present a parallel version of the model that runs on distributed memory computer clusters. The accuracy of the model is validated against the analytical, numerical, and experimental data available in the literature. The validated model is then used to simulate both unsteady- and steady-state flow of an incompressible fluid directly in a representative elementary volume (REV) size micro-CT image of a naturally-occurring sandstone with 3.398 {mu}m resolution. We analyze the quality and consistency of the predicted flow behavior and calculate absolute permeability using the steady-state flow rate.
Soltani, M; Chen, P
2013-01-01
Modeling of interstitial fluid flow involves processes such as fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. To date, majority of microvascular flow modeling has been done at different levels and scales mostly on simple tumor shapes with their capillaries. However, with our proposed numerical model, more complex and realistic tumor shapes and capillary networks can be studied. Both blood flow through a capillary network, which is induced by a solid tumor, and fluid flow in tumor's surrounding tissue are formulated. First, governing equations of angiogenesis are implemented to specify the different domains for the network and interstitium. Then, governing equations for flow modeling are introduced for different domains. The conservation laws for mass and momentum (including continuity equation, Darcy's law for tissue, and simplified Navier-Stokes equation for blood flow through capillaries) are used for simulating interstitial and intravascular flows and Starling's law is used for closing this system of equations and coupling the intravascular and extravascular flows. This is the first study of flow modeling in solid tumors to naturalistically couple intravascular and extravascular flow through a network. This network is generated by sprouting angiogenesis and consisting of one parent vessel connected to the network while taking into account the non-continuous behavior of blood, adaptability of capillary diameter to hemodynamics and metabolic stimuli, non-Newtonian blood flow, and phase separation of blood flow in capillary bifurcation. The incorporation of the outlined components beyond the previous models provides a more realistic prediction of interstitial fluid flow pattern in solid tumors and surrounding tissues. Results predict higher interstitial pressure, almost two times, for realistic model compared to the simplified model. PMID:23840579
Soltani, M.; Chen, P.
2013-01-01
Modeling of interstitial fluid flow involves processes such as fluid diffusion, convective transport in extracellular matrix, and extravasation from blood vessels. To date, majority of microvascular flow modeling has been done at different levels and scales mostly on simple tumor shapes with their capillaries. However, with our proposed numerical model, more complex and realistic tumor shapes and capillary networks can be studied. Both blood flow through a capillary network, which is induced by a solid tumor, and fluid flow in tumor’s surrounding tissue are formulated. First, governing equations of angiogenesis are implemented to specify the different domains for the network and interstitium. Then, governing equations for flow modeling are introduced for different domains. The conservation laws for mass and momentum (including continuity equation, Darcy’s law for tissue, and simplified Navier–Stokes equation for blood flow through capillaries) are used for simulating interstitial and intravascular flows and Starling’s law is used for closing this system of equations and coupling the intravascular and extravascular flows. This is the first study of flow modeling in solid tumors to naturalistically couple intravascular and extravascular flow through a network. This network is generated by sprouting angiogenesis and consisting of one parent vessel connected to the network while taking into account the non-continuous behavior of blood, adaptability of capillary diameter to hemodynamics and metabolic stimuli, non-Newtonian blood flow, and phase separation of blood flow in capillary bifurcation. The incorporation of the outlined components beyond the previous models provides a more realistic prediction of interstitial fluid flow pattern in solid tumors and surrounding tissues. Results predict higher interstitial pressure, almost two times, for realistic model compared to the simplified model. PMID:23840579
(Reactive fluid flow models and applications to diagenesis, mineral deposits and crustal rocks)
Lasaga, A.C.; Rye, D.M.
1991-01-01
The objective of this proposal is to put constraints on fluid and mass flux by combining theoretical and field studies of coupled fluid flow, heat and mass transport, and chemical reaction in hydrothermal and metamorphic systems. We are presently applying the two-dimensional code 2DREACT to the study of reactive flow in hydrothermal system. The code represents a new approach to the modeling of simultaneous reaction and transport which was completed in the first months of this project.
An Image-Based Model of Fluid Flow Through Lymph Nodes.
Cooper, Laura J; Heppell, James P; Clough, Geraldine F; Ganapathisubramani, Bharathram; Roose, Tiina
2016-01-01
The lymphatic system returns fluid to the bloodstream from the tissues to maintain tissue fluid homeostasis. Lymph nodes distributed throughout the system filter the lymphatic fluid. The afferent and efferent lymph flow conditions of lymph nodes can be measured in experiments; however, it is difficult to measure the flow within the nodes. In this paper, we present an image-based modelling approach to investigating how the internal structure of the node affects the fluid flow pathways within the node. Selective plane illumination microscopy images of murine lymph nodes are used to identify the geometry and structure of the tissue within the node and to determine the permeability of the lymph node interstitium to lymphatic fluid. Experimental data are used to determine boundary conditions and optimise the parameters for the model. The numerical simulations conducted within the model are implemented in COMSOL Multiphysics, a commercial finite element analysis software. The parameter fitting resulted in the estimate that the average permeability for lymph node tissue is of the order of magnitude of [Formula: see text]. Our modelling shows that the flow predominantly takes a direct path between the afferent and efferent lymphatics and that fluid is both filtered and absorbed across the blood vessel boundaries. The amount that is absorbed or extravasated in the model is dependent on the efferent lymphatic lumen fluid pressure. PMID:26690921
NASA Technical Reports Server (NTRS)
Cho, Y. I.; Crawford, D. W.; Back, L. H.; Back, M. R.
1987-01-01
A flow visualization study using selective dye injection and frame by frame analysis of a movie provided qualitative and quantitative data on the motion of marked fluid particles in a 60 degree artery branch model for simulation of physiological femoral artery flow. Physical flow features observed included jetting of the branch flow into the main lumen during the brief reverse flow period, flow separation along the main lumen wall during the near zero flow phase of diastole when the core flow was in the downstream direction, and inference of flow separation conditions along the wall opposite the branch later in systole at higher branch flow ratios. There were many similarities between dye particle motions in pulsatile flow and the comparative steady flow observations.
Development of Efficient Real-Fluid Model in Simulating Liquid Rocket Injector Flows
NASA Technical Reports Server (NTRS)
Cheng, Gary; Farmer, Richard
2003-01-01
The characteristics of propellant mixing near the injector have a profound effect on the liquid rocket engine performance. However, the flow features near the injector of liquid rocket engines are extremely complicated, for example supercritical-pressure spray, turbulent mixing, and chemical reactions are present. Previously, a homogeneous spray approach with a real-fluid property model was developed to account for the compressibility and evaporation effects such that thermodynamics properties of a mixture at a wide range of pressures and temperatures can be properly calculated, including liquid-phase, gas- phase, two-phase, and dense fluid regions. The developed homogeneous spray model demonstrated a good success in simulating uni- element shear coaxial injector spray combustion flows. However, the real-fluid model suffered a computational deficiency when applied to a pressure-based computational fluid dynamics (CFD) code. The deficiency is caused by the pressure and enthalpy being the independent variables in the solution procedure of a pressure-based code, whereas the real-fluid model utilizes density and temperature as independent variables. The objective of the present research work is to improve the computational efficiency of the real-fluid property model in computing thermal properties. The proposed approach is called an efficient real-fluid model, and the improvement of computational efficiency is achieved by using a combination of a liquid species and a gaseous species to represent a real-fluid species.
NASA Astrophysics Data System (ADS)
Cartwright, Ian
Advection-dispersion fluid flow models implicitly assume that the infiltrating fluid flows through an already fluid-saturated medium. However, whether rocks contain a fluid depends on their reaction history, and whether any initial fluid escapes. The behaviour of different rocks may be illustrated using hypothetical marble compositions. Marbles with diverse chemistries (e.g. calcite + dolomite + quartz) are relatively reactive, and will generally produce a fluid during heating. By contrast, marbles with more restricted chemistries (e.g. calcite + quartz or calcite-only) may not. If the rock is not fluid bearing when fluid infiltration commences, mineralogical reactions may produce a reaction-enhanced permeability in calcite + dolomite + quartz or calcite + quartz, but not in calcite-only marbles. The permeability production controls the pattern of mineralogical, isotopic, and geochemical resetting during fluid flow. Tracers retarded behind the mineralogical fronts will probably be reset as predicted by the advection-dispersion models; however, tracers that are expected to be reset ahead of the mineralogical fronts cannot progress beyond the permeability generating reaction. In the case of very unreactive lithologies (e.g. pure calcite marbles, cherts, and quartzites), the first reaction to affect the rocks may be a metasomatic one ahead of which there is little pervasive resetting of any tracer. Centimetre-scale layering may lead to the formation of self-perpetuating fluid channels in rocks that are not fluid saturated due to the juxtaposition of reactants. Such layered rocks may show patterns of mineralogical resetting that are not predicted by advection-dispersion models. Patterns of mineralogical and isotopic resetting in marbles from a number of terrains, for example: Chillagoe, Marulan South, Reynolds Range (Australia); Adirondack Mountains, Old Woman Mountains, Notch Peak (USA); and Stephen Cross Quarry (Canada) vary as predicted by these models.
Computational modelling of the flow of viscous fluids in carbon nanotubes
NASA Astrophysics Data System (ADS)
Khosravian, N.; Rafii-Tabar, H.
2007-11-01
Carbon nanotubes will have extensive application in all areas of nano-technology, and in particular in the field of nano-fluidics, wherein they can be used for molecular separation, nano-scale filtering and as nano-pipes for conveying fluids. In the field of nano-medicine, nanotubes can be functionalized with various types of receptors to act as bio-sensors for the detection and elimination of cancer cells, or be used as bypasses and even neural connections. Modelling fluid flow inside nanotubes is a very challenging problem, since there is a complex interplay between the motion of the fluid and the stability of the walls. A critical issue in the design of nano-fluidic devices is the induced vibration of the walls, due to the fluid flow, which can promote structural instability. It has been established that the resonant frequencies depend on the flow velocity. We have studied, for the first time, the flow of viscous fluids through multi-walled carbon nanotubes, using the Euler-Bernoulli classical beam theory to model the nanotube as a continuum structure. Our aim has been to compute the effect of the fluid flow on the structural stability of the nanotubes, without having to consider the details of the fluid-walls interaction. The variations of the resonant frequencies with the flow velocity are obtained for both unembedded nanotubes, and when they are embedded in an elastic medium. It is found that a nanotube conveying a viscous fluid is more stable against vibration-induced buckling than a nanotube conveying a non-viscous fluid, and that the aspect ratio plays the same role in both cases.
Modeling the Effect of Fluid Flow on a Growing Network of Fractures in a Porous Medium
NASA Astrophysics Data System (ADS)
Alhashim, Mohammed; Koch, Donald
2015-11-01
The injection of a viscous fluid at high pressure in a geological formation induces the fracturing of pre-existing joints. Assuming a constant solid-matrix stress field, a weak joint saturated with fluid is fractured when the fluid pressure exceeds a critical value that depends on the joint's orientation. In this work, the formation of a network of fractures in a porous medium is modeled. When the average length of the fractures is much smaller than the radius of a cluster of fractured joints, the fluid flow within the network can be described as Darcy flow in a permeable medium consisting of the fracture network. The permeability and porosity of the medium are functions of the number density of activated joints and consequently depend on the fluid pressure. We demonstrate conditions under which these relationships can be derived from percolation theory. Fluid may also be lost from the fracture network by flowing into the permeable rock matrix. The solution of the model shows that the cluster radius grows as a power law with time in two regimes: (1) an intermediate time regime when the network contains many fractures but fluid loss is negligible; and (2) a long time regime when fluid loss dominates. In both regimes, the power law exponent depends on the Euclidean dimension and the injection rate dependence on time.
Unsteady fluid dynamic model for propeller induced flow fields
NASA Technical Reports Server (NTRS)
Katz, Joseph; Ashby, Dale L.; Yon, Steven
1991-01-01
A potential flow based three-dimensional panel method was modified to treat time dependent flow conditions in which the body's geometry may vary with time. The main objective of this effort was the study of a flow field due to a propeller rotating relative to a nonrotating body which is otherwise moving at a constant forward speed. Calculated surface pressure, thrust and torque coefficient data for a four-bladed marine propeller/body compared favorably with previously published experimental results.
Fault zone architecture and fluid flow: Insights from field data and numerical modeling
NASA Astrophysics Data System (ADS)
Caine, Jonathan Saul; Forster, Craig B.
Fault zones in the upper crust are typically composed of complex fracture networks and discrete zones of comminuted and geochemically altered fault rocks. Determining the patterns and rates of fluid flow in these distinct structural discontinuities is a three-dimensional problem. A series of numerical simulations of fluid flow in a set of three-dimensional discrete fracture network models aids in identifying the primary controlling parameters of fault-related fluid flow, and their interactions, throughout episodic deformation. Four idealized, but geologically realistic, fault zone architectural models are based on fracture data collected along exposures of the Stillwater Fault Zone in Dixie Valley, Nevada and geometric data from a series of normal fault zones in east Greenland. The models are also constrained by an Andersonian model for mechanically compatible fracture networks associated with normal faulting. Fluid flow in individual fault zone components, such as a fault core and damage zone, and full outcrop scale model domains are simulated using a finite element routine. Permeability contrasts between components and permeability anisotropy within components are identified as the major controlling factors in fault-related fluid flow. Additionally, the structural and hydraulic variations in these components are also major controls of flow at the scale of the full model domains. The four models can also be viewed as a set of snapshots in the mechanical evolution of a single fault zone. Changes in the hydraulic parameters within the models mimic the evolution of the permeability structure of each model through a single deformation cycle. The model results demonstrate that small changes in the architecture and hydraulic parameters of individual fault zone components can have very large impacts, up to five orders of magnitude, on the permeability structure of the full model domains. Closure of fracture apertures in each fault zone magnifies the magnitude and
Preece, D.S. Perkins, E.D.
1999-02-10
Techniques for modeling oil well sand production have been developed using the formulations for superquadric discrete elements and Darcy fluid flow. Discrete element models are generated using the new technique of particle cloning. Discrete element sources and sinks allow simulation of sand production from the initial state through the transition to an equilibrium state where particles are created and removed at the same rate.
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.
Fluid simulation of tokamak ion temperature gradient turbulence with zonal flow closure model
NASA Astrophysics Data System (ADS)
Yamagishi, Osamu; Sugama, Hideo
2016-03-01
Nonlinear fluid simulation of turbulence driven by ion temperature gradient modes in the tokamak fluxtube configuration is performed by combining two different closure models. One model is a gyrofluid model by Beer and Hammett [Phys. Plasmas 3, 4046 (1996)], and the other is a closure model to reproduce the kinetic zonal flow response [Sugama et al., Phys. Plasmas 14, 022502 (2007)]. By including the zonal flow closure, generation of zonal flows, significant reduction in energy transport, reproduction of the gyrokinetic transport level, and nonlinear upshift on the critical value of gradient scale length are observed.
Oxygen isotopic transport and exchange during fluid flow: One-dimensional models and applications
Bowman, J.R. ); Willett, S.D. ); Cook, S.J. Environ Corp., Houston, TX )
1994-01-01
In this work the authors investigate the consequences of fluid flow and fluid-rock interaction to the isotopic evolution of fluids and rock with one-dimensional transport models of fluid flow and oxygen isotope exchange. Transport models dealing with stable isotopes are well established in recent geochemical literature. The authors extend previous treatments by presenting the derivation of both analytical and numerical solutions to the transport equations incorporating simultaneously advection, diffusion and hydrodynamic dispersion, and kinetics of isotopic exchange. The increased generality of numerical solutions allows the incorporation of other effects which control the spatial patterns of [delta][sup 18]O values developed in rocks and fluids including multiple reactive species and temperature gradients. The authors discuss the effects of flow parameters, conditions of isotopic exchange, and temperature gradients on the spatial patterns of isotopic shifts produced in rock sequences subjected to fluid flow, and on conventionally calculated W/R ratios for these rock sequences. Finally, the authors examine the implications of oxygen isotope transport for two natural systems where isotopic shifts or gradients could be interpreted in terms of unidirectional fluid infiltration. Solutions of one-dimensional transport equations including the mechanisms of advection, diffusion, hydrodynamic dispersion, and non-equilibrium exchange between water and rock indicate that the time-space evolution of oxygen isotopic compositions of rock and infiltrating fluid is dependent on (1) the rate of fluid infiltration, (2) the diffusive and dispersive properties of the rock matrix, (3) the rate of isotopic exchange, and (4) the rock-water mass oxygen ratio in a unit volume of water-saturated, porous rock. 56 refs., 18 figs., 2 tabs.
A two-phase solid/fluid model for dense granular flows including dilatancy effects
NASA Astrophysics Data System (ADS)
Mangeney, Anne; Bouchut, Francois; Fernandez-Nieto, Enrique; Koné, El-Hadj; Narbona-Reina, Gladys
2016-04-01
Describing grain/fluid interaction in debris flows models is still an open and challenging issue with key impact on hazard assessment [{Iverson et al.}, 2010]. We present here a two-phase two-thin-layer model for fluidized debris flows that takes into account dilatancy effects. It describes the velocity of both the solid and the fluid phases, the compression/dilatation of the granular media and its interaction with the pore fluid pressure [{Bouchut et al.}, 2016]. The model is derived from a 3D two-phase model proposed by {Jackson} [2000] based on the 4 equations of mass and momentum conservation within the two phases. This system has 5 unknowns: the solid and fluid velocities, the solid and fluid pressures and the solid volume fraction. As a result, an additional equation inside the mixture is necessary to close the system. Surprisingly, this issue is inadequately accounted for in the models that have been developed on the basis of Jackson's work [{Bouchut et al.}, 2015]. In particular, {Pitman and Le} [2005] replaced this closure simply by imposing an extra boundary condition at the surface of the flow. When making a shallow expansion, this condition can be considered as a closure condition. However, the corresponding model cannot account for a dissipative energy balance. We propose here an approach to correctly deal with the thermodynamics of Jackson's model by closing the mixture equations by a weak compressibility relation following {Roux and Radjai} [1998]. This relation implies that the occurrence of dilation or contraction of the granular material in the model depends on whether the solid volume fraction is respectively higher or lower than a critical value. When dilation occurs, the fluid is sucked into the granular material, the pore pressure decreases and the friction force on the granular phase increases. On the contrary, in the case of contraction, the fluid is expelled from the mixture, the pore pressure increases and the friction force diminishes. To
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.
A numerical model of deformation and fluid-flow in an evolving thrust wedge
NASA Astrophysics Data System (ADS)
Strayer, Luther M.; Hudleston, Peter J.; Lorig, Loren J.
2001-06-01
To investigate deformation and fluid-flow in an actively deforming tectonic wedge, we model the evolution of a large, two-dimensional (100 km long, 5 km thick), mechanically and hydrologically homogeneous and isotropic pile of sedimentary strata that is deformed to become a thrust wedge. We compare both 'dry' and 'wet' cases, in order to illustrate: (1) the relative importance of fluids on wedge evolution, and (2) the effect of brittle deformation on fluid-flow. We use an elastic-plastic constitutive relation, including a Mohr-Coulomb failure criterion and a non-associated flow rule, and coupled fluid flow, with bulk rock properties that approximate typical foreland sedimentary strata. Simulations are made both with and without dilation. The model is fully dynamic, but inertial forces remain small. Results show that deformation within the wedge is accommodated by reverse-slip movement on shear bands, which migrate in both directions through the wedge as both fore- and back-thrusts. The model has features predicted by the critical-taper theory: (1) overall wedge geometry; (2) crudely self-similar growth during evolution; (3) more intense deformation toward the rear of the wedge. The models show strong overall in-sequence faulting behavior with major thrusts isolating relatively undeformed packages, which are moved in a piggyback manner upon the active thrusts. Intermittent out-of-sequence faulting does however occur, in order to maintain the wedge taper. Fluid-flow in the deforming wedge is dominated by topography, but is also strongly affected by dilational plastic deformation. In all simulations, there is focused fluid flow within fault zones. When mechanical time-stepping is shut off (uncoupled), flow systems evolve to steady-state where inflow equals outflow. By subtracting the two 'states' we isolate the mechanical fluid response from the total coupled system response. The mechanical fluid response is manifest as contours of head and pressure difference and
Morin, Kristen T; Lenz, Michelle S; Labat, Caroline A; Tranquillo, Robert T
2015-05-01
Knowledge is limited about fluid flow in tissues containing engineered microvessels, which can be substantially different in topology than native capillary networks. A need exists for a computational model that allows for flow through tissues dense in nonpercolating and possibly nonperfusable microvessels to be efficiently evaluated. A finite difference (FD) model based on Poiseuille flow through a distribution of straight tubes acting as point sources and sinks, and Darcy flow through the interstitium, was developed to describe fluid flow through a tissue containing engineered microvessels. Accuracy of the FD model was assessed by comparison to a finite element (FE) model for the case of a single tube. Because the case of interest is a tissue with microvessels aligned with the flow, accuracy was also assessed in depth for a corresponding 2D FD model. The potential utility of the 2D FD model was then explored by correlating metrics of flow through the model tissue to microvessel morphometric properties. The results indicate that the model can predict the density of perfused microvessels based on parameters that can be easily measured. PMID:25424905
Chemical vapor deposition fluid flow simulation modelling tool
NASA Technical Reports Server (NTRS)
Bullister, Edward T.
1992-01-01
Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.
Predictions of bubbly flows in vertical pipes using two-fluid models in CFDS-FLOW3D code
Banas, A.O.; Carver, M.B.; Unrau, D.
1995-09-01
This paper reports the results of a preliminary study exploring the performance of two sets of two-fluid closure relationships applied to the simulation of turbulent air-water bubbly upflows through vertical pipes. Predictions obtained with the default CFDS-FLOW3D model for dispersed flows were compared with the predictions of a new model (based on the work of Lee), and with the experimental data of Liu. The new model, implemented in the CFDS-FLOW3D code, included additional source terms in the {open_quotes}standard{close_quotes} {kappa}-{epsilon} transport equations for the liquid phase, as well as modified model coefficients and wall functions. All simulations were carried out in a 2-D axisymmetric format, collapsing the general multifluid framework of CFDS-FLOW3D to the two-fluid (air-water) case. The newly implemented model consistently improved predictions of radial-velocity profiles of both phases, but failed to accurately reproduce the experimental phase-distribution data. This shortcoming was traced to the neglect of anisotropic effects in the modelling of liquid-phase turbulence. In this sense, the present investigation should be considered as the first step toward the ultimate goal of developing a theoretically sound and universal CFD-type two-fluid model for bubbly flows in channels.
NASA Astrophysics Data System (ADS)
Dietterich, H. R.; Lev, E.; Chen, J.; Cashman, K. V.; Honor, C.
2015-12-01
Recent eruptions in Hawai'i, Iceland, and Cape Verde highlight the need for improved lava flow models for forecasting and hazard assessment. Existing models used for lava flow simulation range in assumptions, complexity, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess the capabilities of existing models and test the development of new codes, we conduct a benchmarking study of computational fluid dynamics models for lava flows, including VolcFlow, OpenFOAM, Flow3D, and COMSOL. Using new benchmark scenarios defined in Cordonnier et al. (2015) as a guide, we model Newtonian, Herschel-Bulkley and cooling flows over inclined planes, obstacles, and digital elevation models with a wide range of source conditions. Results are compared to analytical theory, analogue and molten basalt experiments, and measurements from natural lava flows. Our study highlights the strengths and weakness of each code, including accuracy and computational costs, and provides insights regarding code selection. We apply the best-fit codes to simulate the lava flows in Harrat Rahat, a predominately mafic volcanic field in Saudi Arabia. Input parameters are assembled from rheology and volume measurements of past flows using geochemistry, crystallinity, and present-day lidar and photogrammetric digital elevation models. With these data, we use our verified models to reconstruct historic and prehistoric events, in order to assess the hazards posed by lava flows for Harrat Rahat.
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.
Dynamic coupling between fluid flow and vein growth in fractures: a 3D numerical model
NASA Astrophysics Data System (ADS)
Schwarz, J.-O.; Enzmann, F.
2012-04-01
Fluid flow is one of the main mass transport mechanisms in the Earth's crust and abundant mineral vein networks are important indicators for fluid flow and fluid rock interaction. These systems are dynamic and part of the so called RTM processes (reaction-transport-mechanics). Understanding of mineral vein systems requires coupling of these processes. Here we present a conceptional model for dynamic vein growth of syntaxial, posttectonic veins generated by advective fluid flow and show first results of a numerical model for this scenario. Vein generation requires three processes to occur: (i) fracture generation by mechanical stress e.g. hydro-fracturing, (ii) flow of a supersaturated fluid on that fracture and (iii) crystallization of phase(s) on or in the fracture. 3D synthetic fractures are generated with the SynFrac code (Ogilvie, et al. 2006). Subsequently solutions of the Navier-Stokes equation for this fracture are computed by a computational fluid dynamics code called GeoDict (Wiegmann 2007). Transport (advective and diffusive) of chemical species to growth sites in the fracture and vein growth are computed by a self-written MATLAB script. The numerical model discretizes the wall rock and fracture geometry by volumetric pixels (voxels). Based on this representation, the model computes the three basic functions for vein generation: (a) nucleation, (b) fluid flow with transport of chemical species and (c) growth. The following conditions were chosen for these three modules. Nucleation is heterogeneous and occurs instantaneously at the wall rock/fracture interface. Advective and diffusive flow of a supersaturated fluid and related transport of chemical species occurs according to the computed fluid flow field by GeoDict. Concentration of chemical species at the inflow is constant, representing external fluid buffering. Changes/decrease in the concentration of chemical species occurs only due to vein growth. Growth of nuclei is limited either by transport of
NASA Astrophysics Data System (ADS)
Liu, Yu-Liang; Zhu, Jie; Luo, Xiao-Shu
2009-09-01
Based on the fluid flow time-delayed model proposed by Misra et al in internet congestion control, one modified time-delayed model is presented, where the influence of the communication delay on the router queue length is investigated in detail. The main advantage of the new model is that its stability domain is larger even without an extra controller. By linear stability analysis and numerical simulation, the effectiveness and feasibility of the novel model in internet congestion control are verified.
Kanyanta, V; Ivankovic, A; Karac, A
2009-08-01
Fluid-structure interaction (FSI) numerical models are now widely used in predicting blood flow transients. This is because of the importance of the interaction between the flowing blood and the deforming arterial wall to blood flow behaviour. Unfortunately, most of these FSI models lack rigorous validation and, thus, cannot guarantee the accuracy of their predictions. This paper presents the comprehensive validation of a two-way coupled FSI numerical model, developed to predict flow transients in compliant conduits such as arteries. The model is validated using analytical solutions and experiments conducted on polyurethane mock artery. Flow parameters such as pressure and axial stress (and precursor) wave speeds, wall deformations and oscillating frequency, fluid velocity and Poisson coupling effects, were used as the basis of this validation. Results show very good comparison between numerical predictions, analytical solutions and experimental data. The agreement between the three approaches is generally over 95%. The model also shows accurate prediction of Poisson coupling effects in unsteady flows through flexible pipes, which up to this stage have only being predicted analytically. Therefore, this numerical model can accurately predict flow transients in compliant vessels such as arteries. PMID:19482285
A coupled fluid flow, deformation and heat transfer model for a twin roll caster
Bradbury, P.J.; Hunt, J.D.
1995-12-31
A numerical model of coupled fluid flow, solidification and plastic deformation for the twin roll casting process has been developed. Non-orthogonal curvilinear coordinates are used to provide a more realistic description of the roll bite geometry while a variable heat transfer coefficient based on the normal pressure at the strip/roll interface, is incorporated through the roll bite.
Hutnak, M.; Hurwitz, S.; Ingebritsen, S.E.; Hsieh, P.A.
2009-01-01
Ground surface displacement (GSD) in large calderas is often interpreted as resulting from magma intrusion at depth. Recent advances in geodetic measurements of GSD, notably interferometric synthetic aperture radar, reveal complex and multifaceted deformation patterns that often require complex source models to explain the observed GSD. Although hydrothermal fluids have been discussed as a possible deformation agent, very few quantitative studies addressing the effects of multiphase flow on crustal mechanics have been attempted. Recent increases in the power and availability of computing resources allow robust quantitative assessment of the complex time-variant thermal interplay between aqueous fluid flow and crustal deformation. We carry out numerical simulations of multiphase (liquid-gas), multicomponent (H 2O-CO2) hydrothermal fluid flow and poroelastic deformation using a range of realistic physical parameters and processes. Hydrothermal fluid injection, circulation, and gas formation can generate complex, temporally and spatially varying patterns of GSD, with deformation rates, magnitudes, and geometries (including subsidence) similar to those observed in several large calderas. The potential for both rapid and gradual deformation resulting from magma-derived fluids suggests that hydrothermal fluid circulation may help explain deformation episodes at calderas that have not culminated in magmatic eruption.
Fluid friction and wall viscosity of the 1D blood flow model.
Wang, Xiao-Fei; Nishi, Shohei; Matsukawa, Mami; Ghigo, Arthur; Lagrée, Pierre-Yves; Fullana, Jose-Maria
2016-02-29
We study the behavior of the pulse waves of water into a flexible tube for application to blood flow simulations. In pulse waves both fluid friction and wall viscosity are damping factors, and difficult to evaluate separately. In this paper, the coefficients of fluid friction and wall viscosity are estimated by fitting a nonlinear 1D flow model to experimental data. In the experimental setup, a distensible tube is connected to a piston pump at one end and closed at another end. The pressure and wall displacements are measured simultaneously. A good agreement between model predictions and experiments was achieved. For amplitude decrease, the effect of wall viscosity on the pulse wave has been shown as important as that of fluid viscosity. PMID:26862041
NASA Astrophysics Data System (ADS)
Shen, Yi; Diplas, Panayiotis
2008-01-01
SummaryComplex flow patterns generated by irregular channel topography, such as boulders, submerged large woody debris, riprap and spur dikes, provide unique habitat for many aquatic organisms. Numerical modeling of the flow structures surrounding these obstructions is challenging, yet it represents an important tool for aquatic habitat assessment. In this study, the ability of two- (2-D) and three-dimensional (3-D) computational fluid dynamics models to reproduce these localized complex flow features is examined. The 3-D model is validated with laboratory data obtained from the literature for the case of a flow around a hemisphere under emergent and submerged conditions. The performance of the 2-D and 3-D models is then evaluated by comparing the numerical results with field measurements of flow around several boulders located at a reach of the Smith River, a regulated mountainous stream, obtained at base and peak flows. Close agreement between measured values and the velocity profiles predicted by the two models is obtained outside the wakes behind the hemisphere and boulders. However, the results suggest that in the vicinity of these obstructions the 3-D model is better suited for reproducing the circulation flow behavior at both low and high discharges. Application of the 2-D and 3-D models to meso-scale stream flows of ecological significance is furthermore demonstrated by using a recently developed spatial hydraulic metric to quantify flow complexity surrounding a number of brown trout spawning sites. It is concluded that the 3-D model can provide a much more accurate description of the heterogeneous velocity patterns favored by many aquatic species over a broad range of flows, especially under deep flow conditions when the various obstructions are submerged. Issues pertaining to selection of appropriate models for a variety of flow regimes and potential implication of the 3-D model on the development of better habitat suitability criteria are discussed. The
Numerical Modeling on Two phase Fluid flow in a Coupled Fracture-Skin-Matrix System
NASA Astrophysics Data System (ADS)
Valsala Kumari, R.; G, S. K.
2015-12-01
Multiphase flow modeling studies below the ground surface is very essential for designing suitable remediation strategies for contaminated aquifers and for the development of petroleum and geothermal reservoirs. Presence of fractured bedrock beneath the ground surface will make multiphase flow process more complex due to its highly heterogeneous nature. A major challenge in modeling flow within a fractured rock is to capture the interaction between the high permeability fracture and the low permeability rock-matrix. In some instances, weathering and mineral depositions will lead to formation of an additional layer named fracture-skin at the fracture-matrix interface. Porosity and permeability of fracture-skin may significantly vary from the adjacent rock matrix and this variation will result in different flow and transport behavior within the fracture-skin. In the present study, an attempt has been made to model simultaneous flow of two immiscible phases (water and LNAPL) in a saturated coupled fracture-skin-matrix system. A fully-implicit finite difference model has been developed to simulate the variation of pressure and saturation of fluid phases along the fracture and within the rock-matrix. Sensitivity studies have been done to analyze the effect of change of various fracture-skin parameters such as porosity, diffusion coefficient and thickness on pressure and saturation distribution of both wetting and non-wetting fluid phases. It can be concluded from the study that the presence of fracture-skin is significantly affecting the fluid flow at the fracture-matrix interface and it can also be seen from the study that the flow behavior of both fluid phases is sensitive to fracture-skin parameters.
Flow of non-Newtonian blood analog fluids in rigid curved and straight artery models.
Mann, D E; Tarbell, J M
1990-01-01
The influence of non-Newtonian rheology on wall shear rate in steady and oscillatory flow through rigid curved and straight artery models was studied experimentally. Wall shear rates measured by flush mounted hot film anemometry under nearly identical flow conditions are reported for the following four fluids: aqueous glycerin (Newtonian), aqueous polyacrylamide (shear thinning, highly elastic), aqueous Xanthan gum (shear thinning, moderately elastic), and bovine blood. For steady flow conditions there was little difference at any measurement site in the wall shear rate levels measured for the four fluids. However, large differences were apparent for oscillatory flows, particularly at the inner curvature 180 degrees from the entrance of the curved artery model. At that position the peak wall shear rate for polyacrylamide was 5-6 times higher than for glycerin and 2-3 times higher than for bovine blood. It is concluded that polyacylamide is too elastic to provide a good model of blood flow under oscillatory conditions, particularly when there is wall shear reversal. Xanthan gum and glycerin are better analog fluids, but neither is entirely satisfactory. PMID:2271763
Massoudi, M.C.; Tran, P.X.
2007-06-15
After providing a brief review of the constitutive modeling of the stress tensor for granular materials using non-Newtonian fluid models, we study the flow between two horizontal flat plates. It is assumed that the granular media behaves as a non-Newtonian fluid (of the Reiner–Rivlin type); we use the constitutive relation derived by Rajagopal and Massoudi [Rajagopal, K. R. and M. Massoudi, “A Method for measuring material moduli of granular materials: flow in an orthogonal rheometer,” Topical Report, DOE/PETC/TR-90/3, 1990] which can predict the normal stress differences. The lower plate is fixed and heated, and the upper plate (which is at a lower temperature than the lower plate) is set into motion with a constant velocity. The steady fully developed flow and the heat transfer equations are made dimensionless and are solved numerically; the effects of different dimensionless numbers and viscous dissipation are discussed.
Modeling of confined turbulent fluid-particle flows using Eulerian and Lagrangian schemes
NASA Technical Reports Server (NTRS)
Adeniji-Fashola, A.; Chen, C. P.
1990-01-01
Two important aspects of fluid-particulate interaction in dilute gas-particle turbulent flows (the turbulent particle dispersion and the turbulence modulation effects) are addressed, using the Eulerian and Lagrangian modeling approaches to describe the particulate phase. Gradient-diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangina formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. Models proposed for both schemes are used to predict turbulent fully-developed gas-solid vertical pipe flow with reasonable accuracy.
Fluid flow model of the Cerro Prieto Geothermal Field based on well log interpretation
Halfman, S.E.; Lippmann, M.J.; Zelwe, R.; Howard, J.H.
1982-08-10
The subsurface geology of the Cerro Prieto geothermal field was analyzed using geophysical and lithologic logs. The distribution of permeable and relatively impermeable units and the location of faults are shown in a geologic model of the system. By incorporating well completion data and downhole temperature profiles into the geologic model, it was possible to determine the direction of geothermal fluid flow and the role of subsurface geologic features that control this movement.
Fluid flow model of the Cerro Prieto geothermal field based on well log interpretation
Halfman, S.E.; Lippmann, M.J.; Zelwer, R.; Howard, J.H.
1982-10-01
The subsurface geology of the Cerro Prieto geothermal field was analyzed using geophysical and lithologic logs. The distribution of permeable and relatively impermeable units and the location of faults are shown in a geologic model of the system. By incorporating well completion data and downhole temperature profiles into the geologic model, it was possible to determine he direction of geothermal fluid flow and the role of subsurface geologic features that control this movement.
Numerical simulation of cavitating flow of liquid helium in a pipe using multi-fluid model
NASA Astrophysics Data System (ADS)
Ishimoto, J.; Oike, M.; Kamijo, K.
2002-05-01
The two-dimensional characteristics of the cavitating flow of liquid helium in a pipe are numerically investigated to realize the further development and high performance of new cryogenic engineering applications. First, the governing equations of the cavitating flow of liquid helium based on the unsteady thermal nonequilibrium multi-fluid model are presented and several flow characteristics are numerically calculated, taking into account the effect of superfluidity. Based on the numerical results, the two-dimensional structure of the turbulent cavitating flow of liquid helium passing through the orifice is shown in detail, and it is also found that the phase transition of the normal fluid to the superfluid and the generation of superfluid counterflow against normal fluid flow are conspicuous in the large gas phase volume fraction region where the liquid to gas phase change actively occurs. Furthermore, it is clarified that the mechanism of the He I to He II phase transition caused by the temperature decrease is due to the deprivation of latent heat for vaporization from the liquid phase. According to these theoretical results, the fundamental characteristics of the cryogenic cavitating flow are predicted.
Modeling Fluid Flow and Microbial Reactions in the Peru Accretionary Complex
NASA Astrophysics Data System (ADS)
Bekins, B. A.; Matmon, D.
2002-12-01
Accretionary complexes are sites where sediment compaction and deeper reactions drive large-scale flow systems that can affect global solute budgets. Extensive modeling and drilling studies have elucidated the origin of the fluids, pore pressures, duration of flow, and major flow paths in these settings. An important research goal is to quantify the effect of these flow systems on global chemical budgets of reactive solutes such as carbon. The Peru margin represents an end member setting that can serve as a basis to extend the results to other margins. The sediments are relatively high in organic carbon with an average value of 2.6%. The subduction rate at ~9 cm/yr and taper angle at 14-17° are among the largest in the world. Recent microbial studies on Ocean Drilling Program Leg 201 at the Peru accretionary margin provide many key elements needed to quantify the processes affecting organic carbon in an accretionary complex. Pore water chemistry data from Site 1230 located in the Peru accretionary prism indicate that sulfate reduction is important in the top 8 mbsf. Below this depth, methanogenesis is the dominant process and methane concentrations are among the highest measured at any site on Leg 201. The presence of high methane concentrations at shallow depths suggests that methane is transported upward in the prism by fluid flow. Measurements of in-situ pore pressures and temperatures also support the presence of upward fluid flow. A single in-situ pressure measurement at ~100 mbsf indicated an overpressure of 0.14 MPa. For a reasonable formation permeability of ~ 10-16 m2, the measured overpressure is adequate to produce flow at a rate of ~5 mm/yr. This rate is comparable to previous model estimates for flow rates in the Peru accretionary prism. In addition, curvature in the downhole temperature profile can best be explained by upward fluid flow of 1-10 mm/yr. These data are used to constrain a two-dimensional coupled fluid flow and reactive transport model
SPH-DCDEM model for arbitrary geometries in free surface solid-fluid flows
NASA Astrophysics Data System (ADS)
Canelas, Ricardo B.; Crespo, Alejandro J. C.; Domínguez, Jose M.; Ferreira, Rui M. L.; Gómez-Gesteira, Moncho
2016-05-01
A unified discretization of rigid solids and fluids is introduced, allowing for resolved simulations of fluid-solid phases within a meshless framework. The numerical solution, attained by Smoothed Particle Hydrodynamics (SPH) and a variation of Discrete Element Method (DEM), the Distributed Contact Discrete Element Method (DCDEM) discretization, is achieved by directly considering solid-solid and solid-fluid interactions. The novelty of the work is centred on the generalization of the coupling of the DEM and SPH methodologies for resolved simulations, allowing for state-of-the-art contact mechanics theories to be used in arbitrary geometries, while fluid to solid and vice versa momentum transfers are accurately described. The methods are introduced, analysed and discussed. Initial validations on the DCDEM and the fluid coupling are presented, drawing from test cases in the literature. An experimental campaign serves as a validation point for complex, large scale solid-fluid flows, where a set of blocks in several configurations is subjected to a dam-break wave. Blocks are tracked and positions are then compared between experimental data and the numerical solutions. A Particle Image Velocimetry (PIV) technique allows for the quantification of the flow field and direct comparison with numerical data. The results show that the model is accurate and is capable of treating highly complex interactions, such as transport of debris or hydrodynamic actions on structures, if relevant scales are reproduced.
Evaluation of a vortex model of turbulent cavity flow. [for computing noise generation of fluid flow
NASA Technical Reports Server (NTRS)
Hardin, J. C.; Block, P. J. W.
1979-01-01
A two dimensional discrete vortex model is applied to turbulent flow over a cavity at a Reynolds number based on cavity half length of approximately 500,000. These model predictions are compared with experimental data on a time average and on a spectral basis. The comparisons indicate that the vortex model reproduces the mean velocity profiles of this complex flow, while the turbulent intensity and Reynolds stress profiles are overpredicted in magnitude. Spectral analyses show that the vortex model introduces unrealistic low frequency power which becomes less dominant as the model becomes more disordered. The spectral power distribution of the model in this state is nearly in agreement with that of the data, somewhat too low at the higher frequencies and too high at the lower frequencies. The model is shown to be capable of reproducing the cavity feedback oscillation phenomenon.
NASA Astrophysics Data System (ADS)
Leggiero, Michael; Bulusu, Kartik V.; Plesniak, Michael W.
2013-11-01
The main objective of this study was to examine inertial effects in a 180-degree model of curved arteries under pulsatile inflow conditions. Two-component, two-dimensional particle image velocimetery (2C-2D PIV) data were acquired upstream of and at several cross-sectional locations in the curved artery model. A blood-analog fluid comprised of 71% saturated sodium iodide solution, 28% glycerol and 1% distilled water (by volume) was subjected to multi-harmonic pulsatile inflow functions. First, signal time-lag was quantified by cross-correlating the input (voltage-time) supplied to a programmable pump and the output PIV (flow rate-time) measurements. The experiment was then treated as a linear, time-invariant system, and frequency response was estimated for phase shifts across a certain spectrum. Input-output signal dissimilarities were attributable to intrinsic inertial effects of flow. By coupling pressure-time and upstream flow rate-time measurements, the experiment was modeled using system identification methods. Results elucidate the role of inertial effects in fluid flow velocity measurements and the effect of these delays on secondary flow structure detection in a curved artery model. Supported by the NSF Grant No. CBET- 0828903 and GW Center for Biomimetics and Bioinspired Engineering.
NASA Astrophysics Data System (ADS)
Antonellini, Marco; Nella Mollema, Pauline
2016-04-01
Surface outcrops provide natural analogs for aquifers and they offer an opportunity to study the geometry of geologic heterogeneities in three dimensions over a range of scales. We show photographs, maps, quantitative field data of rock fractures and sedimentary features in outcrops exposed in a unique collection of many different settings. These include small-scale sedimentary structures, carbonate nodules, faults, and other fractures as documented in outcrops of porous sandstone (Utah, USA and Italy), tight sandstones (Bolivia), dolomite (Northern Italy), and carbonates (Central Italy). We simulate the geometries observed in outcrops with simple conceptual and numerical models of flow to show how important it is to recognize the appropriate attributes for the description and the process responsible for the formation of geologic heterogeneities. For example, knowing the type of structural heterogeneities (fault, joint, compaction band, stylolite, and vein) and their development mechanics helps to predict the distribution and preferential orientation of these features within an aquifer. This knowledge is particularly important for modeling of fluid flow where geophysical or borehole data are lacking. Geologic heterogeneities of sedimentary, structural or diagenetic (chemical) nature influence the fluid flow properties in many aquifers and reservoirs at scales varying over several orders of magnitude and with a spatial variability ranging from mm to tens of meters. Heterogeneities may enhance or degrade porosity and permeability, they impart anisotropy to permeability and dispersion and affect mass transport-related processes in groundwater. Furthermore, aquifer heterogeneities control aquifer continuity and compartmentalization. In fractured aquifers, geologic and diagenetic heterogeneities may affect connectivity, aperture of the flow channels or the distribution of permeability buffers, barriers and seals. Also variations in layer thickness and lithology within a
Fluid and thermal mixing in a model cold leg and downcomer with loop flow
Rothe, P.H.; Ackerson, M.F.
1982-04-01
This report describes an experimental program of fluid mixing experiments performed at atmospheric pressure in a 1/5-scale transparent model of the cold leg and downcomer of typical Westinghouse and Combustion Engineering pressurized water reactors. The results include transient data from a grid of thermocouples and exensive flow visualization photographs. Substantial mixing of cold injected water with hot primary coolant occurred during many of the tests.
NASA Astrophysics Data System (ADS)
Horiuchi, Shun-suke; Iwamori, Hikaru
2016-05-01
Water plays crucial roles in the subduction zone dynamics affecting the thermal-flow structure through the fluid processes. We aim to understand what controls the dynamics and construct a model to solve consistently fluid generation, fluid transport, its reaction with the solid and resultant viscosity, and thermal-flow structure. We highlight the effect of mechanical weakening of rocks associated with hydration. The viscosity of serpentinite (ηserp) in subduction zones critically controls the flow-thermal structure via extent of mechanical coupling between the subducting slab and overlying mantle wedge. When ηserp is greater than 1021 Pa s, the thermal-flow structure reaches a steady state beneath the volcanic zone, and the melting region expands until Cin (initial water content in the subducting oceanic crust) reaches 3 wt %, and it does not expand from 3 wt %. On the other hand, when ηserp is less than 1019 Pa s, the greater water dependence of viscosity (expressed by a larger fv) confines a hot material to a narrower channel intruding into the wedge corner from a deeper part of the back-arc region. Consequently, the overall heat flux becomes less for a larger fv. When ageba (age of back-arc basin as a rifted lithosphere) = 7.5 Ma, the increase in fv weakens but shifts the melting region toward the trench side because of the narrow channel flow intruding into the wedge corner, where as it shuts down melting when ageba=20 Ma. Several model cases (particularly those with ηserp=1020 to 1021 Pa s and a relatively large fv for Cin=2 to 3 wt %) broadly account for the observations in the Northeast Japan arc (i.e., location and width of volcanic chain, extent of serpentinite, surface heat flow, and seismic tomography), although the large variability of surface heat flow and seismic tomographic images does not allow us to constrain the parameter range tightly.
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.
Reactive fluid flow models and applications to diagenesis, mineral deposits and crustal rocks
Lasaga, A.C.; Rye, D.M.
1992-01-01
This project is obtaining new results and developing new techniques along three directions: (a) experimental studies of water-rock reactions (b) theoretical modeling of coupled fluid flow-chemical reactions and (c) isotopic measurements of both regional isotopic compositions as well as isotopic zoning within individual mineral grains. An important part of the project is the integration of all three approaches into a concerted effort aimed at new understanding of the behavior of fluids and their chemical reactions with minerals in the crust. The experimental work pioneered in our laboratories has produced several startling results on the kinetic rate laws of silicate-water reactions. The approach to equilibrium has been shown to follow a non-linear path in rate constant-free energy space. This behavior is quite distinct from most work done by geochemists on modeling silicate behavior in diagenesis, weathering, hydrothermal systems or environmental models. The work to date has involved albite, kaolinite and gibbsite, which together with silica would comprise a kinetic granite system prototype''. The theoretical modeling has produced a state-of-the-art computer code that can efficiently handle dozens of chemical species, many mineral reactions and variations of fluid flow properties and temperature in both one and two dimensions. In addition the code can now treat oxidation-reduction reactions and isotopic exchange between fluids and minerals. The main thrust of the theoretical modeling has been to develop further the differences between equilibrium, steady state, and non-steady-state behavior of the chemical evolution of open fluid-rock systems. These differences have not been fully appreciated in previous models.
Lasaga, A.C.; Rye, D.M.
1992-10-01
This project is obtaining new results and developing new techniques along three directions: (a) experimental studies of water-rock reactions (b) theoretical modeling of coupled fluid flow-chemical reactions and (c) isotopic measurements of both regional isotopic compositions as well as isotopic zoning within individual mineral grains. An important part of the project is the integration of all three approaches into a concerted effort aimed at new understanding of the behavior of fluids and their chemical reactions with minerals in the crust. The experimental work pioneered in our laboratories has produced several startling results on the kinetic rate laws of silicate-water reactions. The approach to equilibrium has been shown to follow a non-linear path in rate constant-free energy space. This behavior is quite distinct from most work done by geochemists on modeling silicate behavior in diagenesis, weathering, hydrothermal systems or environmental models. The work to date has involved albite, kaolinite and gibbsite, which together with silica would comprise a kinetic ``granite system prototype``. The theoretical modeling has produced a state-of-the-art computer code that can efficiently handle dozens of chemical species, many mineral reactions and variations of fluid flow properties and temperature in both one and two dimensions. In addition the code can now treat oxidation-reduction reactions and isotopic exchange between fluids and minerals. The main thrust of the theoretical modeling has been to develop further the differences between equilibrium, steady state, and non-steady-state behavior of the chemical evolution of open fluid-rock systems. These differences have not been fully appreciated in previous models.
NASA Astrophysics Data System (ADS)
Abdelaziz, Ramadan; Sussumu Komori, Fabio
2015-04-01
Recently, Lattice Boltzmann Modelling (LBM) techniques attract many scientists in various fields of research. This work shows the capability for LBM to simulate the fluid flow and solute transport in porous and fracture media, additionally, how to study behavior of nanofluids submitted to a temperature gradient, which it is an important process in natural aquatic environments, water treatment, and other water related technologies. LBSim is used in this work as Lattice Boltzmann Model simulator software. In this article, a series of cases using the lattice Boltzmann method are presented, showing the capability of the method in simulating phenomena with fluid flow and heat transfer in porous media. Results show that the lattice Boltzmann method delivers reliable and helpful simulations for the analyses of processes in water related technologies. Thus, LBSim is a recommended tool for simulating fluid flow at laminar and turbulent condition, and heat and mass transport under complex geometry and boundary condition. parameter values. Keywords: Lattice Boltzmann Model, LBSim, Fractures Media, Porous Media, nanofluids
Surface flow boundary conditions in modeling land subsidence due to fluid withdrawal.
Baú, Domenico; Ferronato, Massimiliano; Gambolati, Giuseppe; Teatini, Pietro
2004-01-01
Land subsidence due to subsurface fluid (water, gas, oil) withdrawal is often predicted by either finite element or finite difference numerical models based on coupled poroelastic theory, where the soil is represented as a semi-infinite medium bounded by the traction-free (ground) surface. One of the variables playing a most important role on the final outcome is the flow condition used on the traction-free boundary, which may be assumed as either permeable or impermeable. Although occasionally justified, the assumption of no-flow surface seems to be in general rather unrealistic. A permeable boundary where the fluid pressure is fixed to the external atmospheric pressure appears to be more appropriate. This paper addresses the response, in terms of land subsidence, obtained with a coupled poroelastic finite element model that simulates a distributed pumping from a horizontal aquifer confined between two relatively impervious layers, and takes either a permeable boundary surface, i.e., constant hydraulic potential, or an impermeable boundary, i.e., a zero Neumann flow condition. The analysis reveals that land subsidence is rather sensitive to the flow condition implemented on the traction-free boundary. In general, the no-flow condition leads to an overestimate of the predicted ground surface settlement, which could even be 1 order of magnitude larger than that obtained with the permeable boundary. PMID:15318774
Mauldon, A.D.; Karasaki, K.; Martel, S.J.; Long, J.C.S.; Landsfield, M.; Mensch, A. ); Vomvoris, S. )
1993-11-01
One of the characteristics of flow and transport in fractured rock is that the flow may be largely confined to a poorly connected network of fractures. In order to represent this condition, Lawrence Berkeley Laboratory has been developing a new type of fracture hydrology model called an equivalent discontinuum model. In this model we represent the discontinuous nature of the problem through flow on a partially filled lattice. This is done through a statistical inverse technique called [open quotes]simulated annealing.[close quotes] The fracture network model is [open quotes]annealed[close quotes] by continually modifying a base model, or [open quotes]template,[close quotes] so that with each modification, the model behaves more and more like the observed system. This template is constructed using geological and geophysical data to identify the regions that possibly conduct fluid an the probable orientations of channels that conduct fluid. In order to see how the simulated annealing algorithm works, we have developed a synthetic case. In this case, the geometry of the fracture network is completely known, so that the results of annealing to steady state data can be evaluated absolutely. We also analyze field data from the Migration Experiment at the Grimsel Rock Laboratory in Switzerland. 28 refs., 14 figs., 3 tabs.
Otto Laporte Lecture: Fluid Dynamics Prize Talk: Simple Models for Turbulent Flows
NASA Astrophysics Data System (ADS)
Pope, Stephen B.
2009-11-01
We focus on the modeling of two turbulent flows: dispersion from a line source in grid turbulence; and, a lifted non-premixed turbulent jet flame. Stochastic Lagrangian models and PDF methods are described, and are shown to model these flows satisfactorily. For the line source, a Lagrangian approach is taken, with the Langevin equation modeling the velocity following a fluid particle, and with a simple relaxation model for the particle temperature. Comparison with experimental data shows that the resulting model describes accurately the dispersion from single and multiple line sources. These simple stochastic Lagrangian models are then applied to the much more challenging case of a lifted non-premixed jet flame. The stochastic Lagrangian models form the basis for a particle/mesh numerical method for solving a modeled transport equation for the Eulerian joint probability density function (PDF) of velocity and composition. The PDF calculations are in excellent agreement with the experimental data, and exhibit the observed extreme sensitivity of the flame to the temperature of the co-flow. The PDF model calculations presented clearly demonstrate that simple models can be very useful, even though aspects of their behavior may be inaccurate or incomplete. The shortcomings of the Langevin equation are examined, and more advanced models (designed to overcome some of these shortcomings) are described. These include models for fluid-particle acceleration, including the effects of intermittency; models accounting for mean shear, which are correct in the rapid- distortion limit; and models designed for use in conjunction with large-eddy simulations (LES).
NASA Astrophysics Data System (ADS)
Adeniji-Fashola, A. A.
1988-07-01
A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.
NASA Technical Reports Server (NTRS)
Adeniji-Fashola, A. A.
1988-01-01
A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.
An approximate single fluid 3-dimensional magnetohydrodynamic equilibrium model with toroidal flow
NASA Astrophysics Data System (ADS)
Cooper, W. A.; Hirshman, S. P.; Chapman, I. T.; Brunetti, D.; Faustin, J. M.; Graves, J. P.; Pfefferlé, D.; Raghunathan, M.; Sauter, O.; Tran, T. M.; Aiba, N.
2014-09-01
An approximate model for a single fluid three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium with pure isothermal toroidal flow with imposed nested magnetic flux surfaces is proposed. It recovers the rigorous toroidal rotation equilibrium description in the axisymmetric limit. The approximation is valid under conditions of nearly rigid or vanishing toroidal rotation in regions with significant 3D deformation of the equilibrium flux surfaces. Bifurcated helical core equilibrium simulations of long-lived modes in the MAST device demonstrate that the magnetic structure is only weakly affected by the flow but that the 3D pressure distortion is important. The pressure is displaced away from the major axis and therefore is not as noticeably helically deformed as the toroidal magnetic flux under the subsonic flow conditions measured in the experiment. The model invoked fails to predict any significant screening by toroidal plasma rotation of resonant magnetic perturbations in MAST free boundary computations.
Fedosov, Dmitry A.; Karniadakis, George Em; Caswell, Bruce
2010-01-01
Polymer fluids are modeled with dissipative particle dynamics (DPD) as undiluted bead-spring chains and their solutions. The models are assessed by investigating their steady shear-rate properties. Non-Newtonian viscosity and normal stress coefficients, for shear rates from the lower to the upper Newtonian regimes, are calculated from both plane Couette and plane Poiseuille flows. The latter is realized as reverse Poiseuille flow (RPF) generated from two Poiseuille flows driven by uniform body forces in opposite directions along two-halves of a computational domain. Periodic boundary conditions ensure the RPF wall velocity to be zero without density fluctuations. In overlapping shear-rate regimes the RPF properties are confirmed to be in good agreement with those calculated from plane Couette flow with Lees–Edwards periodic boundary conditions (LECs), the standard virtual rheometer for steady shear-rate properties. The concentration and the temperature dependence of the properties of the model fluids are shown to satisfy the principles of concentration and temperature superposition commonly employed in the empirical correlation of real polymer-fluid properties. The thermodynamic validity of the equation of state is found to be a crucial factor for the achievement of time-temperature superposition. With these models, RPF is demonstrated to be an accurate and convenient virtual rheometer for the acquisition of steady shear-rate rheological properties. It complements, confirms, and extends the results obtained with the standard LEC configuration, and it can be used with the output from other particle-based methods, including molecular dynamics, Brownian dynamics, smooth particle hydrodynamics, and the lattice Boltzmann method. PMID:20405981
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.
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)
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.
Optimization of a Two-Fluid Hydrodynamic Model of Churn-Turbulent Flow
Donna Post Guillen
2009-07-01
A hydrodynamic model of two-phase, churn-turbulent flows is being developed using the computational multiphase fluid dynamics (CMFD) code, NPHASE-CMFD. The numerical solutions obtained by this model are compared with experimental data obtained at the TOPFLOW facility of the Institute of Safety Research at the Forschungszentrum Dresden-Rossendorf. The TOPFLOW data is a high quality experimental database of upward, co-current air-water flows in a vertical pipe suitable for validation of computational fluid dynamics (CFD) codes. A five-field CMFD model was developed for the continuous liquid phase and four bubble size groups using mechanistic closure models for the ensemble-averaged Navier-Stokes equations. Mechanistic models for the drag and non-drag interfacial forces are implemented to include the governing physics to describe the hydrodynamic forces controlling the gas distribution. The closure models provide the functional form of the interfacial forces, with user defined coefficients to adjust the force magnitude. An optimization strategy was devised for these coefficients using commercial design optimization software. This paper demonstrates an approach to optimizing CMFD model parameters using a design optimization approach. Computed radial void fraction profiles predicted by the NPHASE-CMFD code are compared to experimental data for four bubble size groups.
An improved model of ER fluids in squeeze-flow through model updating of the estimated yield stress
NASA Astrophysics Data System (ADS)
El Wahed, A. K.; Sproston, J. L.; Stanway, R.; Williams, E. W.
2003-11-01
In the squeeze-flow mode of operation, electrorheological (ER) fluid is placed between two electrodes, which are free to translate in a direction roughly parallel to the direction of the applied electric field. Consequently, the ER fluid is subjected to alternate tensile and compressive strokes and shearing of the fluid also occurs. Available displacements are small but large forces are available from compact devices and there are many potential applications, notably in vibration isolation. The present authors have spent several years developing mathematical models to account for the observed behaviour of ER fluids in squeeze-flow. Previous attempts at modelling squeeze-flow behaviour have been partially successful but there have always been discrepancies. These discrepancies have generally been attributed to the difficulty of estimating the yield stress developed within the ER fluid when an electric field is applied. In the present paper, the authors describe a new approach in which the yield stress is determined iteratively by minimizing the difference between observed and predicted values of the transmitted force. Using this technique, force/displacement and force/velocity plots are predicted and compared with values from an experimental facility. It is shown that agreement between model predictions and experimental observation is excellent and significantly better than those obtained using existing models.
Sanford, Ward E.; Pearson, S.C.P.; Kiyosugi, K.; Lehto, H.L.; Saballos, J.A.; Connor, C.B.
2012-01-01
We investigate geologic controls on circulation in the shallow hydrothermal system of Masaya volcano, Nicaragua, and their relationship to surface diffuse degassing. On a local scale (~250 m), relatively impermeable normal faults dipping at ~60° control the flowpath of water vapor and other gases in the vadose zone. These shallow normal faults are identified by modeling of a NE-SW trending magnetic anomaly of up to 2300 nT that corresponds to a topographic offset. Elevated SP and CO2 to the NW of the faults and an absence of CO2 to the SE suggest that these faults are barriers to flow. TOUGH2 numerical models of fluid circulation show enhanced flow through the footwalls of the faults, and corresponding increased mass flow and temperature at the surface (diffuse degassing zones). On a larger scale, TOUGH2 modeling suggests that groundwater convection may be occurring in a 3-4 km radial fracture zone transecting the entire flank of the volcano. Hot water rising uniformly into the base of the model at 1 x 10-5 kg/m2s results in convection that focuses heat and fluid and can explain the three distinct diffuse degassing zones distributed along the fracture. Our data and models suggest that the unusually active surface degassing zones at Masaya volcano can result purely from uniform heat and fluid flux at depth that is complicated by groundwater convection and permeability variations in the upper few km. Therefore isolating the effects of subsurface geology is vital when trying to interpret diffuse degassing in light of volcanic activity.
Multidimensional Model of Fluid Flow and Heat Transfer in Generation-IV Supercritical Water Reactors
Gallaway, Tara; Antal, Steven P.; Podowski, Michael Z.
2006-07-01
This paper is concerned with the mechanistic modeling and theoretical/computational analysis of flow and heat transfer in future Generation-IV Supercritical Water Cooled Reactors (SCWR). The issues discussed in the paper include: the development of analytical models of the properties of supercritical water, and the application of full three-dimensional computational modeling framework to simulate fluid flow and heat transfer in SCWRs. Several results of calculations are shown, including the evaluation of water properties (density, specific heat, thermal conductivity, viscosity, and Prandtl number) near the pseudo-critical temperature for various supercritical pressures, and the CFD predictions using the NPHASE computer code. It is demonstrated that the proposed approach is very promising for future mechanistic analyses of SCWR thermal-hydraulics and safety. (authors)
Zebrafish models of idiopathic scoliosis link cerebrospinal fluid flow defects to spine curvature.
Grimes, D T; Boswell, C W; Morante, N F C; Henkelman, R M; Burdine, R D; Ciruna, B
2016-06-10
Idiopathic scoliosis (IS) affects 3% of children worldwide, yet the mechanisms underlying this spinal deformity remain unknown. Here we show that ptk7 mutant zebrafish, a faithful developmental model of IS, exhibit defects in ependymal cell cilia development and cerebrospinal fluid (CSF) flow. Transgenic reintroduction of Ptk7 in motile ciliated lineages prevents scoliosis in ptk7 mutants, and mutation of multiple independent cilia motility genes yields IS phenotypes. We define a finite developmental window for motile cilia in zebrafish spine morphogenesis. Notably, restoration of cilia motility after the onset of scoliosis blocks spinal curve progression. Together, our results indicate a critical role for cilia-driven CSF flow in spine development, implicate irregularities in CSF flow as an underlying biological cause of IS, and suggest that noninvasive therapeutic intervention may prevent severe scoliosis. PMID:27284198
Meakin, Paul; Tartakovsky, Alexandre M.
2009-07-14
In the subsurface fluids play a critical role by transporting dissolved minerals, colloids and contaminants (sometimes over long distances), by mediating dissolution and precipitation processes and enabling chemical transformations in solution and at mineral surfaces. Although the complex geometries of fracture apertures, fracture networks and pore spaces may make it difficult to accurately predict fluid flow in saturated (single-phase) subsurface systems, well developed methods are available. The simulation of multiphase fluid flow in the subsurface is much more challenging because of the large density and/or viscosity ratios found in important applications (water/air in the vadose zone, water/oil, water/gas, gas/oil and water/oil/gas in oil reservoirs, water/air/non-aqueous phase liquids (NAPL) in contaminated vadose zone systems and gas/molten rock in volcanic systems, for example). In addition, the complex behavior of fluid-fluid-solid contact lines, and its impact on dynamic contact angles, must also be taken into account, and coupled with the fluid flow. Pore network models and simple statistical physics based models such as the invasion percolation and diffusion-limited aggregation models have been used quite extensively. However, these models for multiphase fluid flow are based on simplified models for pore space geometries and simplified physics. Other methods such a lattice Boltzmann and lattice gas models, molecular dynamics, Monte Carlo methods, and particle methods such as dissipative particle dynamics and smoothed particle hydrodynamics are based more firmly on first principles, and they do not require simplified pore and/or fracture geometries. However, they are less (in some cases very much less) computationally efficient that pore network and statistical physics models. Recently a combination of continuum computation fluid dynamics, fluid-fluid interface tracking or capturing and simple models for the dependence of contact angles on fluid velocity
Paul Meakin; Alexandre Tartakovsky
2009-07-01
In the subsurface fluids play a critical role by transporting dissolved minerals, colloids and contaminants (sometimes over long distances), by mediating dissolution and precipitation processes and enabling chemical transformations in solution and at mineral surfaces. Although the complex geometries of fracture apertures, fracture networks and pore spaces may make it difficult to accurately predict fluid flow in saturated (single-phase) subsurface systems, well developed methods are available. The simulation of multiphase fluid flow in the subsurface is much more challenging because of the large density and/or viscosity ratios found in important applications (water/air in the vadose zone, water/oil, water/gas, gas/oil and water/oil/gas in oil reservoirs, water/air/non-aqueous phase liquids (NAPL) in contaminated vadose zone systems and gas/molten rock in volcanic systems, for example). In addition, the complex behavior of fluid-fluid-solid contact lines, and its impact on dynamic contact angles, must also be taken into account, and coupled with the fluid flow. Pore network models and simple statistical physics based models such as the invasion percolation and diffusion-limited aggregation models have been used quite extensively. However, these models for multiphase fluid flow are based on simplified models for pore space geometries and simplified physics. Other methods such a lattice Boltzmann and lattice gas models, molecular dynamics, Monte Carlo methods, and particle methods such as dissipative particle dynamics and smoothed particle hydrodynamics are based more firmly on first principles, and they do not require simplified pore and/or fracture geometries. However, they are less (in some cases very much less) computationally efficient that pore network and statistical physics models. Recently a combination of continuum computation fluid dynamics, fluid-fluid interface tracking or capturing and simple models for the dependence of contact angles on fluid velocity
NASA Astrophysics Data System (ADS)
Kozlov, A. N.
2009-05-01
This paper reports the results of numerical studies of axisymmetric flows in a coaxial plasma accelerator in the presence of a longitudinal magnetic field. The calculations were performed using a two-dimensional two-fluid magnetohydrodynamic model taking into account the Hall effect and the conductivity tensor of the medium. The numerical experiments confirmed the main features of the plasmadynamic processes found previously using analytical and one-fluid models and made it possible to study plasma flows near the electrodes.
Computational Fluid Dynamics Modeling of a Lithium/Thionyl Chloride Battery with Electrolyte Flow
Gu, W.B.; Jungst, Rudolph G.; Nagasubramanian, Ganesan; Wang, C.Y.; Weidner, John.
1999-06-11
A two-dimensional model is developed to simulate discharge of a lithium/thionyl chloride primary battery. The model accounts for not only transport of species and charge, but also the electrode porosity variations and the electrolyte flow induced by the volume reduction caused by electrochemical reactions. Numerical simulations are performed using a finite volume method of computational fluid dynamics. The predicted discharge curves for various temperatures are compared to the experimental data with excellent agreement. Moreover, the simulation results. in conjunction with computer visualization and animation techniques, confirm that cell utilization in the temperature and current range of interest is limited by pore plugging or clogging of the front side of the cathode as a result of LiCl precipitation. The detailed two-dimensional flow simulation also shows that the electrolyte is replenished from the cell header predominantly through the separator into the front of the cathode during most parts of the discharge, especially for higher cell temperatures.
General Fluid System Simulation Program to Model Secondary Flows in Turbomachinery
NASA Technical Reports Server (NTRS)
Majumdar, Alok K.; Van Hoosier, Katherine P.
1995-01-01
The complexity and variety of turbomachinery flow circuits created a need for a general fluid system simulation program for test data anomaly resolution as well as design review. The objective of the paper is to present a computer program that has been developed to support Marshall Space Flight Center's turbomachinery internal flow analysis efforts. The computer program solves for the mass. energy and species conservation equation at each node and flow rate equation at each branch of the network by a novel numerical procedure which is a combination of both Newton-Ralphson and successive substitution method and uses a thermodynamic property program for computing real gas properties. A generalized, robust, modular, and 'user-friendly' computer program has been developed to model internal flow rates, pressures, temperatures, concentrations of gas mixtures and axial thrusts. The program can be used for any network for compressible and incompressible flows, choked flow, change of phase and gaseous mixturecs. The code has been validated by comparing the predictions with Space Shuttle Main Engine test data.
Incompressible SPH Model for Simulating Violent Free-Surface Fluid Flows
NASA Astrophysics Data System (ADS)
Staroszczyk, Ryszard
2014-06-01
In this paper the problem of transient gravitational wave propagation in a viscous incompressible fluid is considered, with a focus on flows with fast-moving free surfaces. The governing equations of the problem are solved by the smoothed particle hydrodynamics method (SPH). In order to impose the incompressibility constraint on the fluid motion, the so-called projection method is applied in which the discrete SPH equations are integrated in time by using a fractional-step technique. Numerical performance of the proposed model has been assessed by comparing its results with experimental data and with results obtained by a standard (weakly compressible) version of the SPH approach. For this purpose, a plane dam-break flow problem is simulated, in order to investigate the formation and propagation of a wave generated by a sudden collapse of a water column initially contained in a rectangular tank, as well as the impact of such a wave on a rigid vertical wall. The results of simulations show the evolution of the free surface of water, the variation of velocity and pressure fields in the fluid, and the time history of pressures exerted by an impacting wave on a wall.
Analysis Of Residence Time Distribution Of Fluid Flow By Axial Dispersion Model
Sugiharto; Su'ud, Zaki; Kurniadi, Rizal; Waris, Abdul; Abidin, Zainal
2010-12-23
Radioactive tracer {sup 82}Br in the form of KBr-82 with activity {+-} 1 mCi has been injected into steel pipeline to qualify the extent dispersion of water flowing inside it. Internal diameter of the pipe is 3 in. The water source was originated from water tank through which the water flow gravitically into the pipeline. Two collimated sodium iodide detectors were used in this experiment each of which was placed on the top of the pipeline at the distance of 8 and 11 m from injection point respectively. Residence time distribution (RTD) curves obtained from injection of tracer are elaborated numerically to find information of the fluid flow properties. The transit time of tracer calculated from the mean residence time (MRT) of each RTD curves is 14.9 s, therefore the flow velocity of the water is 0.2 m/s. The dispersion number, D/uL, for each RTD curve estimated by using axial dispersion model are 0.055 and 0.06 respectively. These calculations are performed after fitting the simulated axial dispersion model on the experiment curves. These results indicated that the extent of dispersion of water flowing in the pipeline is in the category of intermediate.
Analysis Of Residence Time Distribution Of Fluid Flow By Axial Dispersion Model
NASA Astrophysics Data System (ADS)
Sugiharto, Su'ud, Zaki; Kurniadi, Rizal; Waris, Abdul; Abidin, Zainal
2010-12-01
Radioactive tracer 82Br in the form of KBr-82 with activity ± 1 mCi has been injected into steel pipeline to qualify the extent dispersion of water flowing inside it. Internal diameter of the pipe is 3 in. The water source was originated from water tank through which the water flow gravitically into the pipeline. Two collimated sodium iodide detectors were used in this experiment each of which was placed on the top of the pipeline at the distance of 8 and 11 m from injection point respectively. Residence time distribution (RTD) curves obtained from injection of tracer are elaborated numerically to find information of the fluid flow properties. The transit time of tracer calculated from the mean residence time (MRT) of each RTD curves is 14.9 s, therefore the flow velocity of the water is 0.2 m/s. The dispersion number, D/uL, for each RTD curve estimated by using axial dispersion model are 0.055 and 0.06 respectively. These calculations are performed after fitting the simulated axial dispersion model on the experiment curves. These results indicated that the extent of dispersion of water flowing in the pipeline is in the category of intermediate.
Multicomponent model of deformation and detachment of a biofilm under fluid flow.
Tierra, Giordano; Pavissich, Juan P; Nerenberg, Robert; Xu, Zhiliang; Alber, Mark S
2015-05-01
A novel biofilm model is described which systemically couples bacteria, extracellular polymeric substances (EPS) and solvent phases in biofilm. This enables the study of contributions of rheology of individual phases to deformation of biofilm in response to fluid flow as well as interactions between different phases. The model, which is based on first and second laws of thermodynamics, is derived using an energetic variational approach and phase-field method. Phase-field coupling is used to model structural changes of a biofilm. A newly developed unconditionally energy-stable numerical splitting scheme is implemented for computing the numerical solution of the model efficiently. Model simulations predict biofilm cohesive failure for the flow velocity between [Formula: see text] and [Formula: see text] m s(-1) which is consistent with experiments. Simulations predict biofilm deformation resulting in the formation of streamers for EPS exhibiting a viscous-dominated mechanical response and the viscosity of EPS being less than [Formula: see text]. Higher EPS viscosity provides biofilm with greater resistance to deformation and to removal by the flow. Moreover, simulations show that higher EPS elasticity yields the formation of streamers with complex geometries that are more prone to detachment. These model predictions are shown to be in qualitative agreement with experimental observations. PMID:25808342
Multicomponent model of deformation and detachment of a biofilm under fluid flow
Tierra, Giordano; Pavissich, Juan P.; Nerenberg, Robert; Xu, Zhiliang; Alber, Mark S.
2015-01-01
A novel biofilm model is described which systemically couples bacteria, extracellular polymeric substances (EPS) and solvent phases in biofilm. This enables the study of contributions of rheology of individual phases to deformation of biofilm in response to fluid flow as well as interactions between different phases. The model, which is based on first and second laws of thermodynamics, is derived using an energetic variational approach and phase-field method. Phase-field coupling is used to model structural changes of a biofilm. A newly developed unconditionally energy-stable numerical splitting scheme is implemented for computing the numerical solution of the model efficiently. Model simulations predict biofilm cohesive failure for the flow velocity between and m s−1 which is consistent with experiments. Simulations predict biofilm deformation resulting in the formation of streamers for EPS exhibiting a viscous-dominated mechanical response and the viscosity of EPS being less than . Higher EPS viscosity provides biofilm with greater resistance to deformation and to removal by the flow. Moreover, simulations show that higher EPS elasticity yields the formation of streamers with complex geometries that are more prone to detachment. These model predictions are shown to be in qualitative agreement with experimental observations. PMID:25808342
A fiber matrix model for fluid flow and streaming potentials in the canaliculi of an osteon
NASA Technical Reports Server (NTRS)
Zeng, Y.; Cowin, S. C.; Weinbaum, S.
1994-01-01
A theoretical model is developed to predict the fluid shear stress and streaming potential at the surface of osteocytic processes in the lacunar-canalicular porosity of an osteon when the osteon is subject to mechanical loads that are parallel or perpendicular to its axis. The theory developed in Weinbaum et al. (31) for the flow through a proteoglycan matrix in a canaliculus is employed in a poroelastic model for the osteon. Our formulation is a generalization of that of Petrov et al. (17). Our model predicts that, in order to satisfy the measured frequency dependence of the phase and magnitude of the SGP in macroscopic bone samples, the fiber spacing in the fluid annulus must lie in the narrow range 6-7 nm typical of the spacing of GAG sidechains along a protein monomer. The model predictions for the local SGP profiles in the osteon agree with the experimental observations of Starkebaum et al. (24). The theory predicts that the pore pressure relaxation time, tau d, for a 150-300 microns diameter osteon with the foregoing matrix structure is approximately 0.03-0.13 sec, and that the amplitude of the mean fluid shear stress on the membrane of the osteocytic process at the mean areal radius of the osteon has a maximum at 28 Hz if tau d = 0.06 sec. This maximum, which is independent of the magnitude of the loading, could be important in vivo since the recent experiments of Turner et al. (28) and McLeod et al. (15) have a peak in the strain frequency spectrum between 20 and 30 Hz that also appears to be independent of the type (magnitude) of loading. Numerical predictions for the amplitude of the average fluid shear stress on the osteocytic membrane at the mean areal radius of the osteon show that the fluid shear stress associated with the low amplitude 20-30 Hz spectral strain component is at least as large as the average fluid shear stress associated with the high amplitude 1 Hz stride component, although the latter loading is an order of magnitude larger, and has a
SDEM modelling of deformation associated with a listric fault system and associated fluid flow
NASA Astrophysics Data System (ADS)
Rasmussen, Marie L.; Clausen, Ole R.; Egholm, David L.; Andresen, Katrine J.
2016-04-01
Numerical modelling of geological structures using FEM, DEM and SDEM methods as well as analogue modelling are widely used in order to achieve a better understanding of the kinematics and dynamics during deformation. The methods are furthermore the ultimate source for mapping (observing) the true geometry of geological structures as well as subsurface fluid flow phenomena in 3D seismic data developed for hydrocarbon exploration. Here we use 3D seismic data and SDEM modelling to suggest a dynamic-kinematic evolution of the deformation in the hangingwall of a listric fault overlying an active salt roller. We use the results to obtain a better understanding of the fluid flow in a complex deformed hangingwall. The case study is focused at the D-1 fault trend in the western part of the Norwegian Danish Basin, at the northern slope of the Ringkøbing-Fyn High. The D-1 main fault detaches along the northern flank of a Zechstein salt roller which was active during the Cenozoic. The seismic analysis shows a system of secondary normal antithetic and synthetic faults dipping approximately 50-60dg within the hangingwall. Shallow gas is trapped in the hangingwall and the secondary faults often confine the accumulations i.e. indicating that the secondary faults are sealing. The modelling confirms that the geometry of the secondary faults is highly controlled by the rheology of different layers in the hangingwall but also on the intensity of the salt movement. The modelling also suggests the presence of vertical deformation zones; structures which are not directly observed on the seismic data. The vertical deformation zones are related to the differential vertical movement of the strata due to salt migration. A neural network trained chimney probability cube shows high probabilities for the presence of minor vertical gas chimneys below the gas accumulations suggesting that vertical fluid migration in the hangingwall occurred in areas with significant vertical salt movements. The
Inverse model of fully coupled fluid flow and stress in fractured rock masses
NASA Astrophysics Data System (ADS)
Wang, Y.; Rutqvist, J.
2008-12-01
In order to reflect the real behavior of the seepage field and deformation field during the environment change and construction process£¬the basic equations and FEM methods for fully coupled analysis of fluid flow and stress are developed£¬based on the assumptions of small deformation and incompressible water flow in complicated fractured rock masses. Both the equivalent continuum media model and the discrete media model are adopted. And the modified initial flow method is used to deal with the free surface of unconfined seepage. Due to the difficulty in determining the parameters of water flow field, stress field and their coupling relations, an inverse model is presented for the fully coupled problem in which both the data of water head and displacement are taken into consideration. Objective function is defined based on sensitivity analysis of parameters, and the relative values of water head, displacement on parameters are adopted in the establishment of objective function. A hybrid genetic algorithm is proposed as optimization method. The probability of crossover and mutation is determined according to chromosome fitness and a concept of self- adaptive probability is given. In addition, simplex method is also applied to increase the ability of local search, the operation of accelerated cycle is used in order to decrease optimization time.
Nanoscale Pore Imaging and Pore Scale Fluid Flow Modeling in Chalk
Tomutsa, Liviu; Silin, Dmitriy
2004-08-19
For many rocks of high economic interest such as chalk, diatomite, tight gas sands or coal, nanometer scale resolution is needed to resolve the 3D-pore structure, which controls the flow and trapping of fluids in the rocks. Such resolutions cannot be achieved with existing tomographic technologies. A new 3D imaging method, based on serial sectioning and using the Focused Ion Beam (FIB) technology has been developed. FIB allows for the milling of layers as thin as 10 nanometers by using accelerated Ga+ ions to sputter atoms from the sample surface. After each milling step, as a new surface is exposed, a 2D image of this surface is generated. Next, the 2D images are stacked to reconstruct the 3D pore or grain structure. Resolutions as high as 10 nm are achievable using such a technique. A new robust method of pore-scale fluid flow modeling has been developed and applied to sandstone and chalk samples. The method uses direct morphological analysis of the pore space to characterize the petrophysical properties of diverse formations. Not only petrophysical properties (porosity, permeability, relative permeability and capillary pressures) can be computed but also flow processes, such as those encountered in various IOR approaches, can be simulated. Petrophysical properties computed with the new method using the new FIB data will be presented. Present study is a part of the development of an Electronic Core Laboratory at LBNL/UCB.
A mathematical model of fluid and gas flow in nanoporous media.
Monteiro, Paulo J M; Rycroft, Chris H; Barenblatt, Grigory Isaakovich
2012-12-11
The mathematical modeling of the flow in nanoporous rocks (e.g., shales) becomes an important new branch of subterranean fluid mechanics. The classic approach that was successfully used in the construction of the technology to develop oil and gas deposits in the United States, Canada, and the Union of Soviet Socialist Republics becomes insufficient for deposits in shales. In the present article a mathematical model of the flow in nanoporous rocks is proposed. The model assumes the rock consists of two components: (i) a matrix, which is more or less an ordinary porous or fissurized-porous medium, and (ii) specific organic inclusions composed of kerogen. These inclusions may have substantial porosity but, due to the nanoscale of pores, tubes, and channels, have extremely low permeability on the order of a nanodarcy (~109-²¹ m² ) or less. These inclusions contain the majority of fluid: oil and gas. Our model is based on the hypothesis that the permeability of the inclusions substantially depends on the pressure gradient. At the beginning of the development of the deposit, boundary layers are formed at the boundaries of the low-permeable inclusions, where the permeability is strongly increased and intensive flow from inclusions to the matrix occurs. The resulting formulae for the production rate of the deposit are presented in explicit form. The formulae demonstrate that the production rate of deposits decays with time following a power law whose exponent lies between -1/2 and -1. Processing of experimental data obtained from various oil and gas deposits in shales demonstrated an instructive agreement with the prediction of the model. PMID:23188803
A mathematical model of fluid and gas flow in nanoporous media
Monteiro, Paulo J. M.; Rycroft, Chris H.; Barenblatt, Grigory Isaakovich
2012-01-01
The mathematical modeling of the flow in nanoporous rocks (e.g., shales) becomes an important new branch of subterranean fluid mechanics. The classic approach that was successfully used in the construction of the technology to develop oil and gas deposits in the United States, Canada, and the Union of Soviet Socialist Republics becomes insufficient for deposits in shales. In the present article a mathematical model of the flow in nanoporous rocks is proposed. The model assumes the rock consists of two components: (i) a matrix, which is more or less an ordinary porous or fissurized-porous medium, and (ii) specific organic inclusions composed of kerogen. These inclusions may have substantial porosity but, due to the nanoscale of pores, tubes, and channels, have extremely low permeability on the order of a nanodarcy () or less. These inclusions contain the majority of fluid: oil and gas. Our model is based on the hypothesis that the permeability of the inclusions substantially depends on the pressure gradient. At the beginning of the development of the deposit, boundary layers are formed at the boundaries of the low-permeable inclusions, where the permeability is strongly increased and intensive flow from inclusions to the matrix occurs. The resulting formulae for the production rate of the deposit are presented in explicit form. The formulae demonstrate that the production rate of deposits decays with time following a power law whose exponent lies between and . Processing of experimental data obtained from various oil and gas deposits in shales demonstrated an instructive agreement with the prediction of the model. PMID:23188803
Modeling of flow and heat transfer for fluids at supercritical conditions
NASA Astrophysics Data System (ADS)
Gallaway, Tara
2011-12-01
The Supercritical Water Reactor (SCWR) has been proposed as one of the six Generation IV reactor design concepts under consideration. The key feature of the SCWR is that water at supercritical pressures is used as the reactor coolant. At supercritical pressures, the working fluid does not undergo phase change as it is heated, but rather the fluid properties experience dramatic variations throughout what is known as the pseudo-critical region. Highly nonuniform temperature and uid property distributions are expected in the reactor core, which will have a significant impact on turbulence and heat transfer as well as stability limits for future SCWRs. The goal of this work is to understand and predict the effects of these fluid property variations on turbulence and heat transfer throughout the reactor core and to predict the potential onset of dynamic instabilities. CO2 at supercritical conditions is included in the current study due in some part to its use as a viable simulant fluid in place of water for experimental studies. The use of CO2 at supercritical conditions as a reactor coolant has also gained popularity in recent years. Spline-type property models have been developed for both water and CO2 at supercritical pressures in order to include the property variations into a numerical solver. Turbulence and heat transfer models for fluids at supercritical conditions have been developed and implemented into the NPHASE-CMFD computer code. The results of predictions using the proposed models have been compared to experimental data from the Korea Atomic Energy Research Institute (KAERI) for various heat transfer regimes. While no model is without some deficiency, the Chien Low-Reynolds k -- epsilon model performs best at predicting the experimental data. A stability model has been developed and is presented in this dissertation as well. This model utilizes three different solution methods and tests the effects of inlet temperature, mass flow rate, local loss
Computer modeling of fluid flow and combustion in the ISV (In Situ Vitrification) confinement hood
Johnson, R.W.; Paik, S.
1990-09-01
Safety and suitability objectives for the application of the In Situ Vitrification (ISV) technology at the INEL require that the physical processes involved in ISVV be modeled to determine their operational behavior. The mathematical models that have been determined to address the modeling needs adequately for the ISV analysis package are detailed elsewhere. The present report is concerned with the models required for simulating the reacting flow that occurs in the ISV confinement hood. An experimental code named COYOTE has been secured that appears adequate to model the combustion in the confinement hood. The COYOTE code is a two-dimensional, transient, compressible, Eulerian, gas dynamics code for modeling reactive flows. It recognizes nonuniform Cartesian and cylindrical geometry and is based on the ICE (Implicit Continuous-fluid Eulerian) family of solution methods. It includes models for chemical reactions based on chemical kinetics as well as equilibrium chemistry. The mathematical models contained in COYOTE, their discrete analogs, the solution procedure, code structure and some test problems are presented in the report. 12 refs., 17 figs., 6 tabs.
NASA Astrophysics Data System (ADS)
Yeh, J.-L.; Hwang, W.-S.; Chou, C.-L.
1992-10-01
A three-dimensional mathematical model has been developed based on the incorporation of a computational fluid dynamics technique, called SOLA-SURF, and the K-ɛ turbulence model. Numerical solutions of the three-dimensional turbulent Navier-Stokes equations and the K and ɛ equations together with the free surface treatment are presented to study the turbulent flow behavior of molten steel in tundishes. Computed results describing the three-dimensional flow field, particle path lines, residence time distribution curve during steady-state operation are presented. The values of t min, t peak, and t mean derived from the residence time distribution curve are used to evaluate the effects of using various combinations of flow control devices such as dams, weirs, and dams with a hole in the flow field. The computed results were compared with the experimental data obtained from a full-scale plexiglas/water model of tundish. The comparisons exhibited good consistency.
NASA Astrophysics Data System (ADS)
La Pointe, P. R.
2001-12-01
Many significant problems in rock engineering require consideration of fluid flow through natural fractures in rock or the mechanical response of a fractured rock mass. Accurate prediction of flow volumes, rates and mass transport through natural fracture systems, and their mechanical response, is critical for design and licensing of nuclear waste repositories, optimization of recovery from many petroleum reservoirs, and also in solution mining, groundwater resource development and protection, and hardrock civil engineering projects. Many of these projects require a large-scale, 3D numerical model for flow, transport or mechanical simulation or visualization. A fundamental problem in constructing these models is that fracture data from wells, boreholes, geophysical profiles or surface outcrops represents a small portion of the rock volume under consideration. Not only does the data represent a very small proportion of the reservoir or rock mass, it also represents fracturing in very restricted size scales. Thus scaling analysis is critical to accurately constructing a fracture model from the data. This paper describes, through two case examples, how scaling analysis techniques have been used to develop models of natural fracturing to support the design and licensing of a high level nuclear waste repository in Finland, and for optimization of a tertiary recovery project in an aging oil field in the US. A new technique for scaling fracture sizes is presented. Together, these two examples illustrate the importance of the scaling analyses, pitfalls in carrying out the analyses, and new methods to improve the 3D characterization of naturally fractured rock masses.
Kwon, Ronald Y.; Frangos, John A.
2010-01-01
Skeletal adaptation to mechanical loading has been widely hypothesized to involve the stimulation of osteocytes by interstitial fluid flow (IFF). However, direct investigation of this hypothesis has been difficult due in large part to the inability to directly measure IFF velocities within the lacunar–canalicular system. Measurements of fluorescence recovery after photobleaching (FRAP) within individual lacunae could be used to quantify lacunar–canalicular IFF when combined with mathematical modeling. In this study, we used a computational transport model to characterize the relationship between flow frequency (0.5–10 Hz), peak flow velocity (0–300 μm/s), tracer diffusion coefficient (100–300 μm2/s), and transport enhancement (i.e., (k/k0) − 1, where k and k0 are the transport rates in the presence/absence of flow) during lacunar FRAP investigations. We show that this relationship is well described by a simple power law with frequency-dependent coefficients, and is relatively insensitive to variations in lacunar geometry. Using this power law relationship, we estimated peak IFF velocities in hindlimb mice subjected to intramedullary pressurization using values of k and k0 previously obtained from ex vivo lacunar FRAP investigations. Together, our findings suggest that skeletal adaptation in hindlimb suspended mice subjected to dynamic intramedullary pressure occurred in the presence of IFF at levels associated with physiological loading. PMID:21076644
Mukhopadhyay, S.; Tsang, Y.; Finsterle, S.
2009-01-15
A simple conceptual model has been recently developed for analyzing pressure and temperature data from flowing fluid temperature logging (FFTL) in unsaturated fractured rock. Using this conceptual model, we developed an analytical solution for FFTL pressure response, and a semianalytical solution for FFTL temperature response. We also proposed a method for estimating fracture permeability from FFTL temperature data. The conceptual model was based on some simplifying assumptions, particularly that a single-phase airflow model was used. In this paper, we develop a more comprehensive numerical model of multiphase flow and heat transfer associated with FFTL. Using this numerical model, we perform a number of forward simulations to determine the parameters that have the strongest influence on the pressure and temperature response from FFTL. We then use the iTOUGH2 optimization code to estimate these most sensitive parameters through inverse modeling and to quantify the uncertainties associated with these estimated parameters. We conclude that FFTL can be utilized to determine permeability, porosity, and thermal conductivity of the fracture rock. Two other parameters, which are not properties of the fractured rock, have strong influence on FFTL response. These are pressure and temperature in the borehole that were at equilibrium with the fractured rock formation at the beginning of FFTL. We illustrate how these parameters can also be estimated from FFTL data.
Moller, Nancy; Weare J. H.
2008-05-29
Successful exploitation of the vast amount of heat stored beneath the earth’s surface in hydrothermal and fluid-limited, low permeability geothermal resources would greatly expand the Nation’s domestic energy inventory and thereby promote a more secure energy supply, a stronger economy and a cleaner environment. However, a major factor limiting the expanded development of current hydrothermal resources as well as the production of enhanced geothermal systems (EGS) is insufficient knowledge about the chemical processes controlling subsurface fluid flow. With funding from past grants from the DOE geothermal program and other agencies, we successfully developed advanced equation of state (EOS) and simulation technologies that accurately describe the chemistry of geothermal reservoirs and energy production processes via their free energies for wide XTP ranges. Using the specific interaction equations of Pitzer, we showed that our TEQUIL chemical models can correctly simulate behavior (e.g., mineral scaling and saturation ratios, gas break out, brine mixing effects, down hole temperatures and fluid chemical composition, spent brine incompatibilities) within the compositional range (Na-K-Ca-Cl-SO4-CO3-H2O-SiO2-CO2(g)) and temperature range (T < 350°C) associated with many current geothermal energy production sites that produce brines with temperatures below the critical point of water. The goal of research carried out under DOE grant DE-FG36-04GO14300 (10/1/2004-12/31/2007) was to expand the compositional range of our Pitzer-based TEQUIL fluid/rock interaction models to include the important aluminum and silica interactions (T < 350°C). Aluminum is the third most abundant element in the earth’s crust; and, as a constituent of aluminosilicate minerals, it is found in two thirds of the minerals in the earth’s crust. The ability to accurately characterize effects of temperature, fluid mixing and interactions between major rock-forming minerals and hydrothermal and
Desbonnets, Quentin; Broc, Daniel
2012-07-01
It is well known that a fluid may strongly influence the dynamic behaviour of a structure. Many different physical phenomena may take place, depending on the conditions: fluid flow, fluid at rest, little or high displacements of the structure. Inertial effects can take place, with lower vibration frequencies, dissipative effects also, with damping, instabilities due to the fluid flow (Fluid Induced Vibration). In this last case the structure is excited by the fluid. Tube bundles structures are very common in the nuclear industry. The reactor cores and the steam generators are both structures immersed in a fluid which may be submitted to a seismic excitation or an impact. In this case the structure moves under an external excitation, and the movement is influence by the fluid. The main point in such system is that the geometry is complex, and could lead to very huge sizes for a numerical analysis. Homogenization models have been developed based on the Euler equations for the fluid. Only inertial effects are taken into account. A next step in the modelling is to build models based on the homogenization of the Navier-Stokes equations. The papers presents results on an important step in the development of such model: the analysis of the fluid flow in a oscillating tube bundle. The analysis are made from the results of simulations based on the Navier-Stokes equations for the fluid. Comparisons are made with the case of the oscillations of a single tube, for which a lot of results are available in the literature. Different fluid flow pattern may be found, depending in the Reynolds number (related to the velocity of the bundle) and the Keulegan Carpenter number (related to the displacement of the bundle). A special attention is paid to the quantification of the inertial and dissipative effects, and to the forces exchanges between the bundle and the fluid. The results of such analysis will be used in the building of models based on the homogenization of the Navier
Simulation of the Second Grade Fluid Model for Blood Flow through a Tapered Artery with a Stenosis
NASA Astrophysics Data System (ADS)
Nadeem, S.; Noreen Sher, Akbar
2010-06-01
We analyze the blood flow through a tapered artery, assuming the blood to be a second order fluid model. The resulting nonlinear implicit system of partial differential equations is solved by the perturbation method. The expressions for shear stress, velocity, flow rate, wall shear stress and longitudinal impedance are obtained. The physical behavior of different parameters is also discussed, as are trapping phenomena.
Modeling fluid flow in deformation bands with stabilized localization mixed finite elements
NASA Astrophysics Data System (ADS)
Sun, W.; Ostien, J. T.; Foulk, J. W.; Abdeljawad, F.
2012-12-01
Deformation bands in geological materials refer to narrow zones of inhomogeneous strain. Their onset and propagation may cause significant changes in microstructures and therefore profoundly enhance or suppress fluid flow and induce anisotropy. These changes in hydraulic properties have strong implications in geotechnical engineering, carbon dioxide sequestration and nuclear waste storage. The difficulty in modeling such multiphysics phenomena is threefold. 1. Monolithically coupled promechanics formulation may lead to non-physical oscillation in pore pressure near the undrained limit if identical mesh and basis functions are used for pore pressure and displacement. 2. Onsets of deformation bands may lead to non-converging mesh-dependent results if no length scale is introduced to the finite element formulation. 3. Modeling anisotropy induced by the deformation band may require a very fine mesh to capture the sharp pore pressure gradient and results in a computational intensive system. In this study, we introduce a projection-based technique to stabilize a large deformation finite element model that eliminates the non-physical oscillation in pore pressure. Using a 1D analytical solution as guideline, we introduce a simple scheme that can adaptively update the optimal value for the stabilization parameter that can restore stability without over-diffusing the system. This stabilized model is coupled with a localization element technique used to introduce proper length scale to regularize the governing equations and resolve the fluid flow jumps across the deformation bands. Numerical examples are presented to demonstrate the properties and performance of the proposed localized models. *Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000
North Cascadia heat flux and fluid flow from gas hydrates: Modeling 3-D topographic effects
NASA Astrophysics Data System (ADS)
Li, Hong-lin; He, Tao; Spence, George D.
2014-01-01
The bottom-simulating reflector (BSR) of gas hydrate is well imaged from two perpendicular seismic grids in the region of a large carbonate mound, informally called Cucumber Ridge off Vancouver Island. We use a new method to calculate 3-D heat flow map from the BSR depths, in which we incorporate 3-D topographic corrections after calibrated by the drilling results from nearby (Integrated) Ocean Drilling Program Site 889 and Site U1327. We then estimate the associated fluid flow by relating it to the topographically corrected heat flux anomalies. In the midslope region, a heat flux anomaly of 1 mW/m2 can be associated with an approximate focused fluid flow rate of 0.09 mm/yr. Around Cucumber Ridge, high rates of focused fluid flow were observed at steep slopes with values more than double the average regional diffusive fluid discharge rate of 0.56 mm/yr. As well, in some areas of relatively flat seafloor, the focused fluid flow rates still exceeded 0.5 mm/yr. On the seismic lines the regions of focused fluid flow were commonly associated with seismic blanking zones above the BSR and sometimes with strong reflectors below the BSR, indicating that the faults/fractures provide high-permeability pathways for fluids to carry methane from BSR depths to the seafloor. These high fluid flow regions cover mostly the western portion of our area with gas hydrate concentration estimations of ~6% based on empirical correlations from Hydrate Ridge in south off Oregon, significantly higher than previously recognized values of ~2.5% in the eastern portion determined from Site U1327.
Han, Byunghyun; Scicchitano, Vincent; Imhoff, Paul T
2011-03-01
Laboratory procedures were developed to obtain constitutive relations for fluid flow in refuse. Five different types of experiments were conducted for the same waste sample: a drainage experiment, multi step outflow experiment, total porosity measurement, saturated hydraulic conductivity test, and gas permeability tests. To investigate fundamental processes affecting water movement and moisture retention, samples consisted entirely of newspaper. Samples were prepared in two particle sizes and two compaction pressures and packed in compression cells to replicate stress conditions in landfills. Data were modeled using HYDRUS-1D, which allowed alternative conceptual models of the pore space to be assessed. A dual-permeability model performed significantly better than a single-porosity model for water movement, suggesting that a dual domain description is required to describe water flow in landfills with significant amounts of paper and paperboard. However, a single-porosity model was adequate for describing gas transport. Results indicated that properties of the fracture domain, the large openings between refuse particles, are significantly affected by the size of waste materials and compaction, and may be best studied with field-scale measurements. On the other hand properties of the matrix domain, the smaller pore openings within and between refuse particles, are likely amenable to laboratory study because representative samples sizes should be much smaller. PMID:20970978
MODIFIED CHOKE FLOW CRITERION FOR THE TWO-PHASE TWO-FLUID MODEL
Suneet Singh; Vincent A. Mousseau
2009-05-01
A choked condition exists when mass flow rate becomes independent of the downstream conditions. In other words, no information can propagate in the upstream direction under this condition. The real part of the solution of the characteristic equation for the model represents velocity of the signal propagation and the imaginary part is the growth (or decay) rate of that signal. Therefore, if the real part of these eigenvalues is positive then no signal propagates in the upstream direction (choosing downstream direction to be the positive direction) resulting in the choke flow. In order to develop the choke criterion, a non-dimensional form of the characteristic equation is derived for the standard two-phase two-fluid model. The equation is in the terms of a slip Mach number Ms. It can be shown that the slip Mach number is small for many applications including nuclear reactor safety simulations. The eigenvalues of the characteristic equation are obtained as a power series expansion about the point Ms = 0. These eigenvalues are used to develop a choking criterion for the compressible two-phase flows.
Interstitial Fluid Flow and Drug Delivery in Vascularized Tumors: A Computational Model
Welter, Michael; Rieger, Heiko
2013-01-01
Interstitial fluid is a solution that bathes and surrounds the human cells and provides them with nutrients and a way of waste removal. It is generally believed that elevated tumor interstitial fluid pressure (IFP) is partly responsible for the poor penetration and distribution of therapeutic agents in solid tumors, but the complex interplay of extravasation, permeabilities, vascular heterogeneities and diffusive and convective drug transport remains poorly understood. Here we consider–with the help of a theoretical model–the tumor IFP, interstitial fluid flow (IFF) and its impact upon drug delivery within tumor depending on biophysical determinants such as vessel network morphology, permeabilities and diffusive vs. convective transport. We developed a vascular tumor growth model, including vessel co-option, regression, and angiogenesis, that we extend here by the interstitium (represented by a porous medium obeying Darcy's law) and sources (vessels) and sinks (lymphatics) for IFF. With it we compute the spatial variation of the IFP and IFF and determine its correlation with the vascular network morphology and physiological parameters like vessel wall permeability, tissue conductivity, distribution of lymphatics etc. We find that an increased vascular wall conductivity together with a reduction of lymph function leads to increased tumor IFP, but also that the latter does not necessarily imply a decreased extravasation rate: Generally the IF flow rate is positively correlated with the various conductivities in the system. The IFF field is then used to determine the drug distribution after an injection via a convection diffusion reaction equation for intra- and extracellular concentrations with parameters guided by experimental data for the drug Doxorubicin. We observe that the interplay of convective and diffusive drug transport can lead to quite unexpected effects in the presence of a heterogeneous, compartmentalized vasculature. Finally we discuss various
NASA Astrophysics Data System (ADS)
Belferman, Mariana; Katsman, Regina; Agnon, Amotz
2015-04-01
The Levant has been repeatedly devastated by numerous earthquakes since prehistorical time, as recorded in historical documents, archaeological ruins, and sedimentary archives. In order to understand the role of the dynamics of the water bodies in triggering the deformations in the Levant basin, a new theoretical thermo-mechanical model is constructed and extended by including a fluid flow component. The latter is modeled on a basis of two-way poroelastic coupling with momentum equation. This coupling is essential to capture the fluid flow evolution induced by dynamic water loading and to resolve porosity changes. All the components of the model, namely elasticity, creep, plasticity, fluid flow, etc., have been extensively verified and presented. Results of the initial sensitivity analysis addressing the relative importance of each process in earthquakes triggering are discussed. The rich archives of pre-instrumental destructive earthquakes will set constraints for future modeling under the present formulation.
NASA Astrophysics Data System (ADS)
Song, Yongjia; Hu, Hengshan; Rudnicki, John W.
2016-07-01
Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori-Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.
NASA Astrophysics Data System (ADS)
Meakin, Paul; Tartakovsky, Alexandre M.
2009-07-01
In the subsurface, fluids play a critical role by transporting dissolved minerals, colloids, and contaminants (sometimes over long distances); by mediating dissolution and precipitation processes; and by enabling chemical transformations in solution and at mineral surfaces. Although the complex geometries of fracture apertures, fracture networks, and pore spaces may make it difficult to accurately predict fluid flow in saturated (single-phase) subsurface systems, well-developed methods are available. The simulation of multiphase fluid flow in the subsurface is much more challenging because of the large density and/or viscosity ratios found in important applications (water/air in the vadose zone; water/oil, water/gas, gas/oil, and water/oil/gas in hydrocarbon reservoirs; water/air/nonaqueous phase liquids (nonaqueous phase liquids/dense nonaqueous phase liquids) in contaminated vadose zone systems; and gas/molten rock in volcanic systems, for example). In addition, the complex behavior of fluid-fluid-solid contact lines and their impact on dynamic contact angles must also be taken into account and coupled with the fluid flow. Here we review the methods that are currently being used to simulate pore-scale multiphase fluid flow and reactive transport in fractured and porous media. After the introduction, the review begins with an overview of the fundamental physics of multiphase fluids flow followed by a more detailed discussion of the complex dynamic behavior of contact lines and contact angles, an important barrier to accurate pore-scale modeling and simulation. The main part of the review focuses on five different approaches: pore network models, lattice gas and lattice Boltzmann methods, Monte Carlo methods, particle methods (molecular dynamics, dissipative particle dynamics, and smoothed particle hydrodynamics), and traditional grid-based computational fluid dynamics coupled with interface tracking and a contact angle model. Finally, the review closes with a
Interpretation and modeling of the averaged equations for a fluid-solid flow
NASA Technical Reports Server (NTRS)
Shen, Hayley H.; Hwang, Guann-Jiun
1989-01-01
A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.
A unified approach to fluid-flow, geomechanical, and seismic modelling
NASA Astrophysics Data System (ADS)
Yarushina, Viktoriya; Minakov, Alexander
2016-04-01
The perturbations of pore pressure can generate seismicity. This is supported by observations from human activities that involve fluid injection into rocks at high pressure (hydraulic fracturing, CO2 storage, geothermal energy production) and natural examples such as volcanic earthquakes. Although the seismic signals that emerge during geotechnical operations are small both in amplitude and duration when compared to natural counterparts. A possible explanation for the earthquake source mechanism is based on a number of in situ stress measurements suggesting that the crustal rocks are close to its plastic yield limit. Hence, a rapid increase of the pore pressure decreases the effective normal stress, and, thus, can trigger seismic shear deformation. At the same time, little attention has been paid to the fact that the perturbation of fluid pressure itself represents an acoustic source. Moreover, non-double-couple source mechanisms are frequently reported from the analysis of microseismicity. A consistent formulation of the source mechanism describing microseismic events should include both a shear and isotropic component. Thus, improved understanding of the interaction between fluid flow and seismic deformation is needed. With this study we aim to increase the competence in integrating real-time microseismic monitoring with geomechanical modelling such that there is a feedback loop between monitored deformation and stress field modelling. We propose fully integrated seismic, geomechanical and reservoir modelling. Our mathematical formulation is based on fundamental set of force balance, mass balance, and constitutive poro-elastoplastic equations for two-phase media consisting of deformable solid rock frame and viscous fluid. We consider a simplified 1D modelling setup for consistent acoustic source and wave propagation in poro-elastoplastic media. In this formulation the seismic wave is generated due to local changes of the stress field and pore pressure induced by
Hsieh, Paul A.
2001-01-01
This report serves as a user?s guide for two computer models: TopoDrive and ParticleFlow. These two-dimensional models are designed to simulate two ground-water processes: topography-driven flow and advective transport of fluid particles. To simulate topography-driven flow, the user may specify the shape of the water table, which bounds the top of the vertical flow section. To simulate transport of fluid particles, the model domain is a rectangle with overall flow from left to right. In both cases, the flow is under steady state, and the distribution of hydraulic conductivity may be specified by the user. The models compute hydraulic head, ground-water flow paths, and the movement of fluid particles. An interactive visual interface enables the user to easily and quickly explore model behavior, and thereby better understand ground-water flow processes. In this regard, TopoDrive and ParticleFlow are not intended to be comprehensive modeling tools, but are designed for modeling at the exploratory or conceptual level, for visual demonstration, and for educational purposes.
Vaughan, T J; Mullen, C A; Verbruggen, S W; McNamara, L M
2015-08-01
Load-induced fluid flow acts as an important biophysical signal for bone cell mechanotransduction in vivo, where the mechanical environment is thought to be monitored by integrin and primary cilia mechanoreceptors on the cell body. However, precisely how integrin- and primary cilia-based mechanosensors interact with the surrounding fluid flow stimulus and ultimately contribute to the biochemical response of bone cells within either the in vitro or in vivo environment remains poorly understood. In this study, we developed fluid-structure interaction models to characterise the deformation of integrin- and primary cilia-based mechanosensors in bone cells under fluid flow stimulation. Under in vitro fluid flow stimulation, these models predicted that integrin attachments on the cell-substrate interface were highly stimulated ε(eq) > 200,000 με, while the presence of a primary cilium on the cell also resulted in significant strain amplifications, arising at the ciliary base. As such, these mechanosensors likely play a role in mediating bone mechanotransduction in vitro. Under in vivo fluid flow stimulation, integrin attachments along the canalicular wall were highly stimulated and likely play a role in mediating cellular responses in vivo. The role of the primary cilium as a flow sensor in vivo depended upon its configuration within the lacunar cavity. Specifically, our results showed that a short free-standing primary cilium could not effectively fulfil a flow sensing role in vivo. However, a primary cilium that discretely attaches the lacunar wall can be highly stimulated, due to hydrodynamic pressure in the lacunocanalicular system and, as such, could play a role in mediating bone mechanotransduction in vivo. PMID:25399300
Qayyum, Mubashir; Khan, Hamid; Rahim, M Tariq; Ullah, Inayat
2015-01-01
The aim of this article is to model and analyze an unsteady axisymmetric flow of non-conducting, Newtonian fluid squeezed between two circular plates passing through porous medium channel with slip boundary condition. A single fourth order nonlinear ordinary differential equation is obtained using similarity transformation. The resulting boundary value problem is solved using Homotopy Perturbation Method (HPM) and fourth order Explicit Runge Kutta Method (RK4). Convergence of HPM solution is verified by obtaining various order approximate solutions along with absolute residuals. Validity of HPM solution is confirmed by comparing analytical and numerical solutions. Furthermore, the effects of various dimensionless parameters on the longitudinal and normal velocity profiles are studied graphically. PMID:25738864
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.
A critical evaluation of various turbulence models as applied to internal fluid flows
NASA Technical Reports Server (NTRS)
Nallasamy, M.
1985-01-01
Models employed in the computation of turbulent flows are described and their application to internal flows is evaluated by examining the predictions of various turbulence models in selected flow configurations. The main conclusions are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly; the algebraic stress model performs better in this case; and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well for fully developed flow, for developing flow, the algebraic stress model performance is not good; a Reynolds stress model should be used. False diffusion and inlet boundary conditions are discussed. Countergradient transport and its implications in turbulence modeling is mentioned. Two examples of recirculating flow predictions obtained using PHOENICS code are discussed. The vortex method, large eddy simulation (modeling of subgrid scale Reynolds stresses), and direct simulation, are considered. Some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.
Priyadharshini, S.; Ponalagusamy, R.
2015-01-01
An analysis of blood flow through a tapered artery with stenosis and dilatation has been carried out where the blood is treated as incompressible Herschel-Bulkley fluid. A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that the analytical expression for pressure gradient works well for the values of yield stress till 2.4. The wall shear stress and flow resistance increase significantly with axial distance and the increase is more in the case of converging tapered artery. A comparison study of velocity profiles, wall shear stress, and flow resistance for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids shows that the variation is greater for Herschel-Bulkley fluid than the other fluids. The obtained velocity profiles have been compared with the experimental data and it is observed that blood behaves like a Herschel-Bulkley fluid rather than power law, Bingham, and Newtonian fluids. It is observed that, in the case of a tapered stenosed tube, the streamline pattern follows a convex pattern when we move from r/R = 0 to r/R = 1 and it follows a concave pattern when we move from r/R = 0 to r/R = −1. Further, it is of opposite behaviour in the case of a tapered dilatation tube which forms new information that is, for the first time, added to the literature. PMID:27041979
Priyadharshini, S; Ponalagusamy, R
2015-01-01
An analysis of blood flow through a tapered artery with stenosis and dilatation has been carried out where the blood is treated as incompressible Herschel-Bulkley fluid. A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that the analytical expression for pressure gradient works well for the values of yield stress till 2.4. The wall shear stress and flow resistance increase significantly with axial distance and the increase is more in the case of converging tapered artery. A comparison study of velocity profiles, wall shear stress, and flow resistance for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids shows that the variation is greater for Herschel-Bulkley fluid than the other fluids. The obtained velocity profiles have been compared with the experimental data and it is observed that blood behaves like a Herschel-Bulkley fluid rather than power law, Bingham, and Newtonian fluids. It is observed that, in the case of a tapered stenosed tube, the streamline pattern follows a convex pattern when we move from r/R = 0 to r/R = 1 and it follows a concave pattern when we move from r/R = 0 to r/R = -1. Further, it is of opposite behaviour in the case of a tapered dilatation tube which forms new information that is, for the first time, added to the literature. PMID:27041979
Well-posedness and convergence of cfd two-fluid model for bubbly flows
NASA Astrophysics Data System (ADS)
Vaidheeswaran, Avinash
The current research is focused on developing a well-posed multidimensional CFD two-fluid model (TFM) for bubbly flows. Two-phase flows exhibit a wide range of local flow instabilities such as Kelvin-Helmholtz, Rayleigh-Taylor, plume and jet instabilities. They arise due to the density difference and/or the relative velocity between the two phases. A physically correct TFM is essential to model these instabilities. However, this is not the case with the TFMs in numerical codes, which can be shown to have complex eigenvalues due to incompleteness and hence are ill-posed as initial value problems. A common approach to regularize an incomplete TFM is to add artificial physics or numerically by using a coarse grid or first order methods. However, it eliminates the local physical instabilities along with the undesired high frequency oscillations resulting from the ill-posedness. Thus, the TFM loses the capability to predict the inherent local dynamics of the two-phase flow. The alternative approach followed in the current study is to introduce appropriate physical mechanisms that make the TFM well-posed. First a well-posed 1-D TFM for vertical bubbly flows is analyzed with characteristics, and dispersion analysis. When an incomplete TFM is used, it results in high frequency oscillations in the solution. It is demonstrated through the travelling void wave problem that, by adding the missing short wavelength physics to the numerical TFM, this can be removed by making the model well-posed. To extend the limit of well-posedness beyond the well-known TFM of Pauchon and Banerjee [1], the mechanism of collision is considered, and it is shown by characteristics analysis that the TFM then becomes well-posed for all void fractions of practical interest. The aforementioned ideas are then extended to CFD TFM. The travelling void wave problem is again used to demonstrate that by adding appropriate physics, the problem of ill-posedness is resolved. Furthermore, issues pertaining to
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.
Bhamidipati, J.R.; El-Kaddah, N.
1995-12-31
A mathematical model for the analysis and design of electromagnetic continuous casters (EMC) has been developed. The model addresses all aspects of EMC operation, namely, electromagnetic confinement of the molten pool, and melt solidification, all in terms of basic process variables. A methodology is presented for solving the coupled electromagnetic, heat transfer and fluid flow equations without regridding the solution domain during the search for the equilibrium meniscus shape, which is not known a priori. It is based on the solution of these equations in a fixed computational domain using a variable non-orthogonal coordinate transformation. Within this framework, the electromagnetic field was computed using the technique of mutual inductance, while the temperature and velocity fields in the slab were calculated using the control volume technique. The equilibrium shape of the meniscus was determined from the balance between the magnetic, metallostatic, and surface tension pressures. The model was used to simulate EMC casting of aluminum, and the numerical predictions were found to be in excellent agreement with experimental measurements reported in the open literature. The model presented in this paper provides a theoretical framework for the analysis and optimization of EMC casters.
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
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
Description of fluid dynamics and coupled transports in models of a laminar flow diffusion chamber.
Trávníčková, Tereza; Havlica, Jaromír; Ždímal, Vladimír
2013-08-14
The aim of this study is to assess how much the results of nucleation experiments in a laminar flow diffusion chamber (LFDC) are influenced by the complexity of the model of the transport properties. The effects of the type of fluid dynamic model (the steady state compressible Navier-Stokes system for an ideal gas/parabolic profile approximation) and the contributions of the coupled terms describing the Dufour effects and thermodiffusion on the predicted magnitude of the nucleation maxima and its location were investigated. This study was performed on the model of the homogeneous nucleation of an n-butanol-He vapor mixture in a LFDC. The isothermal dependencies of the nucleation rate on supersaturation were determined at three nucleation temperatures: 265 K, 270 K, and 280 K. For this purpose, the experimental LFDC data measured by A. P. Hyvärinen et al. [J. Chem. Phys. 124, 224304 (2006)] were reevaluated using transport models at different levels of complexity. Our results indicate that the type of fluid dynamical model affects both the position of the nucleation maxima in the LFDC and the maximum value of the nucleation rate. On the other hand, the Dufour effects and thermodiffusion perceptibly influence only the value of the maximal nucleation rate. Its position changes only marginally. The dependence of the maximum experimental nucleation rate on the saturation ratio and nucleation temperature was acquired for each case. Based on this dependence, we presented a method for the comparison and evaluation of the uncertainties of simpler models' solutions for the results, where we assumed that the model with Navier-Stokes equations and both coupled effects taken into account was the basis. From this comparison, it follows that an inappropriate choice of mathematical models could lead to relative errors of the order of several hundred percent in the maximum experimental nucleation rate. In the conclusion of this study, we also provide some general recommendations
A simple model of fluid flow and electrolyte balance in the body
NASA Technical Reports Server (NTRS)
White, R. J.; Neal, L.
1973-01-01
The model is basically a three-compartment model, the three compartments being the plasma, interstitial fluid and cellular fluid. Sodium, potassium, chloride and urea are the only major solutes considered explicitly. The control of body water and electrolyte distribution is affected via drinking and hormone levels. Basically, the model follows the effect of various oral input water loads on solute and water distribution throughout the body.
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.
Kenkeremath, D.
1985-05-01
Numerical simulation models and data bases that were developed for DOE as part of a number of geothermal programs have been assessed with respect to their overall stage of development and usefulness. This report combines three separate studies that focus attention upon: (1) economic models related to geothermal energy; (2) physical geothermal system models pertaining to thermal energy and the fluid medium; and (3) geothermal energy data bases. Computerized numerical models pertaining to the economics of extracting and utilizing geothermal energy have been summarized and catalogued with respect to their availability, utility and function. The 19 models that are discussed in detail were developed for use by geothermal operators, public utilities, and lending institutions who require a means to estimate the value of a given resource, total project costs, and the sensitivity of these values to specific variables. A number of the models are capable of economically assessing engineering aspects of geothermal projects. Computerized simulations of heat distribution and fluid flow have been assessed and are presented for ten models. Five of the models are identified as wellbore simulators and five are described as reservoir simulators. Each model is described in terms of its operational characteristics, input, output, and other pertinent attributes. Geothermal energy data bases are reviewed with respect to their current usefulness and availability. Summaries of eight data bases are provided in catalogue format, and an overall comparison of the elements of each data base is included.
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.
NASA Astrophysics Data System (ADS)
Martin, R. M.; Nicolas, A. N.
2003-04-01
A modeling approach of gas solid flow, taking into account different physical phenomena such as gas turbulence and inter-particle interactions is presented. Moment transport equations are derived for the second order fluctuating velocity tensor which allow to involve practical closures based on single phase turbulence modeling on one hand and kinetic theory of granular media on the other hand. The model is applied to fluid catalytic cracking processes and explosive volcanism. In the industry as well as in the geophysical community, multiphase flows are modeled using a finite volume approach and a multicorrector algorithm in time in order to determine implicitly the pressures, velocities and volume fractions for each phase. Pressures, and velocities are generally determined at mid-half mesh step from each other following the staggered grid approach. This ensures stability and prevents oscillations in pressure. It allows to treat almost all the Reynolds number ranges for all speeds and viscosities. The disadvantages appear when we want to treat more complex geometries or if a generalized curvilinear formulation of the conservation equations is considered. Too many interpolations have to be done and accuracy is then lost. In order to overcome these problems, we use here a similar algorithm in time and a Rhie and Chow interpolation (1983) of the collocated variables and essentially the velocities at the interface. The Rhie and Chow interpolation of the velocities at the finite volume interfaces allows to have no oscillations of the pressure without checkerboard effects and to stabilize all the algorithm. In a first predictor step, fluxes at the interfaces of the finite volumes are then computed using 2nd and 3rd order shock capturing schemes of MUSCL/TVD or Van Leer type, and the orthogonal stress components are treated implicitly while cross viscous/diffusion terms are treated explicitly. Pentadiagonal linear systems are solved in each geometrical direction (the so
Barker, Andrew T. Cai Xiaochuan
2010-02-01
We introduce and study numerically a scalable parallel finite element solver for the simulation of blood flow in compliant arteries. The incompressible Navier-Stokes equations are used to model the fluid and coupled to an incompressible linear elastic model for the blood vessel walls. Our method features an unstructured dynamic mesh capable of modeling complicated geometries, an arbitrary Lagrangian-Eulerian framework that allows for large displacements of the moving fluid domain, monolithic coupling between the fluid and structure equations, and fully implicit time discretization. Simulations based on blood vessel geometries derived from patient-specific clinical data are performed on large supercomputers using scalable Newton-Krylov algorithms preconditioned with an overlapping restricted additive Schwarz method that preconditions the entire fluid-structure system together. The algorithm is shown to be robust and scalable for a variety of physical parameters, scaling to hundreds of processors and millions of unknowns.
NASA Astrophysics Data System (ADS)
Tang, H. S.; Qu, K.; Wu, X. G.
2014-09-01
It is now becoming important to develop our capabilities to simulate coastal ocean flows involved with distinct physical phenomena occurring at a vast range of spatial and temporal scales. This paper presents a hybrid modeling system for such simulation. The system consists of a fully three dimensional (3D) fluid dynamics model and a geophysical fluid dynamics model, which couple with each other in two-way and march in time simultaneously. Particularly, in the hybrid system, the solver for incompressible flow on overset meshes (SIFOM) resolves fully 3D small-scale local flow phenomena, while the unstructured grid finite volume coastal ocean model (FVCOM) captures large-scale background flows. The integration of the two models are realized via domain decomposition implemented with an overset grid method. Numerical experiments on performance of the system in resolving flow patterns and solution convergence rate show that the SIFOM-FVCOM system works as intended, and its solutions compare reasonably with data obtained with measurements and other computational approaches. Its unparalleled capabilities to predict multiphysics and multiscale phenomena with high-fidelity are demonstrated by three typical applications that are beyond the reach of other currently existing models. It is anticipated that the SIFOM-FVCOM system will serve as a new platform to study many emerging coastal ocean problems.
NASA Astrophysics Data System (ADS)
Baumert, Brandon Max; Muller, Susan J.
1997-03-01
Flow visualization of two highly elastic, nonshear-thinning polyisobutylene/polybutene fluids in the gap between concentric cylinders was performed over a range of shear rates and choices of relative cylinder rotations. The observed secondary flows are discussed in terms of destabilizing elastic and centrifugal forces. In the more viscous, more elastic fluid, instabilities are found to be independent of the choice of rotating cylinder and due entirely to elasticity. At the lowest shear rates examined, the first detectable secondary flows are steady counter-rotating vortices forming after a shearing time more than five orders of magnitude greater than the characteristic relaxation time of the fluid. At somewhat higher shear rates, a much more rapidly appearing oscillatory flow is observed to evolve into the steady vortex structure. In the less elastic fluid, the structure first detectable at the lowest shear rates is again steady vortices regardless of the choice of driving cylinder. At all shear rates examined, only elastic stationary vortices are observed in the absence of centrifugal destabilization (outer cylinder rotating). Secondary flows are significantly stronger in the presence of the centrifugal destabilization due to a rotating inner cylinder. Interaction of elasticity and centrifugal forces is found to generate a number of axially translating vortex structures, many of which are described here for the first time. At a shear rate more than five times the critical, another family of instability is observed which closely resembles a purely elastic instability observed by Baumert and Muller (1995). These experimental results are expected to provide a challenging test of numerical simulations of these viscoelastic flows.
Baumert, B.M.; Muller, S.J.
1997-03-01
Flow visualization of two highly elastic, nonshear-thinning polyisobutylene/polybutene fluids in the gap between concentric cylinders was performed over a range of shear rates and choices of relative cylinder rotations. The observed secondary flows are discussed in terms of destabilizing elastic and centrifugal forces. In the more viscous, more elastic fluid, instabilities are found to be independent of the choice of rotating cylinder and due entirely to elasticity. At the lowest shear rates examined, the first detectable secondary flows are steady counter-rotating vortices forming after a shearing time more than five orders of magnitude greater than the characteristic relaxation time of the fluid. At somewhat higher shear rates, a much more rapidly appearing oscillatory flow is observed to evolve into the steady vortex structure. In the less elastic fluid, the structure first detectable at the lowest shear rates is again steady vortices regardless of the choice of driving cylinder. At all shear rates examined, only elastic stationary vortices are observed in the absence of centrifugal destabilization (outer cylinder rotating). Secondary flows are significantly stronger in the presence of the centrifugal destabilization due to a rotating inner cylinder. Interaction of elasticity and centrifugal forces is found to generate a number of axially translating vortex structures, many of which are described here for the first time. At a shear rate more than five times the critical, another family of instability is observed which closely resembles a purely elastic instability observed by Baumert and Muller (1995). These experimental results are expected to provide a challenging test of numerical simulations of these viscoelastic flows. {copyright} {ital 1997 American Institute of Physics.}
Pruess, K.; Karasaki, K.
1982-10-01
Conventional approaches to geothermal reservoir modeling have employed a porous medium approximation, but recently methods have been developed which can take into account the different thermodynamic conditions in rock matrix and fractures. The multiple interacting continua method (MINC) treats the thermal and hydraulic interaction between rock matrix and fractures in terms of a set of geometrical parameters. However, this approach was restricted to idealized fracture distributions with regularly shaped matrix blocks. Fractures in geothermal reservoirs usually occur in nearly parallel sets with a certain scatter in orientation, and a stochastic distribution of spacings and apertures. The MINC-method was extended to realistic fracture systems with stochastic distributions. The interaction between matrix and fractures is parameterized in terms of a proximity function, which represents the volume of matrix rock as a function of distance from the fractures. Monte Carlo techniques were employed to compute proximity functions for a number of two-dimensional systems with regular or stochastic fracture distributions. It is shown how the proximity functions can be used to generate computational grids for modeling fluid and heat flow in fractured reservoirs.
Guidelines for optimizing multilevel ECN using fluid-flow-based TCP model
NASA Astrophysics Data System (ADS)
Quet, Pierre-Francois; Chellappan, Sriram; Durresi, Arjan; Sridharan, Mukundan; Ozbay, Hitay; Jain, Raj
2002-07-01
Congestion avoidance on today's Internet is mainly provided by the combination of the TCP protocol and Active Queue Management (AQM) schemes such as the de facto standard RED (Random Early Detection). When used with ECN (Explicit Congestion Notification), these algorithms can be modeled as a feedback control system in which the feedback information is carried on a single bit. A modification of this scheme called MECN was proposed, where the marking information is carried using 2 bits. MECN conveys more accurate feedback about the network congestion to the source than the current 1-bit ECN. The TCP source reaction was modified so that it takes advantage of the extra information about congestion and adapts faster to the changing congestion scenario leading to a smoother decrease in the sending rates of the sources upon congestion detection and consequently resulting in an increase in the router's throughput. A linearized fluid flow model already developed for ECN is extended to our case. Using control theoretic tools we justify the performance obtained in using the MECN scheme and give guidelines for optimizing its parameters. We use ns simulations to illustrate the performance improvement from the point of better throughput and low level of oscillations in the queue.
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.
Strongly coupled fluid-particle flows in vertical channels. II. Turbulence modeling
NASA Astrophysics Data System (ADS)
Capecelatro, Jesse; Desjardins, Olivier; Fox, Rodney O.
2016-03-01
In Part I, simulations of strongly coupled fluid-particle flow in a vertical channel were performed with the purpose of understanding, in general, the fundamental physics of wall-bounded multiphase turbulence and, in particular, the roles of the spatially correlated and uncorrelated components of the particle velocity. The exact Reynolds-averaged (RA) equations for high-mass-loading suspensions were presented, and the unclosed terms that are retained in the context of fully developed channel flow were evaluated in an Eulerian-Lagrangian (EL) framework. Here, data from the EL simulations are used to validate a multiphase Reynolds-stress model (RSM) that predicts the wall-normal distribution of the two-phase, one-point turbulence statistics up to second order. It is shown that the anisotropy of the Reynolds stresses both near the wall and far away is a crucial component for predicting the distribution of the RA particle-phase volume fraction. Moreover, the decomposition of the phase-average (PA) particle-phase fluctuating energy into the spatially correlated and uncorrelated components is necessary to account for the boundary conditions at the wall. When these factors are properly accounted for in the RSM, the agreement with the EL turbulence statistics is satisfactory at first order (e.g., PA velocities) but less so at second order (e.g., PA turbulent kinetic energy). Finally, an algebraic stress model for the PA particle-phase pressure tensor and the Reynolds stresses is derived from the RSM using the weak-equilibrium assumption.
Quaini, A.; Canic, S.; Glowinski, R.; Igo, S.; Hartley, C.J.; Zoghbi, W.; Little, S.
2011-01-01
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
Quaini, A; Canic, S; Glowinski, R; Igo, S; Hartley, C J; Zoghbi, W; Little, S
2012-01-10
This work presents a validation of a fluid-structure interaction computational model simulating the flow conditions in an in vitro mock heart chamber modeling mitral valve regurgitation during the ejection phase during which the trans-valvular pressure drop and valve displacement are not as large. The mock heart chamber was developed to study the use of 2D and 3D color Doppler techniques in imaging the clinically relevant complex intra-cardiac flow events associated with mitral regurgitation. Computational models are expected to play an important role in supporting, refining, and reinforcing the emerging 3D echocardiographic applications. We have developed a 3D computational fluid-structure interaction algorithm based on a semi-implicit, monolithic method, combined with an arbitrary Lagrangian-Eulerian approach to capture the fluid domain motion. The mock regurgitant mitral valve corresponding to an elastic plate with a geometric orifice, was modeled using 3D elasticity, while the blood flow was modeled using the 3D Navier-Stokes equations for an incompressible, viscous fluid. The two are coupled via the kinematic and dynamic conditions describing the two-way coupling. The pressure, the flow rate, and orifice plate displacement were measured and compared with numerical simulation results. In-line flow meter was used to measure the flow, pressure transducers were used to measure the pressure, and a Doppler method developed by one of the authors was used to measure the axial displacement of the orifice plate. The maximum recorded difference between experiment and numerical simulation for the flow rate was 4%, the pressure 3.6%, and for the orifice displacement 15%, showing excellent agreement between the two. PMID:22138194
An adaptive lattice Boltzmann scheme for modeling two-fluid-phase flow in porous medium systems
NASA Astrophysics Data System (ADS)
Dye, Amanda L.; McClure, James E.; Adalsteinsson, David; Miller, Cass T.
2016-04-01
We formulate a multiple-relaxation-time (MRT) lattice-Boltzmann method (LBM) to simulate two-fluid-phase flow in porous medium systems. The MRT LBM is applied to simulate the displacement of a wetting fluid by a nonwetting fluid in a system corresponding to a microfluidic cell. Analysis of the simulation shows widely varying time scales for the dynamics of fluid pressures, fluid saturations, and interfacial curvatures that are typical characteristics of such systems. Displacement phenomena include Haines jumps, which are relatively short duration isolated events of rapid fluid displacement driven by capillary instability. An adaptive algorithm is advanced using a level-set method to locate interfaces and estimate their rate of advancement. Because the displacement dynamics are confined to the interfacial regions for a majority of the relaxation time, the computational effort is focused on these regions. The proposed algorithm is shown to reduce computational effort by an order of magnitude, while yielding essentially identical solutions to a conventional fully coupled approach. The challenges posed by Haines jumps are also resolved by the adaptive algorithm. Possible extensions to the advanced method are discussed.
Reactive fluid flow models and applications to diagenesis, mineral deposits and crustal rocks
Lasaga, A.C.; Rye, D.M.
1993-08-01
Funds are requested for a combined theoretical and field study of coupled fluid flow, heat and mass transport, and chemical reaction in hydrothermal and metamorphic systems. An existing computer code developed by the applicants which numerically treats multi-component, finite-rate reactions combined with advective and dispersive transport in one and two dimensions and which incorporates isotopic exchange and heat and mass transfer will continue to be developed and applied in a variety of geological settings. The code we have developed simultaneously solves for mass transport and reaction, thus offering a significant improvement in computational efficiency over existing ``batch`` reaction path codes. By coupling flow and chemical reaction in a hydrothermal system, we can explicitly investigate the extent to which characteristic flow-reaction paths govern the chemical evolution of the fluids in a hydrothermal system. The concept of a flow-reaction path is particularly important where certain portions of mature hydrothermal systems may exhaust the buffer capacity of the rock as the primary mineralogy is consumed. In these instances 7 fluids traversing distinct regions within the hydrothermal system may experience very different reaction histories, even where the system can be described as nearly isothermal. The study of paleo-hydrothermal systems can yield some important insights into the chemical dynamics of hydrothermal systems in general. As an example of a paleo-hydrothermal system, we have considered the geochemical evolution of ``porphyry-copper`` type mineralization.
NASA Astrophysics Data System (ADS)
Čanić, Sunčica; Mikelić, Andro; Tambača, Josip
2005-12-01
We derive a closed system of effective equations describing a time-dependent flow of a viscous incompressible Newtonian fluid through a long and narrow elastic tube. The 3D axially symmetric incompressible Navier-Stokes equations are used to model the flow. Two models are used to describe the tube wall: the linear membrane shell model and the linearly elastic membrane and the curved, linearly elastic Koiter shell model. We study the behavior of the coupled fluid-structure interaction problem in the limit when the ratio between the radius and the length of the tube, ɛ, tends to zero. We obtain the reduced equations that are of Biot type with memory. An interesting feature of the reduced equations is that the memory term explicitly captures the viscoelastic nature of the coupled problem. Our model provides significant improvement over the standard 1D approximations of the fluid-structure interaction problem, all of which assume an ad hoc closure assumption for the velocity profile. We performed experimental validation of the reduced model using a mock circulatory flow loop assembled at the Cardiovascular Research Laboratory at the Texas Heart Institute. Experimental results show excellent agreement with the numerically calculated solution. Major applications include blood flow through large human arteries. To cite this article: S. Čanić et al., C. R. Mecanique 333 (2005).
Fluid flow modeling at the Lusi mud eruption, East java, Indonesia.
NASA Astrophysics Data System (ADS)
Collignon, Marine; Schmid, Daniel; Mazzini, Adriano
2016-04-01
The 29th of may 2006, gas water and mud breccia started to erupt at several localities along the Watukosek fault system, in the Sidoarjo Regency in East java, Indonesia. The most prominent eruption, named Lusi, is still active and covering a surface of nearly 7 km2, resulting in the displacement of ~ 30 000 people. Although the origin and the chemical composition of the erupted fluids have been documented, the mechanical and physical properties of the mud are poorly constrained, and many aspects still remain not understood. Very little is known about the internal dynamics of the Lusi conduit(s). In this study, conducted in the framework of the Lusi Lab project (ERC grant n°308126) we use both analytical and numerical methods to better understand the flow dynamics within the main conduit and to try to explain the longevity of the edifice. The 2D numerical model considers a vertical conduit with a reservoir at its base and solves the stokes equations, discretized on a finite element mesh. Although, three phases (solid, liquid and gas) are present in nature, we only consider the liquid phase. The solid phase is treated as rigid particles in suspension in the liquid. The gaseous phase (methane and carbon dioxide) is treated in an analytical manner using the equations of state of the H2O-CO2 and H2O-CH4 systems. Here, we discuss the effects of density, viscosity, gas concentration and clasts concentration and size on the dynamics of the flow in the conduit as well as implications of the conduit stability.
Mathematical model for blood flow through a bifurcated artery using couple stress fluid.
Srinivasacharya, D; Madhava Rao, G
2016-08-01
In this article, the blood flow through a bifurcated artery with mild stenosis is investigated taking blood as couple stress fluid. The artery configuring bifurcation is assumed to be symmetric about the axis of the artery and straight cylinders of finite length. The governing equations are non-dimensionalized and coordinate transformation is used to convert the irregular boundary to a regular boundary. The resulting system of equations is solved numerically using the finite difference method. The variation of shear stress, flow rate and impedance near the apex with pertinent parameters are studied graphically. It has been noticed that shear stress, flow rate and impedance have been changing suddenly with all the parameters on both sides of the apex. This occurs because of the backflow of the streaming blood at the onset of the lateral junction and secondary flow near the apex in the daughter artery. PMID:27235925
Fluid flow modeling of resin transfer molding for composite material wind turbine blade structures.
Cairns, Douglas S.; Rossel, Scott M.
2004-06-01
Resin transfer molding (RTM) is a closed mold process for making composite materials. It has the potential to produce parts more cost effectively than hand lay-up or other methods. However, fluid flow tends to be unpredictable and parts the size of a wind turbine blade are difficult to engineer without some predictive method for resin flow. There were five goals of this study. The first was to determine permeabilities for three fabrics commonly used for RTM over a useful range of fiber volume fractions. Next, relations to estimate permeabilities in mixed fabric lay-ups were evaluated. Flow in blade substructures was analyzed and compared to predictions. Flow in a full-scale blade was predicted and substructure results were used to validate the accuracy of a full-scale blade prediction.
Improved PISO algorithms for modeling density varying flow in conjugate fluid-porous domains
NASA Astrophysics Data System (ADS)
Nordlund, M.; Stanic, M.; Kuczaj, A. K.; Frederix, E. M. A.; Geurts, B. J.
2016-02-01
Two modified segregated PISO algorithms are proposed, which are constructed to avoid the development of spurious oscillations in the computed flow near sharp interfaces of conjugate fluid-porous domains. The new collocated finite volume algorithms modify the Rhie-Chow interpolation to maintain a correct pressure-velocity coupling when large discontinuous momentum sources associated with jumps in the local permeability and porosity are present. The Re-Distributed Resistivity (RDR) algorithm is based on spreading flow resistivity over the grid cells neighboring a discontinuity in material properties of the porous medium. The Face Consistent Pressure (FCP) approach derives an auxiliary pressure value at the fluid-porous interface using momentum balance around the interface. Such derived pressure correction is designed to avoid spurious oscillations as would otherwise arise with a strictly central discretization. The proposed algorithms are successfully compared against published data for the velocity and pressure for two reference cases of viscous flow. The robustness of the proposed algorithms is additionally demonstrated for strongly reduced viscosity, i.e., higher Reynolds number flows and low Darcy numbers, i.e., low permeability of the porous regions in the domain, for which solutions without unphysical oscillations are computed. Both RDR and FCP are found to accurately represent porous media flow near discontinuities in material properties on structured grids.
Inertial Coupling in Two-Phase Flow: a Few Test Cases and Their Impacts on Two-Fluid Modeling.
NASA Astrophysics Data System (ADS)
Cai, Xiaolong
Two-phase flow makes up about one-half of all industrial, biological and environmental fluid mechanics. Analytical methods that are presently useful in these are still mostly ad hoc. Unlike in (single phase) fluid mechanics, there are no universal macroscopic fundamental equations existing yet, although the Navier-Stokes equations describe the microscopic details of the motion of each phase. One promising recent theoretical development in two-phase flow is the two-fluid theory. The modeling of the interfacial interaction force terms are very crucial in this theory since they govern the interfacial momentum transfer and couple the two phases together. One of the interaction forces is the inertial coupling force, or the added (virtual) mass force, which plays a significant role in deciding the well -posedness of the model, especially when dynamic interactions (relative acceleration) between phases are important (Drew, 1983). Attention is focused on the inertial coupling aspect of two-phase modeling in this thesis, following Wallis (1989). Although the concept is classic, there is still much to be known about added mass effects in two -phase flow. The boundary (be it fluid boundary or other surfaces) effect on the added mass coefficient of a sphere is studied, resulting in several new, essentially exact solutions which are valuable contributions to the understanding of hydrodynamic interaction in potential flow. The added masses for isotropic and nonisotropic dispersions of spheres in an unbounded flow are studied. Among the new results are: (i) the "conditional" convergence experienced is physically meaningful; (ii) Geurst's conjecture (1985) is correct, contrary to Smereka & Milton's (1991) finding; (iii) predictions from the "tube model" proposed in Cai & Wallis (1992) are further validated by predictions from the image method. Some new phenomena about the added mass of dispersed flow in finite pipes with different end conditions are also studied. Experimental
On the multidimensional modeling of fluid flow and heat transfer in SCWRS
Gallaway, T.; Antal, S. P.; Podowski, M. Z.
2012-07-01
The Supercritical Water Reactor (SCWR) has been proposed as one of the six Generation IV reactor design concepts under consideration. The key feature of the SCWR is that water at supercritical pressures is used as the reactor coolant. Although at such pressures, fluids do not undergo phase change as they are heated, the fluid properties experience dramatic variations throughout what is known as the pseudo-critical region. Highly nonuniform temperature and fluid property distributions are expected in the reactor core, which will have a significant impact on turbulence and heat transfer in future SCWRs. The goal of the present work has been to understand and predict the effects of these fluid property variations on turbulence and heat transfer throughout the reactor core. Spline-type property models have been formulated for water at supercritical pressures in order to include the dependence of properties on both temperature and pressure into a numerical solver. New models of turbulence and heat transfer for variable-property fluids have been developed and implemented into the NPHASE-CMFD software. The results for these models have been compared to experimental data from the Korea Atomic Energy Research Inst. (KAERI) for various heat transfer regimes. It is found that the Low-Reynolds {kappa}-{epsilon} model performs best at predicting the experimental data. (authors)
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.
Modeling on Fluid Flow and Inclusion Motion in Centrifugal Continuous Casting Strands
NASA Astrophysics Data System (ADS)
Wang, Qiangqiang; Zhang, Lifeng; Sridhar, Seetharaman
2016-05-01
During the centrifugal continuous casting process, unreasonable casting parameters can cause violent level fluctuation, serious gas entrainment, and formation of frozen shell pieces at the meniscus. Thus, in the current study, a three-dimensional multiphase turbulent model was established to study the transport phenomena during centrifugal continuous casting process. The effects of nozzle position, casting and rotational speed on the flow pattern, centrifugal force acting on the molten steel, level fluctuation, gas entrainment, shear stress on mold wall, and motion of inclusions during centrifugal continuous casting process were investigated. Volume of Fluid model was used to simulate the molten steel-air two-phase. The level fluctuation and the gas entrainment during casting were calculated by user-developed subroutines. The trajectory of inclusions in the rotating system was calculated using the Lagrangian approach. The results show that during centrifugal continuous casting, a large amount of gas was entrained into the molten steel, and broken into bubbles of various sizes. The greater the distance to the mold wall, the smaller the centrifugal force. Rotation speed had the most important influence on the centrifugal force distribution at the side region. Angular moving angle of the nozzle with 8° and keeping the rotation speed with 60 revolutions per minute can somehow stabilize the level fluctuation. The increase of angular angle of nozzle from 8 to 18 deg and rotation speed from 40 to 80 revolutions per minute favored to decrease the total volume of entrained bubbles, while the increase of distance of nozzle moving left and casting speed had reverse effects. The trajectories of inclusions in the mold were irregular, and then rotated along the strand length. After penetrating a certain distance, the inclusions gradually moved to the center of billet and gathered there. More work, such as the heat transfer, the solidification, and the inclusions entrapment
Modeling on Fluid Flow and Inclusion Motion in Centrifugal Continuous Casting Strands
NASA Astrophysics Data System (ADS)
Wang, Qiangqiang; Zhang, Lifeng; Sridhar, Seetharaman
2016-08-01
During the centrifugal continuous casting process, unreasonable casting parameters can cause violent level fluctuation, serious gas entrainment, and formation of frozen shell pieces at the meniscus. Thus, in the current study, a three-dimensional multiphase turbulent model was established to study the transport phenomena during centrifugal continuous casting process. The effects of nozzle position, casting and rotational speed on the flow pattern, centrifugal force acting on the molten steel, level fluctuation, gas entrainment, shear stress on mold wall, and motion of inclusions during centrifugal continuous casting process were investigated. Volume of Fluid model was used to simulate the molten steel-air two-phase. The level fluctuation and the gas entrainment during casting were calculated by user-developed subroutines. The trajectory of inclusions in the rotating system was calculated using the Lagrangian approach. The results show that during centrifugal continuous casting, a large amount of gas was entrained into the molten steel, and broken into bubbles of various sizes. The greater the distance to the mold wall, the smaller the centrifugal force. Rotation speed had the most important influence on the centrifugal force distribution at the side region. Angular moving angle of the nozzle with 8° and keeping the rotation speed with 60 revolutions per minute can somehow stabilize the level fluctuation. The increase of angular angle of nozzle from 8 to 18 deg and rotation speed from 40 to 80 revolutions per minute favored to decrease the total volume of entrained bubbles, while the increase of distance of nozzle moving left and casting speed had reverse effects. The trajectories of inclusions in the mold were irregular, and then rotated along the strand length. After penetrating a certain distance, the inclusions gradually moved to the center of billet and gathered there. More work, such as the heat transfer, the solidification, and the inclusions entrapment
NASA Astrophysics Data System (ADS)
Wagner, H. P.; Kaeppeler, H. J.; Auweter-Kurtz, M.
1998-03-01
The investigation of magnetoplasmadynamic (MPD) instabilities in connection with the performance of MPD thrusters has received increased attention for several years. A systematic investigation of macroscopic drift and gradient driven instabilities at the Institut für Raumfahrtsysteme (IRS), carried out since 1987, showed that the occurrence of acoustic wave and space charge instabilities was possible. The calculation of the linear development of instabilities using a three-fluid theory of the plasma flow will be presented. Furthermore, in order to identify the instabilities found in the three-fluid formalism, two-fluid and one-fluid models are investigated and links established. When possible, analytical approximations for the unstable roots of the linear dispersion relation are derived, giving the explicit dependence of the oscillation frequencies and growth rates on the wavevector and the various plasma parameters.
Modelling of fluid flow and heat transfer in a reciprocating compressor
NASA Astrophysics Data System (ADS)
Tuhovcak, J.; Hejcik, J.; Jicha, M.
2015-08-01
Efficiency of reciprocating compressor is strongly dependent on several parameters. The most important are valve behaviour and heat transfer. Valves affect the flow through the suction and discharge line. Heat flow from the walls to working fluid increases the work of the cycle. Understanding of these phenomena inside the compressor is a necessary step in the development process. Commercial CFD tools offer wide range of opportunities how to simulate the flow inside the reciprocating compressor nowadays, however they are too demanding in terms of computational time and mesh creation. Several approaches using various correlation equation exist to describe the heat transfer inside the cylinder, however none of them was validated by measurements due to the complicated settings. The goal of this paper is to show a comparison between these correlations using in-house code based on energy balance through the cycle.
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 Technical Reports Server (NTRS)
Sharp, John R.; Kittredge, Ken; Schunk, Richard G.
2003-01-01
As part of the aero-thermodynamics team supporting the Columbia Accident Investigation Board (CAB), the Marshall Space Flight Center was asked to perform engineering analyses of internal flows in the port wing. The aero-thermodynamics team was split into internal flow and external flow teams with the support being divided between shorter timeframe engineering methods and more complex computational fluid dynamics. In order to gain a rough order of magnitude type of knowledge of the internal flow in the port wing for various breach locations and sizes (as theorized by the CAB to have caused the Columbia re-entry failure), a bulk venting model was required to input boundary flow rates and pressures to the computational fluid dynamics (CFD) analyses. This paper summarizes the modeling that was done by MSFC in Thermal Desktop. A venting model of the entire Orbiter was constructed in FloCAD based on Rockwell International s flight substantiation analyses and the STS-107 reentry trajectory. Chemical equilibrium air thermodynamic properties were generated for SINDA/FLUINT s fluid property routines from a code provided by Langley Research Center. In parallel, a simplified thermal mathematical model of the port wing, including the Thermal Protection System (TPS), was based on more detailed Shuttle re-entry modeling previously done by the Dryden Flight Research Center. Once the venting model was coupled with the thermal model of the wing structure with chemical equilibrium air properties, various breach scenarios were assessed in support of the aero-thermodynamics team. The construction of the coupled model and results are presented herein.
Wetzel, L.R.; Raffensperger, J.P.; Shock, E.L.
2001-01-01
Coordinated geochemical and hydrological calculations guide our understanding of the composition, fluid flow patterns, and thermal structure of near-ridge oceanic crust. The case study presented here illustrates geochemical and thermal changes taking place as oceanic crust ages from 0.2 to 1.0 Myr. Using a finite element code, we model fluid flow and heat transport through the upper few hundred meters of an abyssal hill created at an intermediate spreading rate. We use a reaction path model with a customized database to calculate equilibrium fluid compositions and mineral assemblages of basalt and seawater at 500 bars and temperatures ranging from 150 to 400??C. In one scenario, reaction path calculations suggest that volume increases on the order of 10% may occur within portions of the basaltic basement. If this change in volume occurred, it would be sufficient to fill all primary porosity in some locations, effectively sealing off portions of the oceanic crust. Thermal profiles resulting from fluid flow simulations indicate that volume changes along this possible reaction path occur primarily within the first 0.4 Myr of crustal aging. ?? 2001 Elsevier Science B.V. All rights reserved.
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
Numerical investigation of a space-fractional model of turbulent fluid flow in rectangular ducts
NASA Astrophysics Data System (ADS)
Churbanov, Alexander G.; Vabishchevich, Petr N.
2016-09-01
Models containing fractional derivatives are among the most promising new approaches for description of turbulent flows. In the present work, a steady-state flow in a duct is considered under the condition that the turbulent diffusion is governed by a fractional power of the Laplace operator. To study numerically flows in rectangular channels, finite-difference approximations are employed. The resulting discrete problem is solved by a preconditioned conjugate gradient method. At each iteration, the problem with a fractional power of the grid Laplace operator is solved. Predictions of turbulent flows in ducts at different Reynolds numbers are presented via mean velocity fields.
Modeling of Fluid Flow and Heat Transfer in Nanotube and Nanowire Forests
NASA Astrophysics Data System (ADS)
Martin, Michael
2010-10-01
Bundles of nanotubes, also known as nanotube forests, are under consideration for applications such as chip cooling and pre-concentrators for biodetection. Scaling law analysis shows that the air flow through these forests at atmospheric pressure is in the free-molecular flow regime. Based on the linearized free-molecular flow equations, a model is presented for the pressure drop and heat transfer in these systems. The momentum and energy equations are coupled, requiring that they be solved simultaneously. Results show large pressure drops, and a non-linear pressure distribution, similar to that seen in rarefied micro-channel flows.
Dynamic simulation of wavy-stratified two-phase flow with the one-dimensional two-fluid model
NASA Astrophysics Data System (ADS)
Fullmer, William D.
The one-dimensional two-fluid model is the basis for the description of the transport of mass, momentum and energy in the thermal-hydraulic codes used for nuclear reactor safety analysis. Unlike other physical transport models, the one-dimensional two-fluid model suffers from the possibility of being ill-posed as an initial-boundary value problem depending on the flow conditions and the relevant physical closure laws. Typically, the ill-posedness is dealt with through either excessive numerical damping or the addition of unphysical closure laws designed for the sole purpose of hyperbolization. Unfortunately both methods eliminate the instability along with the problem of ill-posedness causing the model to undoubtedly lose some of its inherent dynamic capability. In this work, a one-dimensional two-fluid model for horizontal or slightly inclined stratified flow is developed. Higher order physical models that are often neglected, such as surface tension and axial viscous stress, are retained for their short-wavelength stability properties. Characteristic, dispersion and nonlinear analyses are performed to demonstrate that the resulting model is linearly well-posed and nonlinearly well-behaved. While it has been known that higher-order differential terms are able to regularize the short-wavelength problem of ill-posedness without removing the long-wavelength instability, the literature is relatively silent on the consequences of using a model under linearly unstable conditions. Using carefully selected conditions in an idealized infinite domain, it is demonstrated for the first time that the one-dimensional two-fluid model exhibits chaotic behavior in addition to limit cycles and asymptotic stability. The chaotic behavior is a consequence of the long-wavelength linear instability (energy source) the nonlinearity (energy transfer) and the short-wavelength dissipation (energy sink). Since the model is chaotic, solutions exhibit a sensitive dependence on initial
NASA Astrophysics Data System (ADS)
Mushtaq, A.; Abbasbandy, S.; Mustafa, M.; Hayat, T.; Alsaedi, A.
2016-01-01
Present work studies the well-known Sakiadis flow of Maxwell fluid along a moving plate in a calm fluid by considering the Cattaneo-Christov heat flux model. This recently developed model has the tendency to describe the characteristics of relaxation time for heat flux. Some numerical local similarity solutions of the associated problem are computed by two approaches namely (i) the shooting method and (ii) the Keller-box method. The solution is dependent on some interesting parameters which include the viscoelastic fluid parameter β, the dimensionless thermal relaxation time γ and the Prandtl number Pr. Our simulations indicate that variation in the temperature distribution with an increase in local Deborah number γ is non-monotonic. The results for the Fourier's heat conduction law can be obtained as special cases of the present study.
A Mathematical Model for the Flow of a Casson Fluid due to Metachronal Beating of Cilia in a Tube
Siddiqui, A. M.; Farooq, A. A.; Rana, M. A.
2015-01-01
A mathematical model is developed to study the transport mechanism of a Casson fluid flow inspired by the metachronal coordination between the beating cilia in a cylindrical tube. A two-dimensional system of nonlinear equations governing the flow problem is formulated by using axisymmetric cylindrical coordinates and then simplified by employing the long wavelength and low Reynolds number assumptions. Exact solutions are derived for the velocity components, the axial pressure gradient, and the stream function. However, the expressions for the pressure rise and the volume flow rate are evaluated numerically. The features of the flow characteristics such as pumping and trapping are illustrated and discussed with the help of graphs. It is observed that the volume flow rate is influenced significantly by the width of plug flow region Hp as well as the cilia length parameter ε. The analysis is also applied and compared with the estimated value of the volume flow rate of epididymal fluid in the ductus efferentes of the human male reproductive tract. PMID:25789334
A mathematical model for the flow of a Casson fluid due to metachronal beating of cilia in a tube.
Siddiqui, A M; Farooq, A A; Rana, M A
2015-01-01
A mathematical model is developed to study the transport mechanism of a Casson fluid flow inspired by the metachronal coordination between the beating cilia in a cylindrical tube. A two-dimensional system of nonlinear equations governing the flow problem is formulated by using axisymmetric cylindrical coordinates and then simplified by employing the long wavelength and low Reynolds number assumptions. Exact solutions are derived for the velocity components, the axial pressure gradient, and the stream function. However, the expressions for the pressure rise and the volume flow rate are evaluated numerically. The features of the flow characteristics such as pumping and trapping are illustrated and discussed with the help of graphs. It is observed that the volume flow rate is influenced significantly by the width of plug flow region H p as well as the cilia length parameter ε. The analysis is also applied and compared with the estimated value of the volume flow rate of epididymal fluid in the ductus efferentes of the human male reproductive tract. PMID:25789334
NASA Astrophysics Data System (ADS)
Chen, Zheng-Shou; Kim, Wu-Joan
2012-03-01
This article presents a numerical investigation concerning the effect of two kinds of axially progressing internal flows (namely, upward and downward) on fluid-structure interaction (FSI) dynamics about a marine riser model which is subject to external shear current. The CAE technology behind the current research is a proposed FSI solution, which combines structural analysis software with CFD technology together. Efficiency validation for the CFD software was carried out first. It has been proved that the result from numerical simulations agrees well with the observation from relating model test cases in which the fluidity of internal flow is ignorable. After verifying the numerical code accuracy, simulations are conducted to study the vibration response that attributes to the internal progressive flow. It is found that the existence of internal flow does play an important role in determining the vibration mode (/dominant frequency) and the magnitude of instantaneous vibration amplitude. Since asymmetric curvature along the riser span emerges in the case of external shear current, the centrifugal and Coriolis accelerations owing to up- and downward internal progressive flows play different roles in determining the fluid-structure interaction response. The discrepancy between them becomes distinct, when the velocity ratio of internal flow against external shear current is relatively high.
NASA Astrophysics Data System (ADS)
Konar, D.
2015-09-01
The uses of ejector for efficient refrigeration are manifold - it has been used, among other applications, in the VCRS to reduce the compression ratio, in the combined ejector absorption cycle to enhance the refrigeration capacity and in the ejector absorber cycle to obtain lower evaporator pressures, higher absorber pressures and pre-absorption of the refrigerant in the ejector. Hence, modeling of flow which may be two phase two fluid as in ejector absorber cycle or two phase single fluid as in VCRS in an ejector assumes utmost importance. However, much work has not been done in this field. The primary objective of the present work is to discuss about the role of ejectors in various refrigeration systems and to model the two phase two fluid flow in the nozzle and the diffuser of an ejector under suitable assumptions. The equations of conservation of mass, momentum and energy have been solved to find the different flow properties like pressure, temperature and velocity of the two phases as function of the length in the diffuser. Different cases pertaining to different flows have been taken care of by appreciating what type of phenomena can actually occur at the interface of the two phases. Higher pressure rise was obtained for a given diffuser length with higher diffuser angles, smaller droplet diameter, higher inlet velocity of the gaseous phase and higher drag coefficients. Among other results, it was also seen that the two phases reached thermal equilibrium faster with higher diffuser angle, smaller droplet diameter and higher heat transfer coefficient.
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.
A model for fluid flow during saturated boiling on a horizontal cylinder
NASA Technical Reports Server (NTRS)
Kheyrandish, K.; Dalton, C.; Lienhard, J. H.
1987-01-01
A model has been developed to represent the vapor removal pattern in the vicinity of a cylinder during nucleate flow boiling across a horizontal cylinder. The model is based on a potential flow representation of the liquid and vapor regions and an estimate of the losses that should occur in the flow. Correlation of the losses shows a weak dependence on the Weber number and a slightly stronger dependence on the saturated liquid-to-vapor density ratio. The vapor jet thickness, which is crucial to the prediction of the burnout heat flux, and the shape of the vapor film are predicted. Both are verified by qualitative experimental observations.
Computational fluid dynamics modeling of two-phase flow in a BWR fuel assembly. Final CRADA Report.
Tentner, A.; Nuclear Engineering Division
2009-10-13
A direct numerical simulation capability for two-phase flows with heat transfer in complex geometries can considerably reduce the hardware development cycle, facilitate the optimization and reduce the costs of testing of various industrial facilities, such as nuclear power plants, steam generators, steam condensers, liquid cooling systems, heat exchangers, distillers, and boilers. Specifically, the phenomena occurring in a two-phase coolant flow in a BWR (Boiling Water Reactor) fuel assembly include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. Traditionally, the best analysis tools for this purpose of two-phase flow phenomena inside the BWR fuel assembly have been the sub-channel codes. However, the resolution of these codes is too coarse for analyzing the detailed intra-assembly flow patterns, such as flow around a spacer element. Advanced CFD (Computational Fluid Dynamics) codes provide a potential for detailed 3D simulations of coolant flow inside a fuel assembly, including flow around a spacer element using more fundamental physical models of flow regimes and phase interactions than sub-channel codes. Such models can extend the code applicability to a wider range of situations, which is highly important for increasing the efficiency and to prevent accidents.
NASA Astrophysics Data System (ADS)
Ovaysi, S.; Piri, M.
2009-12-01
We present a three-dimensional fully dynamic parallel particle-based model for direct pore-level simulation of incompressible viscous fluid flow in disordered porous media. The model was developed from scratch and is capable of simulating flow directly in three-dimensional high-resolution microtomography images of naturally occurring or man-made porous systems. It reads the images as input where the position of the solid walls are given. The entire medium, i.e., solid and fluid, is then discretized using particles. The model is based on Moving Particle Semi-implicit (MPS) technique. We modify this technique in order to improve its stability. The model handles highly irregular fluid-solid boundaries effectively. It takes into account viscous pressure drop in addition to the gravity forces. It conserves mass and can automatically detect any false connectivity with fluid particles in the neighboring pores and throats. It includes a sophisticated algorithm to automatically split and merge particles to maintain hydraulic connectivity of extremely narrow conduits. Furthermore, it uses novel methods to handle particle inconsistencies and open boundaries. To handle the computational load, we present a fully parallel version of the model that runs on distributed memory computer clusters and exhibits excellent scalability. The model is used to simulate unsteady-state flow problems under different conditions starting from straight noncircular capillary tubes with different cross-sectional shapes, i.e., circular/elliptical, square/rectangular and triangular cross-sections. We compare the predicted dimensionless hydraulic conductances with the data available in the literature and observe an excellent agreement. We then test the scalability of our parallel model with two samples of an artificial sandstone, samples A and B, with different volumes and different distributions (non-uniform and uniform) of solid particles among the processors. An excellent linear scalability is
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.
Effect of Compressibility on Hyperbolicity and Choke Flow Criterion of the Two-phase Two-fluid Model
Suneet Singh; Vincent A. Mousseau
2008-09-01
The standard two-phase two-fluid model lacks hyperbolicity which results in oscillations in the numerical solutions. For the incompressible two-phase flows an exact correction term can be derived which when added to the momentum equations makes the model hyperbolic. No such straightforward approach exists for the similar compressible flows. In the current work, the effect of the compressibility on the characteristic equation is analyzed. It is shown that the hyperbolicity of the system depends only on the slip velocity and not on the phasic velocities, independently. Moreover, a slip Mach number is defined and a non-dimensional characteristic equation is derived. It is shown that for the small values of slip Mach number the effect of the compressibility on the hyperbolicity can be ignored. To verify the above analysis, the characteristic equation for the two-phase compressible flows is numerically solved and results compared with the values obtained with the analytical solution for incompressible flows. Numerical solution of the two-phase two-fluid model for the benchmark problem is used to further verify the abovementioned analysis. Furthermore, the eigenvalues of the characteristic equation are obtained as a power series expansion about the point where the slip Mach number is zero. These eigenvalues are used to develop a choking criterion for the compressible two-phase flows.
Fluid and thermal mixing in a model cold leg and downcomer with vent-valve flow. [PWR
Rothe, P.H.; Marscher, W.D.; Block, J.A.
1982-03-01
This report describes an experimental program of fluid mixing experiments performed at atmospheric pressure in a 1/5-scale transparent model of the cold leg and downcomer of typical Babcock and Wilcox pressurized water reactors with vent valves. The results include transient data from a grid of thermocouples and extensive flow visualization photographs. Substantial mixing of cold injected water with hot primary coolant occurred during many of the tests.
NASA Astrophysics Data System (ADS)
Neustupa, T.
2016-06-01
This paper deals with the mathematical model of a steady flow of a heat-conductive incompressible viscous fluid through a spatially periodic plane profile cascade. The corresponding boundary value problem is reduced to one spatial period. We prove the existence of a weak solution of a coupled problem, with various boundary conditions on the parts of the boundary. Particularly, the condition on the outflow is a variant of the so called "do nothing" boundary condition.
NASA Astrophysics Data System (ADS)
Frederick, Jennifer M.; Buffett, Bruce A.
2013-05-01
of cosmogenic iodine, 129I, in the pore fluid of marine sediments often indicate that the pore fluid is much older than the host sediment, even when vertical flow due to sediment compaction is taken into account. Old pore fluid has been used in previous studies to argue for pervasive upward fluid flow and a deep methane source for hydrate deposits. Alternatively, old pore fluid age may reflect more complex flow patterns. We use a two-dimensional numerical transport model to account for the effects of topography and fractures on pore fluid pathlines when sediment permeability is anisotropic. We find that fluid focusing can cause significant lateral migration as well as regions where downward flow reverses direction and returns toward the seafloor. Longer pathlines can produce pore fluid ages much older than that expected with a one-dimensional compaction model. For steady-state models with geometry representative of Blake Ridge (USA), a well-studied hydrate province, we find pore fluid ages beneath regions of topography and within fractured zones that are up to 70 Ma old. Our results suggest that the measurements of 129I/127I reflect a mixture of new and old pore fluid. However, old pore fluid need not originate at great depths. Methane within pore fluids can travel laterally several kilometers, implying an extensive source region around the deposit. This type of focusing should aid hydrate formation beneath topographic highs.
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.
NASA Astrophysics Data System (ADS)
Zhou, J. X.; Shen, X.; Yin, Y. J.; Guo, Z.; Wang, H.
2015-06-01
In this paper, Gas-liquid two phase flow mathematic models of incompressible fluid were proposed to explore the feature of fluid under certain centrifugal force in vertical centrifugal casting (VCC). Modified projection-level-set method was introduced to solve the mathematic models. To validate the simulation results, two methods were used in this study. In the first method, the simulation result of basic VCC flow process was compared with its analytic solution. The relationship between the numerical solution and deterministic analytic solution was presented to verify the correctness of numerical algorithms. In the second method, systematic water simulation experiments were developed. In this initial experiment, special experimental vertical centrifugal device and casting shapes were designed to describe typical mold-filling processes in VCC. High speed camera system and data collection devices were used to capture flow shape during the mold-filling process. Moreover, fluid characteristic at different rotation speed (from 40rpm, 60rpmand 80rpm) was discussed to provide comparative resource for simulation results. As compared with the simulation results, the proposed mathematical models could be proven and the experimental design could help us advance the accuracy of simulation and further studies for VCC.
Pathapati, Subbu-Srikanth; Sansalone, John J
2011-07-01
Computational fluid dynamics (CFD) is emerging as a model for resolving the fate of particulate matter (PM) by unit operations subject to rainfall-runoff loadings. However, compared to steady flow CFD models, there are greater computational requirements for unsteady hydrodynamics and PM loading models. Therefore this study examines if integrating a stepwise steady flow CFD model can reproduce PM separation by common unit operations loaded by unsteady flow and PM loadings, thereby reducing computational effort. Utilizing monitored unit operation data from unsteady events as a metric, this study compares the two CFD modeling approaches for a hydrodynamic separator (HS), a primary clarifier (PC) tank, and a volumetric clarifying filtration system (VCF). Results indicate that while unsteady CFD models reproduce PM separation of each unit operation, stepwise steady CFD models result in significant deviation for HS and PC models as compared to monitored data; overestimating the physical size requirements of each unit required to reproduce monitored PM separation results. In contrast, the stepwise steady flow approach reproduces PM separation by the VCF, a combined gravitational sedimentation and media filtration unit operation that provides attenuation of turbulent energy and flow velocity. PMID:21644537
NASA Astrophysics Data System (ADS)
Wang, Xiao-Lu; Fan, Xiang-Yu; He, Yong-Ming; Nie, Ren-Shi; Huang, Quan-Hua
2013-08-01
Based on material balance and Darcy's law, the governing equation with the quadratic pressure gradient term was deduced. Then the nonlinear model for fluid flow in a multiple-zone composite reservoir including the quadratic gradient term was established and solved using a Laplace transform. A series of standard log-log type curves of 1-zone (homogeneous), 2-zone and 3-zone reservoirs were plotted and nonlinear flow characteristics were analysed. The type curves governed by the coefficient of the quadratic gradient term (β) gradually deviate from those of a linear model with time elapsing. Qualitative and quantitative analyses were implemented to compare the solutions of the linear and nonlinear models. The results showed that differences of pressure transients between the linear and nonlinear models increase with elapsed time and β. At the end, a successful application of the theoretical model data against the field data shows that the nonlinear model will be a good tool to evaluate formation parameters more accurately.
A site-scale model for fluid and heat flow in the unsaturated zone of Yucca Mountain, Nevada
NASA Astrophysics Data System (ADS)
Wu, Yu-Shu; Haukwa, Charles; Bodvarsson, G. S.
1999-05-01
A three-dimensional unsaturated-zone numerical model has been developed to simulate flow and distribution of moisture, gas and heat at Yucca Mountain, Nevada, a potential repository site for high-level radioactive waste. The model takes into account the simultaneous flow dynamics of liquid water, vapor, air and heat in the highly heterogeneous, fractured porous rock in the unsaturated zone (UZ). This model is intended for use in the prediction of the current and future conditions in the UZ so as to aid in the assessment of the system performance of the proposed repository. The modeling approach is based on a mathematical formulation of coupled multiphase, multicomponent fluid and heat flow through porous and fractured rock. Fracture and matrix flow is treated using both dual-permeability and effective-continuum modeling approaches. The model domain covers a total area of approximately 43 km 2, and uses the land surface and the water table as its top and bottom boundaries. In addition, site-specific data, representative surface infiltration, and geothermal conditions are incorporated into the model. The reliability and accuracy of the model have been the subject of a comprehensive model calibration study, in which the model was calibrated against measured data, including liquid saturation, water potential and temperature. It has been found that the model is generally able to reproduce the overall system behavior at Yucca Mountain with respect to moisture profiles, pneumatic pressure variations in different geological units, and ambient geothermal conditions.
NASA Astrophysics Data System (ADS)
Damsgaard, A.; Egholm, D. L.; Piotrowski, J. A.; Tulaczyk, S. M.; Larsen, N. K.
2014-12-01
The coupled mechanical response of ice, water and sediment may control the flow of warm-based glaciers residing on deformable sediment. This is most clearly expressed by the fast flowing ice streams in Greenland and Antarctica, where low levels of basal friction are thought to support the high flow rates. These ice streams are of particular interest since they are large constituents of the polar ice sheet mass balance. The study of these ice streams and their future impact on the ice sheets necessitates a deeper understanding of their basal dynamics, including the rheology of water-saturated sediment. We present the methodology and first results of a coupled numerical model for computational experiments on granular-fluid mixtures under dynamic conditions similar to those in subglacial settings. The granular phase is simulated on a per-particle basis by the soft body discrete element method. The fluid phase is handled as a continuum by solving the incompressible Navier-Stokes equations. The particle and fluid phases are coupled by mass conservation and momentum exchanges. The hydraulic diffusivity and permeability is compared to previous laboratory studies on tills. We demonstrate how the onset and halt of granular deformation is an efficient mechanism to create fluid pressure fluctuations due to local porosity changes. These pressure anomalies are driving transient hydraulic flows, and they influence directly the rheology of granular-fluid mixtures. Our results highlight the nonlinear nature of water saturated granular deformation, and demonstrate how the mechanical behaviour of granular materials may include both brittle and viscous components depending on the rates of deformation and the hydrological properties.
NASA Technical Reports Server (NTRS)
Wu, Jie; Yu, Sheng-Tao; Jiang, Bo-nan
1996-01-01
In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.
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
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.
NASA Astrophysics Data System (ADS)
Boutchko, Rostyslav; Rayz, Vitaliy L.; Vandehey, Nicholas T.; O'Neil, James P.; Budinger, Thomas F.; Nico, Peter S.; Druhan, Jennifer L.; Saloner, David A.; Gullberg, Grant T.; Moses, William W.
2012-01-01
This paper presents experimental and modeling aspects of applying nuclear emission tomography to study fluid flow in laboratory packed porous media columns of the type frequently used in geophysics, geochemistry and hydrology research. Positron emission tomography (PET) and single photon emission computed tomography (SPECT) are used as non-invasive tools to obtain dynamic 3D images of radioactive tracer concentrations. Dynamic sequences obtained using 18F-FDG PET are used to trace flow through a 5 cm diameter × 20 cm tall sand packed column with and without an impermeable obstacle. In addition, a custom-made rotating column setup placed in a clinical two-headed SPECT camera is used to image 99mTc-DTPA tracer propagation in a through-flowing column (10 cm diameter × 30 cm tall) packed with recovered aquifer sediments. A computational fluid dynamics software package FLUENT is used to model the observed flow dynamics. Tracer distributions obtained in the simulations in the smaller column uniformly packed with sand and in the column with an obstacle are remarkably similar to the reconstructed images in the PET experiments. SPECT results demonstrate strongly non-uniform flow patterns for the larger column slurry-packed with sub-surface sediment and slow upward flow. In the numerical simulation of the SPECT study, two symmetric channels with increased permeability are prescribed along the column walls, which result in the emergence of two well-defined preferential flow paths. Methods and results of this work provide new opportunities in hydrologic and biogeochemical research. The primary target application for developed technologies is non-destructive, non-perturbing, quantitative imaging of flow dynamics within laboratory scale porous media systems.
Visualization of Two-Fluid-Phase Flow Dynamics Using Micro-models
NASA Astrophysics Data System (ADS)
Dye, A. L.; McClure, J. E.; Gray, W. G.; Pyrak-Nolte, L. J.; Miller, C. T.
2015-12-01
Microfluidic approaches provide an opportunity to observe phase distributions, interfaces, common curves, and curvatures at high spatial and temporal resolutions in two-fluid-phase porous medium systems. Such observations provide an opportunity to visualize mechanisms, determine closure relations needed for existing and evolving models, and provide a means to validate theoretical predictions. We use a 500-micron cell with a circular solid phase to study two-fluid-phase displacement, including a dense pattern of drainage, imbibition, and scanning curves at very high spatial and temporal distribution. Viscous fingering, Haines jumps, and snap-off processes are visualized. A lattice Boltzmann method (LBM) is used to model the experimental data and shows excellent agreement for equilibrium states and certain aspects of the dynamics. Pressure dynamics, displacement processes, and interface curvature change time scales are shown. The fluid pressure, saturation, curvature, and interfacial area data is used to demonstrate that a unique, non-hysteretic closure relation exists. Further, the data is compared to mechanistic predictions based upon the thermodynamically constrained averaging theory (TCAT), validating key aspects of the theory.
Experimental and numerical modelling of the fluid flow in the continuous casting of steel
NASA Astrophysics Data System (ADS)
Timmel, K.; Miao, X.; Wondrak, T.; Stefani, F.; Lucas, D.; Eckert, S.; Gerbeth, G.
2013-03-01
This article gives an overview of recent research activities with respect to the mold flow in the continuous casting of steel in presence of DC magnetic fields. The magnetic fields appear to be an attractive tool for controlling the melt flow in a contactless way. Various kinds of magnetic systems are already in operation in industrial steel casting, but the actual impact on the melt flow has not been sufficiently verified by experimental studies. The rapid development of innovative diagnostic techniques in low-melting liquid metals over the last two decades enables new possibilities for systematic flow measurements in liquid metal model experiments. A new research program was initiated at HZDR comprising three experimental facilities providing a LIquid Metal Model for continuous CASTing of steel (LIMMCAST). The facilities operate in a temperature range from room temperature up to 400∘C using the low-melting alloys GaInSn and SnBi, respectively. The experimental program is focused on quantitative flow measurements in the mold, the submerged entry nozzle and the tundish. Local potential probes, Ultrasonic Doppler Velocimetry (UDV) and Contactless Inductive Flow Tomography (CIFT) are employed to measure the melt flow. The behavior of two-phase flows in case of argon injection is investigated by means of the Mutual Inductance Tomography (MIT) and X-ray radioscopy. The experimental results provide a substantial data basis for the validation of related numerical simulations. Numerical calculations were performed with the software package ANSYS-CFX with an implemented RANS-SST turbulence model. The non-isotropic nature of MHD turbulence was taken into account by specific modifications of the turbulence model. First results of the LIMMCAST program reveal important findings such as the peculiar, unexpected phenomenon that the application of a DC magnetic field may excite non-steady, non-isotropic large-scale flow oscillations in the mold. Another important result of our
Guidelines in the experimental validation of a 3D heat and fluid flow model of keyhole laser welding
NASA Astrophysics Data System (ADS)
Courtois, Mickael; Carin, Muriel; Le Masson, Philippe; Gaied, Sadok; Balabane, Mikhaël
2016-04-01
During the past few years, numerous sophisticated models have been proposed to predict in a self-consistent way the dynamics of the keyhole, together with the melt pool and vapor jet. However, these models are only partially compared to experimental data, so the reliability of these models is questionable. The present paper aims to propose a more complete experimental set-up in order to validate the most relevant results calculated by these models. A complete heat transfer and fluid flow three-dimensional (3D) model is first proposed in order to describe laser welding in keyhole regimes. The interface is tracked with a level set method and fluid flows are calculated in liquid and gas. The mechanisms of recoil pressure and keyhole creation are highlighted in a fusion line configuration chosen as a reference. Moreover, a complete validation of the model is proposed with guidelines on the variables to observe. Numerous comparisons with dedicated experiments (thermocouples, pyrometry, high-speed camera) are proposed to estimate the validity of the model. In addition to traditional geometric measurements, the main variables calculated, temperatures, and velocities in the melt pool are at the center of this work. The goal is to propose a reference validation for complex 3D models proposed over the last few years.
Nanometer-scale imaging and pore-scale fluid flow modeling inchalk
Tomutsa, Liviu; Silin, Dmitriy; Radmilovich, Velimir
2005-08-23
For many rocks of high economic interest such as chalk,diatomite, tight gas sands or coal, nanometer scale resolution is neededto resolve the 3D-pore structure, which controls the flow and trapping offluids in the rocks. Such resolutions cannot be achieved with existingtomographic technologies. A new 3D imaging method, based on serialsectioning and using the Focused Ion Beam (FIB) technology has beendeveloped. FIB allows for the milling of layers as thin as 10 nanometersby using accelerated Ga+ ions to sputter atoms from the sample surface.After each milling step, as a new surface is exposed, a 2D image of thissurface is generated. Next, the 2D images are stacked to reconstruct the3D pore or grain structure. Resolutions as high as 10 nm are achievableusing this technique. A new image processing method uses directmorphological analysis of the pore space to characterize thepetrophysical properties of diverse formations. In addition to estimationof the petrophysical properties (porosity, permeability, relativepermeability and capillary pressures), the method is used for simulationof fluid displacement processes, such as those encountered in variousimproved oil recovery (IOR) approaches. Computed with the new methodcapillary pressure curves are in good agreement with laboratory data. Themethod has also been applied for visualization of the fluid distributionat various saturations from the new FIB data.
NASA Astrophysics Data System (ADS)
Jacobs, C. T.; Piggott, M. D.
2015-03-01
This model description paper introduces a new finite element model for the simulation of non-linear shallow water flows, called Firedrake-Fluids. Unlike traditional models that are written by hand in static, low-level programming languages such as Fortran or C, Firedrake-Fluids uses the Firedrake framework to automatically generate the model's code from a high-level abstract language called Unified Form Language (UFL). By coupling to the PyOP2 parallel unstructured mesh framework, Firedrake can then target the code towards a desired hardware architecture to enable the efficient parallel execution of the model over an arbitrary computational mesh. The description of the model includes the governing equations, the methods employed to discretise and solve the governing equations, and an outline of the automated solution process. The verification and validation of the model, performed using a set of well-defined test cases, is also presented along with a road map for future developments and the solution of more complex fluid dynamical systems.
Pearson, Natalie C; Shipley, Rebecca J; Waters, Sarah L; Oliver, James M
2014-12-01
A 2D model is developed for fluid flow, mass transport and cell distribution in a hollow fibre membrane bioreactor. The geometry of the modelling region is simplified by excluding the exit ports at either end and focusing on the upper half of the central section of the bioreactor. Cells are seeded on a porous scaffold throughout the extracapillary space (ECS), and fluid pumped through the bioreactor via the lumen inlet and/or exit ports. In the fibre lumen and porous fibre wall, flow is described using Stokes and Darcy governing equations, respectively, while in the ECS porous mixture theory is used to model the cells, culture medium and scaffold. Reaction-advection-diffusion equations govern the concentration of a solute of interest in each region. The governing equations are reduced by exploiting the small aspect ratio of the bioreactor. This yields a coupled system for the cell volume fraction, solute concentration and ECS water pressure which is solved numerically for a variety of experimentally relevant case studies. The model is used to identify different regimes of cell behaviour, and results indicate how the flow rate can be controlled experimentally to generate a uniform cell distribution under regimes relevant to nutrient- and/or chemotactic-driven behaviours. PMID:24036069
Viscoelastic flow modeling in the extrusion of a dough-like fluid
NASA Technical Reports Server (NTRS)
Dhanasekharan, M.; Kokini, J. L.; Janes, H. W. (Principal Investigator)
2000-01-01
This work attempts to investigate the effect of viscoelasticity and three-dimensional geometry in screw channels. The Phan-Thien Tanner (PTT) constitutive equation with simplified model parameters was solved in conjunction with the flow equations. Polyflow, a commercially available finite element code was used to solve the resulting nonlinear partial differential equations. The PTT model predicted one log scale lower pressure buildup compared to the equivalent Newtonian results. However, the velocity profile did not show significant changes for the chosen PTT model parameters. Past Researchers neglected viscoelastic effects and also the three dimensional nature of the flow in extruder channels. The results of this paper provide a starting point for further simulations using more realistic model parameters, which may enable the food engineer to more accurately scale-up and design extrusion processes.
Modeling Coupled Processes for Multiphase Fluid Flow in Mechanically Deforming Faults
NASA Astrophysics Data System (ADS)
McKenna, S. A.; Pike, D. Q.
2011-12-01
Modeling of coupled hydrological-mechanical processes in fault zones is critical for understanding the long-term behavior of fluids within the shallow crust. Here we utilize a previously developed cellular-automata (CA) model to define the evolution of permeability within a 2-D fault zone under compressive stress. At each time step, the CA model calculates the increase in fluid pressure within the fault at every grid cell. Pressure surpassing a critical threshold (e.g., lithostatic stress) causes a rupture in that cell, and pressure is then redistributed across the neighboring cells. The rupture can cascade through the spatial domain and continue across multiple time steps. Stress continues to increase and the size and location of rupture events are recorded until a percolating backbone of ruptured cells exists across the fault. Previous applications of this model consider uncorrelated random fields for the compressibility of the fault material. The prior focus on uncorrelated property fields is consistent with development of a number of statistical physics models including percolation processes and fracture propagation. However, geologic materials typically express spatial correlation and this can have a significant impact on the results of the pressure and permeability distributions. We model correlation of the fault material compressibility as a multiGaussian random field with a correlation length defined as the full-width at half maximum (FWHM) of the kernel used to create the field. The FWHM is varied from < 0.001 to approximately 0.47 of the domain size. The addition of spatial correlation to the compressibility significantly alters the model results including: 1) Accumulation of larger amounts of strain prior to the first rupture event; 2) Initiation of the percolating backbone at lower amounts of cumulative strain; 3) Changes in the event size distribution to a combined power-law and exponential distribution with a smaller power; and 4) Evolution of the
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).
NASA Technical Reports Server (NTRS)
Makuch, Lauren A.
2004-01-01
Humans reach peak bone mass at age 30. After this point, we lose 1 to 2 percent of bone mass each decade. In the microgravity environment of space, astronauts lose bone mass at an accelerated rate of 1 to 2 percent each month. When astronauts travel to Mars, they may be in space for as long as 3 years. During this time, they may lose about half of their bone mass from weight-bearing bones. This loss may be irreversible. The drastic loss in bone that astronauts experience in space makes them much more vulnerable to fractures. In addition, the corresponding removal of calcium from bone results in higher levels of calcium in the blood, which increases the risk of developing kidney stones. Currently, studies are being conducted which investigate factors governing bone adaptation and mechanotransduction. Bone is constantly adapting in response to mechanical stimuli. Increased mechanical loading stimulates bone formation and suppresses bone resorption. Reduction in mechanical loading caused by bedrest, disuse, or microgravity results in decreased bone formation and possibly increased bone resorption. Osteoblasts and osteoclasts are the two main cell types that participate in bone remodeling. Osteoblasts are anabolic (bone-forming) cells and osteoclasts are catabolic (bone-resorbing) cells. In microgravity, the activity of osteoblasts slows down and the activity of osteoclasts may speed up, causing a loss of bone density. Mechanotransduction, the molecular mechanism by which mechanical stimuli are converted to biochemical signals, is not yet understood. Exposure of cells to fluid flow imposes a shear stress on the cells. Several studies have shown that the shear stress that results from fluid flow induces a cellular response similar to that induced by mechanical loading. Thus, fluid flow can be used as an in vitro model to simulate the mechanical stress that bone cells experience in vivo. Previous in vitro studies have shown that fluid flow induces several responses in
Rubab, Khansa; Mustafa, M
2016-01-01
This letter investigates the MHD three-dimensional flow of upper-convected Maxwell (UCM) fluid over a bi-directional stretching surface by considering the Cattaneo-Christov heat flux model. This model has tendency to capture the characteristics of thermal relaxation time. The governing partial differential equations even after employing the boundary layer approximations are non linear. Accurate analytic solutions for velocity and temperature distributions are computed through well-known homotopy analysis method (HAM). It is noticed that velocity decreases and temperature rises when stronger magnetic field strength is accounted. Penetration depth of temperature is a decreasing function of thermal relaxation time. The analysis for classical Fourier heat conduction law can be obtained as a special case of the present work. To our knowledge, the Cattaneo-Christov heat flux model law for three-dimensional viscoelastic flow problem is just introduced here. PMID:27093542
Rubab, Khansa; Mustafa, M.
2016-01-01
This letter investigates the MHD three-dimensional flow of upper-convected Maxwell (UCM) fluid over a bi-directional stretching surface by considering the Cattaneo-Christov heat flux model. This model has tendency to capture the characteristics of thermal relaxation time. The governing partial differential equations even after employing the boundary layer approximations are non linear. Accurate analytic solutions for velocity and temperature distributions are computed through well-known homotopy analysis method (HAM). It is noticed that velocity decreases and temperature rises when stronger magnetic field strength is accounted. Penetration depth of temperature is a decreasing function of thermal relaxation time. The analysis for classical Fourier heat conduction law can be obtained as a special case of the present work. To our knowledge, the Cattaneo-Christov heat flux model law for three-dimensional viscoelastic flow problem is just introduced here. PMID:27093542
Faybishenko, Boris; Doughty, Christine; Stoops, Thomas M.; Wood, thomas R.; Wheatcraft, Stephen W.
1999-12-31
(1) To determine if and when dynamical chaos theory can be used to investigate infiltration of fluid and contaminant transport in heterogeneous soils and fractured rocks. (2) To introduce a new approach to the multiscale characterization of flow and transport in fractured basalt vadose zones and to develop physically based conceptual models on a hierarchy of scales. The following activities are indicative of the success in meeting the project s objectives: A series of ponded infiltration tests, including (1) small-scale infiltration tests (ponded area 0.5 m2) conducted at the Hell s Half Acre site near Shelley, Idaho, and (2) intermediate-scale infiltration tests (ponded area 56 m2) conducted at the Box Canyon site near Arco, Idaho. Laboratory investigations and modeling of flow in a fractured basalt core. A series of small-scale dripping experiments in fracture models. Evaluation of chaotic behavior of flow in laboratory and field experiments using methods from nonlinear dynamics; Evaluation of the impact these dynamics may have on contaminant transport through heterogeneous fractured rocks and soils, and how it can be used to guide remediation efforts; Development of a conceptual model and mathematical and numerical algorithms for flow and transport that incorporate (1) the spatial variability of heterogeneous porous and fractured media, and (2) the description of the temporal dynamics of flow and transport, both of which may be chaotic. Development of appropriate experimental field and laboratory techniques needed to detect diagnostic parameters for chaotic behavior of flow. This approach is based on the assumption that spatial heterogeneity and flow phenomena are affected by nonlinear dynamics, and in particular, by chaotic processes. The scientific and practical value of this approach is that we can predict the range within which the parameters of flow and transport change with time in order to design and manage the remediation, even when we can not predict
Numerical Model of Fluid Flow through Heterogeneous Rock for High Level Radioactive Waste Disposal
NASA Astrophysics Data System (ADS)
Shirai, M.; Chiba, R.; Fomin, S.; Chugunov, V.; Takahashi, T.; Niibori, Y.; Hashida, T.
2007-03-01
An international consensus has emerged that deep geological disposal on land is one of the most appropriate means for high level radioactive wastes (HLW). The fluid transport is slow and radioactive elements are dangerous, so it's impossible to experiment over thousands of years. Instead, numerical model in such natural barrier as fractured underground needs to be considered. Field observations reveal that the equation with fractional derivative is more appropriate for describing physical phenomena than the equation which is based on the Fick's law. Thus, non-Fickian diffusion into inhomogeneous underground appears to be important in the assessment of HLW disposal. A solute transport equation with fractional derivative has been suggested and discussed in literature. However, no attempts were made to apply this equation for modeling of HLW disposal with account for the radioactive decay. In this study, we suggest the use of a novel fractional advection-diffusion equation which accounts for the effect of radioactive disintegration and for interactions between major, macro pores and fractal micro pores. This model is fundamentally different from previous proposed model of HLW, particularly in utilizing fractional derivative. Breakthrough curves numerically obtained by the present model are presented for a variety of rock types with respect to some important nuclides. Results of the calculation showed that for longer distance our model tends to be more conservative than the conventional Fickian model, therefore our model can be said to be safer.
Numerical Model of Fluid Flow through Heterogeneous Rock for High Level Radioactive Waste Disposal
Shirai, M.; Chiba, R.; Takahashi, T.; Hashida, T.; Fomin, S.; Chugunov, V.; Niibori, Y.
2007-03-20
An international consensus has emerged that deep geological disposal on land is one of the most appropriate means for high level radioactive wastes (HLW). The fluid transport is slow and radioactive elements are dangerous, so it's impossible to experiment over thousands of years. Instead, numerical model in such natural barrier as fractured underground needs to be considered. Field observations reveal that the equation with fractional derivative is more appropriate for describing physical phenomena than the equation which is based on the Fick's law. Thus, non-Fickian diffusion into inhomogeneous underground appears to be important in the assessment of HLW disposal. A solute transport equation with fractional derivative has been suggested and discussed in literature. However, no attempts were made to apply this equation for modeling of HLW disposal with account for the radioactive decay. In this study, we suggest the use of a novel fractional advection-diffusion equation which accounts for the effect of radioactive disintegration and for interactions between major, macro pores and fractal micro pores. This model is fundamentally different from previous proposed model of HLW, particularly in utilizing fractional derivative. Breakthrough curves numerically obtained by the present model are presented for a variety of rock types with respect to some important nuclides. Results of the calculation showed that for longer distance our model tends to be more conservative than the conventional Fickian model, therefore our model can be said to be safer.
NASA Astrophysics Data System (ADS)
Chitta, Varun
Modeling of complex flows involving the combined effects of flow transition and streamline curvature using two advanced turbulence models, one in the Reynolds-averaged Navier-Stokes (RANS) category and the other in the hybrid RANS-Large eddy simulation (LES) category is considered in this research effort. In the first part of the research, a new scalar eddy-viscosity model (EVM) is proposed, designed to exhibit physically correct responses to flow transition, streamline curvature, and system rotation effects. The four equation model developed herein is a curvature-sensitized version of a commercially available three-equation transition-sensitive model. The physical effects of rotation and curvature (RC) enter the model through the added transport equation, analogous to a transverse turbulent velocity scale. The eddy-viscosity has been redefined such that the proposed model is constrained to reduce to the original transition-sensitive model definition in nonrotating flows or in regions with negligible RC effects. In the second part of the research, the developed four-equation model is combined with a LES technique using a new hybrid modeling framework, dynamic hybrid RANS-LES. The new framework is highly generalized, allowing coupling of any desired LES model with any given RANS model and addresses several deficiencies inherent in most current hybrid models. In the present research effort, the DHRL model comprises of the proposed four-equation model for RANS component and the MILES scheme for LES component. Both the models were implemented into a commercial computational fluid dynamics (CFD) solver and tested on a number of engineering and generic flow problems. Results from both the RANS and hybrid models show successful resolution of the combined effects of transition and curvature with reasonable engineering accuracy, and for only a small increase in computational cost. In addition, results from the hybrid model indicate significant levels of turbulent
Fluid dynamic modeling and numerical simulation of low-density hypersonic flow
NASA Technical Reports Server (NTRS)
Cheng, H. K.; Wong, Eric Y.
1988-01-01
The concept of a viscous shock-layer and several related versions of continuum theories/methods are examined for their adequacy as a viable framework to study flow physics and aerothermodynamics of relevance to sustained hypersonic flights. Considering the flat plate at angle of attack, or the wedge, as a generic example for the major aerodynamic component of a hypersonic vehicle, the relative importance of the molecular-transport effects behind the shock (in the form of the 'shock slip') and the wall-slip effects are studied. In the flow regime where the shock-transition-zone thickness remains small compared to the shock radius of curvature, a quasi-one-dimensional shock structure under the Burnett/thirteen-moment approximation, as well as particulate/collisional models, can be consistently developed. The fully viscous version of the shock-layer model is shown to provide the crucial boundary condition downstream the shock in this case. The gas-kinetic basis of the continuum description for the flow behind the bow shock, and certain features affecting the non-equilibrium flow chemistry, are also discussed.
NASA Astrophysics Data System (ADS)
Lundgren, P.; Lanari, R.; Manzo, M.; Sansosti, E.; Tizzani, P.; Hutnak, M.; Hurwitz, S.
2008-12-01
Campi Flegrei caldera, Italy, located along the Bay of Naples, has a long history of significant vertical deformation, with the most recent large uplift (>1.5m) occurring in 1983-1984. Each episode of uplift has been followed by a period of subsidence that decreases in rate with time and may be punctuated by brief episodes of lesser uplift. The large amplitude of the major uplifts that occur without volcanic activity, and the subsequent subsidence has been argued as evidence for hydrothermal amplification of any magmatic source. The later subsidence and its temporal decay have been argued as due to diffusion of the pressurized caldera fill material into the less porous surrounding country rock. We present satellite synthetic aperture radar (SAR) interferometry (InSAR) time series analysis of ERS and Envisat data from the European Space Agency, based on exploiting the Small Baseline Subset (SBAS) approach [Berardino et al., 2002]; this allows us to generate maps of relative surface deformation though time, beginning in 1992 through 2007, that are relevant to both ascending and descending satellite orbits. The general temporal behavior is one of subsidence punctuated by several lesser uplift episodes. The spatial pattern of deformation can be modeled through simple inflation/deflation sources in an elastic halfspace. Given the evidence to suggest that fluids may play a significant role in the temporal deformation of Campi Flegrei, rather than a purely magmatic or magma chamber-based interpretation, we model the temporal and spatial evolution of surface deformation as a hydrothermal fluid flow process. We use the TOUGH2-BIOT2 set of numerical codes [Preuss et al., 1999; Hsieh, 1996], which couple multi-phase (liquid-gas) and multi-component (H2O-CO2) fluid flow in a porous or fractured media with plane strain deformation and fluid flow in a linearly elastic porous medium. We explore parameters related to the depth and temporal history of fluid injection, fluid
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.
NASA Technical Reports Server (NTRS)
Habbal, Shadia Rifai; Esser, Ruth; Guhathakurta, Madhulika; Fisher, Richard
1995-01-01
Using the empirical constraints provided by observations in the inner corona and in interplanetary space. we derive the flow properties of the solar wind using a two fluid model. Density and scale height temperatures are derived from White Light coronagraph observations on SPARTAN 201-1 and at Mauna Loa, from 1.16 to 5.5 R, in the two polar coronal holes on 11-12 Apr. 1993. Interplanetary measurements of the flow speed and proton mass flux are taken from the Ulysses south polar passage. By comparing the results of the model computations that fit the empirical constraints in the two coronal hole regions, we show how the effects of the line of sight influence the empirical inferences and subsequently the corresponding numerical results.