ITG sideband coupling models for zonal flows
Stransky, M.
2011-05-15
Four-wave interaction model between ITG mode and zonal flow was derived using fluid equations. In this model, the zonal flow is excited non-linearly by ITG turbulence via Reynolds stress. Numerical simulations show that the system allows for a small range above the ITG threshold where the zonal flow can stabilize an unstable ITG mode, effectively increasing {eta}{sub i} threshold, an effect which has been called the Dimits shift. However, the shift is smaller than in known cases such that in the Cyclone base.
Model of globally coupled Duffing flows
NASA Astrophysics Data System (ADS)
Shimada, Tokuzo; Moriya, Takanobu
2014-03-01
A Duffing oscillator in a certain parameter range shows period-doubling that has the same Feigenbaum ratio as the logistic map, which is an important issue in universality in chaos. In this paper a globally coupled lattice of Duffing flows (GCFL), which is a natural extension of the globally coupled logistic map lattice (GCML), is constructed. It is observed that GCFL inherits various intriguing properties of GCML and that universality at the level of elements is thus lifted to that of systems. Phase diagrams for GCFL are determined, which are essentially the same as those for GCML. Similar to the two-clustered periodic attractor of GCML, the GCFL two-clustered attractor exhibits a successive period-doubling with an increase of population imbalance between the clusters (\\vartheta -bifurcation). A nontrivial distinction between the GCML and GCFL attractors that originates from the symmetry in the Duffing equation is investigated in detail.
Coupling Stokes and Darcy Flow in Melt Migration Modelling
NASA Astrophysics Data System (ADS)
Kaus, B.; Lehmann, R.; Lukáčová-Medvid'ová, M.
2015-12-01
Melt migration can be modelled by coupling variable-viscosity Stokes flow and Darcy flow. Stokes Flow, generally, captures the long-term behavior of the mantle and lithosphere while Darcy flow models the two-phase regime. The major unknowns of the coupled system are solid velocity, fluid pressure and compaction pressure, captured in the so-called three-field formulation of the system. The fluid velocity can be computed in a post-processing step. We present lithosperic-scale results of the fully-coupled system with visco-elasto-plastic rheologies. This comprises elasto-plastic effects from shearing (Mode II) as well as poro-elastic effects and "opening mode" (Mode I) tensile plasticity. The system is solved using the Finite Element Method on triangular or quadrilateral grids in the Matlab-based code MVEP. Triangular meshes are adapted dynamically to better resolve the different deformation modes (diapiric, channeling, diking).
A coupled model of fluid flow in jointed rock
Swenson, Daniel; Martineau, Rick; James, Mark; Brown, Don
1991-01-01
We present a fully coupled model of fluid flow in jointed rock, where the fluid flow depends on the joint openings and the joint openings depend on the fluid pressure. The joints and rock blocks are modeled discretely using the finite element method. Solutions for the fluid and rock are obtained and iteration is performed until both solutions converge. Example applications include an examination of the effects of back-pressure on flow in a geothermal reservoir and transient fluid injection into a reservoir.
A New Equation Solver for Modeling Turbulent Flow in Coupled Matrix-Conduit Flow Models.
Hubinger, Bernhard; Birk, Steffen; Hergarten, Stefan
2016-07-01
Karst aquifers represent dual flow systems consisting of a highly conductive conduit system embedded in a less permeable rock matrix. Hybrid models iteratively coupling both flow systems generally consume much time, especially because of the nonlinearity of turbulent conduit flow. To reduce calculation times compared to those of existing approaches, a new iterative equation solver for the conduit system is developed based on an approximated Newton-Raphson expression and a Gauß-Seidel or successive over-relaxation scheme with a single iteration step at the innermost level. It is implemented and tested in the research code CAVE but should be easily adaptable to similar models such as the Conduit Flow Process for MODFLOW-2005. It substantially reduces the computational effort as demonstrated by steady-state benchmark scenarios as well as by transient karst genesis simulations. Water balance errors are found to be acceptable in most of the test cases. However, the performance and accuracy may deteriorate under unfavorable conditions such as sudden, strong changes of the flow field at some stages of the karst genesis simulations. PMID:26821785
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.
Acoustically-coupled flow-induced vibration of a computational vocal fold model.
Daily, David Jesse; Thomson, Scott L
2013-01-15
The flow-induced vibration of synthetic vocal fold models has been previously observed to be acoustically-coupled with upstream flow supply tubes. This phenomenon was investigated using a finite element model that included flow-structure-acoustic interactions. The length of the upstream duct was varied to explore the coupling between model vibration and subglottal acoustics. Incompressible and slightly compressible flow models were tested. The slightly compressible model exhibited acoustic coupling between fluid and solid domains in a manner consistent with experimental observations, whereas the incompressible model did not, showing the slightly compressible approach to be suitable for simulating acoustically-coupled vocal fold model flow-induced vibration. PMID:23585700
Acoustically-coupled flow-induced vibration of a computational vocal fold model
Daily, David Jesse; Thomson, Scott L.
2012-01-01
The flow-induced vibration of synthetic vocal fold models has been previously observed to be acoustically-coupled with upstream flow supply tubes. This phenomenon was investigated using a finite element model that included flow–structure–acoustic interactions. The length of the upstream duct was varied to explore the coupling between model vibration and subglottal acoustics. Incompressible and slightly compressible flow models were tested. The slightly compressible model exhibited acoustic coupling between fluid and solid domains in a manner consistent with experimental observations, whereas the incompressible model did not, showing the slightly compressible approach to be suitable for simulating acoustically-coupled vocal fold model flow-induced vibration. PMID:23585700
Developing hillslope-based catchment models: coupling Boussinesq and regional scale flow models
NASA Astrophysics Data System (ADS)
Broda, S.; Paniconi, C.; Larocque, M.
2009-04-01
The gaining recognition of hillslopes as fundamental building blocks in watershed hydrology makes them appealing for incorporation into larger scale river basin models. Hillslope processes are commonly simulated by means of the Boussinesq equation and are therefore applicable to single layer flow systems only. Two coupled models are presented to simulate both local hillslope scale and regional scale groundwater flow: 1) the hillslope-storage Boussinesq (hsB) model representing unconfined flow and a steady, analytic element model representing transient regional deep groundwater flow through a succession of steady state stress periods over many hydrological years, and 2) the hsB model and a newly developed analytical solution for 1D transient confined groundwater flow. Recharge zones are defined by means of irregular geometric domains, capturing the plan form geometry of the hillslopes. Lateral flows are calculated in inclined aquifers of homogeneous thickness. Tests are conducted on i) single hillslopes of varying inclination and plan form geometry and ii) a laboratory watershed, and heads and baseflows are compared to the results from a fully coupled 3D Richards equation model. Both approaches presented provide reasonable heads and fluxes for a range of hillslope properties in comparison to the benchmark model, and are promising approaches, applicable to a range of land surface models that lack a detailed description of subsurface flow. However the coupled hsB/1D-analytical model is numerically more stable and computationally more efficient than the coupled hsB/analytic element model.
A coupled stream flow and depth-integrated subsurface flow model for catchment hydrology
NASA Astrophysics Data System (ADS)
Pan, Yi; Weill, Sylvain; Ackerer, Philippe; Delay, Frederick
2015-11-01
Few hydrological models that couple stream flow and subsurface flow in shallow aquifers are based on a compromise between simple and complex depiction of the system, although this compromise could result in tractable tools for various applications. We present a depth-integrated approach in which flows in the vadose and saturated zones are assumed to be parallel to the bottom of the aquifer and thus are integrated in the direction normal to the bottom of the aquifer. The hydrodynamic parameters are also integrated in this direction, and gravity effects are preserved. Stream flow is handled by a diffusive-wave equation that is calculated over a network of one-dimensional bonds. The first-order coupling between the stream and subsurface flows exchanges water fluxes between the stream network and the subsurface compartment according to the hydraulic head differences between the systems. Three synthetic test cases, one including a comparison with a three-dimensional code, are used to evaluate the general behavior of the coupled model. It is shown that the approach reproduces the main hydrological features at the catchment scale, including the generation of runoff, infiltration-exfiltration into (from) the vadose zone, and smooth transient head variations in the aquifer.
Numerical simulation of the two-phase flows in a hydraulic coupling by solving VOF model
NASA Astrophysics Data System (ADS)
Luo, Y.; Zuo, Z. G.; Liu, S. H.; Fan, H. G.; Zhuge, W. L.
2013-12-01
The flow in a partially filled hydraulic coupling is essentially a gas-liquid two-phase flow, in which the distribution of two phases has significant influence on its characteristics. The interfaces between the air and the liquid, and the circulating flows inside the hydraulic coupling can be simulated by solving the VOF two-phase model. In this paper, PISO algorithm and RNG k-ɛ turbulence model were employed to simulate the phase distribution and the flow field in a hydraulic coupling with 80% liquid fill. The results indicate that the flow forms a circulating movement on the torus section with decreasing speed ratio. In the pump impeller, the air phase mostly accumulates on the suction side of the blades, while liquid on the pressure side; in turbine runner, air locates in the middle of the flow passage. Flow separations appear near the blades and the enclosing boundaries of the hydraulic coupling.
A numerical method for a model of two-phase flow in a coupled free flow and porous media system
NASA Astrophysics Data System (ADS)
Chen, Jie; Sun, Shuyu; Wang, Xiao-Ping
2014-07-01
In this article, we study two-phase fluid flow in coupled free flow and porous media regions. The model consists of coupled Cahn-Hilliard and Navier-Stokes equations in the free fluid region and the two-phase Darcy law in the porous medium region. We propose a Robin-Robin domain decomposition method for the coupled Navier-Stokes and Darcy system with the generalized Beavers-Joseph-Saffman condition on the interface between the free flow and the porous media regions. Numerical examples are presented to illustrate the effectiveness of this method.
Vertically integrated models for coupled two-phase flow and geomechanics in porous media
NASA Astrophysics Data System (ADS)
Bjørnarâ, Tore I.; Nordbotten, Jan M.; Park, Joonsang
2016-02-01
Models of reduced dimensionality have been found to be particularly attractive in simulating the fate of injected CO2 in supercritical state in the context of carbon capture and storage. This is motivated by the confluence of three aspects: the strong buoyant segregation of the lighter CO2 phase above water, the relatively long time scales associated with storage, and finally the large aspect ratios that characterize the geometry of typical storage aquifers. However, to date, these models have been confined to considering only the flow problem, as the coupling between reduced dimensionality models for flow and models for geomechanical response has previously not been developed. Herein, we develop a fully coupled, reduced dimension, model for multiphase flow and geomechanics. It is characterized by the aquifer(s) being of lower dimension(s), while the surrounding overburden and underburden being of full dimension. The model allows for general constitutive functions for fluid flow (relative permeability and capillary pressure) and uses the standard Biot coupling between the flow and mechanical equations. The coupled model retains all the simplicities of reduced-dimensional models for flow, including less stiff nonlinear systems of equations (since the upscaled constitutive functions are closer to linear), longer time steps (since the high grid resolution in the vertical direction can be avoided), and less degrees of freedom. We illustrate the applicability of the new coupled model through both a validation study and a practical computational example.
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.
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
NASA Astrophysics Data System (ADS)
Juanes, R.; Jha, B.
2014-12-01
The coupling between subsurface flow and geomechanical deformation is critical in the assessment of the environmental impacts of groundwater use, underground liquid waste disposal, geologic storage of carbon dioxide, and exploitation of shale gas reserves. In particular, seismicity induced by fluid injection and withdrawal has emerged as a central element of the scientific discussion around subsurface technologies that tap into water and energy resources. Here we present a new computational approach to model coupled multiphase flow and geomechanics of faulted reservoirs. We represent faults as surfaces embedded in a three-dimensional medium by using zero-thickness interface elements to accurately model fault slip under dynamically evolving fluid pressure and fault strength. We incorporate the effect of fluid pressures from multiphase flow in the mechanical stability of faults and employ a rigorous formulation of nonlinear multiphase geomechanics that is capable of handling strong capillary effects. We develop a numerical simulation tool by coupling a multiphase flow simulator with a mechanics simulator, using the unconditionally stable fixed-stress scheme for the sequential solution of two-way coupling between flow and geomechanics. We validate our modeling approach using several synthetic, but realistic, test cases that illustrate the onset and evolution of earthquakes from fluid injection and withdrawal. We also present the application of the coupled flow-geomechanics simulation technology to the post mortem analysis of the Mw=5.1, May 2011 Lorca earthquake in south-east Spain, and assess the potential that the earthquake was induced by groundwater extraction.
MODELING COUPLING OF EEL GRASS ZOSTRA MARINA AND WATER FLOW
Ecological effects caused by submerged aquatic vegetation not only depend on the plants and their morphology but also on the flow and transport patterns of dissolved and suspended constituents near the canopy. The height of the canopy is a major parameter in any quantitative an...
Kim, Jeongkon; Scwartz, Franklin W.; Shi, Jianyou; Xu, Tianfu
2003-04-01
A complex pattern of coupling between fluid flow and mass transport develops when heterogeneous reactions occur. For instance, dissolution and precipitation reactions can change the physical properties of a medium, such as permeability and pore geometry. These changes influence fluid flow, which in turn impact the composition of dissolved constituents and solid-phase, and the rate and direction of advective transport. Two-dimensional modeling studies using TOUGHREACT were conducted to investigate the coupling between flow and transport developed as a consequence of difference in density, dissolution/precipitation, and medium heterogeneity. The model includes equilibrium reactions for aqueous species, kinetic reactions between the solid phases and aqueous constituents, and full coupling of porosity and permeability changes resulting from precipitation and dissolution reactions in porous media. Generally, the evolutions in the concentrations of the aqueous phase are intimately related to the reaction-front dynamics. Plugging of the medium contributed to significant transients in patterns of flow and mass transport.
A general kinetic-flow coupling model for FCC riser flow simulation.
Chang, S. L.
1998-05-18
A computational fluid dynamic (CFD) code has been developed for fluid catalytic cracking (FCC) riser flow simulation. Depending on the application of interest, a specific kinetic model is needed for the FCC flow simulation. This paper describes a method to determine a kinetic model based on limited pilot-scale test data. The kinetic model can then be used with the CFD code as a tool to investigate optimum operating condition ranges for a specific FCC unit.
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.
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.
Dynamic coupling of pore-scale and reservoir-scale models for multiphase flow
NASA Astrophysics Data System (ADS)
Sheng, Qiang; Thompson, Karsten
2013-09-01
The concept of coupling pore-scale and continuum-scale models for subsurface flow has long been viewed as beneficial, but implementation has been slow. In this paper, we present an algorithm for direct coupling of a dynamic pore-network model for multiphase flow with a traditional continuum-scale simulator. The ability to run the two models concurrently (exchanging parameters and boundary conditions in real numerical time) is made possible by a new dynamic pore-network model that allows simultaneous injection of immiscible fluids under either transient-state or steady-state conditions. Allowing the pore-scale model to evolve to steady state during each time step provides a unique method for reconciling the dramatically different time and length scales across the coupled models. The model is implemented by embedding networks in selected gridblocks in the reservoir model. The network model predicts continuum-scale parameters such as relative permeability or average capillary pressure from first principles, which are used in the continuum model. In turn, the continuum reservoir simulator provides boundary conditions from the current time step back to the network model to complete the coupling process. The model is tested for variable-rate immiscible displacements under conditions in which relative permeability depends on flow rate, thus demonstrating a situation that cannot be modeled using a traditional approach. The paper discusses numerical challenges with this approach, including the fact that there is not a way to explicitly force pore-scale phase saturation to equal the continuum saturation in the host gridblock without an artificial constraint. Hurdles to implementing this type of modeling in practice are also discussed.
A partially open porous media flow with chaotic advection: towards a model of coupled fields.
Metcalfe, Guy; Lester, Daniel; Ord, Alison; Kulkarni, Pandurang; Trefry, Mike; Hobbs, Bruce E; Regenaur-Lieb, Klaus; Morris, Jeffery
2010-01-13
In nature, dissipative fluxes of fluid, heat and/or reacting species couple to each other and may also couple to deformation of a surrounding porous matrix. We use the well-known analogy of Hele-Shaw flow to Darcy flow to make a model porous medium with porosity proportional to local cell height. Time- and space-varying fluid injection from multiple source/sink wells lets us create many different kinds of chaotic flows and chemical concentration patterns. Results of an initial time-dependent potential flow model illustrate that this is a partially open flow, in which parts of the material transported by the flow remain in the cell forever and parts pass through with residence time and exit time distributions that have self-similar features in the control parameter space of the stirring. We derive analytically the existence boundary in stirring control parameter space between where isolated fluid regions can and cannot remain forever in the open flow. Experiments confirm the predictions. PMID:19948552
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.
NASA Astrophysics Data System (ADS)
Juanes, R.; Jha, B.; Hager, B. H.; Shaw, J. H.; Plesch, A.; Astiz, L.; Dieterich, J. H.; Frohlich, C.
2016-07-01
Seismicity induced by fluid injection and withdrawal has emerged as a central element of the scientific discussion around subsurface technologies that tap into water and energy resources. Here we present the application of coupled flow-geomechanics simulation technology to the post mortem analysis of a sequence of damaging earthquakes (Mw = 6.0 and 5.8) in May 2012 near the Cavone oil field, in northern Italy. This sequence raised the question of whether these earthquakes might have been triggered by activities due to oil and gas production. Our analysis strongly suggests that the combined effects of fluid production and injection from the Cavone field were not a driver for the observed seismicity. More generally, our study illustrates that computational modeling of coupled flow and geomechanics permits the integration of geologic, seismotectonic, well log, fluid pressure and flow rate, and geodetic data and provides a promising approach for assessing and managing hazards associated with induced seismicity.
Coupling of FVCOM and CFD Model for Simulation of Multiscale Coastal Flows
NASA Astrophysics Data System (ADS)
Tang, H.; Qu, K.
2013-12-01
In correspondence with the need to simulate many emerging problems, especially those in nearshore regions such as deepwater oil spill, it is necessary to develop capabilities to predict small-scale, fully 3D phenomena in coastal ocean flows. A feasible as well as effective approach for the development is a hybrid method that couples different models designed for physics at different scales. We have developed a two-way coupling between a fully 3D CFD model and the FVCOM, in which the former captures small-scale 3D flows and the latter predicts large-scale background currents. In this presentation, a few new applications of such approach will be illustrated. The following figures show a result on tidal flow in a bay and past bridge peers. Tests and analysis are made on solution accuracy and computation efficiency, and discussions are presented on how to achieve seamless solution transition at the model interfaces. Computed large-scale background flow Simulated flow past a brigde
NASA Astrophysics Data System (ADS)
Tsai, C. H.; Yeh, G. T.
2015-12-01
In this investigation, a coupled model of multiphase flow, reactive biogeochemical transport, thermal transport and geo-mechanics in subsurface media is presented. It iteratively solves the mass conservation equation for fluid flow, thermal transport equation for temperature, reactive biogeochemical transport equations for concentration distributions, and solid momentum equation for displacement with successive linearization algorithm. With species-based equations of state, density of a phase in the system is obtained by summing up concentrations of all species. This circumvents the problem of having to use empirical functions. Moreover, reaction rates of all species are incorporated in mass conservation equation for fluid flow. Formation enthalpy of all species is included in the law of energy conservation as a source-sink term. Finite element methods are used to discretize the governing equations. Numerical experiments are presented to examine the accuracy and robustness of the proposed model. The results demonstrate the feasibility and capability of present model in subsurface media.
Numerical Modeling of Coupled Water Flow and Heat Transport in Soil and Snow
NASA Astrophysics Data System (ADS)
Kelleners, T.
2015-12-01
A numerical model is developed to calculate coupled water flow and heat transport in seasonally frozen soil and snow. Both liquid water flow and water vapor flow are included. The effect of dissolved ions on soil water freezing point depression is included by combining an expression for osmotic head with the Clapeyron equation and the van Genuchten soil water retention function. The coupled water flow and heat transport equations are solved using the Thomas algorithm and Picard iteration. Ice pressure is always assumed zero and frost heave is neglected. The new model is tested using data from a high-elevation rangeland soil that is subject to significant soil freezing and a mountainous forest soil that is snow-covered for about 8 months of the year. Soil hydraulic parameters are mostly based on measurements and only vegetation parameters are fine-tuned to match measured and calculated soil water content, soil & snow temperature, and snow height. Modeling statistics for both systems show good performance for temperature, intermediate performance for snow height, and relatively low performance for soil water content, in accordance with earlier results with an older version of the model.
Coupled dam-break flow and bed load modelling using HLLC-WAF scheme.
Hosseinzadeh-Tabrizi, Alireza; Ghaeini-Hessaroeyeh, Mahnaz
2015-01-01
A two-dimensional numerical model predicting flow over a mobile bed has been developed. Governing equations consist of the shallow water equations and the Exner equation. The finite volume method on an unstructured triangular grid was deployed to discretize the governing equations. The local Riemann problem is solved by the Harten, Lax and van Leer-contact (HLLC) method in the interface of the cells and the equations are solved using a fully coupled method. Then the flux modelling has been deployed by the total variation diminishing (TVD) version of the weighted average flux (WAF) scheme. The model was verified by comparison of the results and available experimental data for dam-break flow, in a laboratory test, via a channel with sudden enlargement and erodible bed conditions. Comparison of these two sets of results shows that increasing the accuracy of flux modelling caused the model results to have a reasonable agreement with the experimental data. PMID:26398031
Thermochemical Nonequilibrium 2D Modeling of Nitrogen Inductively Coupled Plasma Flow
NASA Astrophysics Data System (ADS)
Yu, Minghao; Yusuke, Takahashi; Hisashi, Kihara; Ken-ichi, Abe; Kazuhiko, Yamada; Takashi, Abe; Satoshi, Miyatani
2015-09-01
Two-dimensional (2D) numerical simulations of thermochemical nonequilibrium inductively coupled plasma (ICP) flows inside a 10-kW inductively coupled plasma wind tunnel (ICPWT) were carried out with nitrogen as the working gas. Compressible axisymmetric Navier-Stokes (N-S) equations coupled with magnetic vector potential equations were solved. A four-temperature model including an improved electron-vibration relaxation time was used to model the internal energy exchange between electron and heavy particles. The third-order accuracy electron transport properties (3rd AETP) were applied to the simulations. A hybrid chemical kinetic model was adopted to model the chemical nonequilibrium process. The flow characteristics such as thermal nonequilibrium, inductive discharge, effects of Lorentz force were made clear through the present study. It was clarified that the thermal nonequilibrium model played an important role in properly predicting the temperature field. The prediction accuracy can be improved by applying the 3rd AETP to the simulation for this ICPWT. supported by Grant-in-Aid for Scientific Research (No. 23560954), sponsored by the Japan Society for the Promotion of Science
Flow-radiation coupling for atmospheric entries using a Hybrid Statistical Narrow Band model
NASA Astrophysics Data System (ADS)
Soucasse, Laurent; Scoggins, James B.; Rivière, Philippe; Magin, Thierry E.; Soufiani, Anouar
2016-09-01
In this study, a Hybrid Statistical Narrow Band (HSNB) model is implemented to make fast and accurate predictions of radiative transfer effects on hypersonic entry flows. The HSNB model combines a Statistical Narrow Band (SNB) model for optically thick molecular systems, a box model for optically thin molecular systems and continua, and a Line-By-Line (LBL) description of atomic radiation. Radiative transfer calculations are coupled to a 1D stagnation-line flow model under thermal and chemical nonequilibrium. Earth entry conditions corresponding to the FIRE 2 experiment, as well as Titan entry conditions corresponding to the Huygens probe, are considered in this work. Thermal nonequilibrium is described by a two temperature model, although non-Boltzmann distributions of electronic levels provided by a Quasi-Steady State model are also considered for radiative transfer. For all the studied configurations, radiative transfer effects on the flow, the plasma chemistry and the total heat flux at the wall are analyzed in detail. The HSNB model is shown to reproduce LBL results with an accuracy better than 5% and a speed up of the computational time around two orders of magnitude. Concerning molecular radiation, the HSNB model provides a significant improvement in accuracy compared to the Smeared-Rotational-Band model, especially for Titan entries dominated by optically thick CN radiation.
Modeling flow in nanoporous, membrane reservoirs and interpretation of coupled fluxes
NASA Astrophysics Data System (ADS)
Geren, Filiz
The average pore size in unconventional, tight-oil reservoirs is estimated to be less than 100 nm. At this pore size, Darcy flow is no longer the dominating flow mechanism and a combination of diffusive flows determines the flow characteristics. Concentration driven self-diffusion has been well known and included in the flow and transport models in porous media. However, when the sizes of the pores and pore-throats decrease down to the size of the hydrocarbon molecules, the porous medium acts like a semi-permeable membrane, and the size of the pore openings dictates the direction of transport between adjacent pores. Accordingly, characterization of flow and transport in tight unconventional plays requires understanding of their membrane properties. This Master of Science thesis first highlights the membrane properties of nanoporous, unconventional reservoirs and then discusses how filtration effects can be incorporated into the models of transport in nanoporous media within the coupled flux concept. The effect of filtration on fluid composition and its impact on black-oil fluid properties like bubble point pressure is also demonstrated. To define filtration and filtration pressure in unconventional, tight-oil reservoirs, analogy to chemical osmosis is applied two pore systems connected with a pore throat, which shows membrane properties. Because the pore throat selectivity permits the passage of fluid molecules by their sizes, given a filtration pressure difference between the two pore systems, the concentration difference between the systems is determined by flash calculations. The results are expressed in the form of filtration (membrane) efficiency, which is essential parameter to define coupled fluxes for porous media flow.
Surface models for coupled modelling of runoff and sewer flow in urban areas.
Ettrich, N; Steiner, K; Thomas, M; Rothe, R
2005-01-01
Traditional methods fail for the purpose of simulating the complete flow process in urban areas as a consequence of heavy rainfall and as required by the European Standard EN-752 since the bi-directional coupling between sewer and surface is not properly handled. The new methodology, developed in the EUREKA-project RisUrSim, solves this problem by carrying out the runoff on the basis of shallow water equations solved on high-resolution surface grids. Exchange nodes between the sewer and the surface, like inlets and manholes, are located in the computational grid and water leaving the sewer in case of surcharge is further distributed on the surface. Dense topographical information is needed to build a model suitable for hydrodynamic runoff calculations; in urban areas, in addition, many line-shaped elements like houses, curbs, etc. guide the runoff of water and require polygonal input. Airborne data collection methods offer a great chance to economically gather densely sampled input data. PMID:16248177
NASA Astrophysics Data System (ADS)
Kalteh, Aman Mohammad
2013-04-01
Reliable and accurate forecasts of river flow is needed in many water resources planning, design development, operation and maintenance activities. In this study, the relative accuracy of artificial neural network (ANN) and support vector regression (SVR) models coupled with wavelet transform in monthly river flow forecasting is investigated, and compared to regular ANN and SVR models, respectively. The relative performance of regular ANN and SVR models is also compared to each other. For this, monthly river flow data of Kharjegil and Ponel stations in Northern Iran are used. The comparison of the results reveals that both ANN and SVR models coupled with wavelet transform, are able to provide more accurate forecasting results than the regular ANN and SVR models. However, it is found that SVR models coupled with wavelet transform provide better forecasting results than ANN models coupled with wavelet transform. The results also indicate that regular SVR models perform slightly better than regular ANN models.
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.
NASA Astrophysics Data System (ADS)
Wissmeier, L. C.; Barry, D. A.
2009-12-01
Computer simulations of water availability and quality play an important role in state-of-the-art water resources management. However, many of the most utilized software programs focus either on physical flow and transport phenomena (e.g., MODFLOW, MT3DMS, FEFLOW, HYDRUS) or on geochemical reactions (e.g., MINTEQ, PHREEQC, CHESS, ORCHESTRA). In recent years, several couplings between both genres of programs evolved in order to consider interactions between flow and biogeochemical reactivity (e.g., HP1, PHWAT). Software coupling procedures can be categorized as ‘close couplings’, where programs pass information via the memory stack at runtime, and ‘remote couplings’, where the information is exchanged at each time step via input/output files. The former generally involves modifications of software codes and therefore expert programming skills are required. We present a generic recipe for remotely coupling the PHREEQC geochemical modeling framework and flow and solute transport (FST) simulators. The iterative scheme relies on operator splitting with continuous re-initialization of PHREEQC and the FST of choice at each time step. Since PHREEQC calculates the geochemistry of aqueous solutions in contact with soil minerals, the procedure is primarily designed for couplings to FST’s for liquid phase flow in natural environments. It requires the accessibility of initial conditions and numerical parameters such as time and space discretization in the input text file for the FST and control of the FST via commands to the operating system (batch on Windows; bash/shell on Unix/Linux). The coupling procedure is based on PHREEQC’s capability to save the state of a simulation with all solid, liquid and gaseous species as a PHREEQC input file by making use of the dump file option in the TRANSPORT keyword. The output from one reaction calculation step is therefore reused as input for the following reaction step where changes in element amounts due to advection
Wagner, B.J.
1992-01-01
Parameter estimation and contaminant source characterization are key steps in the development of a coupled groundwater flow and contaminant transport simulation model. Here a methodologyfor simultaneous model parameter estimation and source characterization is presented. The parameter estimation/source characterization inverse model combines groundwater flow and contaminant transport simulation with non-linear maximum likelihood estimation to determine optimal estimates of the unknown model parameters and source characteristics based on measurements of hydraulic head and contaminant concentration. First-order uncertainty analysis provides a means for assessing the reliability of the maximum likelihood estimates and evaluating the accuracy and reliability of the flow and transport model predictions. A series of hypothetical examples is presented to demonstrate the ability of the inverse model to solve the combined parameter estimation/source characterization inverse problem. Hydraulic conductivities, effective porosity, longitudinal and transverse dispersivities, boundary flux, and contaminant flux at the source are estimated for a two-dimensional groundwater system. In addition, characterization of the history of contaminant disposal or location of the contaminant source is demonstrated. Finally, the problem of estimating the statistical parameters that describe the errors associated with the head and concentration data is addressed. A stage-wise estimation procedure is used to jointly estimate these statistical parameters along with the unknown model parameters and source characteristics. ?? 1992.
NASA Astrophysics Data System (ADS)
Pannetier, R.
2014-12-01
A need for improvements in numerical models for the representation of permafrost and active layer dynamics has been highlighted. Initial efforts to address this have been made, in particular by coupling heat and multiphase flow equations under transient conditions. Implications of honoring the physical representation of such process includes an improved ability to model and analyze effects of projected changes in hydro-climatic variables, as for example the change in thermal regime and the geographic patterns of permafrost degradation. In this contribution the recently developed frozen ground module of PFloTran for fully coupled transient heat and water flow under partially saturated conditions is applied to the subarctic field site of Tarfala in Northern Sweden, where the occurrence of permafrost is sporadic and strongly dependent on snow-depth. Optimizing this model to reproduce subsurface temperature fluctuations in the ground, as recorded in two different boreholes over the last decade, allows to identify the respective roles played by heat diffusion, advection and convection in a changing permafrost environment. The influence of topography is analyzed by applying the model to different domains geometries, where topography is simplified to different degrees. The model has been configured to reproduce temperature observations, providing a estimation of the distribution and thermal regime of permafrost over the entire hill-slope. Based on the configurations that give the most accurate validation against data, changes in permafrost are extrapolated by applying different scenarios of prospective climatic conditions.
Auroral signatures of Bursty Bulk Flows from magnetosphere-ionosphere coupling models
NASA Astrophysics Data System (ADS)
Echim, M.; de Keyser, J. M.; Roth, M. A.
2010-12-01
The relationship between bursty bulk flows (BBFs) in the magnetospheric tail and the activation of auroral forms is well established from satellite and ground-based observations. Starting from a self-consistent description of BBFs based on a Vlasov equilibrium we provide a quantitative evaluation of the associated auroral effects by using a quasi-stationary magnetosphere-ionosphere (MI) coupling model. The self-consistent BBF model is based on a kinetic description of a 1-D plasma slab moving in background plasma and electromagnetic field. The model considers two exact constants of motion and one adiabatic invariant (the magnetic moment). It solves the coupled Vlasov-Maxwell system of equations in one spatial dimension (perpendicular to the BBFs plasma bulk velocity and the main magnetic field) assuming the BBF is a 1D structure elongated in the direction of the background magnetic field. The BBF model provides the self-consistent profile of Φm, the electric potential, showing the formation of convergent electric fields at the dawnward flank of the Earth-ward oriented BBFs. It has been shown that magnetospheric convergent electric fields drive field-aligned (FA) potential drops, FA currents and electron precipitation and acceleration. A stationary MI coupling model developed for discontinuity-like magnetospheric generators with convergent electric fields developed earlier is adapted to describe the coupling between the BBFs and the auroral ionosphere. The kernel of the MI coupling model is the condition of current continuity at the topside ionosphere, from which we compute the electric potential in the ionosphere for a given Φm. The MI coupling model is based on a Knight-type current-voltage relationship and a height-integrated conductivity model that depends on the energy deposited in the ionosphere by precipitating electrons. We show that the convergent electric field formed at the flanks of the BBF drive a FA potential drop and downward electron acceleration
Thermo-osmosis coupled-flow characterization in clay-rocks: experiments and modeling
NASA Astrophysics Data System (ADS)
Tremosa, J.; Goncalves, J.; Matray, J.; Violette, S.
2009-12-01
Water flow in clay-rocks is not only driven by a hydraulic gradient but also by chemical, thermal or electrical gradients. It implies a re-evaluation of the Darcy law by considering all gradients occurring in the clay-rock and their associated coupling coefficients (e.g. the osmotic efficiency to link a chemical gradient to a water flow). The occurrence of such processes in clay-rocks is due to the low hydraulic conductivity of this media and because of electrical charges at the clay minerals surface. Here, we focused on the thermo-osmosis process, a water flow under a temperature gradient, which is poorly characterized in spite of its implications in nuclear waste storage in clay-rocks. A set of thermo-osmotic experiments was performed in an equipped borehole installed in a Toarcian compacted clay at the IRSN’s Underground Research Laboratory in the south of France. The water flow induced by a temperature gradient (from the hotter towards the colder zone) was reproduced by the help of a numerical model, including coupled-flow processes, mass conservation laws and hydro-thermo-mechanical changes (see Figure). A range of thermo-osmotic permeability (kT), between 6.10-12 and 2.10-10 m2.K-1.s-2, was obtained during the experiments depending on the temperature gradient and uncertainties on the model parameters. Values obtained for the Tournemire’s argillite are in the high range of thermo-osmotic permeabilities for argillaceous materials and suggest an effect of pore size on the thermo-osmotic permeability of a clay-rock (kT being higher with little pore size). Another dependence of thermo-osmotic permeability with temperature is observed, with kT decreasing when the temperature increases. These experiments and modeling indicate thermo-osmosis will have an influence on water flow in presence of a temperature gradient and this process is to consider in water flow studies in clay-rocks. Reference: Tremosa et al. Estimating thermo-osmotic coefficients in clay
Numerical simulations of cataclysmic floods: A coupling model of surface and subsurface flow.
NASA Astrophysics Data System (ADS)
Miyamoto, H.; Komatsu, G.; Ito, K.; Tosaka, H.; Tokunaga, T.
1999-09-01
The Martian outflow channels are considered to have been formed by catastrophic water flood processes analogous to the Lake Missoula floods responsible for the formation of the Channeled Scabland in Washington State [e.g., Baker and Milton, 1974]. The estimations of peak discharge rates and total amounts of water play critical roles for the discussion of the Martian water cycle. Therefore, for a more realistic estimation we developed a three-dimensional numerical code of surface flows coupled with subsurface flows. Coupling both surface and subsurface flows in the model is very important because the origin of the outflows is thought to be strongly linked to subsurface aquifers [e.g., Baker, et al., 1991]. Our model can calculate air-water movements on the surface and in the subsurface under various hydrological and geological conditions. We concentrated on the water movement as the first step. We applied our model to the glacial Lake Missoula floods to test our code and to study parameter sensitivities. We followed the glacial lake failure scenario and gave a well-constrained high discharge rate at an area of the lake failure. After the breakout, the flood water flows down-slope and covers a wide area. We compared the calculated areal coverage of the floods and the peak water levels with field observations. We obtained a reasonable result of the water coverage within a DTM containing the Cordilleran Ice Sheet. And also the computed time sequential behaviors of the floods, such as the hydraulic ponding in the Pasco Basin, are consistent with the field data. However, there are significant discrepancies in terms of the water depths between the calculated values and some field observations. This may indicate that the history of the floodings is more complex than our assumption.
MODELING COUPLED PROCESSES OF MULTIPHASE FLOW AND HEAT TRANSFER IN UNSATURATED FRACTURED ROCK
Y. Wu; S. Mukhopadhyay; K. Zhang; G.S. Bodvarsson
2006-02-28
A mountain-scale, thermal-hydrologic (TH) numerical model is developed for investigating unsaturated flow behavior in response to decay heat from the radioactive waste repository at Yucca Mountain, Nevada, USA. The TH model, consisting of three-dimensional (3-D) representations of the unsaturated zone, is based on the current repository design, drift layout, and thermal loading scenario under estimated current and future climate conditions. More specifically, the TH model implements the current geological framework and hydrogeological conceptual models, and incorporates the most updated, best-estimated input parameters. This mountain-scale TH model simulates the coupled TH processes related to mountain-scale multiphase fluid flow, and evaluates the impact of radioactive waste heat on the hydrogeological system, including thermally perturbed liquid saturation, gas- and liquid-phase fluxes, and water and rock temperature elevations, as well as the changes in water flux driven by evaporation/condensation processes and drainage between drifts. For a better description of the ambient geothermal condition of the unsaturated zone system, the TH model is first calibrated against measured borehole temperature data. The ambient temperature calibration provides the necessary surface and water table boundary as well as initial conditions. Then, the TH model is used to obtain scientific understanding of TH processes in the Yucca Mountain unsaturated zone under the designed schedule of repository thermal load.
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)
McDonald, R. R.; Nelson, J. M.
2006-12-01
Understanding the relationship between habitat requirements and the magnitude and duration of flow and sediment supply is an important component of both habitat assessment and restoration strategies. Fish habitat is often defined in terms of velocity, depth or substrate composition; locations where combinations of these factors meet habitat requirements depend on channel morphology, flow magnitude, and, in rivers with mobile beds, time-varying change in channel morphology. Because coupled multi-dimensional flow and sediment transport models provide spatially distributed information on flow and other hydraulic quantities, they permit detailed delineation of habitat. Furthermore, they can be used directly to understand how flow magnitude and duration and sediment supply control channel change and habitat availability. We present an example to illustrate how such models can be used in investigations of fish spawning habitat and availability. In the Kootenai River, Idaho, comparison of observed spawning locations with model derived spatial distributions of depths and velocities suggests that white sturgeon utilize the largest available velocity and depth within an 18-kilometer spawning reach over a range of discharges. This is a somewhat more selective criterion than a simple specification of a range of velocity or depth magnitudes, which illustrates the importance of evaluating habitat over a full range of discharge magnitudes. Observations also suggest that spawning currently occurs over a sandy substrate resulting in suffocation of eggs and little to no recruitment of juvenile sturgeon since closure of Libby Dam in 1974. Extending flow modeling to incorporate sediment- transport and bed evolution suggests that a relatively high magnitude long duration discharge can remove sandy substrate, thereby exposing a coarse gravel lag deposit in some areas and providing needed spawning substrate. These results were qualitatively validated through video surveys of channel
NASA Astrophysics Data System (ADS)
Céline, Longchamp; Irène, Manzella; Abellan, Antonio; Caspar, Olivier; Podladchikov, Yury; Jaboyedoff, Michel
2016-04-01
The objective of this research is to better understand the propagation of dry granular flows by coupling analogue modelling, structural analysis of deposit and numerical modelling, as follows: (a) The analogue modelling use laboratory experiments to investigate both the fluid-like flow of a granular mass falling down a slope and the influence of certain parameters such as the basal roughness, mobilized volume and slope angle. The experimental setup allows unconfined flow which spreads freely on a horizontal surface. Different grainsizes (115, 545 and 2605 μm) and substratum roughness (simulates by aluminium and sandpapers with grainsize from 16 to 425 μm) were used in order to understand their influence on the motion of a granular mass. During our experiments, the runout varied between 4.5cm and 11cm, with an increase of the basal roughness. When the volume varied between 300 and 600cm3 the runout was comprised between 9.2cm and 11.7cm. Finally when the slope angle was increased from 35° to 60°, the observed runout was between 5.3cm and 20cm. (b) Rock avalanche dynamic is analysed by means of a detailed structural analysis of the analogue modelling deposits. A series of 3D measurements were carried out on the deposit and a median filter and a gradient operator along the direction of propagation were applied to the 3D datasets. Treatment yield a more precise mapping of the longitudinal and transversal displacement features observed at the surface of the deposits. (c) The numerical modelling performed during this research is based both on continuum mechanics approach and on solving the shallow water equations. The avalanche was described from an Eulerian point of view within a continuum framework as single phase of incompressible granular material following Mohr-Coulomb friction law. The results obtained with the numerical model resemble those observed with the analogue modelling mentioned above. By coupling these three approaches, we obtained a complete scheme
NASA Astrophysics Data System (ADS)
Bürger, Raimund; Kumar, Sarvesh; Ruiz-Baier, Ricardo
2015-10-01
The sedimentation-consolidation and flow processes of a mixture of small particles dispersed in a viscous fluid at low Reynolds numbers can be described by a nonlinear transport equation for the solids concentration coupled with the Stokes problem written in terms of the mixture flow velocity and the pressure field. Here both the viscosity and the forcing term depend on the local solids concentration. A semi-discrete discontinuous finite volume element (DFVE) scheme is proposed for this model. The numerical method is constructed on a baseline finite element family of linear discontinuous elements for the approximation of velocity components and concentration field, whereas the pressure is approximated by piecewise constant elements. The unique solvability of both the nonlinear continuous problem and the semi-discrete DFVE scheme is discussed, and optimal convergence estimates in several spatial norms are derived. Properties of the model and the predicted space accuracy of the proposed formulation are illustrated by detailed numerical examples, including flows under gravity with changing direction, a secondary settling tank in an axisymmetric setting, and batch sedimentation in a tilted cylindrical vessel.
Coupled Modeling of Groundwater Flow and Land Subsidence with Secular Strain (Creep)
NASA Astrophysics Data System (ADS)
Bakr, M.
2012-12-01
Land subsidence limits sustainable development of many areas around the world. This is especially the case in low lying regions such as deltas which accommodate a significant percentage of the human population. Among the most common human-induced factors for land subsidence, is groundwater extractions. In these cases, groundwater flow and land subsidence are coupled processes, especially in basins with extensive spatial extent of soft soils (e.g. clay, peat). Creep (or secondary consolidation) is a land subsidence component that usually contributes to total land subsidence in soft soils. It leads to a reduction in void ratio at constant effective stress, and consequently, to the development of an apparent pre-consolidation pressure. The creep component has been usually ignored in the analysis of coupled groundwater flow and land subsidence. Here, the focus is the development of a coupled model of groundwater flow and land subsidence in porous media considering secular strain (creep). The Bjerrum method for settlement calculation (Bjerrum, 1967) due to change in effective stresses is coupled with MODFLOW to tackle the problem. In particular, the SUB-WT package of MODFLOW (Leake and Galloway, 2007) is modified where the Bjerrum method is used to calculate the primary and secondary consolidation due to change in effective stresses as a result of groundwater abstraction. The Bjerrum model is based on linear strains relationship. Usage of linear strains means that the model directly supports the common parameters Cr, Cc, Cα (i.e. re-compression, compression, and secondary compression indices; respectively). The Bjerrum model assumes that creep rate will reduce with increasing over-consolidation and that over-consolidation will grow by unloading and by ageing. To verify the coupled model, a hypothetical problem is considered where a simple hydrogeological system consisting of a shallow unconfined aquifer and a deeper confined aquifer separated by a (semi
Samper, Javier; Molinero, Jorg; Changbing, Yang; Zhang, Guoxiang
2003-12-01
The Redox Zone Experiment was carried out at the Aespoe HRL in order to study the redox behavior and the hydrochemistry of an isolated vertical fracture zone disturbed by the excavation of an access tunnel. Overall results and interpretation of the Redox Zone Project were reported by /Banwart et al, 1995/. Later, /Banwart et al, 1999/ presented a summary of the hydrochemistry of the Redox Zone Experiment. Coupled groundwater flow and reactive transport models of this experiment were carried out by /Molinero, 2000/ who proposed a revised conceptual model for the hydrogeology of the Redox Zone Experiment which could explain simultaneously measured drawdown and salinity data. The numerical model was found useful to understand the natural system. Several conclusions were drawn about the redox conditions of recharge waters, cation exchange capacity of the fracture zone and the role of mineral phases such as pyrite, calcite, hematite and goethite. This model could reproduce the measured trends of dissolved species, except for bicarbonate and sulfate which are affected by microbially-mediated processes. In order to explore the role of microbial processes, a coupled numerical model has been constructed which accounts for water flow, reactive transport and microbial processes. The results of this model is presented in this report. This model accounts for groundwater flow and reactive transport in a manner similar to that of /Molinero, 2000/ and extends the preliminary microbial model of /Zhang, 2001/ by accounting for microbially-driven organic matter fermentation and organic matter oxidation. This updated microbial model considers simultaneously the fermentation of particulate organic matter by yeast and the oxidation of dissolved organic matter, a product of fermentation. Dissolved organic matter is produced by yeast and serves also as a substrate for iron-reducing bacteria. Model results reproduce the observed increase in bicarbonate and sulfate concentration, thus
NASA Astrophysics Data System (ADS)
Frampton, Andrew; Destouni, Georgia
2016-04-01
In cold regions, flow in the unsaturated zone is highly dynamic with seasonal variability and changes in temperature, moisture, and heat and water fluxes, all of which affect ground freeze-thaw processes and influence transport of inert and reactive waterborne substances. In arctic permafrost environments, near-surface groundwater flow is further restricted to a relatively shallow and seasonally variable active layer, confined by perennially frozen ground below. The active layer is typically partially saturated with ice, liquid water and air, and is strongly dependent on seasonal temperature fluctuations, thermal forcing and infiltration patterns. Here there is a need for improved understanding of the mechanisms controlling subsurface solute transport in the partially saturated active layer zone. Studying solute transport in cold regions is relevant to improve the understanding of how natural and anthropogenic pollution may change as activities in arctic and sub-arctic regions increase. It is also particularly relevant for understanding how dissolved carbon is transported in coupled surface and subsurface hydrological systems under climate change, in order to better understand the permafrost-hydrological-carbon climate feedback. In this contribution subsurface solute transport under surface warming and degrading permafrost conditions is studied using a physically based model of coupled cryotic and hydrogeological flow processes combined with a particle tracking method. Changes in subsurface water flows and solute transport travel times are analysed for different modelled geological configurations during a 100-year warming period. Results show that for all simulated cases, the minimum and mean travel times increase non-linearly with warming irrespective of geological configuration and heterogeneity structure. The travel time changes are shown to depend on combined warming effects of increase in pathway length due to deepening of the active layer, reduced transport
Modeling transient streaming potentials in coupled saturated-unsaturated zone flow to a pumping well
NASA Astrophysics Data System (ADS)
Malama, B.
2012-12-01
A semi-analytical model for transient response of streaming potentials (SP) to pumping in an unconfined aquifer, taking into account unsaturated zone flow, is presented. Flow in the unsaturated zone is modeled with a linearized Richards' equation with the moisture retention curve and unsaturated hydraulic conduc- tivity assumed to be exponential functions of matric potential. For the case presented here, the same sorption number is assumed for moisture retention and unsaturated hydraulic conductivity. The ratio of the unsaturated to satu- rated electrokinetic coupling coefficient is described by Cℓ,r = kr Sw-(d+1), where kr is relative hydraulic conductivity, Sw is saturation, and d is Archie's second exponent. Hence, based on the assumption of an exponential moisture retention curve, Cℓ,r is also an exponential function of matric potential. Model predicted responses in the saturated and unsaturated zones are compared with measured SP responses to pumping in a bench-scale experimental setup that simulates a radially bounded unconfined aquifer. Potential for using SP responses to esti- mate saturated and unsaturated hydraulic parameters is evaluated. 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. This research is funded by WIPP programs administered by the Office of Environmental Management (EM) of the U.S Department of Energy.
NASA Astrophysics Data System (ADS)
Liu, Qingquan; Cheng, Yuanping; Zhou, Hongxing; Guo, Pinkun; An, Fenghua; Chen, Haidong
2015-05-01
The influence of gas diffusion behavior on gas flow and permeability evolution in coal seams is evaluated in this paper. Coalbed methane (CBM) reservoirs differ from conventional porous media and fractured gas reservoirs due to certain unique features, which lead to two distinct gas pressures: one in fractures and the other in the coal matrix. The latter pressure, also known as the sorption pressure, will be used in calculating sorption-based volume changes. The effective stress laws for single-porosity media is not suitable for CBM reservoirs, and the effective stress laws for multi-porosity media need to be applied. The realization of the above two points is based on the study of the two-phase state of gas migration (involving Fickian diffusion and Darcy flow) in a coal seam. Then, a general porosity and permeability model based on the P-M model is proposed to fit this phenomenon. Moreover, the Klinkenberg effect has been taken into account and set as a reference object. Finally, a coupled gas flow and coal deformation model is proposed and solved by using a finite element method. The numerical results indicate that the effects of gas diffusion behavior and Klinkenberg behavior can have a critical influence on the gas pressure, residual gas content, and permeability evolution during the entire methane degasification period, and the impacts of the two effects are of the same order of magnitude. Without considering the gas diffusion effect, the gas pressure and residual gas content will be underestimated, and the permeability will be overestimated.
Coupled Model for CO2 Leaks from Geological Storage: Geomechanics, Fluid Flow and Phase Transitions
NASA Astrophysics Data System (ADS)
Gor, G.; Prevost, J.
2013-12-01
Deep saline aquifers are considered as a promising option for long-term storage of carbon dioxide. However, risk of CO2 leakage from the aquifers through faults, natural or induced fractures or abandoned wells cannot be disregarded. Therefore, modeling of various leakage scenarios is crucial when selecting a site for CO2 sequestration and choosing proper operational conditions. Carbon dioxide is injected into wells at supercritical conditions (t > 31.04 C, P > 73.82 bar), and these conditions are maintained in the deep aquifers (at 1-2 km depth) due to hydrostatic pressure and geothermal gradient. However, if CO2 and brine start to migrate from the aquifer upward, both pressure and temperature will decrease, and at the depth of 500-750 m, the conditions for CO2 will become subcritical. At subcritical conditions, CO2 starts boiling and the character of the flow changes dramatically due to appearance of the third (vapor) phase and latent heat effects. When modeling CO2 leaks, one needs to couple the multiphase flow in porous media with geomechanics. These capabilities are provided by Dynaflow, a finite element analysis program [1]; Dynaflow has already showed to be efficient for modeling caprock failure causing CO2 leaks [2, 3]. Currently we have extended the capabilities of Dynaflow with the phase transition module, based on two-phase and three-phase isenthalpic flash calculations [4]. We have also developed and implemented an efficient method for solving heat and mass transport with the phase transition using our flash module. Therefore, we have developed a robust tool for modeling CO2 leaks. In the talk we will give a brief overview of our method and illustrate it with the results of simulations for characteristic test cases. References: [1] J.H. Prevost, DYNAFLOW: A Nonlinear Transient Finite Element Analysis Program. Department of Civil and Environmental Engineering, Princeton University, Princeton, NJ. http://www.princeton.edu/~dynaflow/ (last update 2013
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
Approaches to Modeling Coupled Flow and Reaction in a 2-D Cementation Experiment
Steefel, Carl; Cochepin, B.; Trotignon, L.; Bildstein, O.; Steefel, C.; Lagneau, V.; van der Lee, J.
2008-04-01
Porosity evolution at reactive interfaces is a key process that governs the evolution and performances of many engineered systems that have important applications in earth and environmental sciences. This is the case, for example, at the interface between cement structures and clays in deep geological nuclear waste disposals. Although in a different transport regime, similar questions arise for permeable reactive barriers used for biogeochemical remediation in surface environments. The COMEDIE project aims at investigating the coupling between transport, hydrodynamics and chemistry when significant variations of porosity occur. The present work focuses on a numerical benchmark used as a design exercise for the future COMEDIE-2D experiment. The use of reactive transport simulation tools like Hytec and Crunch provides predictions of the physico-chemical evolutions that are expected during the future experiments in laboratory. Focus is given in this paper on the evolution during the simulated experiment of precipitate, permeability and porosity fields. A first case is considered in which the porosity is constant. Results obtained with Crunch and Hytec are in relatively good agreement. Differences are attributable to the models of reactive surface area taken into account for dissolution/precipitation processes. Crunch and Hytec simulations taking into account porosity variations are then presented and compared. Results given by the two codes are in qualitative agreement, with differences attributable in part to the models of reactive surface area for dissolution/precipitation processes. As a consequence, the localization of secondary precipitates predicted by Crunch leads to lower local porosities than for predictions obtained by Hytec and thus to a stronger coupling between flow and chemistry. This benchmark highlights the importance of the surface area model employed to describe systems in which strong porosity variations occur as a result of dissolution
Coupled thermodynamic and two-phase flow modelling of partially melting crust
NASA Astrophysics Data System (ADS)
Riel, Nicolas; Bouilhol, Pierre; Magni, Valentina; van Hunen, Jeroen; Velic, Mirko
2016-04-01
How magmas are formed, transferred and interact in the lower crust to form mid-crust plutonic belts remain a fundamental question to understand the chemical and mechanical evolution of continents. To assess this question we developed a 2-D two-phase flow code using finite volume method. Our formulation takes into account: (i) an extended Darcy's law for fluid flow with first order temperature- and fluid-content dependency for the host-rock viscosity and silica-dependent viscosity for the fluid, (ii) the heat equation assuming thermal equilibrium for both solid and liquid and temperature-dependent diffusivity, (iii) thermodynamic modelling of stable phases via a dynamic coupling with Perple_X, and (iv) chemical advection of both the solid and liquid composition. To model chemical interactions with the host rock during magma transport, the melt is assumed to be either in thermodynamic equilibrium or in thermodynamic disequilibrium, or as function of these two endmembers. We applied our modelling approach to investigate the behaviour and composition of magma during lower crust melting. Our goal is to better understand the formation of felsic crust through melting, segregation and assimilation of lower crustal lithologies, applied to Archaean systems. Our preliminary results show the ascend of silica-rich magmas is slow, occurring on the timescale of millions of years, and is highly controlled by (i) the melting curve of the protolith and (ii) by its chemical degree of interaction with the host rock. The resulting transferred magmas are in good accordance with observed composition forming the grey gneisses of Archean terranes (i.e SiO2-rich > 62%, Mg# = 40-50, Na2O ~6%, MgO = 0.5-1%).
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.
NASA Astrophysics Data System (ADS)
Fiorentini, Marcello; Orlandini, Stefano; Paniconi, Claudio
2015-07-01
A process-based numerical model of integrated surface-subsurface flow is analyzed in order to identify, track, and reduce the mass balance errors affiliated with the model's coupling scheme. The sources of coupling error include a surface-subsurface grid interface that requires node-to-cell and cell-to-node interpolation of exchange fluxes and ponding heads, and a sequential iterative time matching procedure that includes a time lag in these same exchange terms. Based on numerical experiments carried out for two synthetic test cases and for a complex drainage basin in northern Italy, it is shown that the coupling mass balance error increases during the flood recession limb when the rate of change in the fluxes exchanged between the surface and subsurface is highest. A dimensionless index that quantifies the degree of coupling and a saturated area index are introduced to monitor the sensitivity of the model to coupling error. Error reduction is achieved through improvements to the heuristic procedure used to control and adapt the time step interval and to the interpolation algorithm used to pass exchange variables from nodes to cells. The analysis presented illustrates the trade-offs between a flexible description of surface and subsurface flow processes and the numerical errors inherent in sequential iterative coupling with staggered nodal points at the land surface interface, and it reveals mitigation strategies that are applicable to all integrated models sharing this coupling and discretization approach.
A Coupled Finite-Volume Model for 2-D Surface and 3-D Subsurface Flows
Technology Transfer Automated Retrieval System (TEKTRAN)
Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed; therefore, hydrologic modeling tools should consider these interactions to provide reliable predictions, especially during rainfall-runoff processes. This paper presents a fully implicit coupled...
NASA Astrophysics Data System (ADS)
Healy, R. W.
2015-12-01
Water movement through soils is often dominated by preferential flow processes such as fingering and macropore flow. Traditional models of flow in the unsaturated zone are based on the diffusion or Richards equation and assume that diffusive (surface-tension viscous) flow is the only flow process. These models are incapable of accurately simulating preferential flow. Several alternative approaches, including kinematic wave, transfer function, and water-content wave models, have been suggested for simulating water movement through preferential flow paths. The source-responsive model proposed by Nimmo (2010) and Nimmo and Mitchell (2013) is unique among such models in that water transfer from land surface to depth is controlled by the water-application rate at land surface. The source-responsive model has been coupled with a one-dimensional version of the Richards-equation based model of variably saturated flow, VS2DT. The new model, can simulate flow within the preferential (S) domain alone, within the diffuse (D) domain alone, or within both the S and D domains simultaneously. Water exchange between the two domains is treated as a first-order diffusive process, with the exchange coefficient being a function of soil-water content. The new model was used to simulate field and laboratory infiltration experiments described in the literature. Simulations were calibrated against measured soil water contents with the PEST parameter estimation package; values for hydraulic conductivity and 3 van Genuchten and 3 source-responsive parameters were optimized. Although exact matches between measured and simulated water contents were not obtained, the simulation results captured the salient characteristics of the published data sets, including features typical of preferential as well as diffusive flow. Results obtained from simulating flow simultaneously in both the S and D domain provided better matches to measured data than results obtained from simulating flow independently
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
Coupled reactive flow modeling with declining reactivity in fractured geothermal systems
NASA Astrophysics Data System (ADS)
Palguta, J.; Williams, C. F.; Ingebritsen, S.; Hickman, S.; Sonnenthal, E. L.
2011-12-01
Changes in permeability and fluid flow within geothermal systems are driven by geochemical reactions, advective and diffusive transport of heat and solute mass, and evolving thermal and mechanical environments. Representation of these simultaneous processes in numerical models is required for the characterization and simulation of natural geothermal systems. However, identifying and developing mathematical representations for all of the relevant mechanisms that control system behavior presents a major challenge. We have developed two-dimensional simulations of physical and chemical evolution in fractured granite under geothermal conditions with temperatures ranging from 150-300 °C. The goal of this study is to help identify possible sources for existing discrepancies between model results and laboratory-based measurements by adding a new mathematical formulation to the code TOUGHREACT. The revised code is designed to further quantify the link between the progressive evolution of reaction rates and alteration mineralogy. We explicitly couple reaction rates to mineral precipitation/dissolution effects by using an exponential function that defines evolving reactive surface areas in terms of each of the following (i) the accumulated total secondary mineral volume fraction, (ii) the accumulated clay (smectite) portion of the secondary mineral volume fraction, and (iii) the net change in mineral volume fraction (combined effects of dissolution and precipitation). We evaluate the performance of these three modified approaches by comparing simulation results to detailed laboratory measurements of fluid compositions, mineral abundances, and permeability changes in fractured Westerly granite and to previous one-dimensional simulations in which reactive surface areas were adjusted with time to match the observed fracture permeability history. The simulation results offer a potentially useful means of quantifying reactivity loss and of examining the extent to which secondary
NASA Astrophysics Data System (ADS)
Yeh, G. T.; Tsai, C. H.
2015-12-01
This paper presents the development of a THMC (thermal-hydrology-mechanics-chemistry) process model in variably saturated media. The governing equations for variably saturated flow and reactive chemical transport are obtained based on the mass conservation principle of species transport supplemented with Darcy's law, constraint of species concentration, equation of states, and constitutive law of K-S-P (Conductivity-Degree of Saturation-Capillary Pressure). The thermal transport equation is obtained based on the conservation of energy. The geo-mechanic displacement is obtained based on the assumption of equilibrium. Conventionally, these equations have been implicitly coupled via the calculations of secondary variables based on primary variables. The mechanisms of coupling have not been obvious. In this paper, governing equations are explicitly coupled for all primary variables. The coupling is accomplished via the storage coefficients, transporting velocities, and conduction-dispersion-diffusion coefficient tensor; one set each for every primary variable. With this new system of equations, the coupling mechanisms become clear. Physical interpretations of every term in the coupled equations will be discussed. Examples will be employed to demonstrate the intuition and superiority of these explicit coupling approaches. Keywords: Variably Saturated Flow, Thermal Transport, Geo-mechanics, Reactive Transport.
Swain, Eric D.; Wexler, Eliezer J.
1996-01-01
Ground-water and surface-water flow models traditionally have been developed separately, with interaction between subsurface flow and streamflow either not simulated at all or accounted for by simple formulations. In areas with dynamic and hydraulically well-connected ground-water and surface-water systems, stream-aquifer interaction should be simulated using deterministic responses of both systems coupled at the stream-aquifer interface. Accordingly, a new coupled ground-water and surface-water model was developed by combining the U.S. Geological Survey models MODFLOW and BRANCH; the interfacing code is referred to as MODBRANCH. MODFLOW is the widely used modular three-dimensional, finite-difference ground-water model, and BRANCH is a one-dimensional numerical model commonly used to simulate unsteady flow in open- channel networks. MODFLOW was originally written with the River package, which calculates leakage between the aquifer and stream, assuming that the stream's stage remains constant during one model stress period. A simple streamflow routing model has been added to MODFLOW, but is limited to steady flow in rectangular, prismatic channels. To overcome these limitations, the BRANCH model, which simulates unsteady, nonuniform flow by solving the St. Venant equations, was restructured and incorporated into MODFLOW. Terms that describe leakage between stream and aquifer as a function of streambed conductance and differences in aquifer and stream stage were added to the continuity equation in BRANCH. Thus, leakage between the aquifer and stream can be calculated separately in each model, or leakages calculated in BRANCH can be used in MODFLOW. Total mass in the coupled models is accounted for and conserved. The BRANCH model calculates new stream stages for each time interval in a transient simulation based on upstream boundary conditions, stream properties, and initial estimates of aquifer heads. Next, aquifer heads are calculated in MODFLOW based on stream
Implementation and use of direct-flow connections in a coupled ground-water and surface-water model
Swain, Eric D.
1994-01-01
The U.S. Geological Survey's MODFLOW finite-difference ground-water flow model has been coupled with three surface-water packages - the MODBRANCH, River, and Stream packages - to simulate surface water and its interaction with ground water. Prior to the development of the coupling packages, the only interaction between these modeling packages was that leakage values could be passed between MODFLOW and the three surface-water packages. To facilitate wider and more flexible uses of the models, a computer program was developed and added to MODFLOW to allow direct flows or stages to be passed between any of the packages and MODFLOW. The flows or stages calculated in one package can be set as boundary discharges or stages to be used in another package. Several modeling packages can be used in the same simulation depending upon the level of sophistication needed in the various reaches being modeled. This computer program is especially useful when any of the River, Stream, or MODBRANCH packages are used to model a river flowing directly into or out of wetlands in direct connection with the aquifer and represented in the model as an aquifer block. A field case study is shown to illustrate an application.
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.
NASA Astrophysics Data System (ADS)
Ryu, Jaiyoung; Hu, Xiao; Shadden, Shawn C.
2014-11-01
The cerebral circulation is unique in its ability to maintain blood flow to the brain under widely varying physiologic conditions. Incorporating this autoregulatory response is critical to cerebral blood flow modeling, as well as investigations into pathological conditions. We discuss a one-dimensional nonlinear model of blood flow in the cerebral arteries that includes coupling of autoregulatory lumped parameter networks. The model is tested to reproduce a common clinical test to assess autoregulatory function - the carotid artery compression test. The change in the flow velocity at the middle cerebral artery (MCA) during carotid compression and release demonstrated strong agreement with published measurements. The model is then used to investigate vasospasm of the MCA, a common clinical concern following subarachnoid hemorrhage. Vasospasm was modeled by prescribing vessel area reduction in the middle portion of the MCA. Our model showed similar increases in velocity for moderate vasospasms, however, for serious vasospasm (~ 90% area reduction), the blood flow velocity demonstrated decrease due to blood flow rerouting. This demonstrates a potentially important phenomenon, which otherwise would lead to false-negative decisions on clinical vasospasm if not properly anticipated.
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
Ashwin, J.; Ganesh, R.
2010-10-15
Using a generalized hydrodynamic (GH) model, the growth rate spectra of Kelvin-Helmholtz (KH) instability has been obtained analytically for a step shear profile in strongly coupled Yukawa liquids. The class of shear flows studied is assumed to be incompressible in nature. The growth rate spectra calculated exhibit viscous damping at high mode numbers, destabilization at stronger coupling, and in the limit {tau}{sub m} (viscoelastic relaxation time){yields}0, reduce to the regular Navier-Stokes growth rate spectra. A direct comparison is made with previous molecular dynamics (MD) simulations [Ashwin J. and R. Ganesh, Phys. Rev. Lett. 104, 215003 (2010)] of KH instability. We find that for a given value of Reynolds number R and coupling parameter 1<{Gamma}<100, the GH and MD growth rates are in a qualitative agreement. The inclusion of the effect of shear heating as an effective coupling parameter {Gamma}{sub e} appears to improve the quantitative comparison as well.
Zhou, Jianzhong; Song, Lixiang; Kursan, Suncana; Liu, Yi
2015-05-01
A two-dimensional coupled water quality model is developed for modeling the flow-mass transport in shallow water. To simulate shallow flows on complex topography with wetting and drying, an unstructured grid, well-balanced, finite volume algorithm is proposed for numerical resolution of a modified formulation of two-dimensional shallow water equations. The slope-limited linear reconstruction method is used to achieve second-order accuracy in space. The algorithm adopts a HLLC-based integrated solver to compute the flow and mass transport fluxes simultaneously, and uses Hancock's predictor-corrector scheme for efficient time stepping as well as second-order temporal accuracy. The continuity and momentum equations are updated in both wet and dry cells. A new hybrid method, which can preserve the well-balanced property of the algorithm for simulations involving flooding and recession, is proposed for bed slope terms approximation. The effectiveness and robustness of the proposed algorithm are validated by the reasonable good agreement between numerical and reference results of several benchmark test cases. Results show that the proposed coupled flow-mass transport model can simulate complex flows and mass transport in shallow water. PMID:25686488
NASA Astrophysics Data System (ADS)
Liu, Xin; Mohammadian, Abdolmajid; Kurganov, Alexander; Infante Sedano, Julio Angel
2015-11-01
Intense sediment transport and rapid bed evolution are frequently observed under highly-energetic flows, and bed erosion sometimes is of the same magnitude as the flow itself. Simultaneous simulation of multiple physical processes requires a fully coupled system to achieve an accurate hydraulic and morphodynamical prediction. In this paper, we develop a high-order well-balanced finite-volume method for a new fully coupled two-dimensional hyperbolic system consisting of the shallow water equations with friction terms coupled with the equations modeling the sediment transport and bed evolution. The nonequilibrium sediment transport equation is used to predict the sediment concentration variation. Since bed-load, sediment entrainment and deposition have significant effects on the bed evolution, an Exner-based equation is adopted together with the Grass bed-load formula and sediment entrainment and deposition models to calculate the morphological process. The resulting 5 × 5 hyperbolic system of balance laws is numerically solved using a Godunov-type central-upwind scheme on a triangular grid. A computationally expensive process of finding all of the eigenvalues of the Jacobian matrices is avoided: The upper/lower bounds on the largest/smallest local speeds of propagation are estimated using the Lagrange theorem. A special discretization of the bed-slope term is proposed to guarantee the well-balanced property of the designed scheme. The proposed fully coupled model is verified on a number of numerical experiments.
Fu, P; Johnson, S M; Hao, Y; Carrigan, C R
2011-01-18
The primary objective of our current research is to develop a computational test bed for evaluating borehole techniques to enhance fluid flow and heat transfer in enhanced geothermal systems (EGS). Simulating processes resulting in hydraulic fracturing and/or the remobilization of existing fractures, especially the interaction between propagating fractures and existing fractures, represents a critical goal of our project. To this end, we are continuing to develop a hydraulic fracturing simulation capability within the Livermore Distinct Element Code (LDEC), a combined FEM/DEM analysis code with explicit solid-fluid mechanics coupling. LDEC simulations start from an initial fracture distribution which can be stochastically generated or upscaled from the statistics of an actual fracture distribution. During the hydraulic stimulation process, LDEC tracks the propagation of fractures and other modifications to the fracture system. The output is transferred to the Non-isothermal Unsaturated Flow and Transport (NUFT) code to capture heat transfer and flow at the reservoir scale. This approach is intended to offer flexibility in the types of analyses we can perform, including evaluating the effects of different system heterogeneities on the heat extraction rate as well as seismicity associated with geothermal operations. This paper details the basic methodology of our approach. Two numerical examples showing the capability and effectiveness of our simulator are also presented.
Validation of a coupled wave-flow model in a high-energy setting: the mouth of the Columbia River
Elias, Edwin P.L.; Gelfenbaum, Guy R.; Andre van de Wethusen
2012-01-01
A monthlong time series of wave, current, salinity, and suspended-sediment measurements was made at five sites on a transect across the Mouth of Columbia River (MCR). These data were used to calibrate and evaluate the performance of a coupled hydrodynamic and wave model for the MCR based on the Delft3D modeling system. The MCR is a dynamic estuary inlet in which tidal currents, river discharge, and wave-driven currents are all important. Model tuning consisted primarily of spatial adjustments to bottom drag coefficients. In combination with (near-) default parameter settings, the MCR model application is able to simulate the dominant features in the tidal flow, salinity and wavefields observed in field measurements. The wave-orbital averaged method for representing the current velocity profile in the wave model is considered the most realistic for the MCR. The hydrodynamic model is particularly effective in reproducing the observed vertical residual and temporal variations in current structure. Density gradients introduce the observed and modeled reversal of the mean flow at the bed and augment mean and peak flow in the upper half of the water column. This implies that sediment transport during calmer summer conditions is controlled by density stratification and is likely net landward due to the reversal of flow near the bed. The correspondence between observed and modeled hydrodynamics makes this application a tool to investigate hydrodynamics and associated sediment transport.
Resolved granular debris-flow simulations with a coupled SPH-DCDEM model
NASA Astrophysics Data System (ADS)
Birjukovs Canelas, Ricardo; Domínguez, José M.; Crespo, Alejandro J. C.; Gómez-Gesteira, Moncho; Ferreira, Rui M. L.
2016-04-01
Debris flows represent some of the most relevant phenomena in geomorphological events. Due to the potential destructiveness of such flows, they are the target of a vast amount of research (Takahashi, 2007 and references therein). A complete description of the internal processes of a debris-flow is however still an elusive achievement, explained by the difficulty of accurately measuring important quantities in these flows and developing a comprehensive, generalized theoretical framework capable of describing them. This work addresses the need for a numerical model applicable to granular-fluid mixtures featuring high spatial and temporal resolution, thus capable of resolving the motion of individual particles, including all interparticle contacts. This corresponds to a brute-force approach: by applying simple interaction laws at local scales the macro-scale properties of the flow should be recovered by upscaling. This methodology effectively bypasses the complexity of modelling the intermediate scales by resolving them directly. The only caveat is the need of high performance computing, a demanding but engaging research challenge. The DualSPHysics meshless numerical implementation, based on Smoothed Particle Hydrodynamics (SPH), is expanded with a Distributed Contact Discrete Element Method (DCDEM) in order to explicitly solve the fluid and the solid phase. The model numerically solves the Navier-Stokes and continuity equations for the liquid phase and Newton's motion equations for solid bodies. The interactions between solids are modelled with classical DEM approaches (Kruggel-Emden et al, 2007). Among other validation tests, an experimental set-up for stony debris flows in a slit check dam is reproduced numerically, where solid material is introduced trough a hopper assuring a constant solid discharge for the considered time interval. With each sediment particle undergoing tens of possible contacts, several thousand time-evolving contacts are efficiently treated
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.
Coupling fluvial-hydraulic models to predict gravel transport in spatially variable flows
NASA Astrophysics Data System (ADS)
Segura, Catalina; Pitlick, John
2015-05-01
This study investigated spatial-temporal variations of shear stress and bed load transport at three gravel bed river reaches of the Williams Fork River, Colorado. A two-dimensional flow model was used to compute spatial distributions of shear stress (τ) for four discharge levels between one third of bankfull (Qbf) and Qbf. Results indicate that mean τ values are highly variable among sites. However, the properties of the mean-normalized distributions of τ are similar across sites for all flows. The distributions of τ are then used with a transport function to compute bed load transport rates of individual grain size fractions. Probability distributions of the instantaneous unit-width transport rates, qb, indicate that most of the bed load is transported through small portions of the bed with high τ. The mean-normalized probability distributions of qb are different among sites for all flows except at Qbf, when the distributions overlap. We also find that the grain size distribution (GSD) of the bed load adjusts with discharge to resemble the grain size distribution of the subsurface at Qbf. We extend these results to 13 locations in the basin, using the mean-normalized distributions of shear stress and measured subsurface grain sizes to compute bed load transport rates at Qbf. We found a remarkably similar shape of the qb distribution among sites highlighting the basin-wide balance between flow forces and GSD at Qbf and the potential to predict sediment flux at the watershed scale.
NASA Technical Reports Server (NTRS)
Dvorak, Frank A.; Dash, Sanford M.
1987-01-01
Work currently in progress to update an existing transonic circulation control airfoil analysis method is described. Existing methods suffer from two dificiencies: the inability to predict the shock structure of the underexpanded supersonic jets; and the insensitivity of the calculation to small changes in the Coanda surface geometry. A method developed for the analysis of jet exhaust plumes in supersonic flow is being modified for the case of the underexpanded wall jet. In the subsonic case, the same wall jet model was modified to include the calculation of the normal pressure gradient. This model is currently being coupled with the transonic circulation control airfoil analysis.
Markstrom, Steven L.; Niswonger, Richard G.; Regan, R. Steven; Prudic, David E.; Barlow, Paul M.
2008-01-01
The need to assess the effects of variability in climate, biota, geology, and human activities on water availability and flow requires the development of models that couple two or more components of the hydrologic cycle. An integrated hydrologic model called GSFLOW (Ground-water and Surface-water FLOW) was developed to simulate coupled ground-water and surface-water resources. The new model is based on the integration of the U.S. Geological Survey Precipitation-Runoff Modeling System (PRMS) and the U.S. Geological Survey Modular Ground-Water Flow Model (MODFLOW). Additional model components were developed, and existing components were modified, to facilitate integration of the models. Methods were developed to route flow among the PRMS Hydrologic Response Units (HRUs) and between the HRUs and the MODFLOW finite-difference cells. This report describes the organization, concepts, design, and mathematical formulation of all GSFLOW model components. An important aspect of the integrated model design is its ability to conserve water mass and to provide comprehensive water budgets for a location of interest. This report includes descriptions of how water budgets are calculated for the integrated model and for individual model components. GSFLOW provides a robust modeling system for simulating flow through the hydrologic cycle, while allowing for future enhancements to incorporate other simulation techniques.
Coupled reactive mass transport and fluid flow: Issues in model verification
Freedman, Vicky L.; Ibaraki, Motomu
2003-01-03
Model verification and validation are both important steps in the development of reactive transport models. In this paper, a distinction is made between verification and validation, and the focus is on codifying the issues of verification for a numerical, reactive transport flow model. First, the conceptual basis of model verification is reviewed, which shows that verification should be understood as a first step in model development, and be followed by a protocol that assures that the model accurately represents system behavior. Second, commonly used procedures and methods of model verification are presented. In the third part of this paper, an intercomparison of models is used to demonstrate that model verification can be performed despite differences in hydrogeochemical transport code formulations. Results of an example simulation of transport are presented in which the numerical model is tested against other hydrogeochemical codes. Different kinetic formulations between solid and aqueous phases used among numerical models complicates model verification. This test problem involves uranium transport under conditions of varying pH and oxidation potential, with reversible precipitation of calcium uranate and coffinite. Results between the different hydrogeochemical transport codes show differences in oxidation potentials, but similarities in mineral assemblages and aqueous transport patterns. Because model verification can be further complicated by differences in the approach for solving redox problems, a comparison of a fugacity approach to both the external approach (based on hypothetical electron activity) and effective internal approach (based on conservation of electrons) is performed. The comparison demonstrates that the oxygen fugacity approach produces different redox potentials and mineral assemblages than both the effective internal and external approaches.
Freeze, G.A.; Larson, K.W.; Davies, P.B.; Webb, S.W.
1995-10-01
Long-term repository assessment must consider the processes of (1) gas generation, (2) room closure and expansions due to salt creep, and (3) multiphase (brine and gas) fluid flow, as well as the complex coupling between these three processes. The mechanical creep closure code SANCHO was used to simulate the closure of a single, perfectly sealed disposal room filled with water and backfill. SANCHO uses constitutive models to describe salt creep, waste consolidation, and backfill consolidation, Five different gas-generation rate histories were simulated, differentiated by a rate multiplier, f, which ranged from 0.0 (no gas generation) to 1.0 (expected gas generation under brine-dominated conditions). The results of the SANCHO f-series simulations provide a relationship between gas generation, room closure, and room pressure for a perfectly sealed room. Several methods for coupling this relationship with multiphase fluid flow into and out of a room were examined. Two of the methods are described.
NASA Astrophysics Data System (ADS)
Ke, Xinyou; Alexander, J. Iwan D.; Prahl, Joseph M.; Savinell, Robert F.
2015-08-01
A simple analytical model of a layered system comprised of a single passage of a serpentine flow channel and a parallel underlying porous electrode (or porous layer) is proposed. This analytical model is derived from Navier-Stokes motion in the flow channel and Darcy-Brinkman model in the porous layer. The continuities of flow velocity and normal stress are applied at the interface between the flow channel and the porous layer. The effects of the inlet volumetric flow rate, thickness of the flow channel and thickness of a typical carbon fiber paper porous layer on the volumetric flow rate within this porous layer are studied. The maximum current density based on the electrolyte volumetric flow rate is predicted, and found to be consistent with reported numerical simulation. It is found that, for a mean inlet flow velocity of 33.3 cm s-1, the analytical maximum current density is estimated to be 377 mA cm-2, which compares favorably with experimental result reported by others of ∼400 mA cm-2.
Coupling actin flow, adhesion, and morphology in a computational cell motility model
NASA Astrophysics Data System (ADS)
Levine, Herbert
2014-03-01
Eukaryotic cells crawl by means of the coordinated spatiotemporal dynamics of an active polymer gel, consisting of actin, myosin and regulators thereof. Motility is necessarily coupled to shape, as the force generating mechanisms such as polymerization-based protrusions interact with the elasticity of the cell membrane and thereby determine the cell morphology. We have introduced a ``phase-field'' model of crawling cells, utilizing a mathematical approach originally developed for morphology problems arising in the field of liquid-solid phase transitions. Our model can be used to explain the pattern of traction forces applied to the substrate as well as some recent observations concerning oscillatory instabilities of cells moving on one-dimensional fiber tracks.
NASA Astrophysics Data System (ADS)
Seroussi, H. L.; Rignot, E. J.; Morlighem, M.; Larour, E. Y.; Ben Dhia, H.; Aubry, D.
2010-12-01
The recent development of new higher-order, higher-resolution ice sheet models has shown that sophisticated models, such as Full-Stokes, were essential in some parts of the ice sheets, including the grounding line region. These areas are crucial for ice flow projections and can only be rigorously simulated using full 3d models. Higher-order models are well-suited to ice stream dynamics, whereas the shallow-shelf approximation is sufficient for modeling ice shelf flow. Higher-order and full-Stokes model are computationally intensive and prohibitive for large-scale modeling. There is therefore a strong need to combine such different models in order to balance computational cost and physical accuracy for the whole ice sheet. Here we present a new methodology adapted from the Arlequin framework to couple finite element shelfy-stream, higher-order and Full-Stokes models. We achieve this by strongly coupling the different approximations within the same large scale simulation. This technique is applied to the Greenland ice sheet, and compared with single-model approaches. Our new method preserves the conditioning number of the stiffness matrix, and ensures seamless stress regimes across model transition zones, hence improving numerical accuracy compared to existing techniques that use penalties or kinematical constrains. Furthermore, it optimizes the number of degrees of freedom leading to reduced computational cost. This work was performed at the California Institute of Technology's Jet Propulsion Laboratory under a contract with the National Aeronautics and Space Administration's Modeling, Analysis and Prediction (MAP) Program.
SABAE-HW3D: a Meteor-Hydrological Model Coupling the Land Surface to Groundwater Flow
NASA Astrophysics Data System (ADS)
Loukili, Y.; Woodbury, A. D.
2007-12-01
A new coupled model linking the land surface scheme SABAE-HW with the saturated groundwater flow is introduced. SABAE-HW stands for: Soil Atmosphere Boundary, Accurate Evaluations of Heat and Water. It was presented (Loukili et al., 2007) as a parallel extension to the Canadian LAnd Surface Scheme (CLASS) which permits a free choice of soil column depth and layers. This development is important as it allows for the coupling of the atmosphere, land surface and subsurface zones. The generic domain under study is subdivided into soil columns whose surface areas represent GCM grid squares. The depths of columns extend to the saturated zone where a water table lower boundary condition is prescribed. In each column the unsaturated flow is forced by the corresponding meteorological data, and the moisture (liquid and frozen) and temperature profiles are computed by SABAE-HW using the half hour time step. The horizontal flow in the saturated zone is described by the vertically integrated model incorporating a storage term, and discretized using a finite volume scheme operating on the same GCM quadrilateral mesh. Since the saturated flow model affords larger time steps, bottom drainages from soil columns are summed up and input as cell recharges. Meanwhile, as the water table fluctuates, the individual column's mesh is allowed to deform. This physically based coupling strategy was selected for its superiority in providing stable and consistent results. In fact, the convergence to steady state situations in both the unsaturated and saturated zones is realized in many validation numerical tests. Moreover, our SABAE-HW3D code was successfully benchmarked against the finite element code Seep/W (Geo-Slope International, 2002) for steady and transient groundwater flows through different soil types. Even when handling coarse grids, SABAE-HW3D solutions are free of moisture oscillation and under or overshoot near the water table. The model is used to understand and assess properly
NASA Astrophysics Data System (ADS)
Yoshida, Hiroaki; Kinjo, Tomoyuki; Washizu, Hitoshi
2014-10-01
We present a coupled lattice Boltzmann method (LBM) to solve a set of model equations for electrokinetic flows in micro-/nano-channels. The model consists of the Poisson equation for the electrical potential, the Nernst-Planck equation for the ion concentration, and the Navier-Stokes equation for the flows of the electrolyte solution. In the proposed LBM, the electrochemical migration and the convection of the electrolyte solution contributing to the ion flux are incorporated into the collision operator, which maintains the locality of the algorithm inherent to the original LBM. Furthermore, the Neumann-type boundary condition at the solid/liquid interface is then correctly imposed. In order to validate the present LBM, we consider an electro-osmotic flow in a slit between two charged infinite parallel plates, and the results of LBM computation are compared to the analytical solutions. Good agreement is obtained in the parameter range considered herein, including the case in which the nonlinearity of the Poisson equation due to the large potential variation manifests itself. We also apply the method to a two-dimensional problem of a finite-length microchannel with an entry and an exit. The steady state, as well as the transient behavior, of the electro-osmotic flow induced in the microchannel is investigated. It is shown that, although no external pressure difference is imposed, the presence of the entry and exit results in the occurrence of the local pressure gradient that causes a flow resistance reducing the magnitude of the electro-osmotic flow.
NASA Astrophysics Data System (ADS)
Kim, Ho Jun; Lee, Hae June
2016-06-01
The wide applicability of capacitively coupled plasma (CCP) deposition has increased the interest in developing comprehensive numerical models, but CCP imposes a tremendous computational cost when conducting a transient analysis in a three-dimensional (3D) model which reflects the real geometry of reactors. In particular, the detailed flow features of reactive gases induced by 3D geometric effects need to be considered for the precise calculation of radical distribution of reactive species. Thus, an alternative inclusive method for the numerical simulation of CCP deposition is proposed to simulate a two-dimensional (2D) CCP model based on the 3D gas flow results by simulating flow, temperature, and species fields in a 3D space at first without calculating the plasma chemistry. A numerical study of a cylindrical showerhead-electrode CCP reactor was conducted for particular cases of SiH4/NH3/N2/He gas mixture to deposit a hydrogenated silicon nitride (SiN x H y ) film. The proposed methodology produces numerical results for a 300 mm wafer deposition reactor which agree very well with the deposition rate profile measured experimentally along the wafer radius.
Adamovich, Igor V.
2014-04-15
A three-dimensional, nonperturbative, semiclassical analytic model of vibrational energy transfer in collisions between a rotating diatomic molecule and an atom, and between two rotating diatomic molecules (Forced Harmonic Oscillator–Free Rotation model) has been extended to incorporate rotational relaxation and coupling between vibrational, translational, and rotational energy transfer. The model is based on analysis of semiclassical trajectories of rotating molecules interacting by a repulsive exponential atom-to-atom potential. The model predictions are compared with the results of three-dimensional close-coupled semiclassical trajectory calculations using the same potential energy surface. The comparison demonstrates good agreement between analytic and numerical probabilities of rotational and vibrational energy transfer processes, over a wide range of total collision energies, rotational energies, and impact parameter. The model predicts probabilities of single-quantum and multi-quantum vibrational-rotational transitions and is applicable up to very high collision energies and quantum numbers. Closed-form analytic expressions for these transition probabilities lend themselves to straightforward incorporation into DSMC nonequilibrium flow codes.
NASA Astrophysics Data System (ADS)
Smith, M. A.; Bryce, J.; Davis, J.
2007-12-01
2006 show a decrease in 87Sr/86Sr in the downstream direction as the water generally flows from the magmatic series to the gneiss complex, indicating that low 87Sr/86Sr groundwater is infiltrating the river. These results also indicate groundwater inputs from extrabasinal flow. Further refinements of the model, coupled with forthcoming trace element analyses, will provide a stronger means to quantify these contributions along the gradient of the river.
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
Pedestrian flow dynamics in a lattice gas model coupled with an evolutionary game.
Hao, Qing-Yi; Jiang, Rui; Hu, Mao-Bin; Jia, Bin; Wu, Qing-Song
2011-09-01
This paper studies unidirectional pedestrian flow by using a lattice gas model with parallel update rules. Game theory is introduced to deal with conflicts that two or three pedestrians want to move into the same site. Pedestrians are either cooperators or defectors. The cooperators are gentle and the defectors are aggressive. Moreover, pedestrians could change their strategy. The fundamental diagram and the cooperator fraction at different system width W have been investigated in detail. It is found that a two-lane system exhibits a first-order phase transition while a multilane system does not. A microscopic mechanism behind the transition has been provided. Mean-field analysis is carried out to calculate the critical density of the transition as well as the probability of games at large value of W. The spatial distribution of pedestrians is investigated, which is found to be dependent (independent) on the initial cooperator fraction when W is small (large). Finally, the influence of the evolutionary game rule has been discussed. PMID:22060456
NASA Astrophysics Data System (ADS)
Yang, Jubiao; Krane, Michael; Zhang, Lucy
2013-11-01
Vocal fold vibrations and the glottal jet are successfully simulated using the modified Immersed Finite Element method (mIFEM), a fully coupled dynamics approach to model fluid-structure interactions. A self-sustained and steady vocal fold vibration is captured given a constant pressure input at the glottal entrance. The flow rates at different axial locations in the glottis are calculated, showing small variations among them due to the vocal fold motion and deformation. To further facilitate the understanding of the phonation process, two control volume analyses, specifically with Bernoulli's equation and Newton's 2nd law, are carried out for the glottal flow based on the simulation results. A generalized Bernoulli's equation is derived to interpret the correlations between the velocity and pressure temporally and spatially along the center line which is a streamline using a half-space model with symmetry boundary condition. A specialized Newton's 2nd law equation is developed and divided into terms to help understand the driving mechanism of the glottal flow.
LI, MING-HSU; SIEGEL, MALCOLM D.; YEH, GOUR-TSYH
1999-09-20
The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically-based numerical model for simulation of coupled fluid flow and reactive chemical transport including both fast and slow reactions invariably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation-dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.
NASA Astrophysics Data System (ADS)
Muñoz-Cobo, José; Chiva, Sergio; El Aziz Essa, Mohamed; Mendes, Santos
2012-08-01
Two phase flow experiments with different superficial velocities of gas and water were performed in a vertical upward isothermal cocurrent air-water flow column with conditions ranging from bubbly flow, with very low void fraction, to transition flow with some cap and slug bubbles and void fractions around 25%. The superficial velocities of the liquid and the gas phases were varied from 0.5 to 3 m/s and from 0 to 0.6 m/s, respectively. Also to check the effect of changing the surface tension on the previous experiments small amounts of 1-butanol were added to the water. These amounts range from 9 to 75 ppm and change the surface tension. This study is interesting because in real cases the surface tension of the water diminishes with temperature, and with this kind of experiments we can study indirectly the effect of changing the temperature on the void fraction distribution. The following axial and radial distributions were measured in all these experiments: void fraction, interfacial area concentration, interfacial velocity, Sauter mean diameter and turbulence intensity. The range of values of the gas superficial velocities in these experiments covered the range from bubbly flow to the transition to cap/slug flow. Also with transition flow conditions we distinguish two groups of bubbles in the experiments, the small spherical bubbles and the cap/slug bubbles. Special interest was devoted to the transition region from bubbly to cap/slug flow; the goal was to understand the physical phenomena that take place during this transition A set of numerical simulations of some of these experiments for bubbly flow conditions has been performed by coupling a Lagrangian code, that tracks the three dimensional motion of the individual bubbles in cylindrical coordinates inside the field of the carrier liquid, to an Eulerian model that computes the magnitudes of continuous phase and to a 3D random walk model that takes on account the fluctuation in the velocity field of the
Inverse methods-based estimation of plate coupling in a plate motion model governed by mantle flow
NASA Astrophysics Data System (ADS)
Ratnaswamy, V.; Stadler, G.; Gurnis, M.
2013-12-01
Plate motion is primarily controlled by buoyancy (slab pull) which occurs at convergent plate margins where oceanic plates undergo deformation near the seismogenic zone. Yielding within subducting plates, lateral variations in viscosity, and the strength of seismic coupling between plate margins likely have an important control on plate motion. Here, we wish to infer the inter-plate coupling for different subduction zones, and develop a method for inferring it as a PDE-constrained optimization problem, where the cost functional is the misfit in plate velocities and is constrained by the nonlinear Stokes equation. The inverse models have well resolved slabs, plates, and plate margins in addition to a power law rheology with yielding in the upper mantle. Additionally, a Newton method is used to solve the nonlinear Stokes equation with viscosity bounds. We infer plate boundary strength using an inexact Gauss-Newton method with line search for backtracking. Each inverse model is applied to two simple 2-D scenarios (each with three subduction zones), one with back-arc spreading and one without. For each case we examine the sensitivity of the inversion to the amount of surface velocity used: 1) full surface velocity data and 2) surface velocity data simplified using a single scalar average (2-D equivalent to an Euler pole) for each plate. We can recover plate boundary strength in each case, even in the presence of highly nonlinear flow with extreme variations in viscosity. Additionally, we ascribe an uncertainty in each plate's velocity and perform an uncertainty quantification (UQ) through the Hessian of the misfit in plate velocities. We find that as plate boundaries become strongly coupled, the uncertainty in the inferred plate boundary strength decreases. For very weak, uncoupled subduction zones, the uncertainty of inferred plate margin strength increases since there is little sensitivity between plate margin strength and plate velocity. This result is significant
Essaid, H.I.
1986-01-01
A quasi-three dimensional finite difference model which simulates coupled, fresh water and salt water flow, separated by a sharp interface, is used to investigate the effects of storage characteristics, transmissivity, boundary conditions and anisotropy on the transient responses of such flow systems. The magnitude and duration of the departure of aquifer response from the behavior predicted using the Ghyben-Herzberg, one-fluid approach is a function of the ease with which flow can be induced in the salt water region. In many common hydrogeologic settings short-term fresh water head responses, and transitional responses between short-term and long-term, can only be realistically reproduced by including the effects of salt water flow on the dynamics of coastal flow systems. The coupled fresh water-salt water flow modeling approach is able to reproduce the observed annual fresh water head response of the Waialae aquifer of southeastern Oahu, Hawaii. ?? 1986.
Numerical simulation of the stokes wave for the flow around a ship hull coupled with the VOF model
NASA Astrophysics Data System (ADS)
Shengtao, Chen; Jingjun, Zhong; Peng, Sun
2015-06-01
The surface wave generated by flow around a ship hull moving near free surface of water is simulated numerically in this study. The three-dimensional implicit finite volume method (FVM) is applied to solve Reynolds averaged Navier-Stokes (RANS) equation. The realizable k-e turbulence model has been implemented to capture turbulent flow around the ship hull in the free surface zone. The volume of fluid (VOF) method coupled with the Stokes wave theory has been used to determine the free surface effect of water. By using is a six degrees of freedom model, the ship hull's movement is numerically solved with the Stokes wave together. Under the action of Stokes waves on the sea, the interface between the air and water waves at the same regular pattern and so does the pressure and the vertical velocity. The ship hull moves in the same way as the wave. The amplitude of the ship hull's heave is less than the wave height because of the viscosity damping. This method could provide an important reference for the study of ships' movement, wave and hydrodynamics.
Essaid, H.I.
1990-01-01
The model allows for regional simulation of coastal groundwater conditions, including the effects of saltwater dynamics on the freshwater system. Vertically integrated freshwater and saltwater flow equations incorporating the interface boundary condition are solved within each aquifer. Leakage through confining layers is calculated by Darcy's law, accounting for density differences across the layer. The locations of the interface tip and toe, within grid blocks, are tracked by linearly extrapolating the position of the interface. The model has been verified using available analytical solutions and experimental results and applied to the Soquel-Aptos basin, Santa Cruz County, California. -from Author
Xu, Tianfu; Pruess, Karsten
1998-09-01
Coupled modeling of subsurface multiphase fluid and heat flow, solute transport and chemical reactions can be used for the assessment of acid mine drainage remediation, mineral deposition, waste disposal sites, hydrothermal convection, contaminant transport, and groundwater quality. Here they present a numerical simulation model, TOUGHREACT, which considers non-isothermal multi-component chemical transport in both liquid and gas phases. A wide range of subsurface thermo-physical-chemical processes is considered. The model can be applied to one-, two- or three-dimensional porous and fractured media with physical and chemical heterogeneity. The model can accommodate any number of chemical species present in liquid, gas and solid phases. A variety of equilibrium chemical reactions is considered, such as aqueous complexation, gas dissolution/exsolution, cation exchange, and surface complexation. Mineral dissolution/precipitation can proceed either subject to local equilibrium or kinetic conditions. The coupled model employs a sequential iteration approach with reasonable computing efficiency. The development of the governing equations and numerical approach is presented along with the discussion of the model implementation and capabilities. The model is verified for a wide range of subsurface physical and chemical processes. The model is well suited for flow and reactive transport in variably saturated porous and fractured media. In the second of this two-part paper, three applications covering a variety of problems are presented to illustrate the capabilities of the model.
NASA Astrophysics Data System (ADS)
Li, Linmin; Liu, Zhongqiu; Cao, Maoxue; Li, Baokuan
2015-07-01
In the ladle metallurgy process, the bubble movement and slag layer behavior is very important to the refining process and steel quality. For the bubble-liquid flow, bubble movement plays a significant role in the phase structure and causes the unsteady complex turbulent flow pattern. This is one of the most crucial shortcomings of the current two-fluid models. In the current work, a one-third scale water model is established to investigate the bubble movement and the slag open-eye formation. A new mathematical model using the large eddy simulation (LES) is developed for the bubble-liquid-slag-air four-phase flow in the ladle. The Eulerian volume of fluid (VOF) model is used for tracking the liquid-slag-air free surfaces and the Lagrangian discrete phase model (DPM) is used for describing the bubble movement. The turbulent liquid flow is induced by bubble-liquid interactions and is solved by LES. The procedure of bubble coming out of the liquid and getting into the air is modeled using a user-defined function. The results show that the present LES-DPM-VOF coupled model is good at predicting the unsteady bubble movement, slag eye formation, interface fluctuation, and slag entrainment.
A control method applied to mixed traffic flow for the coupled-map car-following model
NASA Astrophysics Data System (ADS)
Cheng, Rong-Jun; Han, Xiang-Lin; Lo, Siu-Ming; Ge, Hong-Xia
2014-03-01
In light of previous work [Phys. Rev. E 60 4000 (1999)], a modified coupled-map car-following model is proposed by considering the headways of two successive vehicles in front of a considered vehicle described by the optimal velocity function. The non-jam conditions are given on the basis of control theory. Through simulation, we find that our model can exhibit a better effect as p = 0.65, which is a parameter in the optimal velocity function. The control scheme, which was proposed by Zhao and Gao, is introduced into the modified model and the feedback gain range is determined. In addition, a modified control method is applied to a mixed traffic system that consists of two types of vehicle. The range of gains is also obtained by theoretical analysis. Comparisons between our method and that of Zhao and Gao are carried out, and the corresponding numerical simulation results demonstrate that the temporal behavior of traffic flow obtained using our method is better than that proposed by Zhao and Gao in mixed traffic systems.
Steinberg, A.M.; Boxx, I.; Stoehr, M.; Meier, W.; Carter, C.D.
2010-12-15
A detailed analysis of the flow-flame interactions associated with acoustically coupled heat-release rate fluctuations was performed for a 10 kW, CH{sub 4}/air, swirl stabilized flame in a gas turbine model combustor exhibiting self-excited thermo-acoustic oscillations at 308 Hz. High-speed stereoscopic particle image velocimetry, OH planar laser induced fluorescence, and OH* chemiluminescence measurements were performed at a sustained repetition rate of 5 kHz, which was sufficient to resolve the relevant combustor dynamics. Using spatio-temporal proper orthogonal decomposition, it was found that the flow-field contained several simultaneous periodic motions: the reactant flux into the combustion chamber periodically oscillated at the thermo-acoustic frequency (308 Hz), a helical precessing vortex core (PVC) circumscribed the burner nozzle at 515 Hz, and the PVC underwent axial contraction and extension at the thermo-acoustic frequency. The global heat release rate fluctuated at the thermo-acoustic frequency, while the heat release centroid circumscribed the combustor at the difference between the thermo-acoustic and PVC frequencies. Hence, the three-dimensional location of the heat release fluctuations depended on the interaction of the PVC with the flame surface. This motivated the compilation of doubly phase resolved statistics based on the phase of both the acoustic and PVC cycles, which showed highly repeatable periodic flow-flame configurations. These include flames stabilized between the inflow and inner recirculation zone, large-scale flame wrap-up by the PVC, radial deflection of the inflow by the PVC, and combustion in the outer recirculation zones. Large oscillations in the flame surface area were observed at the thermo-accoustic frequency that significantly affected the total heat-release oscillations. By filtering the instantaneous reaction layers at different scales, the importance of the various flow-flame interactions affecting the flame area was
Howe, M. S.; McGowan, R. S.
2009-01-01
An analysis is made of the nonlinear interactions between flow in the subglottal vocal tract and glottis, sound waves in the subglottal system and a mechanical model of the vocal folds. The mean flow through the system is produced by a nominally steady contraction of the lungs, and mechanical experiments frequently involve a ‘lung cavity’ coupled to an experimental subglottal tube of arbitrary or ill-defined effective length L, on the basis that the actual value of L has little or no influence on excitation of the vocal folds. A simple, self-exciting single mass mathematical model of the vocal folds is used to investigate the sound generated within the subglottal domain and the unsteady volume flux from the glottis for experiments where it is required to suppress feedback of sound from the supraglottal vocal tract. In experiments where the assumed absorption of sound within the sponge-like interior of the lungs is small, the influence of changes in L can be very significant: when the subglottal tube behaves as an open-ended resonator (when L is as large as half the acoustic wavelength) there is predicted to be a mild increase in volume flux magnitude and a small change in waveform. However, the strong appearance of second harmonics of the acoustic field is predicted at intermediate lengths, when L is roughly one quarter of the acoustic wavelength. In cases of large lung damping, however, only modest changes in the volume flux are predicted to occur with variations in L. PMID:20161450
Hossain, Md Shakhawath; Bergstrom, D J; Chen, X B
2015-11-01
The in vitro chondrocyte cell culture process in a perfusion bioreactor provides enhanced nutrient supply as well as the flow-induced shear stress that may have a positive influence on the cell growth. Mathematical and computational modelling of such a culture process, by solving the coupled flow, mass transfer and cell growth equations simultaneously, can provide important insight into the biomechanical environment of a bioreactor and the related cell growth process. To do this, a two-way coupling between the local flow field and cell growth is required. Notably, most of the computational and mathematical models to date have not taken into account the influence of the cell growth on the local flow field and nutrient concentration. The present research aimed at developing a mathematical model and performing a numerical simulation using the lattice Boltzmann method to predict the chondrocyte cell growth without a scaffold on a flat plate placed inside a perfusion bioreactor. The model considers the two-way coupling between the cell growth and local flow field, and the simulation has been performed for 174 culture days. To incorporate the cell growth into the model, a control-volume-based surface growth modelling approach has been adopted. The simulation results show the variation of local fluid velocity, shear stress and concentration distribution during the culture period due to the growth of the cell phase and also illustrate that the shear stress can increase the cell volume fraction to a certain extent. PMID:25804699
Adjoint-based estimation of plate coupling in a non-linear mantle flow model: theory and examples
NASA Astrophysics Data System (ADS)
Ratnaswamy, Vishagan; Stadler, Georg; Gurnis, Michael
2015-08-01
We develop and validate a systematic approach to infer plate boundary strength and rheological parameters in models of mantle flow from surface velocity observations. Based on a realistic rheological model that includes yielding and strain rate weakening from dislocation creep, we formulate the inverse problem in a Bayesian inference framework. To study the distribution of parameters that are consistent with the observations, we compute the maximum a posteriori (MAP) point, Gaussian approximations of the parameter distribution around that MAP point, and employ Markov Chain Monte Carlo (MCMC) sampling methods. The computation of the MAP point and the Gaussian approximation require first and second derivatives of an objective function subject to non-linear Stokes equations; these derivatives are computed efficiently using adjoint Stokes equations. We set up 2-D numerical experiments with many of the elements expected in a global geophysical inversion. This setup incorporates three subduction zones with slab and weak zone (interplate fault) geometry consistent with average seismic characteristics. With these experiments, we demonstrate that when the temperature field is known, we can recover the strength of plate boundaries, the yield stress and strain rate exponent in the upper mantle. When the number of uncertain parameters increases, there are trade-offs between the inferred parameters. These trade-offs depend on how well the observational data represents the surface velocities, and on the weakness of plate boundaries. As the plate boundary coupling drops below a threshold, the uncertainty of the inferred parameters increases due to insensitivity of plate motion to plate coupling. Comparing the trade-offs between inferred rheological parameters found from the Gaussian approximation of the parameter distribution and from MCMC sampling, we conclude that the Gaussian approximation-which is significantly cheaper to compute-is often a good approximation, in particular
Vibration dissociation coupling in nonequilibrium flows
NASA Technical Reports Server (NTRS)
Candler, Graham V.
1992-01-01
The final report on research between North Carolina State University and the NASA Ames Research Center is presented. The research was aimed at using the Schwartz, Slawsky, Herzfeld (SSH) theory to simulate the vibrational relaxation of nitrogen molecules undergoing dissociation or recombination over a wide range of conditions. The results of these simulations were then treated as exact, and they were used to develop a model for the coupled vibration-dissociation process. This new model is simple enough to be used in computational fluid dynamics codes, but still captures the physics of the complex process. The model is used to simulate the flow over typical geometries to test it and to determine how much impact it has on the flow field. The key elements of this research are summarized.
NASA Astrophysics Data System (ADS)
Gemitzi, Alexandra; Tolikas, Demetrios
A simulation program, which works seamlessly with GIS and simulates flows in coastal aquifers, is presented in the present paper. The model is based on the Galerkin finite element discretization scheme and it simulates both steady and transient freshwater and saltwater flow, assuming that the two fluids are separated by a sharp interface. The model has been verified in simple cases where analytical solutions exist. The simulation program works as a tool of the GIS program, which is the main database that stores and manages all the necessary data. The combined use of the simulation and the GIS program forms an integrated management tool offering a simpler way of simulating and studying saline intrusion in coastal aquifers. Application of the model to the Yermasogia aquifer illustrates the coupled use of modeling and GIS techniques for the examination of regional coastal aquifer systems. Pour étudier un système aquifère côtier, nous avons développé un modèle aux éléments finis en quasi 3-D qui simule les écoulements d'eau douce et d'eau salée en régime aussi bien permanent que transitoire. Les équations qui les régissent sont discrétisées par un schéma de discrétisation de Garlekin aux éléments finis. Le modèle a été vérifié dans des cas simples où il existe des solutions analytiques. Toutes les données nécessaires sont introduites et gérées grâce à un logiciel de gestion de SIG. Le programme de simulation est utilisé comme un outil du logiciel de SIG, constituant ainsi un outil de gestion intégrée dont le but est de simuler et d'étudier l'intrusion saline dans les aquifères côtiers. L'application du modèle à l'aquifère de Yermasogia illustre l'utilisation couplée de la modélisation et des techniques de SIG pour l'étude des systèmes aquifères côtiers régionaux. Se ha desarrollado un modelo casi tridimensional de elementos finitos para simular el flujo de agua dulce y salada, tanto en régimen estacionario como en
NASA Astrophysics Data System (ADS)
Klammler, G.; Rock, G.; Kupfersberger, H.; Fank, J.
2012-04-01
interval of every hydrotope vertical profile until the lowest observed groundwater table is reached. The fluctuation range of the phreatic surface is also discretized in 10 cm intervals and used as outflow boundary condition. By this procedure, the influence of the groundwater table on the water and nitrate mass leaving the unsaturated can be considered taken into account varying soil horizons. To cover saturated flow in the WLF aquifer a 2-dimensional transient horizontal flow and solute transport model is set up. A sequential coupling between the two models is implemented, i.e. a unidirectional transfer of recharge and nitrate mass outflow from the hydrotopes to the saturated compartment. For this purpose, a one-time assignment between the spatial discretization of the hydrotopes and the finite element mesh has to be set up. The resulting groundwater table computed for a given time step with the input from SIMWASER/STOTRASIM is then used to extract the corresponding water and nitrate mass values from the look-up table to be used for the consecutive time step. This process is being repeated until the end of the simulation period. Within this approach there is no direct feedback between the unsaturated and the saturated aquifer compartment, i.e. there is no simultaneous (within the same time step) update of the pressure head - unsaturated head relationship at the soil and the phreatic surface (like is shown e.g. in Walsum and Groedendijk, 2008). For the dominating coarse sand conditions of the WLF aquifer we believe that this simplification is not of further relevance. For higher soil moisture contents (i.e. almost full saturation near the groundwater table) the curve returns to specific retention within a short vertical distance. Thus, there might only be mutual impact between soil and phreatic surface conditions for shallow groundwater tables. However, it should be mentioned here that all other processes in the two compartments (including capillary rise due to clay rich
Siriwardane, Hema J; Gondle, Raj K; Bromhal, Grant S
2013-08-01
Understanding the transport of carbon dioxide (CO{sub 2}) during long-term CO{sub 2} injection into a typical geologic reservoir, such as a saline aquifer, could be complicated because of changes in geochemical, hydrogeological, and hydromechanical behavior. While the caprock layer overlying the target aquifer is intended to provide a tight, impermeable seal in securing injected CO{sub 2}, the presence of geologic uncertainties, such as a caprock fracture or fault, may provide a channel for CO{sub 2} leakage. There could also be a possibility of the activation of a new or existing dormant fault or fracture, which could act as a leakage pathway. Such a leakage event during CO{sub 2} injection may lead to a different pressure and ground response over a period of time. In the present study, multiphase fluid flow simulations in porous media coupled with geomechanics were used to investigate the overburden geologic response and plume behavior during CO{sub 2} injection in the presence of a hypothetical permeable fractured zone in a caprock, existing or activated. Both single-phase and multiphase fluid flow simulations were performed. The CO{sub 2} migration through an existing fractured zone leads to changes in the fluid pressure in the overburden geologic layers and could have a significant impact on ground deformation behavior. Results of the study show that pressure signatures and displacement patterns are significantly different in the presence of a fractured zone in the caprock layer. The variation in pressure and displacement signatures because of the presence of a fractured zone in the caprock at different locations may be useful in identifying the presence of a fault/fractured zone in the caprock. The pressure signatures can also serve as a mechanism to identify the activation of leakage pathways through the caprock during CO{sub 2} injection. Pressure response and ground deformation behavior from sequestration modeling could be useful in the development of
NASA Astrophysics Data System (ADS)
He, Lijuan; Wang, Jiyang
Heat flow measurements show a moderate thermal background (∼61 mW/m 2) in the Bohai Bay Basin, which although experienced multi-phase rifting in the Cenozoic era. In contrast, its surrounding mountain areas are characterized by low heat flow. Constraint by heat flow measurements, the thermal evolution of the Bohai Bay Basin during the Cenozoic era was performed by a numerical basin-mountain model. The model incorporates differential lithosphere stretching and shortening by finite-element method in the Lagrangian frame. The predicted heat flow in the center of the three depressions of the Bohai Bay Basin is calculated to have varied between 51 and 63 mW/m 2 through the Cenozoic evolution, indicating a rather smooth variation of basin thermal state, and a cooling trend from the Oligocene to present-day. Model results also suggest that the Taihang Mountains probably uplifted in the Quaternary, which resulted in low heat flow in the mountain area. Both heat flow constraints and modeling imply that a new phase of rifting in the Pliocene existed in the Huanghua and Bozhong Depressions, which was suggested by tectonic subsidence analysis.
NASA Astrophysics Data System (ADS)
He, L.; Wang, J.
2003-04-01
Bohai Bay Basin is located in the east of North China Craton. Heat flow measurements show a moderate thermal background (~ 61 mW/m2) in the Bohai Bay Basin, which although experienced multi-phase rifting in the Cenozoic era. In contrast, its surrounding mountain areas are characterized by low heat flow. Constraint by heat flow measurements, the thermal evolution of the Bohai Bay Basin during the Cenozoic era was performed by a numerical basin-mountain model. The model incorporates differential lithosphere stretching and shortening by finite-element method in the Lagrangian frame. The predicted heat flow in the center of the three depressions of the Bohai Bay Basin is calculated to have varied between 51 mW/m2 and 63 mW/m2 through the Cenozoic evolution, indicating a rather smooth variation of basin thermal state, and a cooling trend from the Oligocene to present-day. Model results also suggest that the Taihang Mountains probably uplifted in the Quaternary, which resulted in low heat flow in the mountain area. Both heat flow constraints and modeling imply that a new phase of rifting in the Pliocene existed in the Huanghua and Bozhong Depressions, which was suggested by tectonic subsidence analysis.
NASA Astrophysics Data System (ADS)
Francioso, L.; De Pascali, C.; Pescini, E.; De Giorgi, M. G.; Siciliano, P.
2016-06-01
Preventing the flow separation could enhance the performance of propulsion systems and future civil aircraft. To this end, a fast detection of boundary layer separation is mandatory for a sustainable and successful application of active flow control devices, such as plasma actuators. The present work reports on the design, fabrication and functional tests of low-cost capacitive pressure sensors coupled with dielectric barrier discharge (DBD) plasma actuators to detect and then control flow separation. Finite element method (FEM) simulations were used to obtain information on the deflection and the stress distribution in different-shaped floating membranes. The sensor sensitivity as a function of the pressure load was also calculated by experimental tests. The results of the calibration of different capacitive pressure sensors are reported in this work, together with functional tests in a wind tunnel equipped with a curved wall plate on which a DBD plasma actuator was mounted to control the flow separation. The flow behavior was experimentally investigated by particle image velocimetry (PIV) measurements. Statistical and spectral analysis, applied to the output signals of the pressure sensor placed downstream of the profile leading edge, demonstrated that the sensor is able to discriminate different ionic wind velocity and turbulence conditions. The sensor sensitivity in the 0–100 Pa range was experimentally measured and it ranged between 0.0030 and 0.0046 pF Pa‑1 for the best devices.
Technology Transfer Automated Retrieval System (TEKTRAN)
Watershed modeling is a key component of watershed management that involves the simulation of hydrological and fluvial processes for predicting flow and sediment transport within a watershed. For practical purposes, most numerical models have been developed to simulate either runoff and soil erosion...
NASA Astrophysics Data System (ADS)
Carroll, K. C.; Nguyen, B. N.; Fang, Y.; Richmond, M. C.; Murray, C. J.
2011-12-01
Geomechanical alteration of porous media is generally ignored for most shallow subsurface applications, whereas CO2 injection, migration, and trapping in deep saline aquifers will be controlled by coupled multifluid flow, energy transfer, geomechanical, and geochemical processes. The accurate assessment of the risks associated with potential leakage of injected CO2 and the design of effective injection systems requires that we represent these coupled processes within numerical simulators. The objective of this study was to examine the coupling of hydraulic and geomechanical processes for simulation of CO2 injection into the subsurface for carbon sequestration. The impact of nonisothermal multifluid flow and porous media deformation mechanics on CO2 migration and storage was evaluated. We present a sequentially coupled approach for multifluid and geomechanical simulation using STOMP and ABAQUS that has been developed and validated through comparison to the solutions for benchmark problems that were solved with a coupled TOUGH-FLAC simulator. The poroelastic model was implemented with user-subroutines in ABAQUS. We also compare the STOMP-ABAQUS simulator to a new version of STOMP that includes the fully coupled poroelastic simulation within the multifluid flow and transport simulator. The poroelastic model computes stiffness, stresses, and strains using aqueous and gas pressures as well as saturations from STOMP output, and provides STOMP with the updated permeability, porosity, and capillary pressure over time during the simulation. The hydraulic only (uncoupled from mechanics) simulation and the hydrogeomechanical (coupled) simulation results using STOMP-ABAQUS were comparable to the previous results of a TOUGH-FLAC simulator. Results from the STOMP-ABAQUS coupled simulator were essentially identical to the fully coupled STOMP hydrogeomechanical simulator when the sequential coupling occurred at small time steps, and deviations between results increased with
Nuclear-Coupled Flow Instabilities and Their Effects on Dryout
M. Ishii; X. Sunn; S. Kuran
2004-09-27
Nuclear-coupled flow/power oscillations in boiling water reactors (BWRs) are investigated experimentally and analytically. A detailed literature survey is performed to identify and classify instabilities in two-phase flow systems. The classification and the identification of the leading physical mechanisms of the two-phase flow instabilities are important to propose appropriate analytical models and scaling criteria for simulation. For the purpose of scaling and the analysis of the nonlinear aspects of the coupled flow/power oscillations, an extensive analytical modeling strategy is developed and used to derive both frequency and time domain analysis tools.
NASA Astrophysics Data System (ADS)
Khaleel, M. A.; Lin, Z.; Singh, P.; Surdoval, W.; Collin, D.
A 3D simulation tool for modeling solid oxide fuel cells is described. The tool combines the versatility and efficiency of a commercial finite element analysis code, MARC ®, with an in-house developed robust and flexible electrochemical (EC) module. Based upon characteristic parameters obtained experimentally and assigned by the user, the EC module calculates the current density distribution, heat generation, and fuel and oxidant species concentration, taking the temperature profile provided by MARC ® and operating conditions such as the fuel and oxidant flow rate and the total stack output voltage or current as the input. MARC ® performs flow and thermal analyses based on the initial and boundary thermal and flow conditions and the heat generation calculated by the EC module. The main coupling between MARC ® and EC is for MARC ® to supply the temperature field to EC and for EC to give the heat generation profile to MARC ®. The loosely coupled, iterative scheme is advantageous in terms of memory requirement, numerical stability and computational efficiency. The coupling is iterated to self-consistency for a steady-state solution. Sample results for steady states as well as the startup process for stacks with different flow designs are presented to illustrate the modeling capability and numerical performance characteristic of the simulation tool.
George, David L.; Iverson, Richard M.
2011-01-01
Pore-fluid pressure plays a crucial role in debris flows because it counteracts normal stresses at grain contacts and thereby reduces intergranular friction. Pore-pressure feedback accompanying debris deformation is particularly important during the onset of debrisflow motion, when it can dramatically influence the balance of forces governing downslope acceleration. We consider further effects of this feedback by formulating a new, depth-averaged mathematical model that simulates coupled evolution of granular dilatancy, solid and fluid volume fractions, pore-fluid pressure, and flow depth and velocity during all stages of debris-flow motion. To illustrate implications of the model, we use a finite-volume method to compute one-dimensional motion of a debris flow descending a rigid, uniformly inclined slope, and we compare model predictions with data obtained in large-scale experiments at the USGS debris-flow flume. Predictions for the first 1 s of motion show that increasing pore pressures (due to debris contraction) cause liquefaction that enhances flow acceleration. As acceleration continues, however, debris dilation causes dissipation of pore pressures, and this dissipation helps stabilize debris-flow motion. Our numerical predictions of this process match experimental data reasonably well, but predictions might be improved by accounting for the effects of grain-size segregation.
NASA Astrophysics Data System (ADS)
Rodrigues, C. Veiga; Palma, J. M. L. M.; Rodrigues, Á. H.
2016-05-01
The atmospheric flow over a mountainous region has been simulated using a model-chain approach, whereby the flow in a larger region was simulated using a mesoscale model with three nesting levels, down to a 3-km horizontal resolution, within which a fourth nesting level was set with a microscale flow solver and a domain with varying horizontal resolution, around 300 m at the site of interest. Two periods in the summer (July) and autumn (November-December) 2005, each with a duration of two weeks, were selected to test the present approach. Two sites were chosen, comprising a total of seven meteorological masts with wind vanes and anemometers at two heights. The microscale solver improved the wind-speed prediction of the mesoscale model in 10 of the 14 anemometers and replicated the high wind speeds, which were under-predicted in the mesoscale model. The wind conditions in summer varied with the daily cycle, related to regional-scale sea breezes and their interaction with local circulations induced by the topography. Regarding the turbulence intensity, the predicted decay with wind-speed increase was in agreement with the measurements. This study shows the need of both models: the microscale model captures the details of the boundary-layer physics, which would not be possible without the boundary conditions provided by the mesoscale model.
Technology Transfer Automated Retrieval System (TEKTRAN)
Many or most subsurface pollution problems at the field scale involve such simultaneous processes as water flow, multicomponent solute transport, heat transport and biogeochemical processes and reactions. Process-based models that integrate these various processes can be valuable tools for investiga...
Miller, P.J.
1996-07-01
A new reactive flow model for highly non-ideal explosives and propellants is presented. These compositions, which contain large amounts of metal, upon explosion have reaction kinetics that are characteristic of both fast detonation and slow metal combustion chemistry. A reaction model for these systems was incorporated into the two-dimensional, finite element, Lagrangian hydrodynamic code, DYNA2D. A description of how to determine the model parameters is given. The use of the model and variations are applied to AP, Al, and nitramine underwater explosive and propellant systems.
Maxwell, R; Kollet, S; Chow, F; Granvold, P; Duan, Q
2007-02-23
This final report is organized in four sections. Section 1 is the project summary (below), Section 2 is a submitted manuscript that describes the offline, or spinup simulations in detail, Section 3 is also a submitted manuscript that describes the online, or fully-coupled simulations in detail and Section 3, which is report that describes work done via a subcontract with UC Berkeley. The goal of this project was to develop and apply a coupled regional climate, land-surface, groundwater flow model as a means to further understand important mass and energy couplings between regional climate, the land surface, and groundwater. The project involved coupling three distinct submodels that are traditionally used independently with abstracted and potentially oversimplified (inter-model) boundary conditions. This coupled model lead to (1) an improved understanding of the sensitivity and importance of coupled physical processes from the subsurface to the atmosphere; (2) a new tool for predicting hydrologic conditions (rainfall, temperature, snowfall, snowmelt, runoff, infiltration and groundwater flow) at the watershed scale over a range of timeframes; (3) a simulation of hydrologic response of a characteristic watershed that will provide insight into the certainty of hydrologic forecasting, dominance and sensitivity of groundwater dynamics on land-surface fluxes; and (4) a more realistic model representation of weather predictions, precipitation and temperature, at the regional scale. Regional climate models are typically used for the simulation of weather, precipitation and temperature behavior over 10-1000 km domains for weather or climate prediction purposes, and are typically driven by boundary conditions derived from global climate models (GCMs), observations or both. The land or ocean surface typically represents a bottom boundary condition of these models, where important mass (water) and energy fluxes are approximated. The viability and influence of these
De Focatiis, Davide S. A.; Buckley, C. Paul; Embery, John
2008-07-07
This paper investigates the behaviour of a well-characterised monodisperse grade of entangled atactic polystyrene across a very wide temperature and strain rate range through linear and non-linear melt rheology and solid-state deformation. In an effort to construct a constitutive model for large deformations able to describe rheological response right across this wide timescale, two well-established rheological models are combined: the well known RoliePoly (RP) conformational melt model and the Oxford glass-rubber constitutive model for glassy polymers. Comparisons between experimental data and simulations from a numerical implementation of the model illustrate that the model can cope well with the range of deformations in which orientation is limited to length-scales longer than an entanglement length. One approach in which the model can be expanded to incorporate the effects of orientation on shorter length scales using anisotropic viscoplastic flow is briefly discussed.
Cheng, Cheng; Zhang, Xiaobing
2013-05-01
In conventional models for two-phase reactive flow of interior ballistic, the dynamic collision phenomenon of particles is neglected or empirically simplified. However, the particle collision between particles may play an important role in dilute two-phase flow because the distribution of particles is extremely nonuniform. The collision force may be one of the key factors to influence the particle movement. This paper presents the CFD-DEM approach for simulation of interior ballistic two-phase flow considering the dynamic collision process. The gas phase is treated as a Eulerian continuum and described by a computational fluid dynamic method (CFD). The solid phase is modeled by discrete element method (DEM) using a soft sphere approach for the particle collision dynamic. The model takes into account grain combustion, particle-particle collisions, particle-wall collisions, interphase drag and heat transfer between gas and solid phases. The continuous gas phase equations are discretized in finite volume form and solved by the AUSM+-up scheme with the higher order accurate reconstruction method. Translational and rotational motions of discrete particles are solved by explicit time integrations. The direct mapping contact detection algorithm is used. The multigrid method is applied in the void fraction calculation, the contact detection procedure, and CFD solving procedure. Several verification tests demonstrate the accuracy and reliability of this approach. The simulation of an experimental igniter device in open air shows good agreement between the model and experimental measurements. This paper has implications for improving the ability to capture the complex physics phenomena of two-phase flow during the interior ballistic cycle and to predict dynamic collision phenomena at the individual particle scale. PMID:24891728
Hibi, Yoshihiko; Tomigashi, Akira
2015-09-01
Numerical simulations that couple flow in a surface fluid with that in a porous medium are useful for examining problems of pollution that involve interactions among atmosphere, water, and groundwater, including saltwater intrusion along coasts. Coupled numerical simulations of such problems must consider both vertical flow between the surface fluid and the porous medium and complicated boundary conditions at their interface. In this study, a numerical simulation method coupling Navier-Stokes equations for surface fluid flow and Darcy equations for flow in a porous medium was developed. Then, the basic ability of the coupled model to reproduce (1) the drawdown of a surface fluid observed in square-pillar experiments, using pillars filled with only fluid or with fluid and a porous medium and (2) the migration of saltwater (salt concentration 0.5%) in the porous medium using the pillar filled with fluid and a porous medium was evaluated. Simulations that assumed slippery walls reproduced well the results with drawdowns of 10-30 cm when the pillars were filled with packed sand, gas, and water. Moreover, in the simulation of saltwater infiltration by the method developed in this study, velocity was precisely reproduced because the experimental salt concentration in the porous medium after saltwater infiltration was similar to that obtained in the simulation. Furthermore, conditions across the boundary between the porous medium and the surface fluid were satisfied in these numerical simulations of square-pillar experiments in which vertical flow predominated. Similarly, the velocity obtained by the simulation for a system coupling flow in surface fluid with that in a porous medium when horizontal flow predominated satisfied the conditions across the boundary. Finally, it was confirmed that the present simulation method was able to simulate a practical-scale surface fluid and porous medium system. All of these numerical simulations, however, required a great deal of
NASA Technical Reports Server (NTRS)
1995-01-01
The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF
Flow noise source-resonator coupling
Pollack, M.L.
1997-11-01
This paper investigates the coupling mechanism between flow noise sources and acoustic resonators. Analytical solutions are developed for the classical cases of monopole and dipole types of flow noise sources. The effectiveness of the coupling between the acoustic resonator and the noise source is shown to be dependent on the type of noise source as well as its location on the acoustic pressure mode shape. For a monopole source, the maximum coupling occurs when the noise source is most intense near an acoustic pressure antinode (i.e., location of maximum acoustic pressure). A numerical study with the impedance method demonstrates this effect. A dipole source couples most effectively when located near an acoustic pressure node.
NASA Astrophysics Data System (ADS)
Mendicino, Giuseppe; Pedace, Jessica; Senatore, Alfonso
2015-04-01
Cellular Automata are often used for modeling the evolution in time of environmental systems mainly because they are directly compatible with parallel programming. Nevertheless, defining the optimal time step criterion for integrating forward in time numerical processes can further enhance model computational efficiency. To this aim, a numerical stability analysis of an original overland flow model, within the framework of a fully coupled eco-hydrological system based on the Macroscopic Cellular Automata paradigm, is performed. According to the other modules of the system describing soil water flow, soil-surface-atmosphere fluxes and vegetation dynamics, overland flow model equations were derived through a direct discrete formulation (i.e. no differential equations were discretized), adopting the diffusion wave model as an approximation of the full De Saint Venant equations and including the capability of accounting for specific processes, such as the increasing roughness effects due to vegetation growth or surface-soil water exchanges. Suitable formulations of robust tools usually applied in the stability analyses, such as Courant-Friedrichs-Lewy and von Neumann conditions, were initially derived for the CA-based overland flow model. Afterwards, the theoretical stability conditions were compared to experimental time step constraints through several numerical simulations of a 5-h rain event. Specifically, adopting a constant (i.e. not adaptive) time step for simulations, and discretizing head losses in a way that increases model stability, experimental upper limits preventing numerical instability were found for 13 test cases with different slopes, precipitation intensities, vegetation densities and depths of surface depressions. Even though von Neumann condition and experimental values were well positively correlated, the latter were almost always sensibly lower, excluding cases when free surface gradients tended to zero. Therefore, based on the original method
NASA Astrophysics Data System (ADS)
Mikolajewicz, Uwe; Crowley, Thomas J.
1997-06-01
Prior ocean modeling work suggested that an open central American isthmus would cause a collapse of the North Atlantic thermohaline circulation because of free exchange of low salinity water between the Atlantic and the Pacific. Geological data provide some support for this response, but the data also indicate that some North Atlantic Deep Water formation occurred before final closure of the isthmus. We previously postulated that this "early switch on" could reflect a more limited exchange of Atlantic waters with the Pacific. In this study we discuss model sensitivity experiments testing that hypothesis and interpret the response in terms of shifts between multiple steady states of the model. Two simulations are conducted with a version of the Hamburg large-scale geostrophic ocean model that is coupled to an atmospheric energy balance model. Constrictions of throughflow through the central American isthmus is mimicked by locally changing the frictional drag coefficient in the ocean model. Results indicate that modest levels of throughflow can maintain some level of thermohaline circulation. These results support the conjecture in our earlier study. However, the overturning cell is about 300 m shallower than in the control run, with deep water production nearly eliminated in the Labrador Sea. These latter responses should be testable with marine data.
Coupled Groundwater and Heat Flow in the Tahoe Basin Region
NASA Astrophysics Data System (ADS)
Trask, J. C.
2002-12-01
We propose that recent developments in available temperature probe technology and improvements in appropriate modeling software, together with the advent of desktop high-speed computing, have enabled the thermal analysis approach to be an inexpensive, robust, and practical way to constrain groundwater flow fields in a wide variety of environments. We present field measurements and numerical models of coupled heat and groundwater flow in the Tahoe Basin region. In montane regions such as the Tahoe Basin, steep topography provides a driving force for deep groundwater flow. Deep groundwater flow re-routes subsurface heat flow, impacting temperature gradients to depth, including the shallow subsurface (<100m depth). In the Tahoe Basin region, the magnitude of deep groundwater flow on the areal or regional scale has been largely unknown. We present examples of borehole temperature profiles that constrain possible areal groundwater flow patterns, including the magnitude of flow beneath the bottom of boreholes probed.
NASA Astrophysics Data System (ADS)
Wen-Chao, Liu; Yue-Wu, Liu; Cong-Cong, Niu; Guo-Feng, Han; Yi-Zhao, Wan
2016-02-01
The threshold pressure gradient and formation stress-sensitive effect as the two prominent physical phenomena in the development of a low-permeable reservoir are both considered here for building a new coupled moving boundary model of radial flow in porous medium. Moreover, the wellbore storage and skin effect are both incorporated into the inner boundary conditions in the model. It is known that the new coupled moving boundary model has strong nonlinearity. A coordinate transformation based fully implicit finite difference method is adopted to obtain its numerical solutions. The involved coordinate transformation can equivalently transform the dynamic flow region for the moving boundary model into a fixed region as a unit circle, which is very convenient for the model computation by the finite difference method on fixed spatial grids. By comparing the numerical solution obtained from other different numerical method in the existing literature, its validity can be verified. Eventually, the effects of permeability modulus, threshold pressure gradient, wellbore storage coefficient, and skin factor on the transient wellbore pressure, the derivative, and the formation pressure distribution are analyzed respectively. Project supported by the National Natural Science Foundation of China (Grant No. 51404232), the China Postdoctoral Science Foundation (Grant No. 2014M561074), and the National Science and Technology Major Project, China (Grant No. 2011ZX05038003).
Swain, E.D.; Howie, Barbara; Dixon, Joann
1996-01-01
A coupled surface-water model (BRANCH) and ground-water model (MODFLOW) model were tested to simulate the interacting wetlands/surface-water/ ground-water system of southern Dade County. Several options created for the MODFLOW ground- ground-water model were used in representing this field situation. The primary option is the MODBRANCH interfacing software, which allows leakage to be accounted for between the MODFLOW ground-water model and the BRANCH dynamic model for simulation of flow in an interconnected network of open channels. A modification to an existing software routine, which is referred to as BCF2, allows cells in MODFLOW to rewet when dry--a requirement in representing the seasonal wetlands in Dade County. A companion to BCF2 is the modified evapotranspiration routine EVT2. The EVT2 routine changes the cells where evapotranspiration occurs, depending on which cells are wet. The Streamlink package represents direct connections between the canals and wetlands at locations where canals open directly into overland flow. Within the BRANCH model, the capability to represent the numerous hydraulic structures, gated spillways, gated culverts, and pumps was added. The application of these modifications to model surface-water/ground-water interactions in southern Dade County demonstrated the usefulness of the coupled MODFLOW/BRANCH model. Ground-water and surface-water flows are both simulated with dynamic models. Flow exchange between models, intermittent wetting and drying, evapotranspiration, and hydraulic structure operations are all represented appropriately. Comparison was made with a simulation using the RIV1 package instead of MODBRANCH to represent the canals. RIV1 represents the canals by user-defined stages, and computes leakage to the aquifer. Greater accuracy in reproducing measured ground- water heads was achieved with MODBRANCH, which also computes dynamic flow conditions in the canals, unlike RIV1. The surface-water integrated flow and transport
NASA Astrophysics Data System (ADS)
Gregory Forest, M.; Wang, Qi; Bechtel, Stephen E.
Slender asymptotic fiber models are derived from Doi-type 3-D equations for free surface flows of liquid-crystalline polymers. Leading order equations and self-consistent corrections are presented for a variety of physical regimes. We then explore the coupling of orientation effects to slender elongational flow behavior, with particular focus on the interplay between the Rayleigh capillary instability and both stabilizing and destabilizing orientation behavior. In the simple context of constant solutions, we identify physical regimes and precise conditions under which the Rayleigh instability may be completely arrested, as well as other regimes where orientation reduces but does not cancel capillary instability. In addition, we identify sources of additional orientation-dominated instabilities that are evident in both the uniaxial and biaxial nematic liquid crystal order parameters. These models and stability analyses lay the foundation for applications to fiber spinning processes.
Hydraulic jumps in inhomogeneous strongly coupled toroidal dust flows
NASA Astrophysics Data System (ADS)
Piel, Alexander; Wilms, Jochen
2016-07-01
The inhomogeneous flow of strongly coupled dust particles in a toroidal particle trap with harmonic radial confinement is analyzed in the incompressible fluid limit. It is shown that the flow can spontaneously generate shock-like events, which are similar to the hydraulic jump in open channel flows. A definition of the Froude number for this model is given and the critical speed is recovered as the group velocity of surface waves. This hydraulic model is compared with molecular-dynamics simulations, which show that a sudden bifurcation of the flow lines and a localized temperature peak appear just at the point where the critical condition for the hydraulic jump is located.
NASA Astrophysics Data System (ADS)
Munters, Wim; Meneveau, Charles; Meyers, Johan
2014-11-01
In order to incorporate multiple scales of meteorological phenomena in atmospheric simulations, subsequent nesting of meso-scale models is often used. However, the spatial and temporal resolution in such models is too coarse to resolve the three-dimensional turbulent eddies that are characteristic for atmospheric boundary layer flows. This motivates the development of tools to couple meso-scale models to Large-Eddy Simulations (LES), in which turbulent fluctuations are explicitly resolved. A major challenge in this area is the spin-up region near the inlet of the LES in which the flow has to evolve from a RANS-like inflow, originating from the meso-scale model, to a fully turbulent velocity field. We propose a generalized concurrent precursor inflow method capable of imposing boundary conditions for time-varying inflow directions. The method is based on a periodic fully-developed precursor boundary-layer simulation that is dynamically rotated with the wind direction that drives the main LES. In this way realistic turbulent inflow conditions are applied while still retaining flexibility to dynamically adapt to meso-scale variations in wind directions. Applications to wind simulations with varying inflow directions, and comparisons to conventional coupling methods are shown. Work supported by ERC (ActiveWindFarms, Grant No: 306471). CM is supported by NSF (Grant No. 1243482).
NASA Technical Reports Server (NTRS)
Allison, Dennis O.; Cavallo, Peter A.
2003-01-01
An equivalent-plate structural deformation technique was coupled with a steady-state unstructured-grid three-dimensional Euler flow solver and a two-dimensional strip interactive boundary-layer technique. The objective of the research was to assess the extent to which a simple accounting for static model deformations could improve correlations with measured wing pressure distributions and lift coefficients at transonic speeds. Results were computed and compared to test data for a wing-fuselage model of a generic low-wing transonic transport at a transonic cruise condition over a range of Reynolds numbers and dynamic pressures. The deformations significantly improved correlations with measured wing pressure distributions and lift coefficients. This method provided a means of quantifying the role of dynamic pressure in wind-tunnel studies of Reynolds number effects for transonic transport models.
NASA Astrophysics Data System (ADS)
Naz, B. S.; Frans, C. D.; Clarke, G. K.; Nolin, A. W.; Lettenmaier, D. P.; Istanbulluoglu, E.; Burns, P. J.
2011-12-01
Several recent studies have suggested that observed decreases in summer flows in Canada's South Saskatchewan River are partly due to retreat of glaciers in the river's headwaters. Despite the risk posed by declining glaciers to water supply in the high mountain river systems, our ability to accurately predict runoff contribution from partially glacierized basins is limited. Modeling the effect of glacier changes on streamflow response in such basins is complicated due to limited availability of high resolution gridded meteorological data, lack of long term measurements of glaciological parameters and most importantly glacier dynamics are not linked to hydrological processes in many existing physically-based distributed hydrologic models. We investigate the effect of glacier recession on streamflow variations for the Upper Bow River basin, a tributary of the South Saskatchewan, near Lake Louise, Alberta, using the Distributed Hydrology Soil Vegetation Model (DHSVM) coupled with the spatially distributed glacier dynamics model. The coupled model is forced with the North American Regional Reanalysis (NARR) climate data for the period of 1979 - 2010 at a 3-hourly time step. The NARR data are adjusted for spatial variability in precipitation and temperature using the Parameter-elevation Regressions on Independent Slopes Model (PRISM) monthly data at 2.5 arcmin resolution made available through the Climate Western North America (ClimateWNA) database (Wang et al. 2006). Using known subglacial bed topography information, a multidecade spin-up run of the stand alone glacier model is first conducted until the beginning of the simulation period for the coupled model to accurately predict ice thickness confirmed through comparison of modeled ice margins with observed glacier extent. The integrated model initialized with already estimated glacier thickness and ice extent is then run to predict glacier evolution, including spatial extent in combination with other hydrologic
Plasmoelectronics: coupling plasmonic excitation with electron flow.
Warren, Scott C; Walker, David A; Grzybowski, Bartosz A
2012-06-19
Explorations of the coupling of light and charge via localized surface plasmons have led to the discovery that plasmonic excitation can influence macroscopic flows of charge and, conversely, that charging events can change the plasmonic excitation. We discuss recent theory and experiments in the emerging field of plasmoelectronics, with particular emphasis on the application of these materials to challenges in nanotechnology, energy use, and sensing. PMID:22385329
Sill, W.R.; Killpack, T.J.
1982-03-01
The program is applicable to the calculation of self-potential effects due to fluid flow (electrokinetic effects) and heat flow (thermoelectric effects). The geological structure is two dimensional but the sources can be either finite line sources or point sources. The accuracy of the calculated potentials depends on the model discretization and the distance from the source(s). For the default mesh, the accuracy is usually a few percent at a distance of about one unit from the source. Surface boundary conditions for the primary problem require careful consideration as the form of the flow near the air-earth interface can have a profound effect on the resultant electric potentials. For temperature problems the appropriate boundary condition is a constant temperature, which is taken as zero. With this boundary condition there is a normal flux of heat at the surface and there will be induced electrical sources here, if the surface medium has a nonzero coupling coefficient. In the models, zero temperature at the surface is produced by giving the air a very large thermal conductivity.
NASA Astrophysics Data System (ADS)
Guo, L.; Huang, H.; Gaston, D.; Redden, G. D.
2009-12-01
One approach for immobilizing subsurface metal contaminants involves stimulating the in situ production of mineral phases that sequester or isolate contaminants. One example is using calcium carbonate to immobilize strontium. The success of such approaches depends on understanding how various processes of flow, transport, reaction and resulting porosity-permeability change couple in subsurface systems. Reactive transport models are often used for such purpose. Current subsurface reactive transport simulators typically involve a de-coupled solution approach, such as operator-splitting, that solves the transport equations for components and batch chemistry sequentially, which has limited applicability for many biogeochemical processes with fast kinetics and strong medium property-reaction interactions. A massively parallel, fully coupled, fully implicit reactive transport simulator has been developed based on a parallel multi-physics object oriented software environment computing framework (MOOSE) developed at the Idaho National Laboratory. Within this simulator, the system of transport and reaction equations is solved simultaneously in a fully coupled manner using the Jacobian Free Newton-Krylov (JFNK) method with preconditioning. The simulator was applied to model reactive transport in a one-dimensional column where conditions that favor calcium carbonate precipitation are generated by urea hydrolysis that is catalyzed by urease enzyme. Simulation results are compared to both laboratory column experiments and those obtained using the reactive transport simulator STOMP in terms of: the spatial and temporal distributions of precipitates and reaction rates and other major species in the reaction system; the changes in porosity and permeability; and the computing efficiency based on wall clock simulation time.
Blanchard, Solenna; Saillet, Sandrine; Ivanov, Anton; Benquet, Pascal; Bénar, Christian-George; Pélégrini-Issac, Mélanie; Benali, Habib; Wendling, Fabrice
2016-01-01
Developing a clear understanding of the relationship between cerebral blood flow (CBF) response and neuronal activity is of significant importance because CBF increase is essential to the health of neurons, for instance through oxygen supply. This relationship can be investigated by analyzing multimodal (fMRI, PET, laser Doppler…) recordings. However, the important number of intermediate (non-observable) variables involved in the underlying neurovascular coupling makes the discovery of mechanisms all the more difficult from the sole multimodal data. We present a new computational model developed at the population scale (voxel) with physiologically relevant but simple equations to facilitate the interpretation of regional multimodal recordings. This model links neuronal activity to regional CBF dynamics through neuro-glio-vascular coupling. This coupling involves a population of glial cells called astrocytes via their role in neurotransmitter (glutamate and GABA) recycling and their impact on neighboring vessels. In epilepsy, neuronal networks generate epileptiform discharges, leading to variations in astrocytic and CBF dynamics. In this study, we took advantage of these large variations in neuronal activity magnitude to test the capacity of our model to reproduce experimental data. We compared simulations from our model with isolated epileptiform events, which were obtained in vivo by simultaneous local field potential and laser Doppler recordings in rats after local bicuculline injection. We showed a predominant neuronal contribution for low level discharges and a significant astrocytic contribution for higher level discharges. Besides, neuronal contribution to CBF was linear while astrocytic contribution was nonlinear. Results thus indicate that the relationship between neuronal activity and CBF magnitudes can be nonlinear for isolated events and that this nonlinearity is due to astrocytic activity, highlighting the importance of astrocytes in the
Blanchard, Solenna; Saillet, Sandrine; Ivanov, Anton; Benquet, Pascal; Bénar, Christian-George; Pélégrini-Issac, Mélanie; Benali, Habib; Wendling, Fabrice
2016-01-01
Developing a clear understanding of the relationship between cerebral blood flow (CBF) response and neuronal activity is of significant importance because CBF increase is essential to the health of neurons, for instance through oxygen supply. This relationship can be investigated by analyzing multimodal (fMRI, PET, laser Doppler…) recordings. However, the important number of intermediate (non-observable) variables involved in the underlying neurovascular coupling makes the discovery of mechanisms all the more difficult from the sole multimodal data. We present a new computational model developed at the population scale (voxel) with physiologically relevant but simple equations to facilitate the interpretation of regional multimodal recordings. This model links neuronal activity to regional CBF dynamics through neuro-glio-vascular coupling. This coupling involves a population of glial cells called astrocytes via their role in neurotransmitter (glutamate and GABA) recycling and their impact on neighboring vessels. In epilepsy, neuronal networks generate epileptiform discharges, leading to variations in astrocytic and CBF dynamics. In this study, we took advantage of these large variations in neuronal activity magnitude to test the capacity of our model to reproduce experimental data. We compared simulations from our model with isolated epileptiform events, which were obtained in vivo by simultaneous local field potential and laser Doppler recordings in rats after local bicuculline injection. We showed a predominant neuronal contribution for low level discharges and a significant astrocytic contribution for higher level discharges. Besides, neuronal contribution to CBF was linear while astrocytic contribution was nonlinear. Results thus indicate that the relationship between neuronal activity and CBF magnitudes can be nonlinear for isolated events and that this nonlinearity is due to astrocytic activity, highlighting the importance of astrocytes in the
Coupled Flow and Mechanics in Porous and Fractured Media*
NASA Astrophysics Data System (ADS)
Martinez, M. J.; Newell, P.; Bishop, J.
2012-12-01
Numerical models describing subsurface flow through deformable porous materials are important for understanding and enabling energy security and climate security. Some applications of current interest come from such diverse areas as geologic sequestration of anthropogenic CO2, hydro-fracturing for stimulation of hydrocarbon reservoirs, and modeling electrochemistry-induced swelling of fluid-filled porous electrodes. Induced stress fields in any of these applications can lead to structural failure and fracture. The ultimate goal of this research is to model evolving faults and fracture networks and flow within the networks while coupling to flow and mechanics within the intact porous structure. We report here on a new computational capability for coupling of multiphase porous flow with geomechanics including assessment of over-pressure-induced structural damage. The geomechanics is coupled to the flow via the variation in the fluid pore pressures, whereas the flow problem is coupled to mechanics by the concomitant material strains which alter the pore volume (porosity field) and hence the permeability field. For linear elastic solid mechanics a monolithic coupling strategy is utilized. For nonlinear elastic/plastic and fractured media, a segregated coupling is presented. To facilitate coupling with disparate flow and mechanics time scales, the coupling strategy allows for different time steps in the flow solve compared to the mechanics solve. If time steps are synchronized, the controller allows user-specified intra-time-step iterations. The iterative coupling is dynamically controlled based on a norm measuring the degree of variation in the deformed porosity. The model is applied for evaluation of the integrity of jointed caprock systems during CO2 sequestration operations. Creation or reactivation of joints can lead to enhanced pathways for leakage. Similarly, over-pressures can induce flow along faults. Fluid flow rates in fractures are strongly dependent on the
Translation-vibration-dissociation coupling in nonequilibrium hypersonic flows
NASA Technical Reports Server (NTRS)
Candler, Graham
1989-01-01
A new simple and computationally efficient model was developed, describing the evolution of vibrational states during relaxation and dissociation. The model is based on dividing the nitrogen molecules into two types, those in the vibrational states at a lower level, whose vibrational energy is below a cutoff energy, and those in an upper level, with vibrational energy above the cutoff. Dissociation occurs at the upper level, and recombination returns molecules to the lower level. The model was applied to two flows of engineering interest, the flow through a normal Mach 15 shock wave at 60 km, and a supersonic quasi-one-dimensional flow in a nozzle. Results are compared to those obtained by existing translation-vibration-dissociation coupling models, with results indicating significant differences between the models.
Coupled ensemble flow line advection and analysis.
Guo, Hanqi; Yuan, Xiaoru; Huang, Jian; Zhu, Xiaomin
2013-12-01
Ensemble run simulations are becoming increasingly widespread. In this work, we couple particle advection with pathline analysis to visualize and reveal the differences among the flow fields of ensemble runs. Our method first constructs a variation field using a Lagrangian-based distance metric. The variation field characterizes the variation between vector fields of the ensemble runs, by extracting and visualizing the variation of pathlines within ensemble. Parallelism in a MapReduce style is leveraged to handle data processing and computing at scale. Using our prototype system, we demonstrate how scientists can effectively explore and investigate differences within ensemble simulations. PMID:24051840
External iterative coupling strategy for surface-subsurface flow calculations in surface irrigation
Technology Transfer Automated Retrieval System (TEKTRAN)
Coupling the unsteady open-channel flow equations of surface irrigation with the equation of variably saturated porous media flow is a computationally complex problem, because of the dependence of infiltration on flow depths. Several models of this coupled process have been developed, all of which ...
Dynamic coupling of three hydrodynamic models
NASA Astrophysics Data System (ADS)
Hartnack, J. N.; Philip, G. T.; Rungoe, M.; Smith, G.; Johann, G.; Larsen, O.; Gregersen, J.; Butts, M. B.
2008-12-01
The need for integrated modelling is evidently present within the field of flood management and flood forecasting. Engineers, modellers and managers are faced with flood problems which transcend the classical hydrodynamic fields of urban, river and coastal flooding. Historically the modeller has been faced with having to select one hydrodynamic model to cover all the aspects of the potentially complex dynamics occurring in a flooding situation. Such a single hydrodynamic model does not cover all dynamics of flood modelling equally well. Thus the ideal choice may in fact be a combination of models. Models combining two numerical/hydrodynamic models are becoming more standard, typically these models combine a 1D river model with a 2D overland flow model or alternatively a 1D sewer/collection system model with a 2D overland solver. In complex coastal/urban areas the flood dynamics may include rivers/streams, collection/storm water systems along with the overland flow. The dynamics within all three areas is of the same time scale and there is feedback in the system across the couplings. These two aspects dictate a fully dynamic three way coupling as opposed to running the models sequentially. It will be shown that the main challenges of the three way coupling are time step issues related to the difference in numerical schemes used in the three model components and numerical instabilities caused by the linking of the model components. MIKE FLOOD combines the models MIKE 11, MIKE 21 and MOUSE into one modelling framework which makes it possible to couple any combination of river, urban and overland flow fully dynamically. The MIKE FLOOD framework will be presented with an overview of the coupling possibilities. The flood modelling concept will be illustrated through real life cases in Australia and in Germany. The real life cases reflect dynamics and interactions across all three model components which are not possible to reproduce using a two-way coupling alone. The
NASA Astrophysics Data System (ADS)
Ali, M. A.; Ireson, A. M.; Keim, D.
2015-12-01
The Canadian prairies are cold and dry. Surface depressions are ubiquitous, and contain permanent or ephemeral ponds. The ponds are filled by snowmelt and precipitation on the ponds and lose a significant portion of their water to evaporation, but also, depending on their landscape position, may spill to other ponds or channels, recharge groundwater, or received groundwater discharge. Since precipitation and actual evaporation are closely balanced, the pond water balances are very sensitive to change in climate, and the prairies in general have been subject to damaging floods and droughts, in particular in the last decade or two. A 2.25 km2 field site at St Denis, central Saskatchewan, contains over 100 ponds, some permanent, some ephemeral, some saline, some fresh, some recharging groundwater, some receiving groundwater discharge. The site has been extensively studied for almost 50 years, with about one decade of continuous meteorological data, and three years of detailed pond level, soil moisture and temperature, and groundwater data. The objective of this study was to assess the performance of PARFLOW-CLM (a coupled land-atmosphere model) in simulating the pond-groundwater interactions at this site. Our conceptual model of the site includes soil properties that are progressively weathered with depth, and we implement this in a simplified dual permeability mathematical model of the soil hydraulic properties, whereby storage is dominated by the matrix and flow is dominated by macropores. The model performance was surprisingly good, doing quite a good job of capturing the observed groundwater and pond level dynamics. The soil freezing regime is also captured reasonably well, though the timing and pattern of the zero degree isotherm during soil thaw, which is critically important for runoff generation processes, was not captured as well. The model provides credible insights into the spatial patterns of evapotranspiration, and the seasonal dynamics of subsurface
Modeling of coupled geochemical and transport processes: An overview
Carnahan, C.L.
1989-10-01
Early coupled models associated with fluid flow and solute transport have been limited by assumed conditions of constant temperature, fully saturated fluid flow, and constant pore fluid velocity. Developments including coupling of chemical reactions to variable fields of temperature and fluid flow have generated new requirements for experimental data. As the capabilities of coupled models expand, needs are created for experimental data to be used for both input and validation. 25 refs.
Coupling of volatile transport and internal heat flow on Triton
NASA Technical Reports Server (NTRS)
Brown, Robert H.; Kirk, Randolph L.
1994-01-01
Recently Brown et al. (1991) showed that Triton's internal heat source could amount to 5-20% of the absorbed insolation on Triton, thus significantly affecting volatile transport and atmospheric pressure. Subsequently, Kirk and Brown (1991a) used simple analytical models of the effect of internal heat on the distribution of volatiles on Triton's surface, confirming the speculation of Brown et al. that Triton's internal heat flow could strongly couple to the surface volatile distribution. To further explore this idea, we present numerical models of the permanent distribution of nitrogen ice on Triton that include the effects of sunlight, the two-dimensional distribution of internal heat flow, the coupling of internal heat flow to the surface distribution of nitrogen ice, and the finite viscosity of nitrogen ice. From these models we conclude that: (1) The strong vertical thermal gradient induced in Triton's polar caps by internal heat-flow facilitates viscous spreading to lower latitudes, thus opposing the poleward transport of volatiles by sunlight, and, for plausible viscosities and nitrogen inventories, producing permanent caps of considerable latitudinal extent; (2) It is probable that there is a strong coupling between the surface distribution of nitrogen ice on Triton and internal heat flow; (3) Asymmetries in the spatial distribution of Triton's heat flow, possibly driven by large-scale, volcanic activity or convection in Triton's interior, can result in permanent polar caps of unequal latitudinal extent, including the case of only one permanent polar cap; (4) Melting at the base of a permanent polar cap on Triton caused by internal heat flow can significantly enhance viscous spreading, and, as an alternative to the solid-state greenhouse mechanism proposed by Brown et al. (1990), could provide the necessary energy, fluids, and/or gases to drive Triton's geyser-like plumes; (5) The atmospheric collapse predicted to occur on Triton in the next 20 years
NASA Astrophysics Data System (ADS)
Stock, Charles A.; Dunne, John P.; John, Jasmin G.
2014-01-01
Global-scale planktonic ecosystem models exhibit large differences in simulated net primary production (NPP) and assessment of planktonic food web fluxes beyond primary producers has been limited, diminishing confidence in carbon flux estimates from these models. In this study, a global ocean-ice-ecosystem model was assessed against a suite of observation-based planktonic food web flux estimates, many of which were not considered in previous modeling studies. The simulation successfully captured cross-biome differences and similarities in these fluxes after calibration of a limited number of highly uncertain yet influential parameters. The resulting comprehensive carbon budgets suggested that shortened food webs, elevated growth efficiencies, and tight consumer-resource coupling enable oceanic upwelling systems to support 45% of pelagic mesozooplankton production despite accounting for only 22% of ocean area and 34% of NPP. In seasonally stratified regions (42% of ocean area and 40% of NPP), weakened consumer-resource coupling tempers mesozooplankton production to 41% and enhances export below 100 m to 48% of the global total. In oligotrophic systems (36% of ocean area and 26% of NPP), the dominance of small phytoplankton and low consumer growth efficiencies supported only 14% of mesozooplankton production and 17% of export globally. Bacterial production, in contrast, was maintained in nearly constant proportion to primary production across biomes through the compensating effects of increased partitioning of NPP to the microbial food web in oligotrophic ecosystems and increased bacterial growth efficiencies in more productive areas. Cross-biome differences in mesozooplankton trophic level were muted relative to those invoked by previous work such that significant differences in consumer growth efficiencies and the strength of consumer-resource coupling were needed to explain sharp cross-biome differences in mesozooplankton production. Lastly, simultaneous
NASA Astrophysics Data System (ADS)
Garven, G.; Dumoulin, J. A.; Bradley, D. A.; Young, L. E.; Kelley, K. D.; Leach, D. L.
2002-12-01
Crustal heat flow can provide a strong mechanism for driving groundwater flow, particularly in submarine basins where other mechanisms for driving pore fluid flow such as topography, compaction and crustal deformation are too weak or too slow to have a significant effect on disturbing conductive heat flow. Fault zones appear to play a crucial role in focusing fluid migration in basins, as inferred in ancient rocks by many examples of hydrothermal deposits of sediment-hosted ores worldwide. Many rift-hosted deposits of lead, zinc, and barite ore appear to have formed at or near the seafloor by focused venting of hot basinal fluids and modified seawater, although the geophysical nature of these systems is not so well known. For example, the upper Kuna Formation, a finely laminated, black, organic-rich siliceous mudstone and shale in the Western Brooks Range of northwest Alaska, is host to the largest resources of zinc yet discovered in the Earth's crust, containing ore reserves in excess of 175 Mt averaging about 16% Zn and 5% Pb. Although situated today in a highly-deformed series of structural allocthonous plates thrusted during the Jurassic to Cretaceous Brookian Orogeny, the stratiform ores are thought to have formed much earlier in the anoxic, mud-rich Carboniferous-age Kuna Basin when adjacent carbonate platforms were drowned by rifting and tectonic subsidence. Fluid inclusion studies of ore-stage sphalerite and gangue minerals indicate sub-seafloor mineralization temperatures less than 200oC and most likely between 120 to 150 oC, during a period of sediment diagenesis and extensional faulting. We have constructed fully-coupled numerical models of heat and fluid flow to test hydrologic theories for free convection, submarine venting and subsequent ore formation, as constrained by paleoheat flow and petrologic observations. A finite element grid was designed and adapted for a cross section of the Kuna Basin, geologically restored to latest Mississippian time
A coupled heat and water flow apparatus
Mohamed, A.M.O.; Caporouscio, F.; Yong, R.N. ); Cheung, C.H. ); Kjartanson, B.H. )
1993-03-01
Safe and permanent disposal of radioactive waste requires isolation of a number of diverse chemical elements form the environment. The Canadian Nuclear Fuel Waste Management Program is assessing the concept of disposing of waste in a vault excavated at a depth of 500 to 1000 m below the ground surface in plutonic rock of the Canadian Shield. The temperatures and hydraulic potential in the buffer and back fill material were investigated. To study the performance of a compacted buffer material under thermal and isothermal conditions, a coupled heat and water flow apparatus is designed and presented. In the preliminary design, a one-dimensional flow of heat and water was not achieved. however, control of temperature gradient, existence of one-dimensional flow, and uniformity of temperature and volumetric water content distributions at any cross section within the specimen are achieved in the modified design. Experimental results have shown that the temperature stabilizes very rapidly after a period of approximately 0. 107 days. The moisture moves away from the hot end along the longitudinal direction of the specimen due to imposed thermal gradient. The time required for moisture to stabilize is in order of days. 17 refs., 17 figs., 3 tabs.
NASA Astrophysics Data System (ADS)
Beaud, Flavien; Flowers, Gwenn E.; Pimentel, Sam
2014-08-01
Field data and numerical modeling show that glaciations have the potential either to enhance relief or to dampen topography. We aim to model the effect of the subglacial hydraulic system on spatiotemporal patterns of glacial erosion by abrasion and quarrying on time scales commensurate with drainage system fluctuations (e.g., seasonal to annual). We use a numerical model that incorporates a dual-morphology subglacial drainage system coupled to a higher-order ice-flow model and process-specific erosion laws. The subglacial drainage system allows for a dynamic transition between two morphologies: the distributed system, characterized by an increase in basal water pressure with discharge, and the channelized system, which exhibits a decrease in equilibrium water pressure with increasing discharge. We apply the model to a simple synthetic glacier geometry, drive it with prescribed meltwater input variations, and compute sliding and erosion rates over a seasonal cycle. When both distributed and channelized systems are included, abrasion and sliding maxima migrate ~ 20% up-glacier compared to simulations with distributed drainage only. Power-law sliding generally yields to a broader response of abrasion to water pressure changes along the flowline compared to Coulomb-friction sliding. Multi-day variations in meltwater input elicit a stronger abrasion response than either diurnal- or seasonal variations alone for the same total input volume. An increase in water input volume leads to increased abrasion. We find that ice thickness commensurate with ice sheet outlet glaciers can hinder the up-glacier migration of abrasion. Quarrying patterns computed with a recently published law differ markedly from calculated abrasion patterns, with effective pressure being a stronger determinant than sliding speeds of quarrying rates. These variations in calculated patterns of instantaneous erosion as a function of hydrology-, sliding-, and erosion-model formulation, as well as model
NASA Astrophysics Data System (ADS)
Wang, Qiang; He, Zhu; Li, Baokuan; Tsukihashi, Fumitaka
2014-12-01
A transient three-dimensional finite-volume mathematical model has been developed to investigate the coupled physical fields in the electroslag remelting (ESR) process. Through equations solved by the electrical potential method, the electric current, electromagnetic force (EMF), and Joule heating fields are demonstrated. The mold is assumed to be conductive rather than insulated. The volume of fluid approach is implemented for the two-phase flow. Moreover, the EMF and Joule heating, which are the source terms of the momentum and energy sources, are recalculated at each iteration as a function of the phase distribution. The solidification is modeled by an enthalpy-porosity formulation, in which the mushy zone is treated as a porous medium with porosity equal to the liquid fraction. An innovative marking method of the metal pool profile is proposed in the experiment. The effect of the applied current on the ESR process is understood by the model. Good agreement is obtained between the experiment and calculation. The electric current flows to the mold lateral wall especially in the slag layer. A large amount of Joule heating around the metal droplet varies as it falls. The hottest region appears under the outer radius of the electrode tip, close to the slag/metal interface instead of the electrode tip. The metal pool becomes deeper with more power. The maximal temperature increases from 1951 K to 2015 K (1678 °C to 1742 °C), and the maximum metal pool depth increases from 34.0 to 59.5 mm with the applied current ranging from 1000 to 2000 A.
Punjabi, Sangeeta B.; Sahasrabudhe, S. N.; Das, A. K.; Joshi, N. K.; Mangalvedekar, H. A.; Kothari, D. C.
2014-01-15
This paper provides 2D comparative study of results obtained using laminar and turbulent flow model for RF (radio frequency) Inductively Coupled Plasma (ICP) torch. The study was done for the RF-ICP torch operating at 50 kW DC power and 3 MHz frequency located at BARC. The numerical modeling for this RF-ICP torch is done using ANSYS software with the developed User Defined Function. A comparative study is done between laminar and turbulent flow model to investigate how temperature and flow fields change when using different operating conditions such as (a) swirl and no swirl velocity for sheath gas flow rate, (b) variation in sheath gas flow rate, and (c) variation in plasma gas flow rate. These studies will be useful for different material processing applications.
NASA Astrophysics Data System (ADS)
Kihm, J.; Kim, J.
2006-12-01
A series of numerical simulations using a fully coupled hydrogeomechanical numerical model, which is named COWADE123D, is performed to analyze groundwater flow and land deformation in an unsaturated heterogeneous slope and its stability under various rainfall rates. The slope is located along a dam lake in Republic of Korea. The slope consists of the Cretaceous granodiorite and can be subdivided into the four layers such as weathered soil, weathered rock, intermediate rock, and hard rock from its ground surface due to weathering process. The numerical simulation results show that both rainfall rate and heterogeneity play important roles in controlling groundwater flow and land deformation in the unsaturated slope. The slope becomes more saturated, and thus its overall hydrogeomechanical stability deteriorates, especially in the weathered rock and weathered soil layers, as the rainfall increases up to the maximum daily rainfall rate in the return period of one year. However, the slope becomes fully saturated, and thus its hydrogeomechanical responses are almost identical under more than such a critical rainfall rate. From the viewpoint of hydrogeology, the pressure head, and hence the hydraulic head increase as the rainfall rate increases. As a result, the groundwater table rises, the unsaturated zone reduces, the seepage face expands from the slope toe toward the slope crest, and the groundwater flow velocity increases along the seepage face. Particularly, the groundwater flow velocity increases significantly in the weathered soil and weathered rock layers as the rainfall rate increases. This is because their hydraulic conductivity is relatively higher than that of the intermediate rock and hard rock layers. From the viewpoint of geomechanics, the horizontal displacement increases, while the vertical displacement decreases toward the slope toe as the rainfall rate increases. This may result from the buoyancy effect associated with the groundwater table rise as the
Vesselinov, Velimir V.; Pau, George; Finsterle, Stefan
2012-06-28
Agni is open source /community /multi platform framework for model-based analyses (SA, UQ, PE, RA, DS). Agni can be applied independently from other Advanced Simulation Capability for Environmental Management (ASCEM) modules. Agni will be available as Mercurial repository (hg clone https://akuna.labworks.org/hg/Platform) including source code, manual, test and verification examples.
NASA Astrophysics Data System (ADS)
Jha, B.; Plesch, A.; Shaw, J. H.; Hager, B. H.; Juanes, R.
2014-12-01
There has been a recent increase in the number of earthquakes reported in proximity of active oil and gas fields. In particular, the occurrence of a sequence of damaging earthquakes in May 2012 near the Cavone oil field, in Northern Italy, raised the question of whether these earthquakes might have been triggered, or, if not, if future activities might trigger other damaging events. Production and injection of fluids in the underground reservoirs are known to be capable of triggering seismicity by inducing slip on seismogenic faults. However, the effects of injection and production on fault stability in real fields are not always intuitively obvious, and require the development of new-generation coupled flow-geomechanical models that capture the effect of multiphase poromechanics on faults. We study, by way of numerical modeling and simulation, the potential for induced seismicity at the Cavone field. Using a coupled flow and geomechanics model of the field that honors reservoir geology and historical well schedule, we simulate oil production and water injection in the field for a period of three decades leading up to the earthquake sequence. We calculate the change in Coulomb stress on the bounding Mirandola fault, which sourced the May 29, 2012 M 5.8 earthquake. This quantity varies in space and evolves in time with changing pore pressure and total stress in the reservoir. A novel and important aspect of our work is the identification of a potential instability mechanism for a bounding fault at the edge of a reservoir experiencing pressure depletion. The discontinuity in pore pressure across the fault means that there is a discontinuity in effective normal stress and that, therefore, the Coulomb failure criterion must be evaluated locally on both sides of the fault. We track the evolution of the Coulomb stress at the earthquake hypocenter and compare it with the regional tectonic stressing rate to conclude in favor of tectonic origin of the earthquake. In
Flow Behavior Around Coupled, Rotating Turbines in Steady Flow
NASA Astrophysics Data System (ADS)
Fu, Matthew; Dabiri, John
2012-11-01
Counter-rotating vertical axis turbines (VATs) have been shown to yield increased power density in wind farms as compared to typical horizontal axis wind turbine (HAWT) farms. However, the governing physical mechanisms remain poorly understood. Scale model experiments in a free-surface water tunnel were conducted to characterize the effect of parameters such as turbine separation, tip speed ratio, and flow speed on the downstream flow field and the resulting vortex shedding from VATs. The flow field was visualized using particle image velocimetry (PIV) and planar laser induced fluorescence. The results are compared and contrasted with recent studies of counter-rotating circular cylinders to determine if suppression of vortex shedding plays a similarly important role in dictating the overall wake dynamics. This research was made possible through the generosity of Gordon and Betty Moore Foundation and the Caltech SURF Program.
NASA Technical Reports Server (NTRS)
Sinha, Neeraj; Brinckman, Kevin; Jansen, Bernard; Seiner, John
2011-01-01
A method was developed of obtaining propulsive base flow data in both hot and cold jet environments, at Mach numbers and altitude of relevance to NASA launcher designs. The base flow data was used to perform computational fluid dynamics (CFD) turbulence model assessments of base flow predictive capabilities in order to provide increased confidence in base thermal and pressure load predictions obtained from computational modeling efforts. Predictive CFD analyses were used in the design of the experiments, available propulsive models were used to reduce program costs and increase success, and a wind tunnel facility was used. The data obtained allowed assessment of CFD/turbulence models in a complex flow environment, working within a building-block procedure to validation, where cold, non-reacting test data was first used for validation, followed by more complex reacting base flow validation.
Uncertainty Analysis of Model Coupling
NASA Astrophysics Data System (ADS)
Held, H.; Knopf, B.; Schneider von Deimling, T.; Schellnhuber, H.-J.
The Earth System is a highly complex system that is often modelled by coupling sev- eral nonlinear submodules. For predicting the climate with these models, the following uncertainties play an essential role: parameter uncertainty, uncertainty in initial con- ditions or model uncertainty. Here we will address uncertainty in initial conditions as well as model uncertainty. As the process of coupling is an important part of model- ing, the main aspect of this work is the investigation of uncertainties that are due to the coupling process. For this study we use conceptual models that, compared to GCMs, have the advantage that the model itself as well as the output can be treated in a mathematically elabo- rated way. As the time for running the model is much shorter, the investigation is also possible for a longer period, e.g. for paleo runs. In consideration of these facts it is feasible to analyse the whole phase space of the model. The process of coupling is investigated by using different methods of examining low order coupled atmosphere-ocean systems. In the dynamical approach a fully coupled system of the two submodules can be compared to a system where one submodule forces the other. For a particular atmosphere-ocean system, based on the Lorenz model for the atmosphere, there can be shown significant differences in the predictability of a forced system depending whether the subsystems are coupled in a linear or a non- linear way. In [1] it is shown that in the linear case the forcing cannot represent the coupling, but in the nonlinear case, that we investigated in our study, the variability and the statistics of the coupled system can be reproduced by the forcing. Another approach to analyse the coupling is to carry out a bifurcation analysis. Here the bifurcation diagram of a single atmosphere system is compared to that of a cou- pled atmosphere-ocean system. Again it can be seen from the different behaviour of the coupled and the uncoupled system, that the
Chuan Lu; CHI Zhang; Hai Hanag; Timothy C. Johnson
2014-04-01
Successful geological storage and sequestration of carbon dioxide (CO2) require efficient monitoring of the migration of CO2 plume during and after large-scale injection in order to verify the containment of the injected CO2 within the target formation and to evaluate potential leakage risk. Field studies have shown that surface and cross-borehole electrical resistivity tomography (ERT) can be a useful tool in imaging and characterizing solute transport in heterogeneous subsurface. In this synthetic study, we have coupled a 3-D multiphase flow model with a parallel 3-D time-lapse ERT inversion code to explore the feasibility of using time-lapse ERT for simultaneously monitoring the migration of CO2 plume in deep saline formation and potential brine intrusion into shallow fresh water aquifer. Direct comparisons of the inverted CO2 plumes resulting from ERT with multiphase flow simulation results indicate the ERT could be used to delineate the migration of CO2 plume. Detailed comparisons on the locations, sizes and shapes of CO2 plume and intruded brine plumes suggest that ERT inversion tends to underestimate the area review of the CO2 plume, but overestimate the thickness and total volume of the CO2 plume. The total volume of intruded brine plumes is overestimated as well. However, all discrepancies remain within reasonable ranges. Our study suggests that time-lapse ERT is a useful monitoring tool in characterizing the movement of injected CO2 into deep saline aquifer and detecting potential brine intrusion under large-scale field injection conditions.
NASA Astrophysics Data System (ADS)
Smith, P. M.; Baker, L. J.; Asimow, P. D.; Gurnis, M. C.
2007-12-01
Seismic velocity and attenuation studies have shown that 5-20 km thick low velocity layers exist above seismically fast slabs and are associated with broad zones of high attenuation in many subduction zones. These observations are generally interpreted as formation of hydrous phases by dehydration of the slab, although the impact of water in nominally anhydrous minerals (NAM) on seismic wave propagation is largely unknown. Recent petrological experiments on hydrous peridotite at subduction zone conditions suggest that chlorite will be stable adjacent to the subducting slab in sufficient quantities to be a significant water sink. We use a scheme that couples a petrological model (pHMELTS) with a 2-D thermal and variable viscosity flow model (ConMan) to model energy and mass transfer within a subduction zone. By varying input parameters including the convergence rate and slab dip we have developed models for cases in the Costa-Rica and Izu- Bonin-Marianas arc systems and are able to predict major and trace element compositions of primary melts, as well as geophysical observables, such as the topography and geoid. We find that the emergence of a slab- adjacent low-viscosity channel (LVC) is a natural consequence of the thermal and chemical controls on mantle dynamics and feedback between them. In our earlier models, as the LVC is dragged downwards by the subducting slab, hornblende breaks down at about 2.5 GPa and other hydrous phases such as serpentine are secondary in importance to the NAM water reservoir. The spatial limit of the LVC is the water-saturated solidus of the hydrated peridotite; the LVC thickens as the peridotite is progressively depleted by melting and the solidus migrates into the warmer wedge, despite water replenishment at depth. pHMELTS is a hybrid of the pMELTS model of Ghiorso and co-workers and includes amphiboles, serpentines and micas. Chlorite was lacking but we have recently rectified this omission. Following De Capitani and co- workers, we
Modeling blood flow heterogeneity.
King, R B; Raymond, G M; Bassingthwaighte, J B
1996-01-01
It has been known for some time that regional blood flows within an organ are not uniform. Useful measures of heterogeneity of regional blood flows are the standard deviation and coefficient of variation or relative dispersion of the probability density function (PDF) of regional flows obtained from the regional concentrations of tracers that are deposited in proportion to blood flow. When a mathematical model is used to analyze dilution curves after tracer solute administration, for many solutes it is important to account for flow heterogeneity and the wide range of transit times through multiple pathways in parallel. Failure to do so leads to bias in the estimates of volumes of distribution and membrane conductances. Since in practice the number of paths used should be relatively small, the analysis is sensitive to the choice of the individual elements used to approximate the distribution of flows or transit times. Presented here is a method for modeling heterogeneous flow through an organ using a scheme that covers both the high flow and long transit time extremes of the flow distribution. With this method, numerical experiments are performed to determine the errors made in estimating parameters when flow heterogeneity is ignored, in both the absence and presence of noise. The magnitude of the errors in the estimates depends upon the system parameters, the amount of flow heterogeneity present, and whether the shape of the input function is known. In some cases, some parameters may be estimated to within 10% when heterogeneity is ignored (homogeneous model), but errors of 15-20% may result, even when the level of heterogeneity is modest. In repeated trials in the presence of 5% noise, the mean of the estimates was always closer to the true value with the heterogeneous model than when heterogeneity was ignored, but the distributions of the estimates from the homogeneous and heterogeneous models overlapped for some parameters when outflow dilution curves were
Force-coupling method for flows with ellipsoidal particles
NASA Astrophysics Data System (ADS)
Liu, D.; Keaveny, E. E.; Maxey, M. R.; Karniadakis, G. E.
2009-06-01
The force-coupling method, previously developed for spherical particles suspended in a liquid flow, is extended to ellipsoidal particles. In the limit of Stokes flow, there is an exact correspondence with known analytical results for isolated particles. More generally, the method is shown to provide good approximate results for the particle motion and the flow field both in viscous Stokes flow and at finite Reynolds number. This is demonstrated through comparison between fully resolved direct numerical simulations and results from the numerical implementation of the force-coupling method with a spectral/hp element scheme. The motion of settling ellipsoidal particles and neutrally buoyant particles in a Poiseuille flow are discussed.
Kolmogorov flow in two dimensional strongly coupled dusty plasma
Gupta, Akanksha; Ganesh, R. Joy, Ashwin
2014-07-15
Undriven, incompressible Kolmogorov flow in two dimensional doubly periodic strongly coupled dusty plasma is modelled using generalised hydrodynamics, both in linear and nonlinear regime. A complete stability diagram is obtained for low Reynolds numbers R and for a range of viscoelastic relaxation time τ{sub m} [0 < τ{sub m} < 10]. For the system size considered, using a linear stability analysis, similar to Navier Stokes fluid (τ{sub m} = 0), it is found that for Reynolds number beyond a critical R, say R{sub c}, the Kolmogorov flow becomes unstable. Importantly, it is found that R{sub c} is strongly reduced for increasing values of τ{sub m}. A critical τ{sub m}{sup c} is found above which Kolmogorov flow is unconditionally unstable and becomes independent of Reynolds number. For R < R{sub c}, the neutral stability regime found in Navier Stokes fluid (τ{sub m} = 0) is now found to be a damped regime in viscoelastic fluids, thus changing the fundamental nature of transition of Kolmogorov flow as function of Reynolds number R. A new parallelized nonlinear pseudo spectral code has been developed and is benchmarked against eigen values for Kolmogorov flow obtained from linear analysis. Nonlinear states obtained from the pseudo spectral code exhibit cyclicity and pattern formation in vorticity and viscoelastic oscillations in energy.
Unsteady flow analysis of a two-phase hydraulic coupling
NASA Astrophysics Data System (ADS)
Hur, N.; Kwak, M.; Lee, W. J.; Moshfeghi, M.; Chang, C.-S.; Kang, N.-W.
2016-06-01
Hydraulic couplings are being widely used for torque transmitting between separate shafts. A mechanism for controlling the transmitted torque of a hydraulic system is to change the amount of working fluid inside the system. This paper numerically investigates three-dimensional turbulent flow in a real hydraulic coupling with different ratios of charged working fluid. Working fluid is assumed to be water and the Realizable k-ɛ turbulence model together with the VOF method are used to investigate two-phase flow inside the wheels. Unsteady simulations are conducted using the sliding mesh technique. The primary wheel is rotating at a fixed speed of 1780 rpm and the secondary wheel rotates at different speeds for simulating different speed ratios. Results are investigated for different blade angles, speed ratios and also different water volume fractions, and are presented in the form of flow patterns, fluid average velocity and also torques values. According to the results, blade angle severely affects the velocity vector and the transmitted torque. Also in the partially-filled cases, air is accumulated in the center of the wheel forming a toroidal shape wrapped by water and the transmitted torque sensitively depends on the water volume fraction. In addition, in the fully-filled case the transmitted torque decreases as the speed ration increases and the average velocity associated with lower speed ratios are higher.
Strongly coupled turbulent gas-particle flows in vertical channels
NASA Astrophysics Data System (ADS)
Fox, Rodney O.; Capecelatro, Jesse; Desjardins, Olivier
2015-11-01
Eulerian-Lagrangian (EL) simulations of strongly coupled (high mass loading) gas-particle flows in vertical channels are performed with the purpose of exploring the fundamental physics of fully developed, wall-bounded multiphase turbulence. An adaptive spatial filter is developed that accurately decomposes the total granular energy of the particles into correlated and uncorrelated components at each location in the wall-normal direction of the flow. In this manner, Reynolds- and phase-averaged (PA) two-phase turbulence statistics up to second order are reported for both phases and for three values of the PA mean fluid velocity. As expected due to the high mass loading, in all cases the turbulence production due to mean drag dominates production due to mean shear. A multiphase LRR-IP Reynolds-stress turbulence model is developed to predict the turbulent flow statistics as a function of the wall-normal distance. Using a correlation for the vertical drift velocity developed from the EL data, the turbulence model predictions agree satisfactorily with all of one-point EL statistics for the vertical channel flows, as well as for the homogeneous cluster-induced turbulence (CIT) statistics reported previously. Funded by U.S. National Science Foundation (CBET-1437865).
Energy Science and Technology Software Center (ESTSC)
2011-08-24
T2Well/ECO2N is a coupled wellbore and reservoir model for simulating the dynamics of CO2 injection and leakage through wellbores. It can be seen as an extension to standard TOUGH/ECO2N V2.0, and can be applied to situations relevant to geologic CO2 storage involving upward flow (e.g., leakage) and downward flow (injection). The new simulator integrates a wellbore-reservoir system by assigning the wellbore and reservoir to two different sub-domains in which flow is controlled by appropriate physicalmore » laws. In the reservoir, we model flow using a standard multiphase Darcy flow approach. In the wellbores, we use the Drift-Flux Model and related conservation equations for describing transient two-phase non-isothermal wellbore flow of CO2-water mixtures. The mass and thermal energy balance equations are solved numerically by a finite difference scheme with wellbore heat transmission to the surrounding rock handled either semi-analytically or numerically. The momentum balance equation for the flow in the wellbore is solved numerically with a semi-explicit scheme.« less
Radiation-Spray Coupling for Realistic Flow Configurations
NASA Technical Reports Server (NTRS)
El-Asrag, Hossam; Iannetti, Anthony C.
2011-01-01
Three Large Eddy Simulations (LES) for a lean-direct injection (LDI) combustor are performed and compared. In addition to the cold flow simulation, the effect of radiation coupling with the multi-physics reactive flow is analyzed. The flame let progress variable approach is used as a subgrid combustion model combined with a stochastic subgrid model for spray atomization and an optically thin radiation model. For accurate chemistry modeling, a detailed Jet-A surrogate mechanism is utilized. To achieve realistic inflow, a simple recycling technique is performed at the inflow section upstream of the swirler. Good comparison is shown with the experimental data mean and root mean square profiles. The effect of combustion is found to change the shape and size of the central recirculation zone. Radiation is found to change the spray dynamics and atomization by changing the heat release distribution and the local temperature values impacting the evaporation process. The simulation with radiation modeling shows wider range of droplet size distribution by altering the evaporation rate. The current study proves the importance of radiation modeling for accurate prediction in realistic spray combustion configurations, even for low pressure systems.
Linear coupling of acoustic and cyclotron waves in plasma flows
Rogava, Andria; Gogoberidze, Grigol
2005-05-15
It is found that in magnetized electrostatic plasma flows the velocity shear couples ion-acoustic waves with ion-cyclotron waves and leads, under favorable conditions, to their efficient reciprocal transformations. It is shown that in a two-dimensional setup this coupling has a remarkable feature: it is governed by equations that are mathematically equal to the ones describing coupling of sound waves with internal gravity waves [Rogava and Mahajan, Phys. Rev. E 55, 1185 (1997)] in neutral fluids. For flows with low shearing rates a fully analytic, quantitative description of the coupling efficiency, based on a noteworthy quantum-mechanical analogy, is given and transformation coefficients are calculated.
Linear coupling of acoustic and cyclotron waves in plasma flows
NASA Astrophysics Data System (ADS)
Rogava, Andria; Gogoberidze, Grigol
2005-05-01
It is found that in magnetized electrostatic plasma flows the velocity shear couples ion-acoustic waves with ion-cyclotron waves and leads, under favorable conditions, to their efficient reciprocal transformations. It is shown that in a two-dimensional setup this coupling has a remarkable feature: it is governed by equations that are mathematically equal to the ones describing coupling of sound waves with internal gravity waves [Rogava and Mahajan, Phys. Rev. E 55, 1185 (1997)] in neutral fluids. For flows with low shearing rates a fully analytic, quantitative description of the coupling efficiency, based on a noteworthy quantum-mechanical analogy, is given and transformation coefficients are calculated.
Coupled flow, thermal and structural analysis of aerodynamically heated panels
NASA Technical Reports Server (NTRS)
Thornton, Earl A.; Dechaumphai, Pramote
1986-01-01
A finite element approach to coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.
Simplified coupling power model for fibers fusion
NASA Astrophysics Data System (ADS)
Saktioto, J.; Ali, J.; Fadhali, M.
2009-09-01
Fiber coupler fabrication used for an optical waveguide requires lossless power for an optimal application. The previous research coupled fibers were successfully fabricated by injecting hydrogen flow at 1 bar and fused slightly by unstable torch flame in the range of 800-1350°C. Optical parameters may vary significantly over wide range physical properties. Coupling coefficient and refractive index are estimated from the experimental result of the coupling ratio distribution from 1% to 75%. The change of geometrical fiber affects the normalized frequency V even for single mode fibers. V is derived and some parametric variations are performed on the left and right hand side of the coupling region. A partial power is modelled and derived using V, normalized lateral phase constant u, and normalized lateral attenuation constant, w through the second kind of modified Bessel function of the l order, which obeys the normal mode and normalized propagation constant b. Total power is maintained constant in order to comply with the energy conservation law. The power is integrated through V, u, and w over the pulling length of 7500 µm for 1-D. The core radius of a fiber significantly affects V and power partially at coupling region rather than wavelength and refractive index of core and cladding. This model has power phenomena in transmission and reflection for an optical switch and tunable filter.
Geomechanically Coupled Simulation of Flow in Fractured Reservoirs
NASA Astrophysics Data System (ADS)
Barton, C.; Moos, D.; Hartley, L.; Baxter, S.; Foulquier, L.; Holl, H.; Hogarth, R.
2012-12-01
Capturing the necessary and sufficient detail of reservoir hydraulics to accurately evaluate reservoir behavior remains a significant challenge to the exploitation and management of fracture-dominated geothermal reservoirs. In these low matrix permeability reservoirs, stimulation response is controlled largely by the properties of natural and induced fracture networks, which are in turn controlled by the in situ stresses, the fracture distribution and connectivity and the hydraulic behavior of the fractures. This complex interaction of fracture flow systems with the present-day stress field compounds the problem of developing an effective and efficient simulation to characterize, model and predict fractured reservoir performance. We discuss here a case study of the integration of geological, geophysical, geomechanical, and reservoir engineering data to characterize the in situ stresses, the natural fracture network and the controls on fracture permeability in geothermal reservoirs. A 3D geomechanical reservoir model includes constraints on stress magnitudes and orientations, and constraints on mechanical rock properties and the fractures themselves. Such a model is essential to understanding reservoir response to stimulation and production in low matrix permeability, fracture-dominated reservoirs. The geomechanical model for this study was developed using petrophysical, drilling, and wellbore image data along with direct well test measurements and was mapped to a 3D structural grid to facilitate coupled simulation of the fractured reservoir. Wellbore image and stimulation test data were used along with microseismic data acquired during the test to determine the reservoir fracture architecture and to provide control points for a realistic inter-connected discrete fracture network. As most fractures are stress-sensitive, their hydraulic conductivities will change with changes in bottomhole flowing and reservoir pressures, causing variations in production profiles
Session on coupled atmospheric/chemistry coupled models
NASA Technical Reports Server (NTRS)
Thompson, Anne
1993-01-01
The session on coupled atmospheric/chemistry coupled models is reviewed. Current model limitations, current issues and critical unknowns, and modeling activity are addressed. Specific recommendations and experimental strategies on the following are given: multiscale surface layer - planetary boundary layer - chemical flux measurements; Eulerian budget study; and Langrangian experiment. Nonprecipitating cloud studies, organized convective systems, and aerosols - heterogenous chemistry are also discussed.
Dynamics of coupled and uncoupled two-phase flows in a slab mold
NASA Astrophysics Data System (ADS)
Sánchez-Pérez, R.; García-Demedices, L.; Ramos, J. Palafox; Díaz-Cruz, M.; Morales, R. D.
2004-02-01
Two-phase flows in a mold of a slab caster are studied using water modeling, particle-image velocimetry (PIV), and computational fluid-dynamics techniques. Two-way coupled flows are observed in liquidgas systems, because both phases influence each other’s momentum transfer. In addition to this concept, PIV measurements indicate the existence of structurally coupled flows, where the velocity vectors of both phases observe similar orientations. When the drag forces of the liquid, exerted on the bubbles, exceed a certain value of the inertial forces of the liquid phase, at high mass loads of gas (ratio of mass flow rates of the gas phase and the liquid phase), the flow becomes structurally coupled. These types of flows promote large oscillations of the meniscus level. Two jets, liquid and bubble, were identified; the latter always reported larger angles than the first, independent of the gas load. Thus, a gas-rich jet is located closer to the lower edge of the submerged entry nozzle (SEN) port, and the liquid-rich jet is found above this position. The liquid-jet angle approaches that of the SEN port when the flow becomes structurally coupled. Structurally uncoupled flows report gas jets that follow torrent-type patterns which are well explained using a multiphase fluid-dynamics model. Structurally coupled flows yield gas jets with a continuous pattern.
Multiphysics and Multiscale Model Coupling Using Gerris
NASA Astrophysics Data System (ADS)
Keen, T. R.; Dykes, J. D.; Campbell, T. J.
2012-12-01
This work is implementing oceanographic processes encompassing multiple physics and scales using the Gerris Flow Solver (GFS) in order to examine their interdependence and sensitivity to changes in the physical environment. The processes include steady flow due to tides and the wind, phase-averaged wave-forced flow and oscillatory currents, and sediment transport. The 2D steady flow is calculated by the Ocean module contained within GFS. This model solves the Navier-Stokes (N-S) equations using the finite volume method. The model domain is represented by quad-tree adaptive mesh refinement (AMR). A stationary wave field is computed for a specified wave spectrum is uniformly distributed over the domain as a tracer with local wind input parameterized as a source, and dissipation by friction and breaking as a sink. Alongshore flow is included by a radiation stress term; this current is added to the steady flow component from tides and wind. Wave-current interaction is parameterized using a bottom boundary layer model. Sediment transport as suspended and bed load is implemented using tracers that are transported via the advection equations. A bed-conservation equation is implemented to allow changes in seafloor elevation to be used in adjusting the AMR refinement. These processes are being coupled using programming methods that are inherent to GFS and that do not require modification or recompiling of the code. These techniques include passive tracers, C functions that operate as plug-ins, and user-defined C-type macros included with GFS. Our results suggest that the AMR model coupling method is useful for problems where the dynamics are governed by several processes. This study is examining the relative influence of the steady currents, wave field, and sedimentation. Hydrodynamic and sedimentation interaction in nearshore environments is being studied for an idealized beach and for the Sandy Duck storm of Oct. 1998. The potential behavior of muddy sediments on the
Blood flow through sutured and coupled microvascular anastomoses: a comparative computational study.
Wain, Richard A J; Whitty, Justin P M; Dalal, Milind D; Holmes, Michael C; Ahmed, Waqar
2014-07-01
This study uses computational fluid dynamics (CFD) to model blood flow through idealised sutured and coupled arterial anastomoses to investigate the affect of each technique on intravascular blood flow. Local flow phenomena are examined in detail to study characteristics that potentially initiate thrombus formation; for example, changes in velocity profile, wall shear stress (WSS), and shear strain rate (SSR). Idealised geometries of sutured and coupled anastomoses were created with dimensions identical to microvascular suture material and a commercially available coupling device using CFD software. Vessels were modelled as non-compliant 1 mm diameter ducts, and blood was simulated as a Newtonian fluid, in keeping with previous studies. All analyses were steady-state and performed on arteries. The sutured simulation revealed a reduced boundary velocity profile; high WSS; and high SSR at the suture sites. The coupled anastomosis simulation showed a small increase in maximum WSS at the anastomotic region compared to a pristine vessel, however, this was less than half that of the sutured model. The coupled vessel displayed an average WSS equivalent to a pristine vessel simulation. Taken together these observations demonstrate a theoretically more thrombogenic profile in a sutured anastomosis when compared to a coupled vessel. Data from simulations on a coupled anastomosis reveal a profile that is nearly equivalent to that of a pristine vessel. Based purely on the combination of less favourable flow properties shown using these idealised arterial models, the sutured method is potentially more thrombogenic than a coupled anastomosis. PMID:24731801
NASA Astrophysics Data System (ADS)
Macías, D.; Martin, A. P.; García-Lafuente, J.; García, C. M.; Yool, A.; Bruno, M.; Vázquez-Escobar, A.; Izquierdo, A.; Sein, D. V.; Echevarría, F.
2007-08-01
The output of a two-layer hydrodynamic model along a west-east section of the Gibraltar Strait is used to estimate tidal induced mixing between the Mediterranean and Atlantic water layers and to simulate the effects of mixing processes on biogeochemical fluxes and the pelagic community of the area. The hydrodynamic model is used to estimate interfacial mixing and water advection which drive the dynamics of the pelagic community. The model was run for 13 months, in order to analyse the effect of annual modulations in tidal amplitude on mixing. Incorporation of a third intermediate layer leads to a significant improvement in the model results, showing the necessity for a three layer circulation scheme when modelling biogeochemical processes in the Strait of Gibraltar. Pelagic processes are modelled using a simple Nutrient-Phytoplankton-Zooplankton (NPZ) model. The intense physical mixing and advection in the channel are the main influence on plankton dynamics in the area. It is found that residence times within the channel are so short that phytoplankton communities cannot grow appreciably during their transit. As a consequence, the use of a more sophisticated biogeochemical model does not lead to significant changes in the results obtained. According to the model, mixing over the Camarinal Sill causes an average of 16% of the out-flowing nutrients to be returned back to the Mediterranean. This fraction varies between 4% and 35% as a function of the tidal amplitude. The comparison of the model results with field data suggests that in order to obtain an accurate simulation of the plankton ecosystem dynamics in the strait, it is necessary to take into account the full horizontal flow, as recirculation and coast-channel interactions seems to be very important processes in explaining the biological patterns in the area.
Garven, G.; Raffensperger, J.P.; Dumoulin, J.A.; Bradley, D.A.; Young, L.E.; Kelley, K.D.; Leach, D.L.
2003-01-01
The Red Dog deposit is a giant 175 Mton (16% Zn, 5% Pb), shale-hosted Pb-Zn-Ag-Ba ore district situated in the Carboniferous Kuna Basin, Western Brooks Range, Alaska. These SEDEX-type ores are thought to have formed in calcareous turbidites and black mudstone at elevated sub-seafloor temperatures (120-150??C) within a hydrogeologic framework of submarine convection that was structurally organized by large normal faults. The theory for modeling brine migration and heat transport in the Kuna Basin is discussed with application to evaluating flow patterns and heat transport in faulted rift basins and the effects of buoyancy-driven free convection on reactive flow and ore genesis. Finite element simulations show that hydrothermal fluid was discharged into the Red Dog subbasin during a period of basin-wide crustal heat flow of 150-160 mW/m2. Basinal brines circulated to depths as great as 1-3 km along multiple normal faults flowed laterally through thick clastic aquifers acquiring metals and heat, and then rapidly ascended a single discharge fault zone at rates ??? 5 m/year to mix with seafloor sulfur and precipitate massive sulfide ores. ?? 2003 Elsevier Science B.V. All rights reserved.
Multiscale Simulation Framework for Coupled Fluid Flow and Mechanical Deformation
Tchelepi, Hamdi
2014-11-14
A multiscale linear-solver framework for the pressure equation associated with flow in highly heterogeneous porous formations was developed. The multiscale based approach is cast in a general algebraic form, which facilitates integration of the new scalable linear solver in existing flow simulators. The Algebraic Multiscale Solver (AMS) is employed as a preconditioner within a multi-stage strategy. The formulations investigated include the standard MultiScale Finite-Element (MSFE) andMultiScale Finite-Volume (MSFV) methods. The local-stage solvers include incomplete factorization and the so-called Correction Functions (CF) associated with the MSFV approach. Extensive testing of AMS, as an iterative linear solver, indicate excellent convergence rates and computational scalability. AMS compares favorably with advanced Algebraic MultiGrid (AMG) solvers for highly detailed three-dimensional heterogeneous models. Moreover, AMS is expected to be especially beneficial in solving time-dependent problems of coupled multiphase flow and transport in large-scale subsurface formations.
Multiple steady states in coupled flow tank reactors
NASA Astrophysics Data System (ADS)
Hunt, Katharine L. C.; Kottalam, J.; Hatlee, Michael D.; Ross, John
1992-05-01
Coupling between continuous-flow, stirred tank reactors (CSTR's), each having multiple steady states, can produce new steady states with different concentrations of the chemical species in each of the coupled tanks. In this work, we identify a kinetic potential ψ that governs the deterministic time evolution of coupled tank reactors, when the reaction mechanism permits a single-variable description of the states of the individual tanks; examples include the iodate-arsenous acid reaction, a cubic model suggested by Noyes, and two quintic models. Stable steady states correspond to minima of ψ, and unstable steady states to maxima or saddle points; marginally stable states typically correspond to saddle-node points. We illustrate the variation in ψ due to changes in the rate constant for external material intake (k0) and for exchange between tanks (kx). For fixed k0 values, we analyze the changes in numbers and types of steady states as kx increases from zero. We show that steady states disappear by pairwise coalescence; we also show that new steady states may appear with increasing kx, when the reaction mechanism is sufficiently complex. For fixed initial conditions, the steady state ultimately reached in a mixing experiment may depend on the exchange rate constant as a function of time, kx(t) : Adiabatic mixing is obtained in the limit of slow changes in kx(t) and instantaneous mixing in the limit as kx(t)→∞ while t remains small. Analyses based on the potential ψ predict the outcome of mixing experiments for arbitrary kx(t). We show by explicit counterexamples that a prior theory developed by Noyes does not correctly predict the instability points or the transitions between steady states of coupled tanks, to be expected in mixing experiments. We further show that the outcome of such experiments is not connected to the relative stability of steady states in individual tank reactors. We find that coupling may effectively stabilize the tanks. We provide
The Challenges to Coupling Dynamic Geospatial Models
Goldstein, N
2006-06-23
Many applications of modeling spatial dynamic systems focus on a single system and a single process, ignoring the geographic and systemic context of the processes being modeled. A solution to this problem is the coupled modeling of spatial dynamic systems. Coupled modeling is challenging for both technical reasons, as well as conceptual reasons. This paper explores the benefits and challenges to coupling or linking spatial dynamic models, from loose coupling, where information transfer between models is done by hand, to tight coupling, where two (or more) models are merged as one. To illustrate the challenges, a coupled model of Urbanization and Wildfire Risk is presented. This model, called Vesta, was applied to the Santa Barbara, California region (using real geospatial data), where Urbanization and Wildfires occur and recur, respectively. The preliminary results of the model coupling illustrate that coupled modeling can lead to insight into the consequences of processes acting on their own.
Nicole Lautze
2015-01-01
Groundwater flow model for the island of Oahu. Data is from the following sources: Rotzoll, K., A.I. El-Kadi. 2007. Numerical Ground-Water Flow Simulation for Red Hill Fuel Storage Facilities, NAVFAC Pacific, Oahu, Hawaii - Prepared TEC, Inc. Water Resources Research Center, University of Hawaii, Honolulu.; Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume VII – Island of Oahu Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.; and Whittier, R. and A.I. El-Kadi. 2009. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. December 2009.
A novel potential/viscous flow coupling technique for computing helicopter flow fields
NASA Technical Reports Server (NTRS)
Summa, J. Michael; Strash, Daniel J.; Yoo, Sungyul
1990-01-01
Because of the complexity of helicopter flow field, a zonal method of analysis of computational aerodynamics is required. Here, a new procedure for coupling potential and viscous flow is proposed. An overlapping, velocity coupling technique is to be developed with the unique feature that the potential flow surface singularity strengths are obtained directly from the Navier-Stokes at a smoother inner fluid boundary. The closed-loop iteration method proceeds until the velocity field is converged. This coupling should provide the means of more accurate viscous computations of the near-body and rotor flow fields with resultant improved analysis of such important performance parameters as helicopter fuselage drag and rotor airloads.
Integrated Coupling of Surface and Subsurface Flow with HYDRUS-2D
NASA Astrophysics Data System (ADS)
Hartmann, Anne; Šimůnek, Jirka; Wöhling, Thomas; Schütze, Niels
2016-04-01
Describing interactions between surface and subsurface flow processes is important to adequately define water flow in natural systems. Since overland flow generation is highly influenced by rainfall and infiltration, both highly spatially heterogeneous processes, overland flow is unsteady and varies spatially. The prediction of overland flow needs to include an appropriate description of the interactions between the surface and subsurface flow. Coupling surface and subsurface water flow is a challenging task. Different approaches have been developed during the last few years, each having its own advantages and disadvantages. A new approach by Weill et al. (2009) to couple overland flow and subsurface flow based on a generalized Richards equation was implemented into the well-known subsurface flow model HYDRUS-2D (Šimůnek et al., 2011). This approach utilizes the one-dimensional diffusion wave equation to model overland flow. The diffusion wave model is integrated in HYDRUS-2D by replacing the terms of the Richards equation in a pre-defined runoff layer by terms defining the diffusion wave equation. Using this approach, pressure and flux continuity along the interface between both flow domains is provided. This direct coupling approach provides a strong coupling of both systems based on the definition of a single global system matrix to numerically solve the coupled flow problem. The advantage of the direct coupling approach, compared to the loosely coupled approach, is supposed to be a higher robustness, when many convergence problems can be avoided (Takizawa et al., 2014). The HYDRUS-2D implementation was verified using a) different test cases, including a direct comparison with the results of Weill et al. (2009), b) an analytical solution of the kinematic wave equation, and c) the results of a benchmark test of Maxwell et al. (2014), that included several known coupled surface subsurface flow models. Additionally, a sensitivity analysis evaluating the effects
NASA Astrophysics Data System (ADS)
Boadh, R.; Satyanarayana, A. N. V.; Rama Krishna, T. V. B. P. S.; Madala, S.
2014-12-01
Meso-scale atmospheric circulations play important role in air pollution transport and local air quality issues. Rapid industrialization and urbanization are leading to increased air pollution. In most of the metro cities, it is necessary to find the dispersion and ground level concentration (GLCs) of the pollutants emitted from the industries. In the present study AERMOD model has been used for simulation of pollutant dispersion over Nagpur, India. This model requires surface and upper air meteorological observations and various PBL parameters with good temporal resolution in the stand alone mode and mostly not available over India. To fill this gap, in the present study uses WRF-ARW model in getting these boundary layer parameters then will be offline coupled with AERMOD dispersion model. In the study we modified a pre-processor coupler for offline these models. High resolution simulations are conducted with triple nested domains having horizontal resolution of 27, 9 and 3 km; 27 vertical levels by using the 1x1 degree NCEP Final Analysis meteorological fields for initial and boundary conditions. In the study, eight fair weather days in winter and pre-monsoon season (January and April 2009). Sensitivity experiments of WRF-ARW model are conducted with two non-local [Yonsei University (YSU), Asymmetric Convective Model version 2 (ACM2)] and three local turbulence kinetic energy (TKE) closure [Mellor- Yamada Nakanishi and Niino Level 2.5 PBL (MYNN2), Mellor-Yamada-Janjic (MYJ), and quasi-normal scale elimination (QNSE)] turbulence diffusion parameterizations to study the evolution of PBL parameters and thermodynamic structure during the study period for providing the better meteorological parameters to AERMOD model. Sensitivity studies of PBL parameterization schemes after validation and statistical analysis reveal that the non-local PBL scheme YSU and local scheme MYJ could capture the characteristic variations of surface meteorological variables, vertical
Coupled transport in rotor models
NASA Astrophysics Data System (ADS)
Iubini, S.; Lepri, S.; Livi, R.; Politi, A.
2016-08-01
Steady nonequilibrium states are investigated in a one-dimensional setup in the presence of two thermodynamic currents. Two paradigmatic nonlinear oscillators models are investigated: an XY chain and the discrete nonlinear Schrödinger equation. Their distinctive feature is that the relevant variable is an angle in both cases. We point out the importance of clearly distinguishing between energy and heat flux. In fact, even in the presence of a vanishing Seebeck coefficient, a coupling between (angular) momentum and energy arises, mediated by the unavoidable presence of a coherent energy flux. Such a contribution is the result of the ‘advection’ induced by the position-dependent angular velocity. As a result, in the XY model, the knowledge of the two diagonal elements of the Onsager matrix suffices to reconstruct its transport properties. The analysis of the nonequilibrium steady states finally allows to strengthen the connection between the two models.
Inductively coupled plasma torch with laminar flow cooling
Rayson, Gary D.; Shen, Yang
1991-04-30
An improved inductively coupled gas plasma torch. The torch includes inner and outer quartz sleeves and tubular insert snugly fitted between the sleeves. The insert includes outwardly opening longitudinal channels. Gas flowing through the channels of the insert emerges in a laminar flow along the inside surface of the outer sleeve, in the zone of plasma heating. The laminar flow cools the outer sleeve and enables the torch to operate at lower electrical power and gas consumption levels additionally, the laminar flow reduces noise levels in spectroscopic measurements of the gaseous plasma.
NASA Technical Reports Server (NTRS)
Schulman, Richard; Kirk, Daniel; Marsell, Brandon; Roth, Jacob; Schallhorn, Paul
2013-01-01
The SPHERES Slosh Experiment (SSE) is a free floating experimental platform developed for the acquisition of long duration liquid slosh data aboard the International Space Station (ISS). The data sets collected will be used to benchmark numerical models to aid in the design of rocket and spacecraft propulsion systems. Utilizing two SPHERES Satellites, the experiment will be moved through different maneuvers designed to induce liquid slosh in the experiment's internal tank. The SSE has a total of twenty-four thrusters to move the experiment. In order to design slosh generating maneuvers, a parametric study with three maneuvers types was conducted using the General Moving Object (GMO) model in Flow-30. The three types of maneuvers are a translation maneuver, a rotation maneuver and a combined rotation translation maneuver. The effectiveness of each maneuver to generate slosh is determined by the deviation of the experiment's trajectory as compared to a dry mass trajectory. To fully capture the effect of liquid re-distribution on experiment trajectory, each thruster is modeled as an independent force point in the Flow-3D simulation. This is accomplished by modifying the total number of independent forces in the GMO model from the standard five to twenty-four. Results demonstrate that the most effective slosh generating maneuvers for all motions occurs when SSE thrusters are producing the highest changes in SSE acceleration. The results also demonstrate that several centimeters of trajectory deviation between the dry and slosh cases occur during the maneuvers; while these deviations seem small, they are measureable by SSE instrumentation.
Off-Centered Stagnation Point Flow of a Couple Stress Fluid towards a Rotating Disk
Khan, Najeeb Alam; Riaz, Fatima
2014-01-01
An investigation has been made to study the off-centered stagnation flow of a couple stress fluid over a rotating disk. The model developed for the governing problem in the form of partial differential equations has been converted to ordinary differential equations with the use of suitable similarity transformation. The analytical approximation has been made with the most promising analytical approach, homotopy analysis method (HAM). The convergence region of the obtained solution is determined and plotted. The effects of couple stress and nondimensional parameters have been observed on the flows of couple stress fluid. Also comparison has been made with the Newtonian fluid as the special case of considered problem. PMID:24672291
Regan, Robert S.; Niswonger, Richard G.; Markstrom, Steven L.; Barlow, Paul M.
2015-01-01
The spin-up simulation should be run for a sufficient length of time necessary to establish antecedent conditions throughout a model domain. Each GSFLOW application can require different lengths of time to account for the hydrologic stresses to propagate through a coupled groundwater and surface-water system. Typically, groundwater hydrologic processes require many years to come into equilibrium with dynamic climate and other forcing (or stress) data, such as precipitation and well pumping, whereas runoff-dominated surface-water processes respond relatively quickly. Use of a spin-up simulation can substantially reduce execution-time requirements for applications where the time period of interest is small compared to the time for hydrologic memory; thus, use of the restart option can be an efficient strategy for forecast and calibration simulations that require multiple simulations starting from the same day.
Swirl flow turbulence modeling
NASA Technical Reports Server (NTRS)
Abujelala, M. T.; Jackson, T. W.; Lilley, D. G.
1984-01-01
Confined turbulent swirling flow data obtained from a single hot-wire using a six-orientation technique are analyzed numerically. The effects of swirl strength and the presence of a strong contraction nozzle further downstream on deduced parameters is also presented and discussed for the case of chamber-to-inlet diameter ratio D/d = 2. Three swirl strengths are considered with inlet swirl vane angles of 0, 45 and 70 deg. A strong contraction nozzle with an area ratio of 4 is located two chamber-diameters downstream of the inlet to the flowfield. It is found that both the swirl strength and the contraction have strong effects on the turbulence parameters. Generally, the most dramatic effect of increase of swirl strength is the considerable increase in values of all the parameters considered, (rx-viscosity, kinetic energy of turbulence, length scales, and degree of nonisotropy). The presence of a strong contraction nozzle tends to increase the turbulence parameter values in regions of acceleration and to reduce them in deceleration regions. Based on similarity of viscosity and length scale profiles, a C sub mu formulation is deduced which is shown to improve the predictive capability of the standard k-epsilon turbulence model in swirling recirculating flows.
Y. Wu
2004-11-01
The purpose of this report is to document the unsaturated zone (UZ) flow models and submodels, as well as the flow fields that have been generated using the UZ flow model(s) of Yucca Mountain, Nevada. In this report, the term ''UZ model'' refers to the UZ flow model and the several submodels, which include tracer transport, temperature or ambient geothermal, pneumatic or gas flow, and geochemistry (chloride, calcite, and strontium) submodels. The term UZ flow model refers to the three-dimensional models used for calibration and simulation of UZ flow fields. This work was planned in the ''Technical Work Plan (TWP) for: Unsaturated Zone Flow Analysis and Model Report Integration'' (BSC 2004 [DIRS 169654], Section 1.2.7). The table of included Features, Events, and Processes (FEPs), Table 6.2-11, is different from the list of included FEPs assigned to this report in the ''Technical Work Plan for: Unsaturated Zone Flow Analysis and Model Report Integration'' (BSC 2004 [DIRS 169654], Table 2.1.5-1), as discussed in Section 6.2.6. The UZ model has revised, updated, and enhanced the previous UZ model (BSC 2001 [DIRS 158726]) by incorporating the repository design with new grids, recalibration of property sets, and more comprehensive validation effort. The flow fields describe fracture-fracture, matrix-matrix, and fracture-matrix liquid flow rates, and their spatial distributions as well as moisture conditions in the UZ system. These three-dimensional UZ flow fields are used directly by Total System Performance Assessment (TSPA). The model and submodels evaluate important hydrogeologic processes in the UZ as well as geochemistry and geothermal conditions. These provide the necessary framework to test hypotheses of flow and transport at different scales, and predict flow and transport behavior under a variety of climatic conditions. In addition, the limitations of the UZ model are discussed in Section 8.11.
Nicole Lautze
2015-01-01
Groundwater flow model for Kauai. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014.; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume IV – Island of Kauai Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2015.
NASA Astrophysics Data System (ADS)
Tertre, E.; Hubert, F.; Bruzac, S.; Pacreau, M.; Ferrage, E.; Prêt, D.
2013-07-01
The present study aims at testing the validity of using an Na+/Ca2+ ion-exchange model, derived from batch data to interpret experimental Ca2+-for-Na+ exchange breakthrough curves obtained on vermiculite (a common swelling clay mineral in surface environments). The ion-exchange model was constructed considering the multi-site nature of the vermiculite surface as well as the exchange of all aqueous species (Mg2+ derived from the dissolution of the solid and H+). The proposed ion-exchange model was then coupled with a transport model, and the predicted breakthrough curves were compared with the experimental ones obtained using a well stirred flow-through reactor. For a given solute residence time in the reactor (typically 50 min), our thermodynamic model based on instantaneous equilibrium was found to accurately reproduce several of the experimental breakthrough curves, depending on the Na+ and Ca2+ concentrations of the influents pumped through the reactor. However the model failed to reproduce experimental breakthrough curves obtained at high flow rates and low chemical gradient between the exchanger phase and the solution. An alternative model based on a hybrid equilibrium/kinetic approach was thus used and allowed predicting experimental data. Based on these results, we show that a simple parameter can be used to differentiate between thermodynamic and kinetic control of the exchange reaction with water flow. The results of this study are relevant for natural systems where two aquatic environments having contrasted chemistries interact. Indeed, the question regarding the attainment of a full equilibrium in such a system during the contact time of the aqueous phase with the particle/colloid remains most often open. In this context, we show that when a river (a flow of fresh water) encounters marine colloids, a systematic full equilibrium can be assumed (i.e., the absence of kinetic effects) when the residence time of the solute in 1 m3 of the system is ⩾6200 h.
Pan, L.; Oldenburg, C.M.; Wu, Y.-S.; Pruess, K.
2011-02-14
At its most basic level, the injection of CO{sub 2} into geologic CO{sub 2} storage sites involves a system comprising the wellbore and the target reservoir. The wellbore is the only conduit available to emplace CO{sub 2} into reservoirs for long-term storage. At the same time, wellbores in general have been identified as the most likely conduit for CO{sub 2} and brine leakage from geologic carbon sequestration (GCS) sites, especially those in sedimentary basins with historical hydrocarbon production. We have developed a coupled wellbore and reservoir model for simulating the dynamics of CO{sub 2} injection and leakage through wellbores. The model describes the following processes: (1) upward or downward wellbore flow of CO{sub 2} and variable salinity water with transition from supercritical to gaseous CO{sub 2} including Joule-Thomson cooling, (2) exsolution of CO{sub 2} from the aqueous phase as pressure drops, and (3) cross flow into or interaction with layers of surrounding rock (reservoirs). We use the Drift-Flux Model and related conservation equations for describing transient two-phase non-isothermal wellbore flow of CO{sub 2}-water mixtures under different flow regimes and interacting with surrounding rock. The mass and thermal energy balance equations are solved numerically by a finite difference scheme with wellbore heat transmission to the surrounding rock handled either semi-analytically or numerically. The momentum balance equation for the flow in the wellbore is solved numerically with a semi-explicit scheme. This manual provides instructions for compilation and use of the new model, and presents some example problems to demonstrate its use.
Coupled chemo(enzymatic) reactions in continuous flow
Yuryev, Ruslan; Strompen, Simon
2011-01-01
Summary This review highlights the state of the art in the field of coupled chemo(enzymatic) reactions in continuous flow. Three different approaches to such reaction systems are presented herein and discussed in view of their advantages and disadvantages as well as trends for their future development. PMID:22238518
Stochastic power flow modeling
Not Available
1980-06-01
The stochastic nature of customer demand and equipment failure on large interconnected electric power networks has produced a keen interest in the accurate modeling and analysis of the effects of probabilistic behavior on steady state power system operation. The principle avenue of approach has been to obtain a solution to the steady state network flow equations which adhere both to Kirchhoff's Laws and probabilistic laws, using either combinatorial or functional approximation techniques. Clearly the need of the present is to develop sound techniques for producing meaningful data to serve as input. This research has addressed this end and serves to bridge the gap between electric demand modeling, equipment failure analysis, etc., and the area of algorithm development. Therefore, the scope of this work lies squarely on developing an efficient means of producing sensible input information in the form of probability distributions for the many types of solution algorithms that have been developed. Two major areas of development are described in detail: a decomposition of stochastic processes which gives hope of stationarity, ergodicity, and perhaps even normality; and a powerful surrogate probability approach using proportions of time which allows the calculation of joint events from one dimensional probability spaces.
Coupled parametric design of flow control and duct shape
NASA Technical Reports Server (NTRS)
Florea, Razvan (Inventor); Bertuccioli, Luca (Inventor)
2009-01-01
A method for designing gas turbine engine components using a coupled parametric analysis of part geometry and flow control is disclosed. Included are the steps of parametrically defining the geometry of the duct wall shape, parametrically defining one or more flow control actuators in the duct wall, measuring a plurality of performance parameters or metrics (e.g., flow characteristics) of the duct and comparing the results of the measurement with desired or target parameters, and selecting the optimal duct geometry and flow control for at least a portion of the duct, the selection process including evaluating the plurality of performance metrics in a pareto analysis. The use of this method in the design of inter-turbine transition ducts, serpentine ducts, inlets, diffusers, and similar components provides a design which reduces pressure losses and flow profile distortions.
NASA Astrophysics Data System (ADS)
Trévisan, D.; Periáñez, R.
2016-03-01
The residence time of solutes in catchments is longer during low-flow conditions, due to the lengthening of transport routes and the decrease in transfer velocities. In rivers, transient storage depends largely on exchanges with channel storage and the hyporheic zone and reflects the capacity of the river to buffer pollutant loads before they enter the aquatic environment of final receptors. Our objective was to evaluate the fate of solutes along a typical confined river of upland catchments. First, we calculate lateral inflows using a variable-source hydrology approach. Then, water motion and quality in the river channel are predicted by combining hydrodynamics and exchanges with channel storage and the hyporheic zone. The model is mainly parametrized from literature data during baseflow conditions to mimic the fate of adsorptive and non-persistent pollutants. Residence time in surface water, channel storage and the hyporheic zone were found to be sensitive to lateral inflows from groundwater seepage. Channel storage is the main process controlling residence time in upstream conditions, where the riverbed is mainly composed of stones and bedrock. Downstream, along with the formation of sediment deposits and riffle-pool units, hyporheic exchanges also control the lag time in the transfer of solutes. By integrating physically-based processes, the number of parameters is small, but the model still requires a detailed description of stream geometry and morphology. It can be used to evaluate stream restoration or catchment-river management when detailed data of stream geometry and morphology are available.
NASA Astrophysics Data System (ADS)
Arnault, Joel; Wagner, Seven; Rummler, Thomas; Fersch, Benjamin; Bliefernicht, Jan; Andresen, Sabine; Kunstmann, Harald
2016-04-01
The analysis of land-atmosphere feedbacks requires detailed representation of land processes in atmospheric models. Our focus here is on runoff-infiltration partitioning and resolved overland flow. In the standard version of WRF, runoff-infiltration partitioning is described as a purely vertical process. In WRF-Hydro, runoff is enhanced with lateral water flows. The study region is the Sissili catchment (12800 km2) in West Africa and the study period March 2003 - February 2004. Our WRF setup includes an outer and inner domain at 10 and 2 km resolution covering the West African and Sissili region, respectively. In our WRF-Hydro setup the inner domain is coupled with a sub-grid at 500 m resolution to compute overland and river flow. Model results are compared with TRMM precipitation, MTE evapotranspiration, CCI soil moisture, CRU temperature, and streamflow observation. The role of runoff infiltration partitioning and resolved overland flow on land-atmosphere feedbacks is addressed with a sensitivity analysis of WRF results to the runoff-infiltration partitioning parameter and a comparison between WRF and WRF-Hydro results, respectively. In the outer domain, precipitation is sensitive to runoff-infiltration partitioning at the scale of the Sissili area (~100x100 km2), but not of area A (500x2500 km2). In the inner domain, where precipitation patterns are mainly prescribed by lateral boundary conditions, sensitivity is small, but additionally resolved overland flow here clearly increases infiltration and evapotranspiration at the beginning of the wet season when soils are still dry. Our WRF-Hydro setup shows potential for joint atmospheric and terrestrial water balance studies, and reproduces observed daily discharge with a Nash-Sutcliffe model efficiency coefficient of 0.43.
Coupling compositional liquid gas Darcy and free gas flows at porous and free-flow domains interface
NASA Astrophysics Data System (ADS)
Masson, R.; Trenty, L.; Zhang, Y.
2016-09-01
This paper proposes an efficient splitting algorithm to solve coupled liquid gas Darcy and free gas flows at the interface between a porous medium and a free-flow domain. This model is compared to the reduced model introduced in [6] using a 1D approximation of the gas free flow. For that purpose, the gas molar fraction diffusive flux at the interface in the free-flow domain is approximated by a two point flux approximation based on a low-frequency diagonal approximation of a Steklov-Poincaré type operator. The splitting algorithm and the reduced model are applied in particular to the modelling of the mass exchanges at the interface between the storage and the ventilation galleries in radioactive waste deposits.
Reduced Order Modeling Incompressible Flows
NASA Technical Reports Server (NTRS)
Helenbrook, B. T.
2010-01-01
The details: a) Need stable numerical methods; b) Round off error can be considerable; c) Not convinced modes are correct for incompressible flow. Nonetheless, can derive compact and accurate reduced-order models. Can be used to generate actuator models or full flow-field models
P. Dixon
2004-02-11
The purpose of this Model Report is to document the unsaturated zone (UZ) fluid flow and tracer transport models and submodels as well as the flow fields generated utilizing the UZ Flow and Transport Model of Yucca Mountain (UZ Model), Nevada. This work was planned in ''Technical Work Plan (TWP) for: Performance Assessment Unsaturated Zone'' (BSC 2002 [160819], Section 1.10, Work Package AUZM06). The UZ Model has revised, updated, and enhanced the previous UZ Flow Model REV 00 ICN 01 (BSC 2001 [158726]) by incorporation of the conceptual repository design with new grids, recalibration of property sets, and more comprehensive validation effort. The flow fields describe fracture-fracture, matrix-matrix, and fracture-matrix liquid flow rates and their spatial distributions as well as moisture conditions in the UZ system. These 3-D UZ flow fields are used directly by Performance Assessment (PA). The model and submodels evaluate important hydrogeologic processes in the UZ as well as geochemistry and geothermal conditions. These provide the necessary framework to test conceptual hypotheses of flow and transport at different scales and predict flow and transport behavior under a variety of climatic conditions. In addition, this Model Report supports several PA activities, including abstractions, particle-tracking transport simulations, and the UZ Radionuclide Transport Model.
NASA Astrophysics Data System (ADS)
Acharya, Ranadip; Bansal, Rohan; Gambone, Justin J.; Das, Suman
2014-12-01
Scanning laser epitaxy (SLE) is a new laser-based additive manufacturing technology under development at the Georgia Institute of Technology. SLE is aimed at the creation of equiaxed, directionally solidified, and single-crystal deposits of nickel-based superalloys through the melting of alloy powders onto superalloy substrates using a fast scanning Nd:YAG laser beam. The fast galvanometer control movement of the laser (0.2 to 2 m/s) and high-resolution raster scanning (20 to 200 µm line spacing) enables superior thermal control over the solidification process and allows the production of porosity-free, crack-free deposits of more than 1000 µm thickness. Here, we present a combined thermal and fluid flow model of the SLE process applied to alloy CMSX-4 with temperature-dependent thermo-physical properties. With the scanning beam described as a moving line source, the instantaneous melt pool assumes a convex hull shape with distinct leading edge and trailing edge characteristics. Temperature gradients at the leading and trailing edges are of order 2 × 105 and 104 K/m, respectively. Detailed flow analysis provides insights on the flow characteristics of the powder incorporating into the melt pool, showing velocities of order 1 × 10-4 m/s. The Marangoni effect drives this velocity from 10 to 15 times higher depending on the operating parameters. Prediction of the solidification microstructure is based on conditions at the trailing edge of the melt pool. Time tracking of solidification history is incorporated into the model to couple the microstructure prediction model to the thermal-fluid flow model, and to predict the probability of the columnar-to-equiaxed transition. Qualitative agreement is obtained between simulation and experimental result.
A multilingual programming model for coupled systems.
Ong, E. T.; Larson, J. W.; Norris, B.; Tobis, M.; Steder, M.; Jacob, R. L.; Mathematics and Computer Science; Univ. of Wisconsin; Univ. of Chicago; The Australian National Univ.
2008-01-01
Multiphysics and multiscale simulation systems share a common software requirement-infrastructure to implement data exchanges between their constituent parts-often called the coupling problem. On distributed-memory parallel platforms, the coupling problem is complicated by the need to describe, transfer, and transform distributed data, known as the parallel coupling problem. Parallel coupling is emerging as a new grand challenge in computational science as scientists attempt to build multiscale and multiphysics systems on parallel platforms. An additional coupling problem in these systems is language interoperability between their constituent codes. We have created a multilingual parallel coupling programming model based on a successful open-source parallel coupling library, the Model Coupling Toolkit (MCT). This programming model's capabilities reach beyond MCT's native Fortran implementation to include bindings for the C++ and Python programming languages. We describe the method used to generate the interlanguage bindings. This approach enables an object-based programming model for implementing parallel couplings in non-Fortran coupled systems and in systems with language heterogeneity. We describe the C++ and Python versions of the MCT programming model and provide short examples. We report preliminary performance results for the MCT interpolation benchmark. We describe a major Python application that uses the MCT Python bindings, a Python implementation of the control and coupling infrastructure for the community climate system model. We conclude with a discussion of the significance of this work to productivity computing in multidisciplinary computational science.
Coupling entropy of co-processing model on social networks
NASA Astrophysics Data System (ADS)
Zhang, Zhanli
2015-08-01
Coupling entropy of co-processing model on social networks is investigated in this paper. As one crucial factor to determine the processing ability of nodes, the information flow with potential time lag is modeled by co-processing diffusion which couples the continuous time processing and the discrete diffusing dynamics. Exact results on master equation and stationary state are achieved to disclose the formation. In order to understand the evolution of the co-processing and design the optimal routing strategy according to the maximal entropic diffusion on networks, we propose the coupling entropy comprehending the structural characteristics and information propagation on social network. Based on the analysis of the co-processing model, we analyze the coupling impact of the structural factor and information propagating factor on the coupling entropy, where the analytical results fit well with the numerical ones on scale-free social networks.
Modeling of Turbulent Swirling Flows
NASA Technical Reports Server (NTRS)
Shih, Tsan-Hsing; Zhu, Jiang; Liou, William; Chen, Kuo-Huey; Liu, Nan-Suey; Lumley, John L.
1997-01-01
Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.
Effects of coil location and injection flow rate in an inductively coupled RF plasma torch
NASA Astrophysics Data System (ADS)
Wei, D.; Apelian, D.; Farouk, B.
1985-07-01
A numerical model has been developed to investigate the effects of central carrier gas flow rate and coil location in an inductively coupled RF plasma torch. Solution algorithm is based on the primitive variable formulation of the Navier-Stokes equations and includes a pseudo two-dimensional electromagnetic field model. Computational results have shown that with increasing carrier gas flow rate, the plasma plume is penetrated and the back flow due to the magnetic pumping effects is diminished. This facilitates the delivery of powder particles into the discharge region. However, the plasma plume is also disturbed significantly thus enhancing power loss.
Power flow analysis of two coupled plates with arbitrary characteristics
NASA Technical Reports Server (NTRS)
Cuschieri, J. M.
1990-01-01
In the last progress report (Feb. 1988) some results were presented for a parametric analysis on the vibrational power flow between two coupled plate structures using the mobility power flow approach. The results reported then were for changes in the structural parameters of the two plates, but with the two plates identical in their structural characteristics. Herein, limitation is removed. The vibrational power input and output are evaluated for different values of the structural damping loss factor for the source and receiver plates. In performing this parametric analysis, the source plate characteristics are kept constant. The purpose of this parametric analysis is to determine the most critical parameters that influence the flow of vibrational power from the source plate to the receiver plate. In the case of the structural damping parametric analysis, the influence of changes in the source plate damping is also investigated. The results obtained from the mobility power flow approach are compared to results obtained using a statistical energy analysis (SEA) approach. The significance of the power flow results are discussed together with a discussion and a comparison between the SEA results and the mobility power flow results. Furthermore, the benefits derived from using the mobility power flow approach are examined.
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Johnston, Christopher O.; Thompson, Richard A.
2009-01-01
A description of models and boundary conditions required for coupling radiation and ablation physics to a hypersonic flow simulation is provided. Chemical equilibrium routines for varying elemental mass fraction are required in the flow solver to integrate with the equilibrium chemistry assumption employed in the ablation models. The capability also enables an equilibrium catalytic wall boundary condition in the non-ablating case. The paper focuses on numerical implementation issues using FIRE II, Mars return, and Apollo 4 applications to provide context for discussion. Variable relaxation factors applied to the Jacobian elements of partial equilibrium relations required for convergence are defined. Challenges of strong radiation coupling in a shock capturing algorithm are addressed. Results are presented to show how the current suite of models responds to a wide variety of conditions involving coupled radiation and ablation.
Statistical equilibria of the coupled barotropic flow and shallow water flow on a rotating sphere
NASA Astrophysics Data System (ADS)
Ding, Xueru
The motivation of this research is to build equilibrium statistical models that can apply to explain two enigmatic phenomena in the atmospheres of the solar system's planets: (1) the super-rotation of the atmospheres of slowly-rotating terrestrial planets---namely Venus and Titan, and (2) the persistent anticyclonic large vortex storms on the gas giants, such as the Great Red Spot (GRS) on Jupiter. My thesis is composed of two main parts: the first part focuses on the statistical equilibrium of the coupled barotropic vorticity flow (non-divergent) on a rotating sphere; the other one has to do with the divergent shallow water flow rotating sphere system. The statistical equilibria of these two systems are simulated in a wide range of parameter space by Monte Carlo methods based on recent energy-relative enstrophy theory and extended energy-relative enstrophy theory. These kind of models remove the low temperatures defect in the old classical doubly canonical energy-enstrophy theory which cannot support any phase transitions. The other big difference of our research from previous work is that we work on the coupled fluid-sphere system, which consists of a rotating high density rigid sphere, enveloped by a thin shell of fluid. The sphere is considered to have infinite mass and angular momentum; therefore, it can serve as a reservoir of angular momentum. Unlike the fluid sphere system itself, the coupled fluid sphere system allows for the exchange of angular momentum between the atmosphere and the solid planet. This exchange is the key point in any model that is expected to capture coherent structures such as the super-rotation and GRS-like vortices problems in planetary atmospheres. We discovered that slowly-rotating planets can have super-rotation at high energy state. All known slowly-rotating cases in the solar system---Venus and Titan---have super-rotation. Moreover, we showed that the anticyclonicity in the GRS-like structures is closely associated with the
Comparison of Coupled Radiative Flow Solutions with Project Fire 2 Flight Data
NASA Technical Reports Server (NTRS)
Olynick, David R.; Henline, W. D.; Chambers, Lin Hartung; Candler, G. V.
1995-01-01
A nonequilibrium, axisymmetric, Navier-Stokes flow solver with coupled radiation has been developed for use in the design or thermal protection systems for vehicles where radiation effects are important. The present method has been compared with an existing now and radiation solver and with the Project Fire 2 experimental data. Good agreement has been obtained over the entire Fire 2 trajectory with the experimentally determined values of the stagnation radiation intensity in the 0.2-6.2 eV range and with the total stagnation heating. The effects of a number of flow models are examined to determine which combination of physical models produces the best agreement with the experimental data. These models include radiation coupling, multitemperature thermal models, and finite rate chemistry. Finally, the computational efficiency of the present model is evaluated. The radiation properties model developed for this study is shown to offer significant computational savings compared to existing codes.
NASA Astrophysics Data System (ADS)
Ismail, Shuhaida; Shabri, Ani; Abadan, Siti Sarah
2015-02-01
This paper aims to investigate the ability of Empirical Mode Decompositio n (EMD) coupled with Least Square Support Vector Machine (LSSVM) model in order to improve the accuracy of river flow forecasting. To assess the effectiveness of this model, Bernam monthly river flow data, has served as the case study. The proposed model was set at three important stages which are decomposition, component identification and forecasting stages respectively. The first stage is known as decomposition stage where EMD were employed for decomposing the dataset into several numbers of Intrinsic Mode Functions (IMF) and a residue. During on second stage, the meaningful signals are identified using a statistical measure and the new dataset are obtained in this stage. The final stage applied LSSVM as a forecasting tool to perform the river flow forecasting. The performance of the EMD coupled with LSSVM model is compared with the single LSSVM models using various statistics measures of Mean Absolute Error (MAE), Root Mean Square Error (RMSE), correlation-coefficient (R) and Correlation of Efficiency (CE). The comparison results reveal the proposed model of EMD coupled with LSSVM model serves as a useful tool and a promising new method for river flow forecasting.
Mode-coupling effects in anisotropic flow in heavy-ion collisions
NASA Astrophysics Data System (ADS)
Qian, Jing; Heinz, Ulrich; Liu, Jia
2016-06-01
Higher-order anisotropic flows in heavy-ion collisions are affected by nonlinear mode coupling effects. It has been suggested that the associated nonlinear hydrodynamic response coefficients probe the transport properties and are largely insensitive to the spectrum of initial density fluctuations of the medium created in these collisions. To test this suggestion, we explore nonlinear mode coupling effects in event-by-event viscous fluid dynamics, using two different models for the fluctuating initial density profiles, and compare the nonlinear coupling coefficients between the initial eccentricity vectors before hydrodynamic expansion and the final flow vectors after the expansion. For several mode coupling coefficients we find significant sensitivity to the initial fluctuation spectrum. They all exhibit strong sensitivity to the specific shear viscosity at freeze-out, but only weak dependence on the shear viscosity during hydrodynamic evolution.
Centrifuge modelling of granular flows
NASA Astrophysics Data System (ADS)
Cabrera, Miguel Angel; Wu, Wei
2015-04-01
A common characteristic of mass flows like debris flows, rock avalanches and mudflows is that gravity is their main driving force. Gravity defines the intensity and duration of the main interactions between particles and their surrounding media (particle-particle, particle-fluid, fluid-fluid). At the same time, gravity delimits the occurrence of phase separation, inverse segregation, and mass consolidation, among other phenomena. Therefore, in the understanding of the flow physics it is important to account for the scaling of gravity in scaled models. In this research, a centrifuge model is developed to model free surface granular flows down an incline at controlled gravity conditions. Gravity is controlled by the action of an induced inertial acceleration field resulting from the rotation of the model in a geotechnical centrifuge. The characteristics of the induced inertial acceleration field during flow are discussed and validated via experimental data. Flow heights, velocity fields, basal pressure and impact forces are measured for a range of channel inclinations and gravity conditions. Preliminary results enlighten the flow characteristics at variable gravity conditions and open a discussion on the simulation of large scale processes at a laboratory scale. Further analysis on the flow physics brings valuable information for the validation of granular flows rheology.
Simple models for embayment flows
NASA Astrophysics Data System (ADS)
Gibson, F.; Dalziel, S.
2003-04-01
The flow structure in an embayment with a mean external flow has been investigated. The embayment is a relatively quiescent environment separated from the external mean flow by a mixing layer, in a manner analogous to the ventilation of a street canyon in an urban environment. This study aims to improve our knowledge of the exchange between the embayment and the external flow, which is an important mechanism for the transport and dispersion of substances such as nutrients, sediments, heat and pollutants. Understanding of flow in an embayment is therefore vital to the explanation and preservation of its ecology. In an experimental study, a model rectangular embayment was placed in a recirculating flume tank. The aspect ratio and bathymetry of the embayment was varied and the effect on the flow and mixing layer recorded. A neutrally buoyant tracer was added to the flow at various locations to visualise the eddies and the mixing layer. Field experiments in a coastal embayment used an accoustic Doppler current profiler to measure the flow velocities. These measurements demonstrate the existance of a gyre within the bay and support a shear-driven cavity model. In parallel with the experiments and fieldwork, a hierarchy of computer models was used to gain further understanding of the flow. Results from these models are presented alongside the experimental measurements.
Conceptual and Numerical Models for UZ Flow and Transport
H. Liu
2000-03-03
The purpose of this Analysis/Model Report (AMR) is to document the conceptual and numerical models used for modeling of unsaturated zone (UZ) fluid (water and air) flow and solute transport processes. This is in accordance with ''AMR Development Plan for U0030 Conceptual and Numerical Models for Unsaturated Zone (UZ) Flow and Transport Processes, Rev 00''. The conceptual and numerical modeling approaches described in this AMR are used for models of UZ flow and transport in fractured, unsaturated rock under ambient and thermal conditions, which are documented in separate AMRs. This AMR supports the UZ Flow and Transport Process Model Report (PMR), the Near Field Environment PMR, and the following models: Calibrated Properties Model; UZ Flow Models and Submodels; Mountain-Scale Coupled Processes Model; Thermal-Hydrologic-Chemical (THC) Seepage Model; Drift Scale Test (DST) THC Model; Seepage Model for Performance Assessment (PA); and UZ Radionuclide Transport Models.
A Dynamic Coupled Magnetosphere-Ionosphere-Ring Current Model
NASA Astrophysics Data System (ADS)
Pembroke, Asher
In this thesis we describe a coupled model of Earth's magnetosphere that consists of the Lyon-Fedder-Mobarry (LFM) global magnetohydrodynamics (MHD) simulation, the MIX ionosphere solver and the Rice Convection Model (RCM). We report some results of the coupled model using idealized inputs and model parameters. The algorithmic and physical components of the model are described, including the transfer of magnetic field information and plasma boundary conditions to the RCM and the return of ring current plasma properties to the LFM. Crucial aspects of the coupling include the restriction of RCM to regions where field-line averaged plasma-beta ¡=1, the use of a plasmasphere model, and the MIX ionosphere model. Compared to stand-alone MHD, the coupled model produces a substantial increase in ring current pressure and reduction of the magnetic field near the Earth. In the ionosphere, stronger region-1 and region-2 Birkeland currents are seen in the coupled model but with no significant change in the cross polar cap potential drop, while the region-2 currents shielded the low-latitude convection potential. In addition, oscillations in the magnetic field are produced at geosynchronous orbit with the coupled code. The diagnostics of entropy and mass content indicate that these oscillations are associated with low-entropy flow channels moving in from the tail and may be related to bursty bulk flows and bubbles seen in observations. As with most complex numerical models, there is the ongoing challenge of untangling numerical artifacts and physics, and we find that while there is still much room for improvement, the results presented here are encouraging. Finally, we introduce several new methods for magnetospheric visualization and analysis, including a fluid-spatial volume for RCM and a field-aligned analysis mesh for the LFM. The latter allows us to construct novel visualizations of flux tubes, drift surfaces, topological boundaries, and bursty-bulk flows.
Turbulence modeling for separated flow
NASA Technical Reports Server (NTRS)
Durbin, Paul A.
1994-01-01
Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.
Drift-Scale Coupled Processes (DST and THC Seepage) Models
E. Gonnenthal; N. Spyoher
2001-02-05
The purpose of this Analysis/Model Report (AMR) is to document the Near-Field Environment (NFE) and Unsaturated Zone (UZ) models used to evaluate the potential effects of coupled thermal-hydrologic-chemical (THC) processes on unsaturated zone flow and transport. This is in accordance with the ''Technical Work Plan (TWP) for Unsaturated Zone Flow and Transport Process Model Report'', Addendum D, Attachment D-4 (Civilian Radioactive Waste Management System (CRWMS) Management and Operating Contractor (M and O) 2000 [153447]) and ''Technical Work Plan for Nearfield Environment Thermal Analyses and Testing'' (CRWMS M and O 2000 [153309]). These models include the Drift Scale Test (DST) THC Model and several THC seepage models. These models provide the framework to evaluate THC coupled processes at the drift scale, predict flow and transport behavior for specified thermal loading conditions, and predict the chemistry of waters and gases entering potential waste-emplacement drifts. The intended use of this AMR is to provide input for the following: (1) Performance Assessment (PA); (2) Abstraction of Drift-Scale Coupled Processes AMR (ANL-NBS-HS-000029); (3) UZ Flow and Transport Process Model Report (PMR); and (4) Near-Field Environment (NFE) PMR. The work scope for this activity is presented in the TWPs cited above, and summarized as follows: continue development of the repository drift-scale THC seepage model used in support of the TSPA in-drift geochemical model; incorporate heterogeneous fracture property realizations; study sensitivity of results to changes in input data and mineral assemblage; validate the DST model by comparison with field data; perform simulations to predict mineral dissolution and precipitation and their effects on fracture properties and chemistry of water (but not flow rates) that may seep into drifts; submit modeling results to the TDMS and document the models. The model development, input data, sensitivity and validation studies described in
Drift-Scale Coupled Processes (DST and THC Seepage) Models
E. Sonnenthale
2001-04-16
The purpose of this Analysis/Model Report (AMR) is to document the Near-Field Environment (NFE) and Unsaturated Zone (UZ) models used to evaluate the potential effects of coupled thermal-hydrologic-chemical (THC) processes on unsaturated zone flow and transport. This is in accordance with the ''Technical Work Plan (TWP) for Unsaturated Zone Flow and Transport Process Model Report'', Addendum D, Attachment D-4 (Civilian Radioactive Waste Management System (CRWMS) Management and Operating Contractor (M&O) 2000 [1534471]) and ''Technical Work Plan for Nearfield Environment Thermal Analyses and Testing'' (CRWMS M&O 2000 [153309]). These models include the Drift Scale Test (DST) THC Model and several THC seepage models. These models provide the framework to evaluate THC coupled processes at the drift scale, predict flow and transport behavior for specified thermal loading conditions, and predict the chemistry of waters and gases entering potential waste-emplacement drifts. The intended use of this AMR is to provide input for the following: Performance Assessment (PA); Near-Field Environment (NFE) PMR; Abstraction of Drift-Scale Coupled Processes AMR (ANL-NBS-HS-000029); and UZ Flow and Transport Process Model Report (PMR). The work scope for this activity is presented in the TWPs cited above, and summarized as follows: Continue development of the repository drift-scale THC seepage model used in support of the TSPA in-drift geochemical model; incorporate heterogeneous fracture property realizations; study sensitivity of results to changes in input data and mineral assemblage; validate the DST model by comparison with field data; perform simulations to predict mineral dissolution and precipitation and their effects on fracture properties and chemistry of water (but not flow rates) that may seep into drifts; submit modeling results to the TDMS and document the models. The model development, input data, sensitivity and validation studies described in this AMR are required
NASA Astrophysics Data System (ADS)
Györgyi, László; Field, Richard J.
1989-11-01
Deterministic chaos is a well-established phenomenon in continuous-flow, stirred tank reactor (CSTR) experiments with the oscillatory Belousov-Zhabotinskii (BZ) reaction. However, it has not yet been possible to reproduce the experimentally observed, robust chaos in simulations using realistic models of the homogeneous chemical kinetics of the BZ reaction. That it may be necessary to consider spatial inhomogeneities in modeling the BZ chaos is suggested by the existence of strong stirring effects on the aperiodic behavior and by the difficulty of reproducing chaos under the same conditions in reactors of different physical configuration. Such inhomogeneity might result from a lack of perfect mixing in the CSTR, especially near the inlets, or from diffusion of species at low flow rates from the CSTR reaction mixture into the tips of the inlets. The presence of spatial inhomogeneities allows coupling between essentially independent oscillators, a well-known source of chaos. Such a model using a realistic representation of the BZ kinetics leads to an eight-variable set of ordinary differential equations. Numerical analysis of these equations by continuation methods and by numerical integration shows the existence of broad regions of chaos and various hysteresis effects involving limit cycles, a steady state and/or a strange attractor. Tristability was found in calculations in a narrow flow rate range. The computed behavior appears for parameter values closely related to the values used experimentally to obtain similar phenomena, and the visual similarity of the computed and experimental strange attractors is striking. The experimental routes to chaos, period doubling, intermittency, and secondary Hopf bifurcations are all reproduced in the simulations. The computed and experimental structures of periodic windows observed within chaotic regions also are very similar.
Modeling flow and sedimention of slurries
NASA Astrophysics Data System (ADS)
Mondy, L.; Rao, R.; Altobelli, S.; Ingber, M.; Graham, A.
2002-12-01
Many natural processes involve flows of sediments at high particle concentrations. The equations describing such two-phase flows are highly nonlinear. We will give an overview of the performance of a continuum constitutive model of suspensions of particles in liquid for low Reynolds number flows. The diffusive flux model (Leighton and Acrivos, J. Fluid Mech., 1987, and Phillips et al., Phys. Fluids A, 1992) is implemented in a general purpose finite element computational program. This constitutive description couples a Newtonian stress/shear-rate relationship (where the local viscosity of the suspension is dependent on the local volume fraction of solids) with a shear-induced migration model of the suspended particles. The momentum transport, continuity, and diffusive flux equations are solved simultaneously. The formulation is fully three-dimensional and can be run on a parallel computer platform. Recent work introducing a flow-aligned tensor correction to this model has had success in representing the anisotropic force that is seen in curvilinear flows. Gravity effects are added in an approach similar to that of Zhang and Acrivos (Int. J. Multiphase Flow, 1994). The model results are compared with laboratory data obtained with Nuclear Magnetic Resonance (NMR) of evolving particle concentration profiles in complex flows, as well as in batch sedimentation. Interesting secondary flows appear both in the experiment and model. Overall, good agreement is found between the experiments and the simulations. This work was supported by the United States Department of Energy under Contract DE-AC04- 94AL85000. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy. The authors would like to acknowledge support for this work by the U.S. Department of Energy, Division of Engineering and Geosciences, Office of Basic Energy Sciences.
Modeling Size Polydisperse Granular Flows
NASA Astrophysics Data System (ADS)
Lueptow, Richard M.; Schlick, Conor P.; Isner, Austin B.; Umbanhowar, Paul B.; Ottino, Julio M.
2014-11-01
Modeling size segregation of granular materials has important applications in many industrial processes and geophysical phenomena. We have developed a continuum model for granular multi- and polydisperse size segregation based on flow kinematics, which we obtain from discrete element method (DEM) simulations. The segregation depends on dimensionless control parameters that are functions of flow rate, particle sizes, collisional diffusion coefficient, shear rate, and flowing layer depth. To test the theoretical approach, we model segregation in tri-disperse quasi-2D heap flow and log-normally distributed polydisperse quasi-2D chute flow. In both cases, the segregated particle size distributions match results from full-scale DEM simulations and experiments. While the theory was applied to size segregation in steady quasi-2D flows here, the approach can be readily generalized to include additional drivers of segregation such as density and shape as well as other geometries where the flow field can be characterized including rotating tumbler flow and three-dimensional bounded heap flow. Funded by The Dow Chemical Company and NSF Grant CMMI-1000469.
High-resolution reactive transport: A coupled parallel hydrogeochemical model
NASA Astrophysics Data System (ADS)
Beisman, J. J.; Maxwell, R. M.; Steefel, C. I.; Sitchler, A.; Molins, S.
2013-12-01
Subsurface hydrogeochemical systems are an especially complex component of the terrestrial environment and play host to a multitude of interactions. Parameterizations of these interactions are perhaps the least understood component of terrestrial systems, presenting uncertainties in the predictive understanding of biogeochemical cycling and transport. Thorough knowledge of biogeochemical transport processes is critical to the quantification of carbon/nutrient fluxes in the subsurface, and to the development of effective contaminant remediation techniques. Here we present a coupled parallel hydrogeochemical model, ParCrunchFlow, as a tool to further our understanding of governing processes and interactions in natural hydrogeochemical systems. ParCrunchFlow is a coupling of the reactive transport simulator CrunchFlow with the hydrologic model ParFlow. CrunchFlow is a multicomponent reactive flow and transport code that can be used to simulate a range of important processes and environments, including reactive contaminant transport, chemical weathering, carbon sequestration, biogeochemical cycling, and water-rock interaction. ParFlow is a parallel, three-dimensional, variably-saturated, coupled surface-subsurface flow and transport code with the ability to simulate complex topography, geology, and heterogeneity. ParCrunchflow takes advantage of the efficient parallelism built into Parflow, allowing the numerical simulation of reactive transport processes in chemically and physically heterogeneous media at high spatial resolutions. This model provides an ability to further examine the interactions and feedbacks between biogeochemical systems and complex subsurface flow fields. In addition to the details of model construction, results will be presented that show floodplain nutrient cycling and the effects of heterogeneity on small-scale mixing reactions at the Department of Energy's Old Rifle Legacy site.
Convectively coupled Kelvin waves in CMIP5 coupled climate models
NASA Astrophysics Data System (ADS)
Wang, Lu; Li, Tim
2016-04-01
This study provided a quantitative evaluation of convectively coupled Kelvin waves (CCKWs) over the Indian Ocean and the Pacific Ocean simulated by 20 coupled climate models that participated in Coupled Model Intercomparison Project phase 5. The two leading empirical orthogonal function (EOF) modes of filtered daily precipitation anomalies are used to represent the eastward propagating CCKWs in both observations and simulations. The eigenvectors and eigenvalues of the EOF modes represent the spatial patterns and intensity of CCKWs respectively, and the lead-lag relationship between the two EOF principle components describe the phase propagation of CCKWs. A non-dimensional metric was designed in consideration of all the three factors (i.e., pattern, amplitude and phase propagation) for evaluation. The relative rankings of the models based on the skill scores calculated by the metric are conducted for the Indian Ocean and the Pacific Ocean, respectively. Two models (NorESM1-M and MPI-ESM-LR) are ranked among the best 20 % for both the regions. Three models (inmcm4, MRI-CGCM3 and HadGEM2-ES) are ranked among the worst 20 % for both the regions. While the observed CCKW amplitude is greater north of the equator in the Pacific, some models overestimate the CCKW ampliutde in the Southern Hemisphere. This bias is related to the mean state precipitation bias along the south Pacific convergence zone.
Groundwater flow and transport modeling
Konikow, L.F.; Mercer, J.W.
1988-01-01
Deterministic, distributed-parameter, numerical simulation models for analyzing groundwater flow and transport problems have come to be used almost routinely during the past decade. A review of the theoretical basis and practical use of groundwater flow and solute transport models is used to illustrate the state-of-the-art. Because of errors and uncertainty in defining model parameters, models must be calibrated to obtain a best estimate of the parameters. For flow modeling, data generally are sufficient to allow calibration. For solute-transport modeling, lack of data not only limits calibration, but also causes uncertainty in process description. Where data are available, model reliability should be assessed on the basis of sensitivity tests and measures of goodness-of-fit. Some of these concepts are demonstrated by using two case histories. ?? 1988.
HYDROGEN ELECTROLYZER FLOW DISTRIBUTOR MODEL
Shadday, M
2006-09-28
The hybrid sulfur process (HyS) hydrogen electrolyzer consists of a proton exchange membrane (PEM) sandwiched between two porous graphite layers. An aqueous solution of sulfuric acid with dissolved SO{sub 2} gas flows parallel to the PEM through the porous graphite layer on the anode side of the electrolyzer. A flow distributor, consisting of a number of parallel channels acting as headers, promotes uniform flow of the anolyte fluid through the porous graphite layer. A numerical model of the hydraulic behavior of the flow distributor is herein described. This model was developed to be a tool to aid the design of flow distributors. The primary design objective is to minimize spatial variations in the flow through the porous graphite layer. The hydraulic data from electrolyzer tests consists of overall flowrate and pressure drop. Internal pressure and flow distributions are not measured, but these details are provided by the model. The model has been benchmarked against data from tests of the current electrolyzer. The model reasonably predicts the viscosity effect of changing the fluid from water to an aqueous solution of 30 % sulfuric acid. The permeability of the graphite layer was the independent variable used to fit the model to the test data, and the required permeability for a good fit is within the range literature values for carbon paper. The model predicts that reducing the number of parallel channels by 50 % will substantially improve the uniformity of the flow in the porous graphite layer, while maintaining an acceptable pressure drop across the electrolyzer. When the size of the electrolyzer is doubled from 2.75 inches square to 5.5 inches square, the same number of channels as in the current design will be adequate, but it is advisable to increase the channel cross-sectional flow area. This is due to the increased length of the channels.
Full coupling of flow, thermal and mechanical effects in COMSOL for simulation of EGS
NASA Astrophysics Data System (ADS)
Sijacic, D.; Fokker, P. A.; van Wees, J.
2012-12-01
The effective modeling of geothermal reservoirs requires the coupling of geomechanics, fluid flow and thermal processes. Understanding of the complete system with these coupled processes is vital for reservoir stimulation targeted at enhancing reservoir performance and for the understanding, prediction and prevention of induced seismicity. The injection of cold water leads to alterations in strain and in situ stress in the reservoir, which can facilitate fracture initiation, opening, and activation of faults and joints which can lead to induced seismicity. Many commercial EGS sites have reported significant seismicity upon injection and production. However, thermal effects tend to be neglected in models for reservoir stimulation, although there is strong evidence that they can play an important role. We developed a model for EGS in COMSOL with a full coupling of flow, heat transfer and mechanics. Poroelasticity describes the influence of pore pressure on stress and strain, and in turn, the changes in stress and strain will change permeability and porosity, which influences fluid flow. Heat-convecting fluid flow (cold water injected in hot rock) will influence the temperature distribution, which in turn will influence the fluid viscosity (and density), altering again the flow itself. The temperature change will also create thermal stresses, effecting the geomechanics. In the geomechanical model we used Mohr-Coulomb failure and associated shear dilation for fault reactivation. The model is verified and benchmarked against existing models and is currently being related to actual EGS field operations in France and in Iceland. We have especially investigated the role of temperature changes on the stimulation and production stage of geothermal operations. The first step is to estimate temperature profiles around the injection well and the thermal stress effects. In the next step we compare a model that is fully coupled to a model where thermal effects are neglected and
Benchmark initiative on coupled multiphase flow and geomechanical processes during CO2 injection
NASA Astrophysics Data System (ADS)
Benisch, K.; Annewandter, R.; Olden, P.; Mackay, E.; Bauer, S.; Geiger, S.
2012-12-01
CO2 injection into deep saline aquifers involves multiple strongly interacting processes such as multiphase flow and geomechanical deformation, which threat to the seal integrity of CO2 repositories. Coupled simulation codes are required to establish realistic prognoses of the coupled process during CO2 injection operations. International benchmark initiatives help to evaluate, to compare and to validate coupled simulation results. However, there is no published code comparison study so far focusing on the impact of coupled multiphase flow and geomechanics on the long-term integrity of repositories, which is required to obtain confidence in the predictive capabilities of reservoir simulators. We address this gap by proposing a benchmark study. A wide participation from academic and industrial institutions is sought, as the aim of building confidence in coupled simulators become more plausible with many participants. Most published benchmark studies on coupled multiphase flow and geomechanical processes have been performed within the field of nuclear waste disposal (e.g. the DECOVALEX project), using single-phase formulation only. As regards CO2 injection scenarios, international benchmark studies have been published comparing isothermal and non-isothermal multiphase flow processes such as the code intercomparison by LBNL, the Stuttgart Benchmark study, the CLEAN benchmark approach and other initiatives. Recently, several codes have been developed or extended to simulate the coupling of hydraulic and geomechanical processes (OpenGeoSys, ELIPSE-Visage, GEM, DuMuX and others), which now enables a comprehensive code comparison. We propose four benchmark tests of increasing complexity, addressing the coupling between multiphase flow and geomechanical processes during CO2 injection. In the first case, a horizontal non-faulted 2D model consisting of one reservoir and one cap rock is considered, focusing on stress and strain regime changes in the storage formation and the
Modeling Reactive Transport in Coupled Groundwater-Conduit Systems
NASA Astrophysics Data System (ADS)
Spiessl, S. M.; Sauter, M.; Zheng, C.; Viswanathan, H. S.
2002-05-01
Modeling reactive transport in coupled groundwater-conduit systems requires consideration of two transport time scales in the flow and transport models. Consider for example a subsurface mine consisting of a network of highly conductive shafts, drifts or ventilation raises (i.e., conduits) within the considerably less permeable ore material (i.e., matrix). In the conduits, potential contaminants can travel much more rapidly than in the background aquifer (matrix). Since conduits cannot necessarily be regarded as a continuum, double continuum models are only of limited use for simulation of contaminant transport in such coupled groundwater-conduit systems. This study utilizes a "hybrid" flow and transport model in which contaminants can in essence be transported at a slower time scale in the matrix and at a faster time scale in the conduits. The hybrid flow model uses an approach developed by Clemens et al. (1996), which is based on the modelling of flow in a discrete pipe network, coupled to a continuum representing the low-permeability inter-conduit matrix blocks. Laminar or turbulent flow can be simulated in the different pipes depending on the flow conditions in the model domain. The three-dimensional finite-difference groundwater flow model MODFLOW (Harbaugh and McDonald, 1996) is used to simulate flow in the continuum. Contaminant transport within the matrix is simulated with a continuum approach using the three-dimensional multi-species solute transport model MT3DMS (Zheng and Wang, 1999), while that in the conduit system is simulated with a one-dimensional advective transport model. As a first step for reactive transport modeling in such systems, only equilibrium reactions among multiple species are considered by coupling the hybrid transport model to a geochemical speciation package. An idealized mine network developed by Viswanathan and Sauter (2001) is used as a test problem in this study. The numerical experiment is based on reference date collected from
Turbulence modeling for hypersonic flows
NASA Technical Reports Server (NTRS)
Marvin, J. G.; Coakley, T. J.
1989-01-01
Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.
Air-structure coupling features analysis of mining contra-rotating axial flow fan cascade
NASA Astrophysics Data System (ADS)
Chen, Q. G.; Sun, W.; Li, F.; Zhang, Y. J.
2013-12-01
The interaction between contra-rotating axial flow fan blade and working gas has been studied by means of establishing air-structure coupling control equation and combining Computational Fluid Dynamics (CFD) and Computational solid mechanics (CSM). Based on the single flow channel model, the Finite Volume Method was used to make the field discrete. Additionally, the SIMPLE algorithm, the Standard k-ε model and the Arbitrary Lagrangian-Eulerian dynamic grids technology were utilized to get the airflow motion by solving the discrete governing equations. At the same time, the Finite Element Method was used to make the field discrete to solve dynamic response characteristics of blade. Based on weak coupling method, data exchange from the fluid solver and the solid solver was processed on the coupling interface. Then interpolation was used to obtain the coupling characteristics. The results showed that the blade's maximum amplitude was on the tip of the last-stage blade and aerodynamic force signal could reflect the blade working conditions to some extent. By analyzing the flow regime in contra-rotating axial flow fan, it could be found that the vortex core region was mainly in the blade surface, the hub and the blade clearance. In those regions, the turbulence intensity was very high. The last-stage blade's operating life is shorter than that of the pre-stage blade due to the fatigue fracture occurs much more easily on the last-stage blade which bears more stress.
Internal thermal coupling in direct-flow coaxial vacuum tube collectors
Glembin, J.; Rockendorf, G.; Scheuren, J.
2010-07-15
This investigation covers the impact of low flow rates on the efficiency of coaxial vacuum tube collectors. Measurements show an efficiency reduction of 10% if reducing the flow rate from 78 kg/m{sup 2} h to 31 kg/m{sup 2} h for a collector group with 60 parallel vacuum tubes with a coaxial flow conduit at one-sided connection. For a more profound understanding a model of the coaxial tube was developed which defines the main energy fluxes including the internal thermal coupling. The tube simulations show a non-linear temperature profile along the tube with the maximum temperature in the outer pipe. Due to heat transfer to the entering flow this maximum is not located at the fluid outlet. The non-linearity increases with decreasing flow rates. The experimentally determined flow distribution allows simulating the measured collector array. The simulation results confirm the efficiency decrease at low flow rates. The flow distribution has a further impact on efficiency reduction, but even at an ideal uniform flow, a considerable efficiency reduction at low flow rates is to be expected. As a consequence, low flow rates should be prevented for coaxial tube collectors, thus restricting the possible operation conditions. The effect of constructional modifications like diameter or material variations is presented. Finally the additional impact of a coaxial manifold design is discussed. (author)
NASA Technical Reports Server (NTRS)
Demuren, A. O.
1994-01-01
Various approaches to the modeling of jets in cross flow are reviewed. These are grouped into four classes, namely: empirical models, integral models, perturbation models, and numerical models. Empirical models depend largely on the correlation of experimental data and are mostly useful for first-order estimates of global properties such as jet trajectory and velocity and temperature decay rates. Integral models are based on some ordinary-differential form of the conservation laws, but require substantial empirical calibration. They allow more details of the flow field to be obtained; simpler versions have to assume similarity of velocity and temperature profiles, but more sophisticated ones can actually calculate these profiles. Perturbation models require little empirical input, but the need for small parameters to ensure convergent expansions limits their application to either the near-field or the far-field. Therefore, they are mostly useful for the study of flow physics. Numerical models are based on conservation laws in partial-differential form. They require little empirical input and have the widest range of applicability. They also require the most computational resources. Although many qualitative and quantitative features of jets in cross flow have been predicted with numerical models, many issues affecting accuracy such as grid resolution and turbulence model are not completely resolved.
A strongly conservative finite element method for the coupling of Stokes and Darcy flow
NASA Astrophysics Data System (ADS)
Kanschat, G.; Rivière, B.
2010-08-01
We consider a model of coupled free and porous media flow governed by Stokes and Darcy equations with the Beavers-Joseph-Saffman interface condition. This model is discretized using divergence-conforming finite elements for the velocities in the whole domain. Discontinuous Galerkin techniques and mixed methods are used in the Stokes and Darcy subdomains, respectively. This discretization is strongly conservative in Hdiv( Ω) and we show convergence. Numerical results validate our findings and indicate optimal convergence orders.
Testing the global flow reconstruction method on coupled chaotic oscillators
NASA Astrophysics Data System (ADS)
Plachy, Emese; Kolláth, Zoltán
2010-03-01
Irregular behaviour of pulsating variable stars may occur due to low dimensional chaos. To determine the quantitative properties of the dynamics in such systems, we apply a suitable time series analysis, the global flow reconstruction method. The robustness of the reconstruction can be tested through the resultant quantities, like Lyapunov dimension and Fourier frequencies. The latter is specially important as it is directly derivable from the observed light curves. We have performed tests using coupled Rossler oscillators to investigate the possible connection between those quantities. In this paper we present our test results.
Coupling surface and subsurface flows with curved interfaces
NASA Astrophysics Data System (ADS)
Song, Pu; Yotov, Ivan
2013-11-01
A mortar multiscale method is developed for the coupled Stokes andDarcy flows with the Beavers-Joseph-Saffman interface condition in irregular domains. Conforming Stokes elements and multipoint flux mixed finite elements in Darcy are used to discretize the subdomains on the fine scale. A coarse scale mortar finite element space is used to approximate interface stresses and pressures and impose weakly continuity of velocities and fluxes. Matching conditions on curved interfaces are imposed by mapping the physical grids to reference grids with flat interfaces.
On the coupled evolution of oceanic internal waves and quasi-geostrophic flow
NASA Astrophysics Data System (ADS)
Wagner, Gregory LeClaire
Oceanic motion outside thin boundary layers is primarily a mixture of quasi-geostrophic flow and internal waves with either near-inertial frequencies or the frequency of the semidiurnal lunar tide. This dissertation seeks a deeper understanding of waves and flow through reduced models that isolate their nonlinear and coupled evolution from the Boussinesq equations. Three physical-space models are developed: an equation that describes quasi-geostrophic evolution in an arbitrary and prescribed field of hydrostatic internal waves; a three-component model that couples quasi-geostrophic flow to both near-inertial waves and the near-inertial second harmonic; and a model for the slow evolution of hydrostatic internal tides in quasi-geostrophic flow of near-arbitrary scale. This slow internal tide equation opens the path to a coupled model for the energetic interaction of quasi-geostrophic flow and oceanic internal tides. Four results emerge. First, the wave-averaged quasi-geostrophic equation reveals that finite-amplitude waves give rise to a mean flow that advects quasi-geostrophic potential vorticity. Second is the definition of a new material invariant: Available Potential Vorticity, or APV. APV isolates the part of Ertel potential vorticity available for balanced-flow evolution in Eulerian frames and proves necessary in the separating waves and quasi-geostrophic flow. The third result, hashed out for near-inertial waves and quasi-geostrophic flow, is that wave-flow interaction leads to energy exchange even under conditions of weak nonlinearity. For storm-forced oceanic near-inertial waves the interaction often energizes waves at the expense of flow. We call this extraction of balanced quasi-geostrophic energy 'stimulated generation' since it requires externally-forced rather than spontaneously-generated waves. The fourth result is that quasi-geostrophic flow can encourage or 'catalyze' a nonlinear interaction between a near-inertial wave field and its second harmonic
Effects of dynamically variable saturation and matrix-conduit coupling of flow in karst aquifers
Reimann, T.; Geyer, T.; Shoemaker, W.B.; Liedl, R.; Sauter, M.
2011-01-01
Well-developed karst aquifers consist of highly conductive conduits and a relatively low permeability fractured and/or porous rock matrix and therefore behave as a dual-hydraulic system. Groundwater flow within highly permeable strata is rapid and transient and depends on local flow conditions, i.e., pressurized or nonpressurized flow. The characterization of karst aquifers is a necessary and challenging task because information about hydraulic and spatial conduit properties is poorly defined or unknown. To investigate karst aquifers, hydraulic stresses such as large recharge events can be simulated with hybrid (coupled discrete continuum) models. Since existing hybrid models are simplifications of the system dynamics, a new karst model (ModBraC) is presented that accounts for unsteady and nonuniform discrete flow in variably saturated conduits employing the Saint-Venant equations. Model performance tests indicate that ModBraC is able to simulate (1) unsteady and nonuniform flow in variably filled conduits, (2) draining and refilling of conduits with stable transition between free-surface and pressurized flow and correct storage representation, (3) water exchange between matrix and variably filled conduits, and (4) discharge routing through branched and intermeshed conduit networks. Subsequently, ModBraC is applied to an idealized catchment to investigate the significance of free-surface flow representation. A parameter study is conducted with two different initial conditions: (1) pressurized flow and (2) free-surface flow. If free-surface flow prevails, the systems is characterized by (1) a time lag for signal transmission, (2) a typical spring discharge pattern representing the transition from pressurized to free-surface flow, and (3) a reduced conduit-matrix interaction during free-surface flow. Copyright 2011 by the American Geophysical Union.
Coupling between pre-onset flows and substorm onset waves
NASA Astrophysics Data System (ADS)
Nishimura, T.; Lyons, L. R.; Angelopoulos, V.; Donovan, E.; Mende, S. B.
2015-12-01
A critical, long-standing problem in substorm research is identification of the sequence of events leading to substorm expansion phase onset. Recent THEMIS all-sky imager (ASI) array observations have shown a repeatable pre-onset sequence, which is initiated by a poleward boundary intensification (PBI) and is followed by auroral streamers moving equatorward (earthward flow in the plasma sheet) and then by substorm onset. On the other hand, substorm onset is also preceded by azimuthally propagating waves, indicating a possible importance of wave instability for triggering substorm onset. However, it has been difficult to identify the link between fast flows and waves. We have found an isolated substorm event that was well-instrumented with the Poker Flat incoherent scatter radar (PFISR), THEMIS white-light ASI, and multi-spectral ASI, where the auroral onset occurred within the PFISR and ASI fields-of-view. This substorm onset was preceded by a PBI, and ionospheric flows propagated equatorward from the polar cap, crossed the PBI and reached the growth phase arc. This sequence provides evidence that flows from open magnetic field lines propagate across the open-closed boundary and reach the near-Earth plasma sheet prior to the onset. Quasi-stable oscillations in auroral luminosity and ionospheric density are found along the growth phase arc. These pre-onset auroral waves amplified abruptly at the onset time, soon after the equatorward flows reached the onset region. This sequence suggests a coupling process where pre-existing stable waves in the near-Earth plasma sheet interact with flows from further downtail and then evolve to onset instability.
Energy structure of MHD flow coupling with outer resistance circuit
NASA Astrophysics Data System (ADS)
Huang, Z. Y.; Liu, Y. J.; Chen, Y. Q.; Peng, Z. L.
2015-08-01
Energy structure of MHD flow coupling with outer resistance circuit is studied to illuminate qualitatively and quantitatively the energy relation of this basic MHD flow system with energy input and output. Energy structure are analytically derived based on the Navier-Stocks equations for two-dimensional fully-developed flow and generalized Ohm's Law. The influences of applied magnetic field, Hall parameter and conductivity on energy structure are discussed based on the analytical results. Associated energies in MHD flow are deduced and validated by energy conservation. These results reveal that energy structure consists of two sub structures: electrical energy structure and internal energy structure. Energy structure and its sub structures provide an integrated theoretical energy path of the MHD system. Applied magnetic field and conductivity decrease the input energy, dissipation by fluid viscosity and internal energy but increase the ratio of electrical energy to input energy, while Hall parameter has the opposite effects. These are caused by their different effects on Bulk velocity, velocity profiles, voltage and current in outer circuit. Understanding energy structure helps MHD application designers to actively adjust the allocation of different parts of energy so that it is more reasonable and desirable.
Coupling environmental models and geospatial data processing
NASA Astrophysics Data System (ADS)
Brandmeyer, Jo Ellen
2000-10-01
This research investigated geospatial functions for solving environmental problems from the perspective of the environmental modeler. Its purpose is to better understand the different approaches to coupling complex models and geospatial data processing, plus the implications for the coupled system. To this end, various coupling methodologies were systematically explored using a geographic information system (GIS) and an emissions processor (SMOKE) for air quality models (AQMs). SMOKE converts an emissions inventory into the format required by an AQM. A GIS creates a file describing the spatial distribution of emissions among the cells in a modeling domain. To demonstrate advantages of a coupled GIS---environmental model system, two methods of spatially distributing on-road mobile emissions to cells were examined. The existing method calculates emissions for each road class, but distributes emissions to the cells using population density. For the new method a GIS builds road density by class and then distributes the emissions using road density. Comparing these methods reveals a significantly different spatial pattern of emissions. Next, various model-coupling methodologies were analyzed, revealing numerous coupling approaches, some of which were categorized in the literature. Critiquing these categorizations while comparing them with documented implementations led to the development of a new coupling hierarchy. The properties of each hierarchical level are discussed with the advantages and limitations of each design. To successfully couple models, the spatial and temporal scales of all models in the coupled system and the spatiotemporal extents of the data must be reconciled. Finally, a case study demonstrated methodologies for coupling SMOKE and a GIS. One methodology required a new approach utilizing dynamically linked libraries. Consequently, emissions were processed using SMOKE from a GIS. Also, a new method of converting data from netCDF files into a database
Stationary spiral flow in polytropic stellar models
Pekeris, C.L.
1980-06-01
It is shown that, in addition to the static Emden solution, a self-gravitating polytropic gas has a dynamic option in which there is stationary flow along spiral trajectories wound around the surfaces of concentric tori. The motion is obtained as a solution of a partial differential equation which is satisfied by the meridional stream function, coupled with Poisson's equation and a Bernoulli-type equation for the pressure (density). The pressure is affected by the whole of the Bernoulli term rather than by the centrifugal part only, which acts for a rotating model, and it may be reduced down to zero at the center. The spiral type of flow is illustrated for an incompressible fluid (n = 0), for which an exact solution is obtained. The features of the dynamic constant-density model are discussed as a basis for future comparison with the solution for compressible models.
MOUNTAIN-SCALE COUPLED PROCESSES (TH/THC/THM)MODELS
Y.S. Wu
2005-08-24
This report documents the development and validation of the mountain-scale thermal-hydrologic (TH), thermal-hydrologic-chemical (THC), and thermal-hydrologic-mechanical (THM) models. These models provide technical support for screening of features, events, and processes (FEPs) related to the effects of coupled TH/THC/THM processes on mountain-scale unsaturated zone (UZ) and saturated zone (SZ) flow at Yucca Mountain, Nevada (BSC 2005 [DIRS 174842], Section 2.1.1.1). The purpose and validation criteria for these models are specified in ''Technical Work Plan for: Near-Field Environment and Transport: Coupled Processes (Mountain-Scale TH/THC/THM, Drift-Scale THC Seepage, and Drift-Scale Abstraction) Model Report Integration'' (BSC 2005 [DIRS 174842]). Model results are used to support exclusion of certain FEPs from the total system performance assessment for the license application (TSPA-LA) model on the basis of low consequence, consistent with the requirements of 10 CFR 63.342 [DIRS 173273]. Outputs from this report are not direct feeds to the TSPA-LA. All the FEPs related to the effects of coupled TH/THC/THM processes on mountain-scale UZ and SZ flow are discussed in Sections 6 and 7 of this report. The mountain-scale coupled TH/THC/THM processes models numerically simulate the impact of nuclear waste heat release on the natural hydrogeological system, including a representation of heat-driven processes occurring in the far field. The mountain-scale TH simulations provide predictions for thermally affected liquid saturation, gas- and liquid-phase fluxes, and water and rock temperature (together called the flow fields). The main focus of the TH model is to predict the changes in water flux driven by evaporation/condensation processes, and drainage between drifts. The TH model captures mountain-scale three-dimensional flow effects, including lateral diversion and mountain-scale flow patterns. The mountain-scale THC model evaluates TH effects on water and gas
NASA Astrophysics Data System (ADS)
Kim, Jihoon; Moridis, George J.
2013-10-01
We developed a hydraulic fracturing simulator by coupling a flow simulator to a geomechanics code, namely T+M simulator. Modeling of the vertical fracture development involves continuous updating of the boundary conditions and of the data connectivity, based on the finite element method for geomechanics. The T+M simulator can model the initial fracture development during the hydraulic fracturing operations, after which the domain description changes from single continuum to double or multiple continua in order to rigorously model both flow and geomechanics for fracture-rock matrix systems. The T+H simulator provides two-way coupling between fluid-heat flow and geomechanics, accounting for thermo-poro-mechanics, treats nonlinear permeability and geomechanical moduli explicitly, and dynamically tracks changes in the fracture(s) and in the pore volume. We also fully account for leak-off in all directions during hydraulic fracturing. We first test the T+M simulator, matching numerical solutions with the analytical solutions for poromechanical effects, static fractures, and fracture propagations. Then, from numerical simulation of various cases of the planar fracture propagation, shear failure can limit the vertical fracture propagation of tensile failure, because of leak-off into the reservoirs. Slow injection causes more leak-off, compared with fast injection, when the same amount of fluid is injected. Changes in initial total stress and contributions of shear effective stress to tensile failure can also affect formation of the fractured areas, and the geomechanical responses are still well-posed.
Modeling Combustion in Supersonic Flows
NASA Technical Reports Server (NTRS)
Drummond, J. Philip; Danehy, Paul M.; Bivolaru, Daniel; Gaffney, Richard L.; Tedder, Sarah A.; Cutler, Andrew D.
2007-01-01
This paper discusses the progress of work to model high-speed supersonic reacting flow. The purpose of the work is to improve the state of the art of CFD capabilities for predicting the flow in high-speed propulsion systems, particularly combustor flow-paths. The program has several components including the development of advanced algorithms and models for simulating engine flowpaths as well as a fundamental experimental and diagnostic development effort to support the formulation and validation of the mathematical models. The paper will provide details of current work on experiments that will provide data for the modeling efforts along with with the associated nonintrusive diagnostics used to collect the data from the experimental flowfield. Simulation of a recent experiment to partially validate the accuracy of a combustion code is also described.
Optimizing information flow in small genetic networks. IV. Spatial coupling
NASA Astrophysics Data System (ADS)
Sokolowski, Thomas R.; Tkačik, Gašper
2015-06-01
We typically think of cells as responding to external signals independently by regulating their gene expression levels, yet they often locally exchange information and coordinate. Can such spatial coupling be of benefit for conveying signals subject to gene regulatory noise? Here we extend our information-theoretic framework for gene regulation to spatially extended systems. As an example, we consider a lattice of nuclei responding to a concentration field of a transcriptional regulator (the input) by expressing a single diffusible target gene. When input concentrations are low, diffusive coupling markedly improves information transmission; optimal gene activation functions also systematically change. A qualitatively different regulatory strategy emerges where individual cells respond to the input in a nearly steplike fashion that is subsequently averaged out by strong diffusion. While motivated by early patterning events in the Drosophila embryo, our framework is generically applicable to spatially coupled stochastic gene expression models.