Comment on ''Critical behavior of a traffic flow model''
Chowdhury, D.; Kertesz, J.; Nagel, K.; Santen, L.; Schadschneider, A.
2000-03-01
We show that the dynamical structure factor investigated by Roters et al. [Phys. Rev. E 59, 2672 (1999)] does not allow the determination of the precise nature of the transition in the Nagel-Schreckenberg cellular automata model for traffic flow. We provide evidence for the existence of a crossover instead of a critical point. (c) 2000 The American Physical Society.
Effect of turbulence models on criticality conditions in swirling flows
NASA Technical Reports Server (NTRS)
Gatski, Thomas B.; Spall, Robert E.
1995-01-01
The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.
Critical assessment of Reynolds stress turbulence models using homogeneous flows
NASA Technical Reports Server (NTRS)
Shabbir, Aamir; Shih, Tsan-Hsing
1992-01-01
In modeling the rapid part of the pressure correlation term in the Reynolds stress transport equations, extensive use has been made of its exact properties which were first suggested by Rotta. These, for example, have been employed in obtaining the widely used Launder, Reece and Rodi (LRR) model. Some recent proposals have dropped one of these properties to obtain new models. We demonstrate, by computing some simple homogeneous flows, that doing so does not lead to any significant improvements over the LRR model and it is not the right direction in improving the performance of existing models. The reason for this, in our opinion, is that violation of one of the exact properties can not bring in any new physics into the model. We compute thirteen homogeneous flows using LRR (with a recalibrated rapid term constant), IP and SSG models. The flows computed include the flow through axisymmetric contraction; axisymmetric expansion; distortion by plane strain; and homogeneous shear flows with and without rotation. Results show that for most general representation for a model linear in the anisotropic tensor, performs either better or as good as the other two models of the same level.
NASA Astrophysics Data System (ADS)
Tutar, M.; Holdø, A. E.
2001-04-01
The numerical simulation of transitional flow around a two-dimensional stationary circular cylinder is presented using two groups of turbulence models in a sub-critical flow regime. In the first group, enhanced two-equation turbulence models based on the eddy viscosity concept are used. They include the non-linear k- model with extended models, such as renormalization group (RNG) and the anisotropic model. In the second group, flow simulation is carried out using the large eddy simulation (LES) method, which is based on a standard sub-grid scale (SGS) model with a near-wall approach. This near-wall model, without using the law of wall, is achieved in a finite element code. The numerical results extracted from these simulations are compared with each other and with the experimental data in order to determine the relative performance of these turbulence models and to find the best model for the flow of interest. Although most of the LES simulations have been previously carried out using finite volume methods, results from using the present model show that the finite element method (FEM) can also be used with confidence. Copyright
NASA Astrophysics Data System (ADS)
Fourrate, K.; Loulidi, M.
2006-01-01
We suggest a disordered traffic flow model that captures many features of traffic flow. It is an extension of the Nagel-Schreckenberg (NaSch) stochastic cellular automata for single line vehicular traffic model. It incorporates random acceleration and deceleration terms that may be greater than one unit. Our model leads under its intrinsic dynamics, for high values of braking probability pr, to a constant flow at intermediate densities without introducing any spatial inhomogeneities. For a system of fast drivers pr→0, the model exhibits a density wave behavior that was observed in car following models with optimal velocity. The gap of the disordered model we present exhibits, for high values of pr and random deceleration, at a critical density, a power law distribution which is a hall mark of a self organized criticality phenomena.
A. Alsaed
2004-09-14
The ''Disposal Criticality Analysis Methodology Topical Report'' (YMP 2003) presents the methodology for evaluating potential criticality situations in the monitored geologic repository. As stated in the referenced Topical Report, the detailed methodology for performing the disposal criticality analyses will be documented in model reports. Many of the models developed in support of the Topical Report differ from the definition of models as given in the Office of Civilian Radioactive Waste Management procedure AP-SIII.10Q, ''Models'', in that they are procedural, rather than mathematical. These model reports document the detailed methodology necessary to implement the approach presented in the Disposal Criticality Analysis Methodology Topical Report and provide calculations utilizing the methodology. Thus, the governing procedure for this type of report is AP-3.12Q, ''Design Calculations and Analyses''. The ''Criticality Model'' is of this latter type, providing a process evaluating the criticality potential of in-package and external configurations. The purpose of this analysis is to layout the process for calculating the criticality potential for various in-package and external configurations and to calculate lower-bound tolerance limit (LBTL) values and determine range of applicability (ROA) parameters. The LBTL calculations and the ROA determinations are performed using selected benchmark experiments that are applicable to various waste forms and various in-package and external configurations. The waste forms considered in this calculation are pressurized water reactor (PWR), boiling water reactor (BWR), Fast Flux Test Facility (FFTF), Training Research Isotope General Atomic (TRIGA), Enrico Fermi, Shippingport pressurized water reactor, Shippingport light water breeder reactor (LWBR), N-Reactor, Melt and Dilute, and Fort Saint Vrain Reactor spent nuclear fuel (SNF). The scope of this analysis is to document the criticality computational method. The criticality
NASA Technical Reports Server (NTRS)
Bertsimas, Dimitris; Odoni, Amedeo
1997-01-01
This document presents a critical review of the principal existing optimization models that have been applied to Air Traffic Flow Management (TFM). Emphasis will be placed on two problems, the Generalized Tactical Flow Management Problem (GTFMP) and the Ground Holding Problem (GHP), as well as on some of their variations. To perform this task, we have carried out an extensive literature review that has covered more than 40 references, most of them very recent. Based on the review of this emerging field our objectives were to: (i) identify the best available models; (ii) describe typical contexts for applications of the models; (iii) provide illustrative model formulations; and (iv) identify the methodologies that can be used to solve the models. We shall begin our presentation below by providing a brief context for the models that we are reviewing. In Section 3 we shall offer a taxonomy and identify four classes of models for review. In Sections 4, 5, and 6 we shall then review, respectively, models for the Single-Airport Ground Holding Problem, the Generalized Tactical FM P and the Multi-Airport Ground Holding Problem (for the definition of these problems see Section 3 below). In each section, we identify the best available models and discuss briefly their computational performance and applications, if any, to date. Section 7 summarizes our conclusions about the state of the art.
A critical evaluation of various turbulence models as applied to internal fluid flows
NASA Technical Reports Server (NTRS)
Nallasamy, M.
1985-01-01
Models employed in the computation of turbulent flows are described and their application to internal flows is evaluated by examining the predictions of various turbulence models in selected flow configurations. The main conclusions are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly; the algebraic stress model performs better in this case; and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well for fully developed flow, for developing flow, the algebraic stress model performance is not good; a Reynolds stress model should be used. False diffusion and inlet boundary conditions are discussed. Countergradient transport and its implications in turbulence modeling is mentioned. Two examples of recirculating flow predictions obtained using PHOENICS code are discussed. The vortex method, large eddy simulation (modeling of subgrid scale Reynolds stresses), and direct simulation, are considered. Some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.
Modeling nutrient dynamics under critical flow conditions in three tributaries of St. Louis Bay.
Liu, Zhijun; Kingery, William L; Huddleston, David H; Hossain, Faisal; Chen, Wei; Hashim, Noor B; Kieffer, Janna M
2008-05-01
Previous research results indicated that dry weather condition has complicated impacts on nitrogen dynamics; monitored and modeling data showed both increased and decreased levels. In order to facilitate the total maximum daily loads (TMDLs) development at three tributaries of St. Louis Bay estuary, the nitrogen dynamics were investigated for two designed critical flow conditions by integrating Hydrological Simulation Program Fortran (HSPF), Environmental Fluid Dynamics Code (EFDC), and Water Quality Analysis Simulation Program (WASP). The total amount of precipitation during the dry year corresponded to a flow condition with return period of 50 years, and 10-year return period for wet year. The dry year contributed more total nitrogen (TN) loads per unit flow volume. At the upstream tributaries, the computed peak reach-averaged TN concentrations were significantly higher for dry weather simulation than wet conditions, whereas at the near-bay tributary, there were no significant differences in the peak TN concentrations. Hence, for the upstream tributaries, the nitrogen TMDL calculation should be based on dry weather condition since the decision-makers are more concerned about the worse scenario. PMID:18393072
Uncertainties in modelling and scaling of critical flows and pump model in TRAC-PF1/MOD1
Rohatgi, U.S.; Yu, Wen-Shi
1987-01-01
The USNRC has established a Code Scalability, Applicability and Uncertainty (CSAU) evaluation methodology to quantify the uncertainty in the prediction of safety parameters by the best estimate codes. These codes can then be applied to evaluate the Emergency Core Cooling System (ECCS). The TRAC-PF1/MOD1 version was selected as the first code to undergo the CSAU analysis for LBLOCA applications. It was established through this methodology that break flow and pump models are among the top ranked models in the code affecting the peak clad temperature (PCT) prediction for LBLOCA. The break flow model bias or discrepancy and the uncertainty were determined by modelling the test section near the break for 12 Marviken tests. It was observed that the TRAC-PF1/MOD1 code consistently underpredicts the break flow rate and that the prediction improved with increasing pipe length (larger L/D). This is true for both subcooled and two-phase critical flows. A pump model was developed from Westinghouse (1/3 scale) data. The data represent the largest available test pump relevant to Westinghouse PWRs. It was then shown through the analysis of CE and CREARE pump data that larger pumps degrade less and also that pumps degrade less at higher pressures. Since the model developed here is based on the 1/3 scale pump and on low pressure data, it is conservative and will overpredict the degradation when applied to PWRs.
Managing critical materials with a technology-specific stocks and flows model.
Busch, Jonathan; Steinberger, Julia K; Dawson, David A; Purnell, Phil; Roelich, Katy
2014-01-21
The transition to low carbon infrastructure systems required to meet climate change mitigation targets will involve an unprecedented roll-out of technologies reliant upon materials not previously widespread in infrastructure. Many of these materials (including lithium and rare earth metals) are at risk of supply disruption. To ensure the future sustainability and resilience of infrastructure, circular economy policies must be crafted to manage these critical materials effectively. These policies can only be effective if supported by an understanding of the material demands of infrastructure transition and what reuse and recycling options are possible given the future availability of end-of-life stocks. This Article presents a novel, enhanced stocks and flows model for the dynamic assessment of material demands resulting from infrastructure transitions. By including a hierarchical, nested description of infrastructure technologies, their components, and the materials they contain, this model can be used to quantify the effectiveness of recovery at both a technology remanufacturing and reuse level and a material recycling level. The model's potential is demonstrated on a case study on the roll-out of electric vehicles in the UK forecast by UK Department of Energy and Climate Change scenarios. The results suggest policy action should be taken to ensure Li-ion battery recycling infrastructure is in place by 2025 and NdFeB motor magnets should be designed for reuse. This could result in a reduction in primary demand for lithium of 40% and neodymium of 70%. PMID:24328245
Managing Critical Materials with a Technology-Specific Stocks and Flows Model
2013-01-01
The transition to low carbon infrastructure systems required to meet climate change mitigation targets will involve an unprecedented roll-out of technologies reliant upon materials not previously widespread in infrastructure. Many of these materials (including lithium and rare earth metals) are at risk of supply disruption. To ensure the future sustainability and resilience of infrastructure, circular economy policies must be crafted to manage these critical materials effectively. These policies can only be effective if supported by an understanding of the material demands of infrastructure transition and what reuse and recycling options are possible given the future availability of end-of-life stocks. This Article presents a novel, enhanced stocks and flows model for the dynamic assessment of material demands resulting from infrastructure transitions. By including a hierarchical, nested description of infrastructure technologies, their components, and the materials they contain, this model can be used to quantify the effectiveness of recovery at both a technology remanufacturing and reuse level and a material recycling level. The model’s potential is demonstrated on a case study on the roll-out of electric vehicles in the UK forecast by UK Department of Energy and Climate Change scenarios. The results suggest policy action should be taken to ensure Li-ion battery recycling infrastructure is in place by 2025 and NdFeB motor magnets should be designed for reuse. This could result in a reduction in primary demand for lithium of 40% and neodymium of 70%. PMID:24328245
Numerical modeling of continuous flow microwave heating: a critical comparison of COMSOL and ANSYS.
Salvi, D; Boldor, Dorin; Ortego, J; Aita, G M; Sabliov, C M
2010-01-01
Numerical models were developed to simulate temperature profiles in Newtonian fluids during continuous flow microwave heating by one way coupling electromagnetism, fluid flow, and heat transport in ANSYS 8.0 and COMSOL Multiphysics v3.4. Comparison of the results from the COMSOL model with the results from a pre-developed and validated ANSYS model ensured accuracy of the COMSOL model. Prediction of power Loss by both models was in close agreement (5-13% variation) and the predicted temperature profiles were similar. COMSOL provided a flexible model setup whereas ANSYS required coupling incompatible elements to transfer load between electromagnetic, fluid flow, and heat transport modules. Overall, both software packages provided the ability to solve multiphysics phenomena accurately. PMID:21721467
Critical evaluation of three hemodynamic models for the numerical simulation of intra-stent flows.
Chabi, Fatiha; Champmartin, Stéphane; Sarraf, Christophe; Noguera, Ricardo
2015-07-16
We evaluate here three hemodynamic models used for the numerical simulation of bare and stented artery flows. We focus on two flow features responsible for intra-stent restenosis: the wall shear stress and the re-circulation lengths around a stent. The studied models are the Poiseuille profile, the simplified pulsatile profile and the complete pulsatile profile based on the analysis of Womersley. The flow rate of blood in a human left coronary artery is considered to compute the velocity profiles. "Ansys Fluent 14.5" is used to solve the Navier-Stokes and continuity equations. As expected our results show that the Poiseuille profile is questionable to simulate the complex flow dynamics involved in intra-stent restenosis. Both pulsatile models give similar results close to the strut but diverge far from it. However, the computational time for the complete pulsatile model is five times that of the simplified pulsatile model. Considering the additional "cost" for the complete model, we recommend using the simplified pulsatile model for future intra-stent flow simulations. PMID:26044195
A critical assessment of viscous models of trench topography and corner flow
NASA Technical Reports Server (NTRS)
Zhang, J.; Hager, B. H.; Raefsky, A.
1984-01-01
Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.
Turbulent flow past a backward-facing step - A critical evaluation of two-equation models
NASA Technical Reports Server (NTRS)
Thangam, S.; Speziale, C. G.
1992-01-01
The ability of two-equation models to accurately predict separated flows is analyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward facing step is chosen as a test case in an effort to resolve the variety of conflicting results that were published during the past decade concerning the performance of two-equation models. It is found that the errors in the reported predictions of the k-epsilon model have two major origins: (1) numerical problems arising from inadequate resolution, and (2) inaccurate predictions for normal Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in near wall modelling play a substantially smaller role. Detailed calculations are presented which strongly indicate the standard k-epsilon model - when modified with an independently calibrated anisotropic eddy viscosity - can yield surprisingly good predictions for the backstep problem.
Knoeri, Christof; Wäger, Patrick A; Stamp, Anna; Althaus, Hans-Joerg; Weil, Marcel
2013-09-01
Emerging technologies such as information and communication-, photovoltaic- or battery technologies are expected to increase significantly the demand for scarce metals in the near future. The recently developed methods to evaluate the criticality of mineral raw materials typically provide a 'snapshot' of the criticality of a certain material at one point in time by using static indicators both for supply risk and for the impacts of supply restrictions. While allowing for insights into the mechanisms behind the criticality of raw materials, these methods cannot account for dynamic changes in products and/or activities over time. In this paper we propose a conceptual framework intended to overcome these limitations by including the dynamic interactions between different possible demand and supply configurations. The framework integrates an agent-based behaviour model, where demand emerges from individual agent decisions and interaction, into a dynamic material flow model, representing the materials' stocks and flows. Within the framework, the environmental implications of substitution decisions are evaluated by applying life-cycle assessment methodology. The approach makes a first step towards a dynamic criticality assessment and will enhance the understanding of industrial substitution decisions and environmental implications related to critical metals. We discuss the potential and limitation of such an approach in contrast to state-of-the-art methods and how it might lead to criticality assessments tailored to the specific circumstances of single industrial sectors or individual companies. PMID:23453658
Lee, S.Y. ); Schrock, V.E. )
1990-11-01
Analysis of loss of coolant accident (LOCA) scenarios in nuclear reactor safety evaluation depends on knowledge of many complex phenomena. A primary phenomenon controlling the sequence of events, by determining the residual coolant mass inventory within the primary system, is the critical flow process. Critical flow of a flashing liquid is complicated by marked departure from thermal equilbrium. Several complex models have been proposed to represent the non-equilibrium effects, including six-equation two-fluid models. In the present paper a new cavity flooding model is used for the evaluation of pressure undershoot at flashing inception. This model is similar to the one developed by Fabic (1964) for the evaluation of liquid superheat required for boiling on a surface subjected to transient heating. The model contains an experimentally deduced factor, which is correlated against stagnation subcooling using the experimental data of Amos and Schrock (1983, 1984), Jeandey et al. (1981), and the Marviken tests (Anon., 1979). The model was then tested against seven additional data sets and shown to be very accurate in predicted mass flux (standard deviation of 10.9% for all data). The cavity flooding model is thought to represent the true physics more correctly than does the earlier model, which had its origin in molecular fluctuation theory.
NASA Astrophysics Data System (ADS)
Dănilă, B.; Harko, T.; Mocanu, G.
2015-11-01
We investigate the transition to self-organized criticality in a two-dimensional model of a flux tube with a background flow. The magnetic induction equation, represented by a partial differential equation with a stochastic source term, is discretized and implemented on a two-dimensional cellular automaton. The energy released by the automaton during one relaxation event is the magnetic energy. As a result of the simulations, we obtain the time evolution of the energy release, of the system control parameter, of the event lifetime distribution and of the event size distribution, respectively, and we establish that a self-organized critical state is indeed reached by the system. Moreover, energetic initial impulses in the magnetohydrodynamic flow can lead to one-dimensional signatures in the magnetic two-dimensional system, once the self-organized critical regime is established. The applications of the model for the study of gamma-ray bursts (GRBs) is briefly considered, and it is shown that some astrophysical parameters of the bursts, like the light curves, the maximum released energy and the number of peaks in the light curve can be reproduced and explained, at least on a qualitative level, by working in a framework in which the systems settles in a self-organized critical state via magnetic reconnection processes in the magnetized GRB fireball.
Kanellopoulos, Giorgos; van der Weele, Ko
2012-06-01
We study the transport of granular matter through a staircaselike array of K vertically vibrated compartments. Given a constant inflow rate Qflow establishes itself along the entire length of the system. However, as soon as Q grows beyond the critical value Q{cr}(K) the particles form a cluster and the flow comes to a halt. Interestingly, this clustering is preceded by a subcritical warning signal: for Q values just below Q{cr}(K) the density profile along the conveyor belt spontaneously develops a pattern in which the compartments are alternatingly densely and sparsely populated. In a previous paper [Kanellopoulos and van der Weele, Int. J. Bifurcation Chaos 21, 2305 (2011)] this pattern was shown to be the result of a period-doubling bifurcation. The present paper aims at unravelling the physical mechanism that lies at the basis of the pattern formation. To this end we study the continuum version of the same system, replacing the compartment number k=1,...,K by a continuous variable x. The dynamics of the system is now described (instead of by K coupled ordinary differential equations) by a single partial differential equation of the Fokker-Planck type, with a drift and a diffusive term that both depend on the density. The drift term turns out to be responsible for the subcritical density oscillations, thereby paving the way for the eventual clustering which sets in when the diffusion coefficient becomes negative. The observed sequence of events in the granular transport system is thus explained as an interplay between drift and (anti)diffusion. PMID:23005082
Caballero-Guzman, Alejandro; Nowack, Bernd
2016-06-01
Material flow analysis (MFA) is a useful tool to predict the flows of engineered nanomaterials (ENM) to the environment. The quantification of release factors is a crucial part of MFA modeling. In the last years an increasing amount of literature on release of ENM from materials and products has been published. The purpose of this review is to analyze the strategies implemented by MFA models to include these release data, in particular to derive transfer coefficients (TC). Our scope was focused on those articles that analyzed the release from applications readily available in the market in settings that resemble average use conditions. Current MFA studies rely to a large extent on extrapolations, authors' assumptions, expert opinions and other informal sources of data to parameterize the models. We were able to qualitatively assess the following aspects of the release literature: (i) the initial characterization of ENM provided, (ii) quantitative information on the mass of ENM released and its characterization, (iii) description of transformation reactions and (iv) assessment of the factors determining release. Although the literature on ENM release is growing, coverage of exposure scenarios is still limited; only 20% of the ENMs used industrially and 36% of the product categories involved have been investigated in release studies and only few relevant release scenarios have been described. Furthermore, the information provided is rather incomplete concerning descriptions and characterizations of ENMs and the released materials. Our results show that both the development of methods to define the TCs and of protocols to enhance assessment of ENM release from nano-applications will contribute to increase the exploitability of the data provided for MFA models. The suggestions we provide in this article will likely contribute to an improved exposure modeling by providing ENM release estimates closer to reality. PMID:26970875
Natural gas flow through critical nozzles
NASA Technical Reports Server (NTRS)
Johnson, R. C.
1969-01-01
Empirical method for calculating both the mass flow rate and upstream volume flow rate through critical flow nozzles is determined. Method requires knowledge of the composition of natural gas, and of the upstream pressure and temperature.
Fraser, D.W.H.; Abdelmessih, A.H.
1995-09-01
A general unified model is developed to predict one-component critical two-phase pipe flow. Modelling of the two-phase flow is accomplished by describing the evolution of the flow between the location of flashing inception and the exit (critical) plane. The model approximates the nonequilibrium phase change process via thermodynamic equilibrium paths. Included are the relative effects of varying the location of flashing inception, pipe geometry, fluid properties and length to diameter ratio. The model predicts that a range of critical mass fluxes exist and is bound by a maximum and minimum value for a given thermodynamic state. This range is more pronounced at lower subcooled stagnation states and can be attributed to the variation in the location of flashing inception. The model is based on the results of an experimental study of the critical two-phase flow of saturated and subcooled water through long tubes. In that study, the location of flashing inception was accurately controlled and adjusted through the use of a new device. The data obtained revealed that for fixed stagnation conditions, the maximum critical mass flux occurred with flashing inception located near the pipe exit; while minimum critical mass fluxes occurred with the flashing front located further upstream. Available data since 1970 for both short and long tubes over a wide range of conditions are compared with the model predictions. This includes test section L/D ratios from 25 to 300 and covers a temperature and pressure range of 110 to 280{degrees}C and 0.16 to 6.9 MPa. respectively. The predicted maximum and minimum critical mass fluxes show an excellent agreement with the range observed in the experimental data.
NASA Technical Reports Server (NTRS)
Thangam, S.; Speziale, C. G.
1991-01-01
The ability of two-equation models to accurately predict separated flows is analyzed from a combined theoretical and computational standpoint. Turbulent flow past a backward facing step is chosen as a test case in an effort to resolve the variety of conflicting results that were published during the past decade concerning the performance of two-equation models. It is found that the errors in the reported predictions of the k-epsilon model have two major origins: (1) numerical problems arising from inadequate resolution, and (2) inaccurate predictions for normal Reynolds stress differences arising from the use of an isotropic eddy viscosity. Inadequacies in near wall modelling play a substantially smaller role. Detailed calculations are presented which strongly indicate the standard k-epsilon model - when modified with an independently calibrated anisotropic eddy viscosity - can yield surprisingly good predictions for the backstep problem.
NASA Technical Reports Server (NTRS)
Steffen, C. J., Jr.
1993-01-01
Turbulent backward-facing step flow was examined using four low turbulent Reynolds number k-epsilon models and one standard high Reynolds number technique. A tunnel configuration of 1:9 (step height: exit tunnel height) was used. The models tested include: the original Jones and Launder; Chien; Launder and Sharma; and the recent Shih and Lumley formulation. The experimental reference of Driver and Seegmiller was used to make detailed comparisons between reattachment length, velocity, pressure, turbulent kinetic energy, Reynolds shear stress, and skin friction predictions. The results indicated that the use of a wall function for the standard k-epsilon technique did not reduce the calculation accuracy for this separated flow when compared to the low turbulent Reynolds number techniques.
Critical heat flux in subcooled flow boiling
NASA Astrophysics Data System (ADS)
Hall, David Douglas
The critical heat flux (CHF) phenomenon was investigated for water flow in tubes with particular emphasis on the development of methods for predicting CHF in the subcooled flow boiling regime. The Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) CHF database for water flow in a uniformly heated tube was compiled from the world literature dating back to 1949 and represents the largest CHF database ever assembled with 32,544 data points from over 100 sources. The superiority of this database was proven via a detailed examination of previous databases. The PU-BTPFL CHF database is an invaluable tool for the development of CHF correlations and mechanistic models that are superior to existing ones developed with smaller, less comprehensive CHF databases. In response to the many inaccurate and inordinately complex correlations, two nondimensional, subcooled CHF correlations were formulated, containing only five adjustable constants and whose unique functional forms were determined without using a statistical analysis but rather using the parametric trends observed in less than 10% of the subcooled CHF data. The correlation based on inlet conditions (diameter, heated length, mass velocity, pressure, inlet quality) was by far the most accurate of all known subcooled CHF correlations, having mean absolute and root-mean-square (RMS) errors of 10.3% and 14.3%, respectively. The outlet (local) conditions correlation was the most accurate correlation based on local CHF conditions (diameter, mass velocity, pressure, outlet quality) and may be used with a nonuniform axial heat flux. Both correlations proved more accurate than a recent CHF look-up table commonly employed in nuclear reactor thermal hydraulic computer codes. An interfacial lift-off, subcooled CHF model was developed from a consideration of the instability of the vapor-liquid interface and the fraction of heat required for liquid-vapor conversion as opposed to that for bulk liquid heating. Severe
Critical velocity of superfluid helium flow in narrow pore filters
NASA Astrophysics Data System (ADS)
Hofmann, A.
1990-05-01
The flow rates of superfluid helium passing through a spongelike plug made from a stack of high-porosity cellulose-nitrate membrane filters were measured, comparing three grades of filters: 10, 50, and 200 nm. The critical flow rate was analyzed from the onset of nonzero difference in the chemical potential across the plug. The flow, driven by activating a heater, was measured with an acoustic flowmeter. It is shown that the critical flow rate increases with decreasing filtration grade and that the upper limit of the superfluid flow velocity can be predicted from the quoted filtration grade of the filter and from a shape factor derived from elementary geometric conditions of densely packed spheres. This model yields correct temperature dependence of the critical velocity.
Critical Infrastructure Modeling System
Energy Science and Technology Software Center (ESTSC)
2004-10-01
The Critical Infrastructure Modeling System (CIMS) is a 3D modeling and simulation environment designed to assist users in the analysis of dependencies within individual infrastructure and also interdependencies between multiple infrastructures. Through visual cuing and textual displays, a use can evaluate the effect of system perturbation and identify the emergent patterns that evolve. These patterns include possible outage areas from a loss of power, denial of service or access, and disruption of operations. Method ofmore » Solution: CIMS allows the user to model a system, create an overlay of information, and create 3D representative images to illustrate key infrastructure elements. A geo-referenced scene, satellite, aerial images or technical drawings can be incorporated into the scene. Scenarios of events can be scripted, and the user can also interact during run time to alter system characteristics. CIMS operates as a discrete event simulation engine feeding a 3D visualization.« less
Critical Infrastructure Modeling System
2004-10-01
The Critical Infrastructure Modeling System (CIMS) is a 3D modeling and simulation environment designed to assist users in the analysis of dependencies within individual infrastructure and also interdependencies between multiple infrastructures. Through visual cuing and textual displays, a use can evaluate the effect of system perturbation and identify the emergent patterns that evolve. These patterns include possible outage areas from a loss of power, denial of service or access, and disruption of operations. Method of Solution: CIMS allows the user to model a system, create an overlay of information, and create 3D representative images to illustrate key infrastructure elements. A geo-referenced scene, satellite, aerial images or technical drawings can be incorporated into the scene. Scenarios of events can be scripted, and the user can also interact during run time to alter system characteristics. CIMS operates as a discrete event simulation engine feeding a 3D visualization.
Roaming form factors for the tricritical to critical Ising flow
NASA Astrophysics Data System (ADS)
Horváth, D. X.; Dorey, P. E.; Takács, G.
2016-07-01
We study the massless flows described by the staircase model introduced by Al.B. Zamolodchikov through the analytic continuation of the sinh-Gordon S-matrix, focusing on the renormalisation group flow from the tricritical to the critical Ising model. We show that the properly defined roaming limits of certain sinh-Gordon form factors are identical to the form factors of the order and disorder operators for the massless flow. As a by-product, we also construct form factors for a semi-local field in the sinh-Gordon model, which can be associated with the twist field in the ultraviolet limiting free massless bosonic theory.
Critical Velocity in Open Capillary Channel Flows
NASA Technical Reports Server (NTRS)
Rosendahl, Uwe; Dreyer, Michael E.; Rath, Hans J.; Motil, Brian; Singh, Bhim S. (Technical Monitor)
2001-01-01
We investigate forced liquid flows through open capillary channels with free surfaces experimentally. The experiments were performed under low gravity conditions in the Bremen Drop Tower and on board the sounding rocket TEXUS-37. Open capillary channels (vanes) are used in surface tension tanks to transport the propellant and to provide a flow path for the bubble-free liquid supply to the thrusters. Since the free surfaces can only withstand a certain pressure differential between the liquid and ambient, the flow rate in the channel is limited. The maximum flow rate is achieved when the surfaces collapse and gas is ingested into the outlet. Since experimental and theoretical data of this flow rate limitation is lacking, the safety factors for the application of vanes in surface tension tanks must be unnecessary high. The aim of the investigation is to determine the maximum liquid flow rate and the corresponding critical flow velocity. The characteristic nondimensional parameters, OHNESORGE number, and gap ratio, cover a wide range of usual vanes. For the theoretical approach a one-dimensional momentum balance was set up. The numerical solution yields the maximum volume flux and the position of the free surface in good agreement with the experiments.
Critical Velocities in Open Capillary Flow
NASA Technical Reports Server (NTRS)
Dreyer, Michael; Langbein, Dieter; Rath, Hans J.
1996-01-01
This paper describes the proposed research program on open capillary flow and the preliminary work performed theoretically and in drop tower experiments. The work focuses on the fundamental physical understanding of the flow through capillary bound geometries, where the circumference of the cross section of the flow path contains free surfaces. Examples for such a flow configuration are capillary vanes in surface tension tanks, flow along edges and corners and flow through liquid bridges. The geometries may be classified by their cross section areas, wetted circumferences and the radii of curvature of the free surfaces. In the streaming float zone the flow path is bound by a free surface only. The ribbon vane is a model for vane types used in surface tension tanks, where a structure in proximity to the tank wall forms a capillary gap. A groove is used in heat pipes for the transportation of the condensed working fluid to the heat source and a wedge may occur in a spaceborne experiment where fluid has to be transported by the means of surface tension. The research objectives are the determination of the maximum volume flux, the observation of the free surfaces and the liquid flow inside the flow path as well as the evaluation of the limiting capillary wave speed. The restriction of the maximum volume flux is due to convective forces (flow velocity exceeding the capillary wave speed) and/or viscous forces, i.e. the viscous head loss along the flow path must be compensated by the capillary pressure due to the curved free surface. Exceeding the maximum volume flux leads to the choking of the flow path, thus the free surface collapses and.gas ingestion occurs at the outlet. The means are ground-based experimental work with plateau tanks and in a drop tower, a sounding rocket flight, and theoretical analysis with integral balances as well as full three dimensional CFD solutions for flow with free surfaces.
Characterization of non equilibrium effects on high quality critical flows
Camelo, E.; Lemonnier, H.; Ochterbeck, J.
1995-09-01
The appropriate design of various pieces of safety equipment such as relief systems, relies on the accurate description of critical flow phenomena. Most of the systems of industrial interest are willing to be described by one-dimensional area-averaged models and a large fraction of them involves multi-component high gas quality flows. Within these circumstances, the flow is very likely to be of an annular dispersed nature and its description by two-fluid models requires various closure relations. Among the most sensitive closures, there is the interfacial area and the liquid entrained fraction. The critical flowrate depends tremendously on the accurate description of the non equilibrium which results from the correctness of the closure equations. In this study, two-component flows are emphasized and non equilibrium results mainly form the differences in the phase velocities. It is therefore of the utmost importance to have reliable data to characterize non equilibrium phenomena and to assess the validity of the closure models. A comprehensive description of air-water nozzle flows, with emphasis on the effect of the nozzle geometry, has been undertaken and some of the results are presented here which helps understanding the overall flow dynamics. Besides the critical flowrate, the presented material includes pressure profiles, droplet size and velocity, liquid film flowrate and liquid film thickness.
Numerical analysis of critical two-phase flow in a convergent-divergent nozzle
Romstedt, P.; Werner, W.
1986-01-01
The numerical calculation of critical two-phase flow in a convergent-divergent nozzle is complicated by a singularity of the fluid flow equations at the unknown critical point. A method of calculating critical state and its location without any additional assumptions is described. The critical state is identified by its mathematical properties: characteristics and solvability of linear systems with a singular matrix. Because the numerically estimable mathematical properties are the only necessary conditions for the existence of critical flow, some physical ''compatibility criteria'' (flow velocity equals model-consistent two-phase sonic velocity; critical flow is independent of downstream flow state variations) are used as substitutes for mathematically sufficient conditions. Numerical results are shown for the critical flow through LOBI nozzles and for the Super Moby Dick experiment. The two-phase flow is described by a model with equal phase velocities and thermodynamic nonequilibrium.
Gardel, Ls; Afonso, M; Frias, C; Gomes, Me; Reis, Rl
2014-01-01
This work evaluated in vivo performance of a tissue-engineered bone-like matrix obtained by culturing cell-scaffold constructs in a flow perfusion bioreactor, designed to enable culture of large constructs, envisioning the regeneration of critical-sized defects. A blend of starch with polycaprolactone scaffolds was seeded with goat bone marrow stromal cells (GBMSCs) cultured in the perfusion bioreactor for 14 days using osteogenic medium. Cell seeded scaffolds cultured in static conditions acted as controls. After 14 days, constructs (42 mm length and 16 mm in diameter) were implanted in critical size defects performed in the tibial bone of six adult goats from which the bone marrow had been collected previously. Explants were retrieved after six and 12 weeks of implantation and characterized using scanning electron microscopy, energy-dispersive spectroscopy, micro-computed tomography and radiographic analysis to assess tissue morphology and calcification. Explants were histologically analyzed, using Hematoxylin & Eosin and Masson Trichrome staining. Results provided relevant information about the performance and functionality of starch with polycaprolactone-goat bone marrow stromal cell constructs in a critical size orthotopic defect performed in a large animal model and demonstrated that culture of the starch with polycaprolactone scaffolds with the autologous cells in perfusion culture provide a good therapy for the healing and regenerative process of bone defects. PMID:24413026
Experimental investigation of critical flow of supercritical carbon dioxide
NASA Astrophysics Data System (ADS)
Mignot, Guillaume Paul H.
A blowdown facility (0.125 m3) has been built to perform measurements of the critical flow rate of carbon dioxide over a wide range of conditions up to a supercritical pressure of 240 bars and up to a supercritical temperature of 260°C, i.e. three times the critical pressure and two times the critical temperature. The influence of the rupture geometry was investigated using a set of exit pipes with varying entrance shape, roughness and length to diameter ratio ranging from 3.7 to 168. The study showed that a rough sharp edge entrance tube had a lower critical mass flow rate compared to a smooth round entrance tube. For length to diameter ratios larger than 14.7, although two-phase effects were observed, the fluid behavior could be accurately modeled using a homogeneous equilibrium model with friction. For length to diameter ratio smaller than 14.7, the critical mass flux results exhibited a plateau, indicating that the critical mass flow rate was governed by the vena contracta. Stagnation pressure, stagnation temperature and mass time traces were scaled successfully using the initial mass and the initial mass flow rate. An exception was observed for the high density low temperature case due to non equilibrium effects occurring within the vessel. The compressibility of the flow in association with the contraction induced multidimensional and repetitive shock structures within the tube. These have been predicted with computational fluid dynamics modeling for perfect gas conditions. To measure experimentally the fluid state within the tube, an optical absorption technique has been developed, calibrated and tested in two geometries and during an integral blowdown test. Results showed that this new technique lead to the correct qualitative trends in the pressure measurements but that it needed to be calibrated against a more accurate high pressure database obtained for carbon dioxide.
Critical Dimensionless Shields Values for Bankfull Flow
NASA Astrophysics Data System (ADS)
Bunte, K.; Abt, S. R.; Swingle, K. W.
2009-12-01
The critical dimensionless shear stress τ*c (depicted in the Shields curve as a function of the Reynolds particle number Rep) quantifies the slope - flow depth product at which particles from a relatively well-sorted bed with a mean particle size Dm visually start to move. Contrary to many applications of the Shields curve, τ*c was not designed to predict the bed particle size becoming mobile at bankfull flow. This study sketches a bankfull and a critical bankfull Shields-type curve needed to predict the bankfull mobile particle size. Most studies on bedmaterial entrainment are performed in streams that are wadeable at low flow, but exceed wadeability around bankfull flow. Within the bounds of these stream dimensions, the study drafts a relationship of τ*bf versus Rep (bankfull Shields curve) for a sequence of stream types (in sensu Montgomery and Buffington 1997) ranging from steep cobble headwater streams to plane-bed and pool-riffle gravel-beds to low gradient valley streams with sand and silt beds. Probable values of stream gradient, bed D50 size, and a (roughness corrected) hydraulic radius can be assigned to each stream type. The resulting bankfull curve takes values of τ*bf near 10 for silt and 1 for sand-bedded streams, drops to around 0.05 - 0.02 for mobile gravel-bed streams and then increases towards 0.1 and 0.2 for the steepest streams. Transforming τ*bf into the critical bankfull curve τ*cbf from which to predict the bankfull entrainable particle size requires information on the bedload particle size that becomes mobile at bankfull flow. To estimate τ*cbf over a variety of coarse gravel- and cobble-bed mountain streams, the authors used flow competence curves measured with bedload traps at 10 sites. Bedload traps have a 0.3 m by 0.2 m opening, and 1-1.6 m long trailing net with a 4 mm mesh; mounting traps onto ground plates anchored to the stream bottom permits 1-hr sampling times; 4 to 6 traps are typically installed across the stream width
Mathematical modeling of near-critical convection
Cox, B.L.; Pruess, K.; McKibbin, R.
1988-01-01
Fluid and heat flow at temperatures approaching or exceeding that at the critical point (374ºC for pure water, higher for saline fluids) may be encountered in deep zones of geothermal systems and above cooling intrusives. Laboratory experiments have demonstrated strong enhancements in heat transfer at near-critical conditions (Dunn and Hardee, 1981). We have developed special numerical techniques for modeling porous flow at near-critical conditions, which can handle the extreme non-linearities in water properties near the critical point. Our numerical experiments show strong enhancements of convective heat transfer at near-critical conditions; however, the heat transfer rates obtained in the numerical simulations are considerably smaller than those seen in the laboratory experiments by Dunn and Hardee. We discuss possible reasons for this discrepancy and develop suggestions for additional laboratory experiments.
Mathematical modeling of near-critical convection
Cox, B.L.; Pruess, K.; McKibbin, R.
1988-01-01
Fluid and heat flow ar temperatures approaching or exceeding that at the critical point (374)degree)C for pure water, higher for saline fluids) may be encountered in deep zones of geotehrmal systems and above cooling intrusives. Laboratory experiments have demonstrated strong enhancements in heat transfer at near-critial conditions. We have developed special numerical techniques for modeling porous flow at near-critical conditions, which can handle the extreme non-linearities in water properties near the critical point. Our numerical exeperiments show strong enhancements of coventive heat transfer at near-critical conditions;however, the heat transfer rates obtained in the numerical simulations are considerably smaller than those seen in the laboratory experiments by Dunn and Hardee. We discuss possible reasons for this discrepancy and develop suggestions for additional laboratory experiments. 11 refs., 7 figs., 2 tabs.
Numerical Simulation of Two-Phase Critical Flow with the Phase Change in the Nozzle Tube
NASA Astrophysics Data System (ADS)
Ishigaki, Masahiro; Watanabe, Tadashi; Nakamura, Hideo
Two-phase critical flow in the nozzle tube is analyzed numerically by the best estimate code TRACE and the CFD code FLUENT, and the performance of the mass flow rate estimation by the numerical codes is discussed. For the best estimate analysis by the TRACE code, the critical flow option is turned on. The mixture model is used with the cavitation model and the evaporation-condensation model for the numerical simulation by the FLUENT code. Two test cases of the two-phase critical flow are analyzed. One case is the critical flashing flow in a convergent-divergent nozzle (Super Moby Dick experiment), and the other case is the break nozzle flow for a steam generator tube rupture experiment of pressurized water reactors at Large Scale Test Facility of Japan Atomic Energy Agency. The calculation results of the mass flow rates by the numerical simulations show good agreements with the experimental results.
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.
Choked flow of fluid nitrogen with emphasis on the thermodynamic critical region
NASA Technical Reports Server (NTRS)
Hendricks, R. C.; Simoneau, R. J.; Ehlers, R. C.
1972-01-01
Experimental measurements of critical flow rate and pressure ratio for nitrogen flowing through a nozzle are presented. Data for selected stagnation isotherms from 87.5 to 234 K with pressures to 9.3 MN/m2 are compared to an equilibrium model with real fluid properties and also a nonequilibrium model. Critical flow pressure ratio along an isotherm tends to peak while the flow rate indicates an inflection. The point is closely associated with the transposed critical temperature and represents a change in the fluid structure.
NASA Astrophysics Data System (ADS)
Eslamian, Morteza; Saghir, M. Ziad
2009-06-01
In this paper thermodiffusion models developed to estimate the thermal diffusion factor in nonideal liquid mixtures are reviewed; the merits and shortcomings of each model are discussed in detail. Most of these models are multicomponent in principle; however our focus here is on binary mixtures. Two rather different groups of models are identified: models needing a matching parameter to be obtained usually from the outside of thermodynamics, and the self-contained or independent models. Derivation of the matching parameter models using linear non-equilibrium thermodynamics and the details of how to find the matching parameters are investigated. The physical meaning of parameters such as the net heat of transport and the activation energy of viscous flow is elucidated, as the literature is overwhelmed with confusing and misleading information. The so-called dynamic and static models and their relations to the matching and non-matching parameter models are also discussed. We conclude that modeling the net heat of transport by the activation energy of self-diffusion may provide better results than approximating it by the activation energy of viscous flow. Nonetheless, the matching parameter models, which use the activation energy of viscous flow, are more dynamic and predict the thermal diffusion factor better than the non-matching parameter or static models, such as those of Kempers and Haase.
The choking pressure ratio of a critical flow venturi
NASA Technical Reports Server (NTRS)
Hillbrath, H. S.; Dill, W. P.; Wacker, W. A.
1973-01-01
The critical flow venturi has many important applications in the measurement and control of gas flow. In many of these applications, it is desirable to minimize the pressure loss required to maintain critical flow conditions. The performance of the venturi may be characterized by the ratio of outlet static pressure to inlet total pressure just sufficiently small to produce critical flow. This ratio is called choking pressure ratio. The optimization of diffusers for critical flow venturis is discussed and suggestions for designs practice are presented. Test results are given for six different diffuser configurations, and a comparison is made with data on 11 configurations from other investigators. It is shown that, for the small divergence angles considered, a simply defined diffuser effectiveness parameter is approximately independent of flow conditions and may be used to predict choking pressure ratio.
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
Investigation on critical heat flux of flow in pipes
NASA Astrophysics Data System (ADS)
Zhu, Senyuan
1990-08-01
This paper experimentally and theoretically investigates the critical heat flux of flow in pipes. From the analysis of the boiling mechanism and processing by means of the analogy principle of two-phase flow, a criterion equation to express critical heat flux has been derived. Correlated with six different coolants, 355 experimental data, the constant A and exponents K, m, and n are obtained. With a dimensionless correction term to calculate the effect on the varying slotted height of the cooling jacket, the previous equation will be a general equation to calculate the critical heat flux of flow in pipes.
Two critical issues in Langevin simulation of gas flows
NASA Astrophysics Data System (ADS)
Zhang, Jun; Fan, Jing
2014-12-01
A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-Bénard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.
Two critical issues in Langevin simulation of gas flows
Zhang, Jun; Fan, Jing
2014-12-09
A stochastic algorithm based on the Langevin equation has been recently proposed to simulate rarefied gas flows. Compared with the direct simulation Monte Carlo (DSMC) method, the Langevin method is more efficient in simulating small Knudsen number flows. While it is well-known that the cell sizes and time steps should be smaller than the mean free path and the mean collision time, respectively, in DSMC simulations, the Langevin equation uses a drift term and a diffusion term to describe molecule movements, so no direct molecular collisions have to be modeled. This enables the Langevin simulation to proceed with a much larger time step than that in the DSMC method. Two critical issues in Langevin simulation are addressed in this paper. The first issue is how to reproduce the transport properties as that described by kinetic theory. Transport coefficients predicted by Langevin equation are obtained by using Green-Kubo formulae. The second issue is numerical scheme with boundary conditions. We present two schemes corresponding to small time step and large time step, respectively. For small time step, the scheme is similar to DSMC method as the update of positions and velocities are uncoupled; for large time step, we present an analytical solution of the hitting time, which is the crucial factor for accurate simulation. Velocity-Couette flow, thermal-Couette flow, Rayleigh-Bénard flow and wall-confined problem are simulated by using these two schemes. Our study shows that Langevin simulation is a promising tool to investigate small Knudsen number flows.
NASA Technical Reports Server (NTRS)
Steffen, Christopher J., Jr.
1993-01-01
Turbulent backward-facing step flow was examined using four low turbulent Reynolds number k-epsilon models and one standard high Reynolds number technique. A tunnel configuration of 1:9 (step height: exit tunnel height) was used. The models tested include: the original Jones and Launder; Chien; Launder and Sharma; and the recent Shih and Lumley formulation. The experimental reference of Driver and Seegmiller was used to make detailed comparisons between reattachment length, velocity, pressure, turbulent kinetic energy, Reynolds shear stress, and skin friction predictions. The results indicated that the use of a wall function for the standard k-epsilon technique did not reduce the calculation accuracy for this separated flow when compared to the low turbulent Reynolds number techniques.
NASA Astrophysics Data System (ADS)
Steffen, Christopher J., Jr.
1993-06-01
Turbulent backward-facing step flow was examined using four low turbulent Reynolds number k-epsilon models and one standard high Reynolds number technique. A tunnel configuration of 1:9 (step height: exit tunnel height) was used. The models tested include: the original Jones and Launder; Chien; Launder and Sharma; and the recent Shih and Lumley formulation. The experimental reference of Driver and Seegmiller was used to make detailed comparisons between reattachment length, velocity, pressure, turbulent kinetic energy, Reynolds shear stress, and skin friction predictions. The results indicated that the use of a wall function for the standard k-epsilon technique did not reduce the calculation accuracy for this separated flow when compared to the low turbulent Reynolds number techniques.
Disappearance of criticality in thermal explosion for reactive viscous flows
Shonhiwa, T.; Zaturska, M.B.
1987-02-01
The thermal stability of a reactive viscous flow was first investigated by Adler, who considered the steady development flow between symmetrically heated parallel walls. He used a power series in a defined viscous heating parameter and, assuming the Frank-Kamenetskii exponential approximation to the Arrhenius term, obtained an expression for the critical Frank-Kamenetskii parameter in series form. Subsequently Zaturska used the same approach for four other simple flows, each an exact solution of the Navier-Stokes equations, namely, Poiseuille pipe flow, axial flow between concentric circular cylinders, rotating flow between concentric circular cylinders, and plane Couette flow, and also applied the step-function reaction-rate approximation and constant reaction-rate approximation to the five flows.
The critical thinking curriculum model
NASA Astrophysics Data System (ADS)
Robertson, William Haviland
The Critical Thinking Curriculum Model (CTCM) utilizes a multidisciplinary approach that integrates effective learning and teaching practices with computer technology. The model is designed to be flexible within a curriculum, an example for teachers to follow, where they can plug in their own critical issue. This process engages students in collaborative research that can be shared in the classroom, across the country or around the globe. The CTCM features open-ended and collaborative activities that deal with current, real world issues which leaders are attempting to solve. As implemented in the Critical Issues Forum (CIF), an educational program administered by Los Alamos National Laboratory (LANL), the CTCM encompasses the political, social/cultural, economic, and scientific realms in the context of a current global issue. In this way, students realize the importance of their schooling by applying their efforts to an endeavor that ultimately will affect their future. This study measures student attitudes toward science and technology and the changes that result from immersion in the CTCM. It also assesses the differences in student learning in science content and problem solving for students involved in the CTCM. A sample of 24 students participated in classrooms at two separate high schools in New Mexico. The evaluation results were analyzed using SPSS in a MANOVA format in order to determine the significance of the between and within-subjects effects. A comparison ANOVA was done for each two-way MANOVA to see if the comparison groups were equal. Significant findings were validated using the Scheffe test in a Post Hoc analysis. Demographic information for the sample population was recorded and tracked, including self-assessments of computer use and availability. Overall, the results indicated that the CTCM did help to increase science content understanding and problem-solving skills for students, thereby positively effecting critical thinking. No matter if the
Near-Critical CO2 Flow Measurements and Visualization
NASA Astrophysics Data System (ADS)
Kazemifar, Farzan; Kyritsis, Dimitrios
2012-11-01
Carbon dioxide capturing and sequestration is one of the proposed solutions for reducing greenhouse gas emission. This technique will be used in big industrial plants with very high CO2 emissions. Handling such large flow rates requires high pressure and low temperature (in order to maximize density and minimize volumetric flow rate) which brings us close to the critical point of CO2 at approximately 74 bar and 31°C. This necessitates studying near-critical CO2 flows. In our experiment setup CO2 is compressed to supercritical pressures using a hydraulic accumulator. Pressurized CO2 then flows through the test section, which is a 2-ft long stainless steel tube with ID = 0.084 in. The flow rate is controlled by a needle valve downstream of the test section and the mass flow rate is measured using a coriolis mass flow meter. Temperature and pressure are monitored using two K-type thermocouples and pressure transducers at the inlet and exit of the test section. The pressure difference across the pipe is measured separately using a differential pressure transducer. In another set of experiments, the aforementioned test section is replaced with an optically accessible test section. In this setup high-speed imaging is used to visualize the flow inside the test section. We studied the recorded data in order to identify distinct flow regimes based on pressure drop as a function of pressure, temperature and mass flow rate. Acknowledgements: International Institute for Carbon-Neutral Energy Research (I2CNER).
Experimental study of the critical mode of steam-water flow
NASA Astrophysics Data System (ADS)
Shulyupin, A. N.
2011-12-01
This paper presents the results of measurements of static pressure in the critical mode of steam-water flow from long pipes of diameter 50.2 and 100.1 mm at various distances from their outlet section. The measurements were performed at flow rates and enthalpies of the mixture characteristic of geothermal coolants. It is found that the mixture flow rates obtained experimentally significantly exceed the estimation results obtained using well-known models.
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.
Numerical investigation of transition critical Reynolds number of channel flow.
NASA Astrophysics Data System (ADS)
Zhang, Yongming
2015-11-01
Two critical Reynolds numbers are mentioned in investigation of laminar-turbulent transition. One is instability critical Reynolds number from linear stability theory (LST). The other is transition critical Reynolds number at which transition occurs in reality, which is significantly lower than the former in general. The determination of transition critical Reynolds number is of important practical significance in some engineering problems. Theoretical method has not been proposed for its determination, so it has to depend on experiments. However, for some flows with important practical significance, such as hypersonic boundary layer, transition critical Reynolds number cannot be determined by experiments in current situation. In this paper, transition critical Reynolds number of incompressible channel flow is determined by direct numerical simulations (DNS). It is found as Re =1114, which agrees with experimental data. In subsequent paper, transition critical Reynolds number of boundary layer will be investigation by the similar method. Project supported by the National Natural Science Foundation of China (Nos. 11202147, 11332007, 11172203, and 91216111) and the Specialized Research Fund (New Teacher Class) for the Doctoral Program of Higher Education (No. 20120032120007).
The critical layer in pipe flow at high Reynolds number.
Viswanath, D
2009-02-13
We report the computation of a family of travelling wave solutions of pipe flow up to Re=75000. As in all lower branch solutions, streaks and rolls feature prominently in these solutions. For large Re, these solutions develop a critical layer away from the wall. Although the solutions are linearly unstable, the two unstable eigenvalues approach 0 as Re-->infinity at rates given by Re-0.41 and Re-0.87; surprisingly, the solutions become more stable as the flow becomes less viscous. The formation of the critical layer and other aspects of the Re-->infinity limit could be universal to lower branch solutions of shear flows. We give implementation details of the GMRES-hookstep and Arnoldi iterations used for computing these solutions and their spectra, while pointing out the new aspects of our method. PMID:18990661
Critical velocity of superfluid helium flow in fine pore filters
NASA Astrophysics Data System (ADS)
Hofmann, A.
1990-03-01
Membrane filters with porosities of more than 70 pct have been proven to be useful for the fabrication of thermomechanical pumps (fountain effect pumps) with mass fluxes up to 1.7 g/sq cm at about 0.25 bar head of pressure. Different pumps made of 14 mm thick stacks of commercially available Sartorius cellulose nitrate membranes with filtration grades of 200, 50 and 10 nm have been investigated at various temperatures between 1.6 K at the inlet and 2.1 K at the outlet. By analysis of pressure and temperature measurements, it is shown that the chemical potential difference between inlet and outlet is zero up to critical flow rates very close to the peak flow and that it increases steeply when the critical flow is exceeded.
Upchurch, E.R. )
1989-08-01
The critical flow rate of a gas well is the minimum flow rate required to prevent accumulation of liquids in the tubing. Theoretical models currently available for estimating critical flow rates are restricted to wells with water/gas ratios less than 150bbl/MMcf (0.84 X 10/sup -3/ m/sup 3//m/sup 3/). For wells producing at higher water/gas ratios from normally pressured waterdrive reservoirs, a method of estimating critical flow rates is derived through use of an empirical multiphase-flow correlation.
Parametric analyses of planned flowing uranium hexafluoride critical experiments
NASA Technical Reports Server (NTRS)
Rodgers, R. J.; Latham, T. S.
1976-01-01
Analytical investigations were conducted to determine preliminary design and operating characteristics of flowing uranium hexafluoride (UF6) gaseous nuclear reactor experiments in which a hybrid core configuration comprised of UF6 gas and a region of solid fuel will be employed. The investigations are part of a planned program to perform a series of experiments of increasing performance, culminating in an approximately 5 MW fissioning uranium plasma experiment. A preliminary design is described for an argon buffer gas confined, UF6 flow loop system for future use in flowing critical experiments. Initial calculations to estimate the operating characteristics of the gaseous fissioning UF6 in a confined flow test at a pressure of 4 atm, indicate temperature increases of approximately 100 and 1000 K in the UF6 may be obtained for total test power levels of 100 kW and 1 MW for test times of 320 and 32 sec, respectively.
Critical assessment of automated flow cytometry data analysis techniques.
Aghaeepour, Nima; Finak, Greg; Hoos, Holger; Mosmann, Tim R; Brinkman, Ryan; Gottardo, Raphael; Scheuermann, Richard H
2013-03-01
Traditional methods for flow cytometry (FCM) data processing rely on subjective manual gating. Recently, several groups have developed computational methods for identifying cell populations in multidimensional FCM data. The Flow Cytometry: Critical Assessment of Population Identification Methods (FlowCAP) challenges were established to compare the performance of these methods on two tasks: (i) mammalian cell population identification, to determine whether automated algorithms can reproduce expert manual gating and (ii) sample classification, to determine whether analysis pipelines can identify characteristics that correlate with external variables (such as clinical outcome). This analysis presents the results of the first FlowCAP challenges. Several methods performed well as compared to manual gating or external variables using statistical performance measures, which suggests that automated methods have reached a sufficient level of maturity and accuracy for reliable use in FCM data analysis. PMID:23396282
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.
The secondary flow near a baroclinic planetary wave critical line
NASA Technical Reports Server (NTRS)
Schoeberl, M. R.
1981-01-01
A critical line (CL) is the surface where the phase speed of a wave in a fluid is equal to the speed of the background flow. The considered investigation is concerned with one aspect of the simplest model of a CL in which the CL is assumed to totally absorb energy from steady, stationary, planetary waves. The aspect of interest is the secondary mean circulation near a CL in a baroclinic atmosphere. The motivation for this study is the observation of the nearly vertical CL by O'Neill and Taylor (1979) which appeared during the sudden warming of 1976/77. Even though the treatment of the CL is highly idealized in the investigation, there is evidence which indicates a very large rate of change in the zonally averaged temperature along a CL may occur. The Lagrangian-mean properties of an idealized baroclinic CL are also examined. It is found that the Lagrangian jets may provide an important transport process for exchange of stratospheric and tropospheric air.
The Biomantle-Critical Zone Model
NASA Astrophysics Data System (ADS)
Johnson, D. L.; Lin, H.
2006-12-01
It is a fact that established fields, like geomorphology, soil science, and pedology, which treat near surface and surface processes, are undergoing conceptual changes. Disciplinary self examinations are rife. New practitioners are joining these fields, bringing novel and interdisciplinary ideas. Such new names as "Earth's critical zone," "near surface geophysics," and "weathering engine" are being coined for research groups. Their agendas reflect an effort to integrate and reenergize established fields and break new ground. The new discipline "hydropedology" integrates soil science with hydrologic principles, and recent biodynamic investigations have spawned "biomantle" concepts and principles. One force behind these sea shifts may be retrospectives whereby disciplines periodically re-invent themselves to meet new challenges. Such retrospectives may be manifest in the recent Science issue on "Soils, The Final Frontier" (11 June, 2004), and in recent National Research Council reports that have set challenges to science for the next three decades (Basic Research Opportunities in Earth Science, and Grand Challenges for the Environmental Sciences, both published in 2001). In keeping with such changes, we advocate the integration of biomantle and critical zone concepts into a general model of Earth's soil. (The scope of the model automatically includes the domain of hydropedology.) Our justification is that the integration makes for a more appealing holistic, and realistic, model for the domain of Earth's soil at any scale. The focus is on the biodynamics of the biomantle and water flow within the critical zone. In this general model the biomantle is the epidermis of the critical zone, which extends to the base of the aquifer. We define soil as the outer layer of landforms on planets and similar bodies altered by biological, chemical, and/or physical agents. Because Earth is the only planet with biological agents, as far as we know, it is the only one that has all
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.
Soil organic phosphorus flows to water via critical and non-critical hydrological source areas
NASA Astrophysics Data System (ADS)
Wang, Ying; Surridge, Ben; Haygarth, Phil
2015-04-01
Soil organic phosphorus flows to water via critical and non-critical hydrological source areas Ying Wang, Ben W.J. Surridge, Philip M. Haygarth Lancaster Environment Centre, Lancaster University, LA1 4YQ, UK Critical source areas (CSAs) are zones in the landscape where easily connected hydrology coincides with a phosphorus (P) sources in the soil. The P export risks in CSAs are hypothesised to be higher compared with non-critical source areas (Non-CSAs) and specifically that the magnitudes of P forms in CSA areas were higher than Non-CSAs. Past research on CSAs has often neglected forms of organic P, such as DNA and phospholipids which are among the most potentially biodegradable organic P compounds. The objectives of this study were i) to quantify the magnitude of organic P compounds in agricultural soils and specifically determine whether these magnitudes differed significantly between CSAs and Non-CSAs; ii) determine the variation of P magnitude between and within individual fields; iii) identify the P delivery concentrations in soil solution after raining events in CSAs. The study focussed on soils collected from the Morland sub-catchment of the River Eden catchment in Cumbria, northern England. CSA and Non-CSA pairs were identified using the SCIMAP modelling and field assessment providing five CSA - Non-CSA pairs in total. The results showed that there are significant differences in the total P (TP) concentrations, the proportions of DNA-P, WETP (water extractable total P), WERP (water extractable reactive P) and WEUP (water extractable unreactive P) between CSA and Non-CSA. We also found that the concentrations of all the P forms showed distribution variation between fields or even within the same field. Liable organic P such as DNA-P and PLD-P was presented considerable proportions of total P in soil, especially DNA-P which had a good correlation with TP. DNA-P in the ten areas accounted for a considerable proportion of soil TP (4.9 to 16.6%). Given the
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.
A critical review on photochemical conversions in flow analysis.
Rocha, Diogo L; Kamogawa, Marcos Y; Rocha, Fábio R P
2015-10-01
Photochemical conversions are cost-effective and environmental friendly processes that require mild experimental conditions and avoid generation of highly acidic wastes. Treated samples are then compatible with most of the analytical techniques. These characteristics become more relevant when the photoconversions are accomplished to flow analysis, thus allowing exploitation of incomplete reactions, the effective use of the photogenerated unstable radicals and in-line sample treatment. Decreasing of reagent consumption and waste generation, sample processing in a closed environment, and improvement of efficiency of the photochemical processes are other inherent advantages. These aspects are critically reviewed in this article, which emphasizes applications to fractionation and speciation analysis, photo-induced luminescence, miniaturization, and in-line waste treatment. Design of flow-through photochemical cells, use of auxiliary reagents in homogeneous and heterogeneous media, and configurations of flow manifolds are also discussed. PMID:26481985
CRITICAL ASSESSMENT OF AUTOMATED FLOW CYTOMETRY DATA ANALYSIS TECHNIQUES
Aghaeepour, Nima; Finak, Greg; Hoos, Holger; Mosmann, Tim R.; Gottardo, Raphael; Brinkman, Ryan; Scheuermann, Richard H.
2013-01-01
Traditional methods for flow cytometry (FCM) data processing rely on subjective manual gating. Recently, several groups have developed computational methods for identifying cell populations in multidimensional FCM data. The Flow Cytometry: Critical Assessment of Population Identification Methods (FlowCAP) challenges were established to compare the performance of these methods on two tasks – mammalian cell population identification to determine if automated algorithms can reproduce expert manual gating, and sample classification to determine if analysis pipelines can identify characteristics that correlate with external variables (e.g., clinical outcome). This analysis presents the results of the first of these challenges. Several methods performed well compared to manual gating or external variables using statistical performance measures, suggesting that automated methods have reached a sufficient level of maturity and accuracy for reliable use in FCM data analysis. PMID:23396282
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.
Enhancement of critical heat flux in tubes using staged tangential flow injection: (Progress report)
Dhir, V.K.
1987-01-01
Experimental studies of the enhancement in single and two phase heat transfer from tubes subjected to tangential flow injection have been continuing. Investigations using water as the test liquid have been focused on: single phase heat transfer coefficients; two phase heat transfer coefficients under subcooled boiling conditions; subcooled critical heat fluxes; and modeling of the enhancement under swirl flow conditions. With tangential injection up to four fold increase in the average heat transfer coefficient has been observed. During subcooled boiling the enhancement is relatively small. However swirl induced centripetal force increases vapor escape velocity and as a result higher critical heat fluxes can be accommodated. In the range of flow parameters studied up to 40% enhancement in critical heat flux has been observed with single stage injection. This enhancement is slightly less than that obtained with Freon-113. The mechanistic reasons for this observation are currently being investigated.
Enhancement of critical heat flux in tubes using staged tangential flow injection
NASA Astrophysics Data System (ADS)
Dhir, V. K.
Experimental studies of the enhancement in single and two phase heat transfer from tubes subjected to tangential flow injection have been continuing. Investigations using water as the test liquid have been focused on: single phase heat transfer coefficients; two phase heat transfer coefficients under subcooled boiling conditions; subcooled critical heat fluxes; and modeling of the enhancement under swirl flow conditions. With tangential injection up to four fold increase in the average heat transfer coefficient has been observed. During subcooled boiling the enhancement is relatively small. However swirl induced centripetal force increases vapor escape velocity and as a result higher critical heat fluxes can be accommodated. In the range of flow parameters studied up to 40% enhancement in critical heat flux has been observed with single stage injection. This enhancement is slightly less than that obtained with Freon-113. The mechanistic reasons for this observation are currently being investigated.
Critical ignition in rapidly expanding self-similar flows
NASA Astrophysics Data System (ADS)
Radulescu, Matei I.; Maxwell, Brian M.
2010-06-01
The generic problem of ignition of a particle undergoing an expansion given by a power law rate of decay behind a decaying shock is addressed in the present study. It is demonstrated, using a one-step Arrhenius irreversible reaction, that a sufficiently strong expansion wave can quench the reaction. The critical conditions for extinction are obtained in closed form in terms of the time scale for the expansion process and the thermochemical properties of the gas, yielding a critical Damkohler number, i.e., the ratio of the expansion time scale to the homogeneous ignition time scale, given by (γ -1)(Ea/RT)-1/n, where n is the power law exponent of the self-similar expansion. The critical ignition criteria, which are valid in the asymptotic limit n(γ -1)(Ea/RT)=O(1), were found in excellent agreement with numerical results. The applicability of the results obtained are discussed for ignition in rapidly expanding flows which occur behind decaying shock waves, as encountered in problems of detonation initiation by a Taylor-Sedov blast wave, and reacting jet startup, and for reactions in steady hypersonic flows around projectiles.
Siphon flows in isolated magnetic flux tubes. IV - Critical flows with standing tube shocks
NASA Technical Reports Server (NTRS)
Thomas, John H.; Montesinos, Benjamin
1991-01-01
Critical siphon flows in arched, isolated magnetic flux tubes are studied within the thin flux tube approximation, with a view toward applications to intense magnetic flux concentrations in the solar photosphere. The results of calculations of the strength and position of the standing tube shock in the supercritical downstream branch of a critical siphon flow are presented, as are calculations of the flow variables all along the flux tube and the equilibrium path of the flux tube in the surrounding atmosphere. It is suggested that arched magnetic flux tubes, with magnetic field strength increased by a siphon flow, may be associated with some of the intense, discrete magnetic elements observed in the solar photosphere.
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.
Influence of two-phase flow characteristic on critical heat flux in low pressure
Inoue, Akira; Lee, Sangryoul
1996-08-01
Estimation of the critical heat flux (CHF) in a boiling two-phase flow is one of the important subjects for the safety of water-cooled reactors and other energy systems. In the case of a boiling two-phase flow at low pressure, flow pattern and void fraction are easy to change by the power input and the flow becomes more complex due to low density of gas phase. The CHF is affected by the flow pattern. In this study, the CHFs were measured over wide quality range from the subcooled boiling to the annular-mist flow. By using Pyrex glass tube as a test channel, the two-phase flow situation was observed. Graphite rod or stainless steel tube was used as a heater rod and installed at the center of the glass tube. Two-phase flow was formed by steam injection to circulating water at an upstream region of the test section. The flow pattern was kept nearly constant over the length of test section due to the low input power density into the fluid. Then, the characteristics of CHF could be investigated at each flow patterns of bubbly, slug, annular and annular-mist flow. In the subcooled boiling region of bubbly flow, the CHF decreased with increase of quality and was less sensitive to flow rate. In the slug flow region, the CHF showed a minimum value. With more increase of quality in the annular flow, the CHF increased and reached a peak value at a certain quality depending on a flow rate. The peak of CHF occurred almost at a constant vapor mass velocity. In the annular-mist flow region, the CHF decreased with increase of quality. In the region, the effect of heated length on the CHF was systematically measured and validity of an analytical model considering dryout of liquid film based on formation of a dry patch was investigated.
Criticality in the Burridge-Knopoff model
Clancy, Ian; Corcoran, David
2005-04-01
Criticality is a potential origin of the scale invariance observed in the Gutenberg-Richter law for earthquakes. In support of this hypothesis, the Burridge-Knopoff (BK) model of an earthquake fault system is known to exhibit a dynamic phase transition, but the critical nature of the transition is uncertain. Here it is shown that the BK model exhibits a dynamic transition from large-scale stick-slip to small-scale creep motion and through a finite size scaling analysis the critical nature of this transition is established. The order parameter describing the critical transition suggests that the Olami-Feder-Christensen model may be tuned to criticality through its assumptions describing the relaxation of the system.
River Discharge Estimation Using Imaged Critical Flow Phenomena
NASA Astrophysics Data System (ADS)
Fonstad, M. A.; Grant, G.
2015-12-01
A wide variety of river science applications require remote estimation of discharge. Flow in steep rivers often approach critical flow (Froude number equal to one), as evidenced by trains of standing waves either perpendicular or at a high angle relative to the flow direction. Previous work has demonstrated that where such waves are present, water depth can be estimated without the need for a roughness coefficient such as Manning's n. We extend this prior work to remotely calculate river discharge. The wavelength of standing waves can be measured using high-resolution remote sensing imagery. Velocity can then be calculated from wavelength using the Kennedy wave equation. Assuming critical flow allows depth to be calculated using the modified Froude number equation, and with an additional measurement of imaged river width, river discharge can be computed directly as the product of these three values. We test this approach using high-resolution Google Earth imagery of rivers with standing waves near existing stream gages. We also demonstrate the utility of this approach by extracting a drainage area-to-discharge relationship for a large watershed. There are certain challenges with this approach, the greatest being the need for high resolution (meter-scale or better) imagery to see and measure standing waves. Such waves are also easily confused with wind waves, turbulence or other surface effects. Nevertheless, this approach offers promise for both estimating discharge in places without extensive gage networks, and also estimating discharge in archival imagery. Repeat imaging of the same areas might also be used to construct at-a-station and downstream hydraulic geometry relationships.
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.
Flow Boiling Critical Heat Flux in Reduced Gravity
NASA Technical Reports Server (NTRS)
Mudawar, Issam; Zhang, Hui; Hasan, Mohammad M.
2004-01-01
This study provides systematic method for reducing power consumption in reduced gravity systems by adopting minimum velocity required to provide adequate CHF and preclude detrimental effects of reduced gravity . This study proves it is possible to use existing 1 ge flow boiling and CHF correlations and models to design reduced gravity systems provided minimum velocity criteria are met
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.
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.
Hydrodynamic instabilities of near-critical CO2 flow in microchannels: Lattice Boltzmann simulation
NASA Astrophysics Data System (ADS)
Holdych, D. J.; Georgiadis, J. G.; Buckius, R. O.
2004-05-01
Motivated by systematic CO2 evaporation experiments which recently became available (J. Pettersen, "Flow vaporization of CO2 in microchannel tubes," Doctor technicae thesis, Norwegian University of Science and Technology, 2002), the present work constitutes an exploratory investigation of isothermal flow of CO2 near its liquid-vapor critical point through a long 5 μm diameter microchannel. A modified van der Waals constitutive model—with properties closely approximating those of "real" near-critical CO2—is incorporated in a two-dimensional lattice Boltzmann hydrodynamics model by embedding a dimensionless parameter X, with X→1 denoting the "real" fluid. The hydrodynamic phenomena resulting by imposing a constant pressure gradient along a periodic channel are investigated by considering two regimes in tandem: (1) transition from bubbly to annular flow with a liquid film formed at the channel walls and (2) destabilization of the liquid film by the Kelvin-Helmholtz instability. Due to numerical constraints, intrinsic modeling errors are introduced and are shown to be associated with discrepancies in the relative vapor-liquid interfacial thickness, which is expressed by X. The effects of these errors are investigated both theoretically and numerically in the physical limit X→1. Numerically determined flow patterns compare qualitatively well with direct visualization results obtained by Pettersen. Overall, the characteristics of isothermal near-critical two-phase flow in microchannels can be reproduced by the appropriate modification of the thermophysical properties of CO2.
Phenomenology of a flow around a circular cylinder at sub-critical and critical Reynolds numbers
NASA Astrophysics Data System (ADS)
Capone, Alessandro; Klein, Christian; Di Felice, Fabio; Miozzi, Massimo
2016-07-01
In this work, the flow around a circular cylinder is investigated at Reynolds numbers ranging from 79 000 up to 238 000 by means of a combined acquisition system based on Temperature Sensitive Paint (TSP) and particle velocimetry. The proposed setup allows simultaneous and time-resolved measurement of absolute temperature and relative skin friction fields onto the cylinder surface and near-wake velocity field. Combination of time-resolved surface measurements and planar near-field velocity data allows the investigation of the profound modifications undergone by the wall shear stress topology and its connections to the near-field structure as the flow regime travels from the sub-critical to the critical regime. Laminar boundary-layer separation, transition, and re-attachment are analyzed in the light of temperature, relative skin friction maps, and Reynolds stress fields bringing about a new perspective on the relationship between boundary layer development and shear layer evolution. The fast-responding TSP employed allows high acquisition frequency and calculation of power spectral density from surface data. Correlation maps of surface and near-wake data provide insight into the relationship between boundary-layer evolution and vortex shedding. We find that as the Reynolds number approaches the critical state, the separation line oscillations feature an increasingly weaker spectrum peak compared to the near-wake velocity spectrum. In the critical regime, separation line oscillations are strongly reduced and the correlation to the local vorticity undergoes an overall decrease giving evidence of modifications in the vortex shedding mechanism.
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.
Nonlocal modeling of granular flows down inclines.
Kamrin, Ken; Henann, David L
2015-01-01
Flows of granular media down a rough inclined plane demonstrate a number of nonlocal phenomena. We apply the recently proposed nonlocal granular fluidity model to this geometry and find that the model captures many of these effects. Utilizing the model's dynamical form, we obtain a formula for the critical stopping height of a layer of grains on an inclined surface. Using an existing parameter calibration for glass beads, the theoretical result compares quantitatively to existing experimental data for glass beads. This provides a stringent test of the model, whose previous validations focused on driven steady-flow problems. For layers thicker than the stopping height, the theoretical flow profiles display a thickness-dependent shape whose features are in agreement with previous discrete particle simulations. We also address the issue of the Froude number of the flows, which has been shown experimentally to collapse as a function of the ratio of layer thickness to stopping height. While the collapse is not obvious, two explanations emerge leading to a revisiting of the history of inertial rheology, which the nonlocal model references for its homogeneous flow response. PMID:25376561
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.
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.
Potts-model critical manifolds revisited
NASA Astrophysics Data System (ADS)
Scullard, Christian R.; Lykke Jacobsen, Jesper
2016-03-01
We compute critical polynomials for the q-state Potts model on the Archimedean lattices, using a parallel implementation of the algorithm of Jacobsen (2014 J. Phys. A: Math. Theor 47 135001) that gives us access to larger sizes than previously possible. The exact polynomials are computed for bases of size 6 × 6 unit cells, and the root in the temperature variable v={{{e}}}K-1 is determined numerically at q = 1 for bases of size 8 × 8. This leads to improved results for bond percolation thresholds, and for the Potts-model critical manifolds in the real (q, v) plane. In the two most favourable cases, we find now the kagome-lattice threshold to eleven digits and that of the (3,{12}2) lattice to thirteen. Our critical manifolds reveal many interesting features in the antiferromagnetic region of the Potts model, and determine accurately the extent of the Berker-Kadanoff phase for the lattices studied.
Model flocks in a steady vortical flow.
Baggaley, A W
2015-05-01
We modify the standard Vicsek model to clearly distinguish between intrinsic noise due to imperfect alignment between organisms and extrinsic noise due to fluid motion. We then consider the effect of a steady vortical flow, the Taylor-Green vortex, on the dynamics of the flock, for various flow speeds, with a fixed intrinsic particle speed. We pay particular attention to the morphology of the flow, and quantify its filamentarity. Strikingly, above a critical flow speed there is a pronounced increase in the filamentarity of the flock, when compared to the zero-flow case. This is due to the fact that particles appear confined to areas of low vorticity; a familiar phenomena, commonly seen in the clustering of inertial particles in vortical flows. Hence, the cooperative motion of the particles gives them an effective inertia, which is seen to have a profound effect on the morphology of the flock, in the presence of external fluid motion. Finally, we investigate the angle between the flow and the particles direction of movement and find it follows a power-law distribution. PMID:26066260
Model flocks in a steady vortical flow
NASA Astrophysics Data System (ADS)
Baggaley, A. W.
2015-05-01
We modify the standard Vicsek model to clearly distinguish between intrinsic noise due to imperfect alignment between organisms and extrinsic noise due to fluid motion. We then consider the effect of a steady vortical flow, the Taylor-Green vortex, on the dynamics of the flock, for various flow speeds, with a fixed intrinsic particle speed. We pay particular attention to the morphology of the flow, and quantify its filamentarity. Strikingly, above a critical flow speed there is a pronounced increase in the filamentarity of the flock, when compared to the zero-flow case. This is due to the fact that particles appear confined to areas of low vorticity; a familiar phenomena, commonly seen in the clustering of inertial particles in vortical flows. Hence, the cooperative motion of the particles gives them an effective inertia, which is seen to have a profound effect on the morphology of the flock, in the presence of external fluid motion. Finally, we investigate the angle between the flow and the particles direction of movement and find it follows a power-law distribution.
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.
PHYSICAL MODELING OF CONTRACTED FLOW.
Lee, Jonathan K.
1987-01-01
Experiments on steady flow over uniform grass roughness through centered single-opening contractions were conducted in the Flood Plain Simulation Facility at the U. S. Geological Survey's Gulf Coast Hydroscience Center near Bay St. Louis, Miss. The experimental series was designed to provide data for calibrating and verifying two-dimensional, vertically averaged surface-water flow models used to simulate flow through openings in highway embankments across inundated flood plains. Water-surface elevations, point velocities, and vertical velocity profiles were obtained at selected locations for design discharges ranging from 50 to 210 cfs. Examples of observed water-surface elevations and velocity magnitudes at basin cross-sections are presented.
Quantum Criticality in the Biased Dicke Model.
Zhu, Hanjie; Zhang, Guofeng; Fan, Heng
2016-01-01
The biased Dicke model describes a system of biased two-level atoms coupled to a bosonic field, and is expected to produce new phenomena that are not present in the original Dicke model. In this paper, we study the critical properties of the biased Dicke model in the classical oscillator limits. For the finite-biased case in this limit, We present analytical results demonstrating that the excitation energy does not vanish for arbitrary coupling. This indicates that the second order phase transition is avoided in the biased Dicke model, which contrasts to the original Dicke model. We also analyze the squeezing and the entanglement in the ground state, and find that a finite bias will strongly modify their behaviors in the vicinity of the critical coupling point. PMID:26786239
Quantum Criticality in the Biased Dicke Model
Zhu, Hanjie; Zhang, Guofeng; Fan, Heng
2016-01-01
The biased Dicke model describes a system of biased two-level atoms coupled to a bosonic field, and is expected to produce new phenomena that are not present in the original Dicke model. In this paper, we study the critical properties of the biased Dicke model in the classical oscillator limits. For the finite-biased case in this limit, We present analytical results demonstrating that the excitation energy does not vanish for arbitrary coupling. This indicates that the second order phase transition is avoided in the biased Dicke model, which contrasts to the original Dicke model. We also analyze the squeezing and the entanglement in the ground state, and find that a finite bias will strongly modify their behaviors in the vicinity of the critical coupling point. PMID:26786239
Debris flows: Experiments and modelling
NASA Astrophysics Data System (ADS)
Turnbull, Barbara; Bowman, Elisabeth T.; McElwaine, Jim N.
2015-01-01
Debris flows and debris avalanches are complex, gravity-driven currents of rock, water and sediments that can be highly mobile. This combination of component materials leads to a rich morphology and unusual dynamics, exhibiting features of both granular materials and viscous gravity currents. Although extreme events such as those at Kolka Karmadon in North Ossetia (2002) [1] and Huascarán (1970) [2] strongly motivate us to understand how such high levels of mobility can occur, smaller events are ubiquitous and capable of endangering infrastructure and life, requiring mitigation. Recent progress in modelling debris flows has seen the development of multiphase models that can start to provide clues of the origins of the unique phenomenology of debris flows. However, the spatial and temporal variations that debris flows exhibit make this task challenging and laboratory experiments, where boundary and initial conditions can be controlled and reproduced, are crucial both to validate models and to inspire new modelling approaches. This paper discusses recent laboratory experiments on debris flows and the state of the art in numerical models.
Preserving Flow Variability in Watershed Model Calibrations
Background/Question/Methods Although watershed modeling flow calibration techniques often emphasize a specific flow mode, ecological conditions that depend on flow-ecology relationships often emphasize a range of flow conditions. We used informal likelihood methods to investig...
Single point modeling of rotating turbulent flows
NASA Technical Reports Server (NTRS)
Hadid, A. H.; Mansour, N. N.; Zeman, O.
1994-01-01
A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.
Single point modeling of rotating turbulent flows
NASA Astrophysics Data System (ADS)
Hadid, A. H.; Mansour, N. N.; Zeman, O.
1994-12-01
A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.
Experimental Modelling of Debris Flows
NASA Astrophysics Data System (ADS)
Paleo Cageao, P.; Turnbull, B.; Bartelt, P.
2012-04-01
Debris flows are gravity-driven mass movements typically containing water, sediments, soil and rocks. These elements combine to give a flow complex phenomenology that exhibits characteristics common to diverse geophysical flows from dry granular media (e.g. levee formation) to viscous gravity currents (viscous fingering and surge instabilities). The exceptional speeds and range debris flows can achieve motivate the need for a co-ordinated modelling approach that can provide insight into the key physical processes that dictate the hazard associated with the flows. There has been recent progress in theoretical modelling approaches that capture the details of the multi-component nature of debris flows. The promise of such models is underlined by their qualitatively successful comparison with field-scale experimental data. The aim of the present work is to address the technical difficulties in achieving a controlled and repeatable laboratory-scale experiment for robust testing of these multi-component models. A laboratory experiment has been designed and tested that can provide detailed information of the internal structure of debris flows. This constitutes a narrow Perspex chute that can be tilted to any angle between 0° and ≈ 60°. A mixture of glycerine and glass balls was initially held behind a lock-gate, before being released down the chute. The evolving flow was captured through high speed video, analysed with a Particle Image Velocimetry algorithm to provide the changing velocity field. A wide parameter space has been tested, allowing variations in particle size, dispersity, surface roughness, fluid viscosity, slope angle and solid volume fraction. While matching key similarity criteria, such as Froude number, with a typical field event, these experiments allow close examination of a wide range of physical scenarios for the robust testing of new multi-component flow models. Further diagnostics include force plate and pore pressure measurements, with a view
On the criticality of inferred models
NASA Astrophysics Data System (ADS)
Mastromatteo, Iacopo; Marsili, Matteo
2011-10-01
Advanced inference techniques allow one to reconstruct a pattern of interaction from high dimensional data sets, from probing simultaneously thousands of units of extended systems—such as cells, neural tissues and financial markets. We focus here on the statistical properties of inferred models and argue that inference procedures are likely to yield models which are close to singular values of parameters, akin to critical points in physics where phase transitions occur. These are points where the response of physical systems to external perturbations, as measured by the susceptibility, is very large and diverges in the limit of infinite size. We show that the reparameterization invariant metrics in the space of probability distributions of these models (the Fisher information) are directly related to the susceptibility of the inferred model. As a result, distinguishable models tend to accumulate close to critical points, where the susceptibility diverges in infinite systems. This region is the one where the estimate of inferred parameters is most stable. In order to illustrate these points, we discuss inference of interacting point processes with application to financial data and show that sensible choices of observation time scales naturally yield models which are close to criticality.
NASA Astrophysics Data System (ADS)
Hasan, Nusair; Farouk, Bakhtier
2013-11-01
Forced convective thermal transport characteristics of supercritical carbon dioxide in vertical flow are numerically investigated. A tube with a circular cross-section and heated side-wall is considered. A real-fluid model for representing the thermo-physical properties of the supercritical fluid along with the fully compressible form of the Navier-Stokes equations and an implicit time-marching scheme is used to solve the problem. Thermo-physical properties of near-critical supercritical fluids show diverging characteristics. Large variations of density of near-critical supercritical fluid in forced convective flow can induce thermo-hydraulic instability similar to density wave oscillations. The developed numerical model is used for studying the effect of geometrical parameters of the tube, wall heat flux and pressure on steady-state convective thermal transport as well as the stability behavior of the supercritical fluid near its critical point. The enhancement or deterioration of heat transfer caused by the temperature-induced variation of physical properties (especially specific heat) is also investigated, as well as the effect of buoyancy on the forced convective flow.
Traffic flow forecasting: Comparison of modeling approaches
Smith, B.L.; Demetsky, M.J.
1997-08-01
The capability to forecast traffic volume in an operational setting has been identified as a critical need for intelligent transportation systems (ITS). In particular, traffic volume forecasts will support proactive, dynamic traffic control. However, previous attempts to develop traffic volume forecasting models have met with limited success. This research effort focused on developing traffic volume forecasting models for two sites on Northern Virginia`s Capital Beltway. Four models were developed and tested for the freeway traffic flow forecasting problem, which is defined as estimating traffic flow 15 min into the future. They were the historical average, time-series, neural network, and nonparametric regression models. The nonparametric regression model significantly outperformed the other models. A Wilcoxon signed-rank test revealed that the nonparametric regression model experienced significantly lower errors than the other models. In addition, the nonparametric regression model was easy to implement, and proved to be portable, performing well at two distinct sites. Based on its success, research is ongoing to refine the nonparametric regression model and to extend it to produce multiple interval forecasts.
Critical Velocity of a Superfluid Bose Gas Flowing in a Random Potential
NASA Astrophysics Data System (ADS)
Haga, Taiki
2016-05-01
We investigate the critical velocity of a weakly interacting Bose gas flowing in a random potential. By applying the Bogoliubov theory to a disordered Bose system with a steady flow, we determine the critical velocity for weak and moderate disorder. We also calculate the superfluid density and the condensate density as a function of the disorder strength and the flow velocity, and their behaviors near the critical velocity are discussed.
Modeling groundwater flow on MPPs
Ashby, S.F.; Falgout, R.D.; Smith, S.G.; Tompson, A.F.B.
1993-10-01
The numerical simulation of groundwater flow in three-dimensional heterogeneous porous media is examined. To enable detailed modeling of large contaminated sites, preconditioned iterative methods and massively parallel computing power are combined in a simulator called PARFLOW. After describing this portable and modular code, some numerical results are given, including one that demonstrates the code`s scalability.
Turbulence modeling for compressible flows
NASA Technical Reports Server (NTRS)
Marvin, J. G.
1977-01-01
Material prepared for a course on Applications and Fundamentals of Turbulence given at the University of Tennessee Space Institute, January 10 and 11, 1977, is presented. A complete concept of turbulence modeling is described, and examples of progess for its use in computational aerodynimics are given. Modeling concepts, experiments, and computations using the concepts are reviewed in a manner that provides an up-to-date statement on the status of this problem for compressible flows.
NASA Astrophysics Data System (ADS)
Dillard, Seth; Mousel, John; Buchholz, James; Udaykumar, H. S.
2009-11-01
A preliminary method has been developed to model complex moving boundaries interacting with fluids in two dimensions using video files. Image segmentation techniques are employed to generate sharp object interfaces which are cast as level sets embedded in a Cartesian flow domain. In this way, boundary evolution is effected directly through imagery rather than by way of functional approximation. Videos of an American eel swimming in a water tunnel apparatus and a guinea pig duodenum undergoing peristaltic contractions in vitro serve as external and internal flow examples, which are evaluated for wake structure and mixing efficacy, respectively.
Higgs inflation from standard model criticality
NASA Astrophysics Data System (ADS)
Hamada, Yuta; Kawai, Hikaru; Oda, Kin-ya; Park, Seong Chan
2015-03-01
The observed Higgs mass MH=125.9 ±0.4 GeV leads to the criticality of the standard model, that is, the Higgs potential becomes flat around the scale 1 017 - 18 GeV for the top mass 171.3 GeV. Earlier we proposed a Higgs inflation scenario in which this criticality plays a crucial role. In this paper, we investigate the detailed cosmological predictions of this scenario in light of the latest Planck and BICEP2 results. We also consider the Higgs portal scalar dark matter model, and compute the Higgs one-loop effective potential with the two-loop renormalization group improvement. We find a constraint on the coupling between the Higgs boson and dark matter which depends on the inflationary parameters.
A critical evaluation of numerical algorithms and flow physics in complex supersonic flows
NASA Astrophysics Data System (ADS)
Aradag, Selin
In this research, two different complex supersonic flows are selected to apply CFD to Navier-Stokes simulations. First test case is "Supersonic Flow over an Open Rectangular Cavity". Open cavity flow fields are remarkably complicated with internal and external regions that are coupled via self-sustained shear layer oscillations. Supersonic flow past a cavity has numerous applications in store carriage and release. Internal carriage of stores, which can be modeled using a cavity configuration, is used for supersonic aircraft in order to reduce radar cross section, aerodynamic drag and aerodynamic heating. Supersonic, turbulent, three-dimensional unsteady flow past an open rectangular cavity is simulated, to understand the physics and three-dimensional nature of the cavity flow oscillations. Influences of numerical parameters such as numerical flux scheme, computation time and flux limiter on the computed flow are determined. Two dimensional simulations are also performed for comparison purposes. The next test case is "The Computational Design of Boeing/AFOSR Mach 6 Wind Tunnel". Due to huge differences between geometrical scales, this problem is both challenging and computationally intensive. It is believed that most of the experimental data obtained from conventional ground testing facilities are not reliable due to high levels of noise associated with the acoustic fluctuations from the turbulent boundary layers on the wind tunnel walls. Therefore, it is very important to have quiet testing facilities for hypersonic flow research. The Boeing/AFOSR Mach 6 Wind tunnel in Purdue University has been designed as a quiet tunnel for which the noise level is an order of magnitude lower than that in conventional wind tunnels. However, quiet flow is achieved in the Purdue Mach 6 tunnel for only low Reynolds numbers. Early transition of the nozzle wall boundary layer has been identified as the cause of the test section noise. Separation bubbles on the bleed lip and associated
A holographic model for quantum critical responses
NASA Astrophysics Data System (ADS)
Myers, Robert C.; Sierens, Todd; Witczak-Krempa, William
2016-05-01
We analyze the dynamical response functions of strongly interacting quantum critical states described by conformal field theories (CFTs). We construct a self-consistent holographic model that incorporates the relevant scalar operator driving the quantum critical phase transition. Focusing on the finite temperature dynamical conductivity σ( ω, T ), we study its dependence on our model parameters, notably the scaling dimension of the relevant operator. It is found that the conductivity is well-approximated by a simple ansatz proposed in [1] for a wide range of parameters. We further dissect the conductivity at large frequencies ω ≫ T using the operator product expansion, and show how it reveals the spectrum of our model CFT. Our results provide a physically-constrained framework to study the analytic continuation of quantum Monte Carlo data, as we illustrate using the O(2) Wilson-Fisher CFT. Finally, we comment on the variation of the conductivity as we tune away from the quantum critical point, setting the stage for a comprehensive analysis of the phase diagram near the transition.
Critical Surface of the Hexagonal Polygon Model
NASA Astrophysics Data System (ADS)
Grimmett, Geoffrey R.; Li, Zhongyang
2016-05-01
The hexagonal polygon model arises in a natural way via a transformation of the 1-2 model on the hexagonal lattice, and it is related to the high temperature expansion of the Ising model. There are three types of edge, and three corresponding parameters α ,β ,γ >0. By studying the long-range order of a certain two-edge correlation function, it is shown that the parameter space (0,∞)^3 may be divided into subcritical and supercritical regions, separated by critical surfaces satisfying an explicitly known formula. This result complements earlier work on the Ising model and the 1-2 model. The proof uses the Pfaffian representation of Fisher, Kasteleyn, and Temperley for the counts of dimers on planar graphs.
Flows In Model Human Femoral Arteries
NASA Technical Reports Server (NTRS)
Back, Lloyd H.; Kwack, Eug Y.; Crawford, Donald W.
1990-01-01
Flow is visualized with dye traces, and pressure measurements made. Report describes experimental study of flow in models of human femoral artery. Conducted to examine effect of slight curvature of artery on flow paths and distribution of pressure.
Modeling shrouded stator cavity flows in axial-flow compressors
Wellborn, S.R.; Tolchinsky, I.; Okiishi, T.H.
2000-01-01
Experiments and computational analyses were completed to understand the nature of shrouded stator cavity flows. From this understanding, a one-dimensional model of the flow through shrouded stator cavities was developed. This model estimates the leakage mass flow, temperature rise, and angular momentum increase through the cavity, given geometry parameters and the flow conditions at the interface between the cavity and primary flow path. This cavity model consists of two components, one that estimates the flow characteristics through the labyrinth seals and the other that predicts the transfer of momentum due to windage. A description of the one-dimensional model is given. The incorporation and use of the one-dimensional model in a multistage compressor primary flow analysis tool is described. The combination of this model and the primary flow solver was used to reliably simulate the significant impact on performance of the increase of hub seal leakage in a twelve-stage axial-flow compressor. Observed higher temperatures of the hub region fluid, different stage matching, and lower overall efficiencies and core flow than expected could be correctly linked to increased hub seal clearance with this new technique. The importance of including these leakage flows in compressor simulations is shown.
Magnetic resonance measurement of fluid dynamics and transport in tube flow of a near-critical fluid
NASA Astrophysics Data System (ADS)
Bray, Joshua M.; Rassi, Erik M.; Seymour, Joseph D.; Codd, Sarah L.
2014-07-01
An ability to predict fluid dynamics and transport in supercritical fluids is essential for optimization of applications such as carbon sequestration, enhanced oil recovery, "green" solvents, and supercritical coolant systems. While much has been done to model supercritical velocity distributions, experimental characterization is sparse, owing in part to a high sensitivity to perturbation by measurement probes. Magnetic resonance (MR) techniques, however, detect signal noninvasively from the fluid molecules and thereby overcome this obstacle to measurement. MR velocity maps and propagators (i.e., probability density functions of displacement) were acquired of a flowing fluid in several regimes about the critical point, providing quantitative data on the transport and fluid dynamics in the system. Hexafluoroethane (C2F6) was pumped at 0.5 ml/min in a cylindrical tube through an MR system, and propagators as well as velocity maps were measured at temperatures and pressures below, near, and above the critical values. It was observed that flow of C2F6 with thermodynamic properties far above or below the critical point had the Poiseuille flow distribution of an incompressible Newtonian fluid. Flows with thermodynamic properties near the critical point exhibit complex flow distributions impacted by buoyancy and viscous forces. The approach to steady state was also observed and found to take the longest near the critical point, but once it was reached, the dynamics were stable and reproducible. These data provide insight into the interplay between the critical phase transition thermodynamics and the fluid dynamics, which control transport processes.
Stochastic models for turbulent reacting flows
Kerstein, A.
1993-12-01
The goal of this program is to develop and apply stochastic models of various processes occurring within turbulent reacting flows in order to identify the fundamental mechanisms governing these flows, to support experimental studies of these flows, and to further the development of comprehensive turbulent reacting flow models.
QCD critical point in a quasiparticle model
Srivastava, P. K.; Tiwari, S. K.; Singh, C. P.
2010-07-01
Recent theoretical investigations have unveiled a rich structure in the quantum chromodynamics phase diagram, which consists of quark-gluon plasma and the hadronic phases but also supports the existence of a crossover transition ending at a critical end point (CEP). We find a too large variation in the determination of the coordinates of the CEP in the temperature (T) baryon chemical potential ({mu}{sub B}) plane; and, therefore, its identification in the current heavy-ion experiments becomes debatable. Here we use an equation of state for a deconfined quark-gluon plasma using a thermodynamically-consistent quasiparticle model involving noninteracting quarks and gluons having thermal masses. We further use a thermodynamically-consistent excluded-volume model for the hadron gas, which was recently proposed by us. Using these equations of state, a first-order deconfining phase transition is constructed using Gibbs's criteria. This leads to an interesting finding that the phase transition line ends at a critical end point (CEP) beyond which a crossover region exists. Using our thermal hadron gas model, we obtain a chemical freeze out curve, and we find that the CEP lies in close proximity to this curve as proposed by some authors. The coordinates of CEP are found to lie within the reach of Relativistic heavy-ion collider experiment.
Catchment organisation, free energy dynamics and network control on critical zone water flows
NASA Astrophysics Data System (ADS)
Zehe, E.; Ehret, U.; Kleidon, A.; Jackisch, C.; Scherer, U.; Blume, T.
2012-04-01
as that these flow structures organize and dominate flows of water, dissolved matter and sediments during rainfall driven conditions at various scales: - Surface connected vertical flow structures of anecic worm burrows or soil cracks organize and dominated vertical flows at the plot scale - this is usually referred to as preferential flow; - Rill networks at the soil surface organise and dominate hillslope scale overland flow response and sediment yields; - Subsurface pipe networks at the bedrock interface organize and dominate hillslope scale lateral subsurface water and tracer flows; - The river net organizes and dominates flows of water, dissolved matter and sediments to the catchment outlet and finally across continental gradients to the sea. Fundamental progress with respect to the parameterization of hydrological models, subscale flow networks and to understand the adaptation of hydro-geo ecosystems to change could be achieved by discovering principles that govern the organization of catchments flow networks in particular at least during steady state conditions. This insight has inspired various scientists to suggest principles for organization of ecosystems, landscapes and flow networks; as Bejans constructural law, Minimum Energy Expenditure , Maximum Entropy Production. In line with these studies we suggest that a thermodynamic/energetic treatment of the catchment is might be a key for understanding the underlying principles that govern organisation of flow and transport. Our approach is to employ a) physically based hydrological model that address at least all the relevant hydrological processes in the critical zone in a coupled way, behavioural representations of the observed organisation of flow structures and textural elements, that are consistent with observations in two well investigated research catchments and have been tested against distributed observations of soil moisture and catchment scale discharge; to simulate the full concert of hydrological
Criticality in a model of banking crises
NASA Astrophysics Data System (ADS)
Iori, Giulia; Jafarey, Saqib
2001-10-01
An interbank market lets participants pool the risk arising from the combination of illiquid investments and random withdrawals by depositors. But it also creates the potential for one bank's failure to trigger off avalanches of further failures. We simulate a model of interbank lending to study the interplay of these two effects. We show that when banks are similar in size and exposure to risk, avalanche effects are small so that widening the interbank market leads to more stability. But as heterogeneity increases, avalanche effects become more important. By varying the heterogeneity and connectivity across banks, the system enters a critical regime with a power law distribution of avalanche sizes.
Global scale groundwater flow model
NASA Astrophysics Data System (ADS)
Sutanudjaja, Edwin; de Graaf, Inge; van Beek, Ludovicus; Bierkens, Marc
2013-04-01
As the world's largest accessible source of freshwater, groundwater plays vital role in satisfying the basic needs of human society. It serves as a primary source of drinking water and supplies water for agricultural and industrial activities. During times of drought, groundwater sustains water flows in streams, rivers, lakes and wetlands, and thus supports ecosystem habitat and biodiversity, while its large natural storage provides a buffer against water shortages. Yet, the current generation of global scale hydrological models does not include a groundwater flow component that is a crucial part of the hydrological cycle and allows the simulation of groundwater head dynamics. In this study we present a steady-state MODFLOW (McDonald and Harbaugh, 1988) groundwater model on the global scale at 5 arc-minutes resolution. Aquifer schematization and properties of this groundwater model were developed from available global lithological model (e.g. Dürr et al., 2005; Gleeson et al., 2010; Hartmann and Moorsdorff, in press). We force the groundwtaer model with the output from the large-scale hydrological model PCR-GLOBWB (van Beek et al., 2011), specifically the long term net groundwater recharge and average surface water levels derived from routed channel discharge. We validated calculated groundwater heads and depths with available head observations, from different regions, including the North and South America and Western Europe. Our results show that it is feasible to build a relatively simple global scale groundwater model using existing information, and estimate water table depths within acceptable accuracy in many parts of the world.
Modeling of turbulent separated flows for aerodynamic applications
NASA Technical Reports Server (NTRS)
Marvin, J. G.
1983-01-01
Steady, high speed, compressible separated flows modeled through numerical simulations resulting from solutions of the mass-averaged Navier-Stokes equations are reviewed. Emphasis is placed on benchmark flows that represent simplified (but realistic) aerodynamic phenomena. These include impinging shock waves, compression corners, glancing shock waves, trailing edge regions, and supersonic high angle of attack flows. A critical assessment of modeling capabilities is provided by comparing the numerical simulations with experiment. The importance of combining experiment, numerical algorithm, grid, and turbulence model to effectively develop this potentially powerful simulation technique is stressed.
A Geophysical Flow Experiment in a Compressible Critical Fluid
NASA Technical Reports Server (NTRS)
Hegseth, John; Garcia, Laudelino
1996-01-01
The first objective of this experiment is to build an experimental system in which, in analogy to a geophysical system, a compressible fluid in a spherical annulus becomes radially stratified in density through an A.C. electric field. When this density gradient is demonstrated, the system will be augmented so that the fluid can be driven by heating and rotation and tested in preparation for a microgravity experiment. This apparatus consists of a spherical capacitor filled with critical fluid in a temperature controlled environment. To make the fluid critical, the apparatus will be operated near the critical pressure, critical density, and critical temperature of the fluid. This will result in a highly compressible fluid because of the properties of the fluid near its critical point. A high voltage A.C. source applied across the capacitor will create a spherically symmetric central force because of the dielectric properties of the fluid in an electric field gradient. This central force will induce a spherically symmetric density gradient that is analogous to a geophysical fluid system. To generate such a density gradient the system must be small (approx. 1 inch diameter). This small cell will also be capable of driving the critical fluid by heating and rotation. Since a spherically symmetric density gradient can only be made in microgravity, another small cell, of the same geometry, will be built that uses incompressible fluid. The driving of the fluid by rotation and heating in these small cells will be developed. The resulting instabilities from the driving in these two systems will then be studied. The second objective is to study the pattern forming instabilities (bifurcations) resulting from the well controlled experimental conditions in the critical fluid cell. This experiment will come close to producing conditions that are geophysically similar and will be studied as the driving parameters are changed.
Particle and flow field holography: A critical survey
NASA Technical Reports Server (NTRS)
Trolinger, James D.
1987-01-01
A brief background is provided for the fields of particle and flow visualization holography. A summary of methods currently in use is given, followed by a discussion of more recent and unique applications. The problem of data reduction is discussed. A state of the art summary is then provided with a prognosis of the future of the field. Particle and flow visualization holography are characterized as powerful tools currently in wide use and with significant untapped potential.
Chen, Lin; Zhao, Yan; Zhang, Xin-Rong
2015-04-01
Near-critical/supercritical fluids have been widely proposed in material process, energy conversion and chemical engineering, etc. The present study is focused on the near-critical CO2 Poiseuille Rayleigh-Benard convective flow in microchannels. Careful numerical procedures are carried out by compressible Navier-Stokes equations, coupled energy and near-critical CO2 fluid state equations. In the physical model, sudden application of boundary heat fluxes in the boundaries is assumed. The flow and heat transfer characteristics of such Poiseuille Rayleigh-Benard configuration in microscales are systematically explored. For the convection onset, strong near-critical vortex flows are found for a relative wide range of initial and input conditions in microchannels. It is found that typical near-critical thin, hot boundary layer (HBL) plays critical role in the basic stability evolution process. The hot boundary layer formation process and the characteristics of the transition phenomena, convection structure, heat transfer behaviors as well as future development are also presented in this paper. It is hoped that this study can contribute to near-critical hydrodynamics in microscales. PMID:26353532
Ma, Rui; Zheng, Chunmiao; Tonkin, Matt; Zachara, John M
2011-04-01
Correct interpretation of tracer test data is critical for understanding transport processes in the subsurface. This task can be greatly complicated by the presence of intraborehole flows in a highly dynamic flow environment. At a new tracer test site (Hanford IFRC) a dynamic flow field created by changes in the stage of the adjacent Columbia River, coupled with a heterogeneous hydraulic conductivity distribution, leads to considerable variations in vertical hydraulic gradients. These variations, in turn, create intraborehole flows in fully-screened (6.5m) observation wells with frequently alternating upward and downward movement. This phenomenon, in conjunction with a highly permeable aquifer formation and small horizontal hydraulic gradients, makes modeling analysis and model calibration a formidable challenge. Groundwater head data alone were insufficient to define the flow model boundary conditions, and the movement of the tracer was highly sensitive to the dynamics of the flow field. This study shows that model calibration can be significantly improved by explicitly considering (a) dynamic flow model boundary conditions and (b) intraborehole flow. The findings from this study underscore the difficulties in interpreting tracer tests and understanding solute transport under highly dynamic flow conditions. PMID:21216023
Velocity and temperature profiles in near-critical nitrogen flowing past a horizontal flat plate
NASA Technical Reports Server (NTRS)
Simoneau, R. J.
1977-01-01
Boundary layer velocity and temperature profiles were measured for nitrogen near its thermodynamic critical point flowing past a horizontal flat plate. The results were compared measurements made for vertically upward flow. The boundary layer temperatures ranged from below to above the thermodynamic critical temperature. For wall temperatures below the thermodynamic critical temperature there was little variation between the velocity and temperature profiles in three orientations. In all three orientations the point of crossing into the critical temperature region is marked by a significant flattening of the velocity and temperature profiles and also a decrease in heat transfer coefficient.
Post Saltation Sediment Transport by Shallow Super Critical Flows
Technology Transfer Automated Retrieval System (TEKTRAN)
Laboratory flume experiments of shallow overland flow show that the sediment transport is not a random phenomenon but occurs in a highly organized manner. The organization reveals evolution of bed covered by sediments with several scales ranging from saltation of particles at very low concentrations...
Turbulence modeling for complex hypersonic flows
NASA Technical Reports Server (NTRS)
Huang, P. G.; Coakley, T. J.
1993-01-01
The paper presents results of calculations for a range of 2D turbulent hypersonic flows using two-equation models. The baseline models and the model corrections required for good hypersonic-flow predictions will be illustrated. Three experimental data sets were chosen for comparison. They are: (1) the hypersonic flare flows of Kussoy and Horstman, (2) a 2D hypersonic compression corner flow of Coleman and Stollery, and (3) the ogive-cylinder impinging shock-expansion flows of Kussoy and Horstman. Comparisons with the experimental data have shown that baseline models under-predict the extent of flow separation but over-predict the heat transfer rate near flow reattachment. Modifications to the models are described which remove the above-mentioned deficiencies. Although we have restricted the discussion only to the selected baseline models in this paper, the modifications proposed are universal and can in principle be transferred to any existing two-equation model formulation.
Miniaturizing free-flow electrophoresis - a critical review.
Kohlheyer, Dietrich; Eijkel, Jan C T; van den Berg, Albert; Schasfoort, Richard B M
2008-03-01
Free-flow electrophoresis (FFE) separation methods have been developed and investigated for around 50 years and have been applied not only to many types of analytes for various biomedical applications, but also for the separation of inorganic and organic substances. Its continuous sample preparation and mild separation conditions make it also interesting for online monitoring and detection applications. Since 1994 several microfluidic, miniaturized FFE devices were developed and experimentally characterized. In contrast to their large-scale counterparts microfluidic FFE (mu-FFE) devices offer new possibilities due to the very rapid separations within several seconds or below and the requirement for sample volumes in the microliter range. Eventually, these mu-FFE systems might find application in so-called lab-on-a-chip devices for real-time monitoring and separation applications. This review gives detailed information on the results so far published on mu-FFE chips, comprising its four main modes, namely free-flow zone electrophoresis (FFZE), free-flow IEF (FFIEF), free-flow ITP (FFITP), and free-flow field-step electrophoresis (FFFSE). The principles of the different FFE modes and the basic underlying theory are given and discussed with special emphasis on miniaturization. Different designs as well as fabrication methods and applied materials are discussed and evaluated. Furthermore, the separation results shown indicate that similar separation quality with respect to conventional FFE systems, as defined by the resolution and peak capacity, can be achieved with mu-FFE separations when applying much lower electrical voltages. Furthermore, innovations still occur and several approaches for hyphenated, more integrated systems have been proposed so far, some of which are discussed here. This review is intended as an introduction and early compendium for research and development within this field. PMID:18232029
ERIC Educational Resources Information Center
Subramaniam, Maithreyi; Hanafi, Jaffri; Putih, Abu Talib
2016-01-01
This study adopted 30 first year graphic design students' artwork, with critical analysis using Feldman's model of art criticism. Data were analyzed quantitatively; descriptive statistical techniques were employed. The scores were viewed in the form of mean score and frequencies to determine students' performances in their critical ability.…
West Maui Groundwater Flow Model
Nicole Lautze
2015-01-01
Groundwater flow model for West Maui. 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 V – Island of Maui 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.
Hawaii Island Groundwater Flow Model
Nicole Lautze
2015-01-01
Groundwater flow model for Hawaii Island. Data is from the following sources: 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 II – Island of Hawaii 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. 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.
East Maui Groundwater Flow Model
Nicole Lautze
2015-01-01
Groundwater flow model for East Maui. 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 V – Island of Maui 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.
Development of a Water Based, Critical Flow, Non-Vapor Compression cooling Cycle
Hosni, Mohammad H.
2014-03-30
Expansion of a high-pressure liquid refrigerant through the use of a thermostatic expansion valve or other device is commonplace in vapor-compression cycles to regulate the quality and flow rate of the refrigerant entering the evaporator. In vapor-compression systems, as the condensed refrigerant undergoes this expansion, its pressure and temperature drop, and part of the liquid evaporates. We (researchers at Kansas State University) are developing a cooling cycle that instead pumps a high-pressure refrigerant through a supersonic converging-diverging nozzle. As the liquid refrigerant passes through the nozzle, its velocity reaches supersonic (or critical-flow) conditions, substantially decreasing the refrigerant’s pressure. This sharp pressure change vaporizes some of the refrigerant and absorbs heat from the surrounding conditions during this phase change. Due to the design of the nozzle, a shockwave trips the supersonic two-phase refrigerant back to the starting conditions, condensing the remaining vapor. The critical-flow refrigeration cycle would provide space cooling, similar to a chiller, by running a secondary fluid such as water or glycol over one or more nozzles. Rather than utilizing a compressor to raise the pressure of the refrigerant, as in a vapor-cycle system, the critical-flow cycle utilizes a high-pressure pump to drive refrigerant liquid through the cooling cycle. Additionally, the design of the nozzle can be tailored for a given refrigerant, such that environmentally benign substances can act as the working fluid. This refrigeration cycle is still in early-stage development with prototype development several years away. The complex multi-phase flow at supersonic conditions presents numerous challenges to fully understanding and modeling the cycle. With the support of DOE and venture-capital investors, initial research was conducted at PAX Streamline, and later, at Caitin. We (researchers at Kansas State University) have continued development
Characterization of Turbulent Flows for Turbulence Modeling
NASA Astrophysics Data System (ADS)
Reynolds, W. C.; Haire, S. L.
1998-11-01
A diagram for the characterization of turbulent flows using the invariants of the mean velocity gradient tensor is introduced. All mean flows, from irrotationally strained flows to shearing flows, to purely rotational flows, can be identified on this diagram. Different flow fields which occupy the same region on the diagram are said to be comprised of the same topological features. The current state of turbulence modeling can be identified on the diagram based on the type of mean flow fields which can be accurately computed. Regions on the diagram can be shown for which current capabilities in turbulence modeling fail to accurately resolve the turbulent structures. Relevant mean field topology is identified for future work in turbulence modeling. Using this analysis, we suggest a number of flows to be computed by DNS or LES and used as testing cases for new models.
A toy terrestrial carbon flow model
NASA Technical Reports Server (NTRS)
Parton, William J.; Running, Steven W.; Walker, Brian
1992-01-01
A generalized carbon flow model for the major terrestrial ecosystems of the world is reported. The model is a simplification of the Century model and the Forest-Biogeochemical model. Topics covered include plant production, decomposition and nutrient cycling, biomes, the utility of the carbon flow model for predicting carbon dynamics under global change, and possible applications to state-and-transition models and environmentally driven global vegetation models.
Observing and Modeling Earth's Energy Flows
NASA Astrophysics Data System (ADS)
Stevens, Bjorn; Schwartz, Stephen E.
2012-07-01
importantly to this adjustment and thus contribute both to uncertainty in estimates of radiative forcing and to uncertainty in the response. Models are indispensable to calculation of the adjustment of the system to a compositional change but are known to be flawed in their representation of clouds. Advances in tracking Earth's energy flows and compositional changes on daily through decadal timescales are shown to provide both a critical and constructive framework for advancing model development and evaluation.
Deconfined Quantum Criticality, Scaling Violations, and Classical Loop Models
NASA Astrophysics Data System (ADS)
Nahum, Adam; Chalker, J. T.; Serna, P.; Ortuño, M.; Somoza, A. M.
2015-10-01
Numerical studies of the transition between Néel and valence bond solid phases in two-dimensional quantum antiferromagnets give strong evidence for the remarkable scenario of deconfined criticality, but display strong violations of finite-size scaling that are not yet understood. We show how to realize the universal physics of the Néel-valence-bond-solid (VBS) transition in a three-dimensional classical loop model (this model includes the subtle interference effect that suppresses hedgehog defects in the Néel order parameter). We use the loop model for simulations of unprecedentedly large systems (up to linear size L =512 ). Our results are compatible with a continuous transition at which both Néel and VBS order parameters are critical, and we do not see conventional signs of first-order behavior. However, we show that the scaling violations are stronger than previously realized and are incompatible with conventional finite-size scaling, even if allowance is made for a weakly or marginally irrelevant scaling variable. In particular, different approaches to determining the anomalous dimensions ηVBS and ηN é el yield very different results. The assumption of conventional finite-size scaling leads to estimates that drift to negative values at large sizes, in violation of the unitarity bounds. In contrast, the decay with distance of critical correlators on scales much smaller than system size is consistent with large positive anomalous dimensions. Barring an unexpected reversal in behavior at still larger sizes, this implies that the transition, if continuous, must show unconventional finite-size scaling, for example, from an additional dangerously irrelevant scaling variable. Another possibility is an anomalously weak first-order transition. By analyzing the renormalization group flows for the noncompact CP n -1 field theory (the n -component Abelian Higgs model) between two and four dimensions, we give the simplest scenario by which an anomalously weak first
Calculation of the Minimum Critical Reynolds Number for Circular Pipe Flows
NASA Astrophysics Data System (ADS)
Kanda, Hidesada
1999-11-01
A conceptual model was constructed for the problem of determining in circular pipes the conditions under which the transition from laminar to turbulent flow occurs. From many previous experimental investigations, it became clear that (i) plots of the transition length versus the Reynolds number (Re) show that the transition occurs in the entrance region under the condition of a natural disturbance, and (ii) plots of the critical Re versus the ratio of bellmouth diameter to the pipe diameter show that for the case of a straight pipe the critical Re takes a minimum value of about 2000. In the entrance region, the velocity profile changes from uniform at the inlet to parabolic at the entrance length. We found that the radial component of the curl of vorticity multiplied by (2/Re), which we call the normal wall strength, works as an acceleration force and decreases inversely as Re increases. Hence, the onset of the transition should depend on whether or not the acceleration power provided by the normal wall strength exceeds a required value. In this study we calculated the acceleration power via finite difference calculations, and thus obtained the minimum critical Reynolds number of 2040 when using J0 = 101 radial grid points.
VISCOPLASTIC FLUID MODEL FOR DEBRIS FLOW ROUTING.
Chen, Cheng-lung
1986-01-01
This paper describes how a generalized viscoplastic fluid model, which was developed based on non-Newtonian fluid mechanics, can be successfully applied to routing a debris flow down a channel. The one-dimensional dynamic equations developed for unsteady clear-water flow can be used for debris flow routing if the flow parameters, such as the momentum (or energy) correction factor and the resistance coefficient, can be accurately evaluated. The writer's generalized viscoplastic fluid model can be used to express such flow parameters in terms of the rheological parameters for debris flow in wide channels. A preliminary analysis of the theoretical solutions reveals the importance of the flow behavior index and the so-called modified Froude number for uniformly progressive flow in snout profile modeling.
Models for water steam condensing flows
NASA Astrophysics Data System (ADS)
Wróblewski, Włodzimierz; Dykas, Sławomir; Chmielniak, Tadeusz
2012-08-01
The paper presents a description of selected models dedicated to steam condensing flow modelling. The models are implemented into an in-house computational fluid dynamics code that has been successfully applied to wet steam flow calculation for many years now. All models use the same condensation model that has been validated against the majority of available experimental data. The state equations for vapour and liquid water, the physical model as well as the numerical techniques of solution to flow governing equations have been presented. For the single-fluid model, the Reynolds-averaged Navier-Stokes equations for vapour/liquid mixture are solved, whereas the two-fluid model solves separate flow governing equations for the compressible, viscous and turbulent vapour phase and for the compressible and inviscid liquid phase. All described models have been compared with relation to the flow through the Laval nozzle.
Critical Thinking: Frameworks and Models for Teaching
ERIC Educational Resources Information Center
Fahim, Mansoor; Eslamdoost, Samaneh
2014-01-01
Developing critical thinking since the educational revolution gave rise to flourishing movements toward embedding critical thinking (CT henceforth) stimulating classroom activities in educational settings. Nevertheless the process faced with complications such as teachability potentiality, lack of practical frameworks concerning actualization of…
A Model for Teaching Critical Thinking
ERIC Educational Resources Information Center
Emerson, Marnice K.
2013-01-01
In an age in which information is available almost instantly and in quantities unimagined just a few decades ago, most educators would agree that teaching adult learners to think critically about what they are reading, seeing, and hearing has never been more important. But just what is critical thinking? Do adult learners agree with educators that…
Critical behavior of the Widom--Rowlinson lattice model
Dickman, R.; Stell, G.
1995-06-01
We report extensive Monte Carlo simulations of the Widom--Rowlinson lattice model in two and three dimensions. Our results yield precise values for the critical activities and densities, and clearly place the critical behavior in the Ising universality class.
Catchment organisation, free energy dynamics and network control on critical zone water flows
NASA Astrophysics Data System (ADS)
Zehe, E.; Ehret, U.; Kleidon, A.; Jackisch, C.; Scherer, U.; Blume, T.
2012-04-01
as that these flow structures organize and dominate flows of water, dissolved matter and sediments during rainfall driven conditions at various scales: - Surface connected vertical flow structures of anecic worm burrows or soil cracks organize and dominated vertical flows at the plot scale - this is usually referred to as preferential flow; - Rill networks at the soil surface organise and dominate hillslope scale overland flow response and sediment yields; - Subsurface pipe networks at the bedrock interface organize and dominate hillslope scale lateral subsurface water and tracer flows; - The river net organizes and dominates flows of water, dissolved matter and sediments to the catchment outlet and finally across continental gradients to the sea. Fundamental progress with respect to the parameterization of hydrological models, subscale flow networks and to understand the adaptation of hydro-geo ecosystems to change could be achieved by discovering principles that govern the organization of catchments flow networks in particular at least during steady state conditions. This insight has inspired various scientists to suggest principles for organization of ecosystems, landscapes and flow networks; as Bejans constructural law, Minimum Energy Expenditure , Maximum Entropy Production. In line with these studies we suggest that a thermodynamic/energetic treatment of the catchment is might be a key for understanding the underlying principles that govern organisation of flow and transport. Our approach is to employ a) physically based hydrological model that address at least all the relevant hydrological processes in the critical zone in a coupled way, behavioural representations of the observed organisation of flow structures and textural elements, that are consistent with observations in two well investigated research catchments and have been tested against distributed observations of soil moisture and catchment scale discharge; to simulate the full concert of hydrological
A stochastic index flow model of flow duration curves
NASA Astrophysics Data System (ADS)
Castellarin, Attilio; Vogel, Richard M.; Brath, Armando
2004-03-01
Annual flow duration curves (AFDCs) are used increasingly because unlike traditional period of record flow duration curves (FDCs), they provide confidence intervals for the median AFDC, they enable one to assign return periods to individual AFDCs, and they offer opportunities for developing a generalized stochastic model of daily streamflow. Previous stochastic models of FDCs and AFDCs were unable to reproduce the variance of AFDCs. We introduce an index flow approach to modeling the relationship between an FDC and AFDCs of daily streamflow series, which is able to reproduce the FDC, as well as the mean, median, and variance of the AFDCs without resorting to assumptions regarding the seasonal or persistence structure of daily streamflow series. Our approach offers additional opportunities for the regionalization of flow duration curves and for the generation of time series of daily streamflows at ungauged sites. Our approach is tested on three river basins in eastern central Italy.
Goh, Segun; Lee, Keumsook; Choi, Moo Young; Fortin, Jean-Yves
2014-01-01
Social systems have recently attracted much attention, with attempts to understand social behavior with the aid of statistical mechanics applied to complex systems. Collective properties of such systems emerge from couplings between components, for example, individual persons, transportation nodes such as airports or subway stations, and administrative districts. Among various collective properties, criticality is known as a characteristic property of a complex system, which helps the systems to respond flexibly to external perturbations. This work considers the criticality of the urban transportation system entailed in the massive smart card data on the Seoul transportation network. Analyzing the passenger flow on the Seoul bus system during one week, we find explicit power-law correlations in the system, that is, power-law behavior of the strength correlation function of bus stops and verify scale invariance of the strength fluctuations. Such criticality is probed by means of the scaling and renormalization analysis of the modified gravity model applied to the system. Here a group of nearby (bare) bus stops are transformed into a (renormalized) "block stop" and the scaling relations of the network density turn out to be closely related to the fractal dimensions of the system, revealing the underlying structure. Specifically, the resulting renormalized values of the gravity exponent and of the Hill coefficient give a good description of the Seoul bus system: The former measures the characteristic dimensionality of the network whereas the latter reflects the coupling between distinct transportation modes. It is thus demonstrated that such ideas of physics as scaling and renormalization can be applied successfully to social phenomena exemplified by the passenger flow. PMID:24599221
Goh, Segun; Lee, Keumsook; Choi, MooYoung; Fortin, Jean-Yves
2014-01-01
Social systems have recently attracted much attention, with attempts to understand social behavior with the aid of statistical mechanics applied to complex systems. Collective properties of such systems emerge from couplings between components, for example, individual persons, transportation nodes such as airports or subway stations, and administrative districts. Among various collective properties, criticality is known as a characteristic property of a complex system, which helps the systems to respond flexibly to external perturbations. This work considers the criticality of the urban transportation system entailed in the massive smart card data on the Seoul transportation network. Analyzing the passenger flow on the Seoul bus system during one week, we find explicit power-law correlations in the system, that is, power-law behavior of the strength correlation function of bus stops and verify scale invariance of the strength fluctuations. Such criticality is probed by means of the scaling and renormalization analysis of the modified gravity model applied to the system. Here a group of nearby (bare) bus stops are transformed into a (renormalized) “block stop” and the scaling relations of the network density turn out to be closely related to the fractal dimensions of the system, revealing the underlying structure. Specifically, the resulting renormalized values of the gravity exponent and of the Hill coefficient give a good description of the Seoul bus system: The former measures the characteristic dimensionality of the network whereas the latter reflects the coupling between distinct transportation modes. It is thus demonstrated that such ideas of physics as scaling and renormalization can be applied successfully to social phenomena exemplified by the passenger flow. PMID:24599221
Site-Scale Saturated Zone Flow Model
G. Zyvoloski
2003-12-17
The purpose of this model report is to document the components of the site-scale saturated-zone flow model at Yucca Mountain, Nevada, in accordance with administrative procedure (AP)-SIII.lOQ, ''Models''. This report provides validation and confidence in the flow model that was developed for site recommendation (SR) and will be used to provide flow fields in support of the Total Systems Performance Assessment (TSPA) for the License Application. The output from this report provides the flow model used in the ''Site-Scale Saturated Zone Transport'', MDL-NBS-HS-000010 Rev 01 (BSC 2003 [162419]). The Site-Scale Saturated Zone Transport model then provides output to the SZ Transport Abstraction Model (BSC 2003 [164870]). In particular, the output from the SZ site-scale flow model is used to simulate the groundwater flow pathways and radionuclide transport to the accessible environment for use in the TSPA calculations. Since the development and calibration of the saturated-zone flow model, more data have been gathered for use in model validation and confidence building, including new water-level data from Nye County wells, single- and multiple-well hydraulic testing data, and new hydrochemistry data. In addition, a new hydrogeologic framework model (HFM), which incorporates Nye County wells lithology, also provides geologic data for corroboration and confidence in the flow model. The intended use of this work is to provide a flow model that generates flow fields to simulate radionuclide transport in saturated porous rock and alluvium under natural or forced gradient flow conditions. The flow model simulations are completed using the three-dimensional (3-D), finite-element, flow, heat, and transport computer code, FEHM Version (V) 2.20 (software tracking number (STN): 10086-2.20-00; LANL 2003 [161725]). Concurrently, process-level transport model and methodology for calculating radionuclide transport in the saturated zone at Yucca Mountain using FEHM V 2.20 are being
Investigation of Body Force Effects on Flow Boiling Critical Heat Flux
NASA Technical Reports Server (NTRS)
Zhang, Hui; Mudawar, Issam; Hasan, Mohammad M.
2002-01-01
The bubble coalescence and interfacial instabilities that are important to modeling critical heat flux (CHF) in reduced-gravity systems can be sensitive to even minute body forces. Understanding these complex phenomena is vital to the design and safe implementation of two-phase thermal management loops proposed for space and planetary-based thermal systems. While reduced gravity conditions cannot be accurately simulated in 1g ground-based experiments, such experiments can help isolate the effects of the various forces (body force, surface tension force and inertia) which influence flow boiling CHF. In this project, the effects of the component of body force perpendicular to a heated wall were examined by conducting 1g flow boiling experiments at different orientations. FC-72 liquid was boiled along one wall of a transparent rectangular flow channel that permitted photographic study of the vapor-liquid interface at conditions approaching CHF. High-speed video imaging was employed to capture dominant CHF mechanisms. Six different CHF regimes were identified: Wavy Vapor Layer, Pool Boiling, Stratification, Vapor Counterflow, Vapor Stagnation, and Separated Concurrent Vapor Flow. CHF showed great sensitivity to orientation for flow velocities below 0.2 m/s, where very small CHF values where measured, especially with downflow and downward-facing heated wall orientations. High flow velocities dampened the effects of orientation considerably. Figure I shows representative images for the different CHF regimes. The Wavy Vapor Layer regime was dominant for all high velocities and most orientations, while all other regimes were encountered at low velocities, in the downflow and/or downward-facing heated wall orientations. The Interfacial Lift-off model was modified to predict the effects of orientation on CHF for the dominant Wavy Vapor Layer regime. The photographic study captured a fairly continuous wavy vapor layer travelling along the heated wall while permitting liquid
Technology Transfer Automated Retrieval System (TEKTRAN)
A flow pattern is characterized by aquifer features and the number, type, and distribution of stagnation points (locations where the discharge is zero). This article identifies a condition for transition of flow patterns in two-dimensional groundwater flow obeying Darcy's law by examining changes in...
NASA Astrophysics Data System (ADS)
Liu, Richeng; Li, Bo; Jiang, Yujing
2016-02-01
Transition of fluid flow from the linear to the nonlinear regime has been confirmed in single rock fractures when the Reynolds number (Re) exceeds some critical values, yet the criterion for such a transition in discrete fracture networks (DFNs) has received little attention. This study conducted flow tests on crossed fracture models with a single intersection and performed numerical simulations on fluid flow through DFNs of various geometric characteristics. The roles of aperture, surface roughness, and number of intersections of fractures on the variation of the critical hydraulic gradient (Jc) for the onset of nonlinear flow through DFNs were systematically investigated. The results showed that the relationship between hydraulic gradient (J) and flow rate can be well quantified by Forchheimer's law; when J drops below Jc, it reduces to the widely used cubic law, by diminishing the nonlinear term. Larger apertures, rougher fracture surfaces, and a greater number of intersections in a DFN would result in the onset of nonlinear flow at a lower Jc. Mathematical expressions of Jc and the coefficients involved in Forchheimer's law were developed based on multi-variable regressions of simulation results, which can help to choose proper governing equations when solving problems associated with fluid flow in fracture networks.
Nonextensive critical effects in the NJL model
NASA Astrophysics Data System (ADS)
Rozynek, Jacek
2015-05-01
It is shown how nonextensive Nambu-Jona-Lasinio (NJL) calculations in the critical region in the vicinity of phase transition to quark matter allow to determine a nonextensive parameter q from a specific heat capacity Cq.
Model of Transition from Laminar to Turbulent Flow
NASA Astrophysics Data System (ADS)
Kanda, Hidesada
2001-11-01
For circular pipe flows, a model of transition from laminar to turbulent flow has already been proposed and the minimum critical Reynolds number of approximately 2040 was obtained (Kanda, 1999). In order to prove the validity of the model, another verification is required. Thus, for plane Poiseuille flow, results of previous investigations were studied, focusing on experimental data on the critical Reynolds number Rc, the entrance length, and the transition length. Consequently, concerning the natural transition, it was confirmed from the experimental data that (i) the transition occurs in the entrance region, (ii) Rc increases as the contraction ratio in the inlet section increases, and (iii) the minimum Rc is obtained when the contraction ratio is the smallest or one, and there is no-bellshaped entrance or straight parallel plates. Its value exists in the neighborhood of 1300, based on the channel height and the average velocity. Although, for Hagen-Poiseuille flow, the minimum Rc is approximately 2000, based on the pipe diameter and the average velocity, there seems to be no significant difference in the transition from laminar to turbulent flow between Hagen-Poiseuille flow and plane Poiseuille flow (Kanda, 2001). Rc is determined by the shape of the inlet. Kanda, H., 1999, Proc. of ASME Fluids Engineering Division - 1999, FED-Vol. 250, pp. 197-204. Kanda, H., 2001, Proc. of ASME Fluids Engineering Division - 2001.
Numerical modeling of fluidic flow meters
NASA Astrophysics Data System (ADS)
Choudhury, D.; Patel, B. R.
1992-05-01
The transient fluid flow in fluidic flow meters has been modeled using Creare.x's flow modeling computer program FLUENT/BFC that solves the Navier-Stokes equations in general curvilinear coordinates. The numerical predictions of fluid flow in a fluidic flow meter have been compared with the available experimental results for a particular design, termed the PC-4 design. Overall flow structures such as main jet bending, and primary and secondary vortices predicted by FLUENT/BFC are in excellent agreement with flow visualization results. The oscillation frequencies of the PC-4 design have been predicted for a range of flow rates encompassing laminar and turbulent flow and the results are in good agreement with experiments. The details of the flow field predictions reveal that an important factor that determines the onset of oscillations in the fluidic flow meter is the feedback jet momentum relative to the main jet momentum. The insights provided by the analysis of the PC-4 fluidic flow meter design have led to an improved design. The improved design has sustained oscillations at lower flow rates compared with the PC-4 design and has a larger rangeability.
Performance testing of the Silo Flow Model
Stadler, S.P.; O`Connor, D.; Gould, A.F.
1994-12-31
Several instruments are commercially available for on-line analysis of coal properties such as total moisture, ash, sulfur, and mineral matter content. These instruments have found use in coal cleaning and coal-fired utility applications. However, in many instances, the coal is stored in large bunkers or silos after on-line analysis, making the data gathered from on-line analysis a poor predictor of short-term coal quality due to the flow pattern and mixing within the silo. A computerized model, the Silo Flow Model, has been developed to model the flow of coal through a silo or bunker thus providing a prediction of the output coal quality based on on-line measurements of the quality of coal entering the silo. A test procedure was developed and demonstrated to test the performance of the Silo Flow Model. The testing was performed using controlled addition of silver nitrate to the coal, in conjunction with surface profile measurements using an array of ultrasonic gauges and data acquired from plant instrumentation. Results obtained from initial testing provided estimates of flow-related parameters used in the Silo flow Model. Similar test techniques are also used to compare predicted and actual silver content at the silo outlet as a measure of model performance. This paper describes test procedures used to validate the Silo Flow Model, the testing program, and the results obtained to data. The Silo Flow Model performance is discussed and compared against other modeling approaches.
Experimental Flow Models for SSME Flowfield Characterization
NASA Technical Reports Server (NTRS)
Abel, L. C.; Ramsey, P. E.
1989-01-01
Full scale flow models with extensive instrumentation were designed and manufactured to provide data necessary for flow field characterization in rocket engines of the Space Shuttle Main Engine (SSME) type. These models include accurate flow path geometries from the pre-burner outlet through the throat of the main combustion chamber. The turbines are simulated with static models designed to provide the correct pressure drop and swirl for specific power levels. The correct turbopump-hot gas manifold interfaces were designed into the flow models to permit parametric/integration studies for new turbine designs. These experimental flow models provide a vehicle for understanding the fluid dynamics associated with specific engine issues and also fill the more general need for establishing a more detailed fluid dynamic base to support development and verification of advanced math models.
Finite element modeling of nonisothermal polymer flows
NASA Technical Reports Server (NTRS)
Roylance, D.
1981-01-01
A finite element formulation designed to simulate polymer melt flows in which both conductive and convective heat transfer are important is described, and the numerical model is illustrated by means of computer experiments using extruder drag flow and entry flow as trial problems. Fluid incompressibility is enforced by a penalty treatment of the element pressures, and the thermal convective transport is modeled by conventional Galerkin and optimal upwind treatments.
Pursuing parameters for critical-density dark matter models
NASA Astrophysics Data System (ADS)
Liddle, Andrew R.; Lyth, David H.; Schaefer, R. K.; Shafi, Q.; Viana, Pedro T. P.
1996-07-01
We present an extensive comparison of models of structure formation with observations, based on linear and quasi-linear theory. We assume a critical matter density, and study both cold dark matter models and cold plus hot dark matter models. We explore a wide range of parameters, by varying the fraction of hot dark matter , the Hubble parameter h and the spectral index of density perturbations n, and allowing for the possibility of gravitational waves from inflation influencing large-angle microwave background anisotropies. New calculations are made of the transfer functions describing the linear power spectrum, with special emphasis on improving the accuracy on short scales where there are strong constraints. For assessing early object formation, the transfer functions are explicitly evaluated at the appropriate redshift. The observations considered are the four-year COBE observations of microwave background anisotropies, peculiar velocity flows, the galaxy correlation function, and the abundances of galaxy clusters, quasars and damped Lyman alpha systems. Each observation is interpreted in terms of the power spectrum filtered by a top-hat window function. We find that there remains a viable region of parameter space for critical-density models when all the dark matter is cold, though h must be less than 0.5 before any fit is found and n significantly below unity is preferred. Once a hot dark matter component is invoked, a wide parameter space is acceptable, including n 1. The allowed region is characterized by 0.35 and 0.60 n 1.25, at 95 per cent confidence on at least one piece of data. There is no useful lower bound on h, and for curious combinations of the other parameters it is possible to fit the data with h as high as 0.65.
Critical behavior of the Ising model on random fractals
NASA Astrophysics Data System (ADS)
Monceau, Pascal
2011-11-01
We study the critical behavior of the Ising model in the case of quenched disorder constrained by fractality on random Sierpinski fractals with a Hausdorff dimension df≃1.8928. This is a first attempt to study a situation between the borderline cases of deterministic self-similarity and quenched randomness. Intensive Monte Carlo simulations were carried out. Scaling corrections are much weaker than in the deterministic cases, so that our results enable us to ensure that finite-size scaling holds, and that the critical behavior is described by a new universality class. The hyperscaling relation is compatible with an effective dimension equal to the Hausdorff one; moreover the two eigenvalues exponents of the renormalization flows are shown to be different from the ones calculated from ɛ expansions, and from the ones obtained for fourfold symmetric deterministic fractals. Although the space dimensionality is not integer, lack of self-averaging properties exhibits some features very close to the ones of a random fixed point associated with a relevant disorder.
Approximate Model for Turbulent Stagnation Point Flow.
Dechant, Lawrence
2016-01-01
Here we derive an approximate turbulent self-similar model for a class of favorable pressure gradient wedge-like flows, focusing on the stagnation point limit. While the self-similar model provides a useful gross flow field estimate this approach must be combined with a near wall model is to determine skin friction and by Reynolds analogy the heat transfer coefficient. The combined approach is developed in detail for the stagnation point flow problem where turbulent skin friction and Nusselt number results are obtained. Comparison to the classical Van Driest (1958) result suggests overall reasonable agreement. Though the model is only valid near the stagnation region of cylinders and spheres it nonetheless provides a reasonable model for overall cylinder and sphere heat transfer. The enhancement effect of free stream turbulence upon the laminar flow is used to derive a similar expression which is valid for turbulent flow. Examination of free stream enhanced laminar flow suggests that the rather than enhancement of a laminar flow behavior free stream disturbance results in early transition to turbulent stagnation point behavior. Excellent agreement is shown between enhanced laminar flow and turbulent flow behavior for high levels, e.g. 5% of free stream turbulence. Finally the blunt body turbulent stagnation results are shown to provide realistic heat transfer results for turbulent jet impingement problems.
Kinetic model for dilute traffic flow
NASA Astrophysics Data System (ADS)
Balouchi, Ashkan; Browne, Dana A.
The flow of traffic represents a many-particle non-equilibrium problem with important practical consequences. Traffic behavior has been studied using a variety of approaches, including fluid dynamics models, Boltzmann equation, and recently cellular automata (CA). The CA model for traffic flow that Nagel and Schreckenberg (NS) introduced can successfully mimic many of the known features of the traffic flow. We show that in the dilute limit of the NS model, where vehicles exhibit free flow, cars show significant nearest neighbor correlation primarily via a short-range repulsion. introduce an approximate analytic model to describe this dilute limit. We show that the distribution of the distance between consecutive vehicles obeys a drift-diffusion equation. We compared this model with direct simulations. The steady state solution and relaxation of this model agrees well with direct simulations. We explore how this model breaks down as the transition to jams occurs.
Modeling of Time Varying Slag Flow in Coal Gasifiers
Pilli, Siva Prasad; Johnson, Kenneth I.; Williford, Ralph E.; Sundaram, S. K.; Korolev, Vladimir N.; Crum, Jarrod V.
2008-08-30
There is considerable interest within government agencies and the energy industries across the globe to further advance the clean and economical conversion of coal into liquid fuels to reduce our dependency on imported oil. To date, advances in these areas have been largely based on experimental work. Although there are some detailed systems level performance models, little work has been done on numerical modeling of the component level processes. If accurate models are developed, then significant R&D time might be saved, new insights into the process might be gained, and some good predictions of process or performance can be made. One such area is the characterization of slag deposition and flow on the gasifier walls. Understanding slag rheology and slag-refractory interactions is critical to design and operation of gasifiers with extended refractory lifetimes and also to better control of operating parameters so that the overall gasifier performance with extended service life can be optimized. In the present work, the literature on slag flow modeling was reviewed and a model similar to Seggiani’s was developed to simulate the time varying slag accumulation and flow on the walls of a Prenflo coal gasifier. This model was further extended and modified to simulate a refractory wall gasifier including heat transfer through the refractory wall with flowing slag in contact with the refractory. The model was used to simulate temperature dependent slag flow using rheology data from our experimental slag testing program. These modeling results as well as experimental validation are presented.
Mathematical and computational models of plasma flows
NASA Astrophysics Data System (ADS)
Brushlinsky, K. V.
Investigations of plasma flows are of interest, firstly, due to numerous applications, and secondly, because of their general principles, which form a special branch of physics: the plasma dynamics. Numerical simulation and computation, together with theoretic and experimental methods, play an important part in these investigations. Speaking on flows, a relatively dense plasma is mentioned, so its mathematical models appertain to the fluid mechanics, i.e., they are based on the magnetohydrodynamic description of plasma. Time dependent two dimensional models of plasma flows of two wide-spread types are considered: the flows across the magnetic field and those in the magnetic field plane.
Estimation of Critical Flow Velocity for Collapse of Gas Test Loop Booster Fuel Assembly
Guillen; Mark J. Russell
2006-07-01
This paper presents calculations performed to determine the critical flow velocity for plate collapse due to static instability for the Gas Test Loop booster fuel assembly. Long, slender plates arranged in a parallel configuration can experience static divergence and collapse at sufficiently high coolant flow rates. Such collapse was exhibited by the Oak Ridge High Flux Reactor in the 1940s and the Engineering Test Reactor at the Idaho National Laboratory in the 1950s. Theoretical formulas outlined by Miller, based upon wide-beam theory and Bernoulli’s equation, were used for the analysis. Calculations based upon Miller’s theory show that the actual coolant flow velocity is only 6% of the predicted critical flow velocity. Since there is a considerable margin between the theoretically predicted plate collapse velocity and the design velocity, the phenomena of plate collapse due to static instability is unlikely.
Analytical models for complex swirling flows
NASA Astrophysics Data System (ADS)
Borissov, A.; Hussain, V.
1996-11-01
We develops a new class of analytical solutions of the Navier-Stokes equations for swirling flows, and suggests ways to predict and control such flows occurring in various technological applications. We view momentum accumulation on the axis as a key feature of swirling flows and consider vortex-sink flows on curved axisymmetric surfaces with an axial flow. We show that these solutions model swirling flows in a cylindrical can, whirlpools, tornadoes, and cosmic swirling jets. The singularity of these solutions on the flow axis is removed by matching them with near-axis Schlichting and Long's swirling jets. The matched solutions model flows with very complex patterns, consisting of up to seven separation regions with recirculatory 'bubbles' and vortex rings. We apply the matched solutions for computing flows in the Ranque-Hilsch tube, in the meniscus of electrosprays, in vortex breakdown, and in an industrial vortex burner. The simple analytical solutions allow a clear understanding of how different control parameters affect the flow and guide selection of optimal parameter values for desired flow features. These solutions permit extension to other problems (such as heat transfer and chemical reaction) and have the potential of being significantly useful for further detailed investigation by direct or large-eddy numerical simulations as well as laboratory experimentation.
Seshasayanan, Kannabiran; Alexakis, Alexandros
2016-01-01
We investigate the critical transition from an inverse cascade of energy to a forward energy cascade in a two-dimensional magnetohydrodynamic flow as the ratio of magnetic to mechanical forcing amplitude is varied. It is found that the critical transition is the result of two competing processes. The first process is due to hydrodynamic interactions and cascades the energy to the large scales. The second process couples small-scale magnetic fields to large-scale flows, transferring the energy back to the small scales via a nonlocal mechanism. At marginality the two cascades are both present and cancel each other. The phase space diagram of the transition is sketched. PMID:26871152
SRMAFTE facility checkout model flow field analysis
NASA Technical Reports Server (NTRS)
Dill, Richard A.; Whitesides, Harold R.
1992-01-01
The Solid Rocket Motor Air Flow Equipment (SRMAFTE) facility was constructed for the purpose of evaluating the internal propellant, insulation, and nozzle configurations of solid propellant rocket motor designs. This makes the characterization of the facility internal flow field very important in assuring that no facility induced flow field features exist which would corrupt the model related measurements. In order to verify the design and operation of the facility, a three-dimensional computational flow field analysis was performed on the facility checkout model setup. The checkout model measurement data, one-dimensional and three-dimensional estimates were compared, and the design and proper operation of the facility was verified. The proper operation of the metering nozzles, adapter chamber transition, model nozzle, and diffuser were verified. The one-dimensional and three-dimensional flow field estimates along with the available measurement data are compared.
Regression modeling of ground-water flow
Cooley, R.L.; Naff, R.L.
1985-01-01
Nonlinear multiple regression methods are developed to model and analyze groundwater flow systems. Complete descriptions of regression methodology as applied to groundwater flow models allow scientists and engineers engaged in flow modeling to apply the methods to a wide range of problems. Organization of the text proceeds from an introduction that discusses the general topic of groundwater flow modeling, to a review of basic statistics necessary to properly apply regression techniques, and then to the main topic: exposition and use of linear and nonlinear regression to model groundwater flow. Statistical procedures are given to analyze and use the regression models. A number of exercises and answers are included to exercise the student on nearly all the methods that are presented for modeling and statistical analysis. Three computer programs implement the more complex methods. These three are a general two-dimensional, steady-state regression model for flow in an anisotropic, heterogeneous porous medium, a program to calculate a measure of model nonlinearity with respect to the regression parameters, and a program to analyze model errors in computed dependent variables such as hydraulic head. (USGS)
A Critical Information Literacy Model: Library Leadership within the Curriculum
ERIC Educational Resources Information Center
Swanson, Troy
2011-01-01
It is a time for a new model for teaching students to find, evaluate, and use information by drawing on critical pedagogy theory in the education literature. This critical information literacy model views the information world as a dynamic place where authors create knowledge for many reasons; it seeks to understand students as information users,…
Laminar-Turbulent Transition: A Hysteresis Curve of Two Critical Reynolds Numbers in Pipe Flow
NASA Astrophysics Data System (ADS)
Kanda, Hidesada
2006-11-01
A laminar-turbulent transition model (DFD 2004) has been constructed for pipe flows: (1) Natural transition occurs in the entrance region, and (2) Entrance shape determines a critical Reynolds number Rc. To verify the model, we have carried out experiments similar to Reynolds's color-dye experiment with 5 bellmouth entrances and a straight pipe. Then, we observed the following: (i) two different types of Rc exist, Rc1 from laminar to turbulent and Rc2 from turbulent to laminar, and (ii) the ratio of bellmouth diameter BD to pipe diameter D affects the values of Rc1 and Rc2. For each entrance, Rc1 has a maximum value Rc1(max) and Rc2 has a minimum value Rc2(min). When overlapping the two curves of Rc1(max) and Rc2(min) against BD/D, a hysteresis curve is confirmed. All Rc values exist inside this hysteresis curve. Consequently, Rc takes a minimum value Rc(min) of approximately 2000 when BD/D is at a minimum, i.e., at BD/D = 1, Rc(min) = Rc1(max) = Rc2(min) = 2000. Regarding Reynolds's Rc of 12,830, we observed Rc1(max) of approximately 13,000 at BD/D above 1.54. Therefore, the model has been partly verified.
Setting up the critical rainfall line for debris flows via support vector machines
NASA Astrophysics Data System (ADS)
Tsai, Y. F.; Chan, C. H.; Chang, C. H.
2015-10-01
The Chi-Chi earthquake in 1999 caused tremendous landslides which triggered many debris flows and resulted in significant loss of public lives and property. To prevent the disaster of debris flow, setting a critical rainfall line for each debris-flow stream is necessary. Firstly, 8 predisposing factors of debris flow were used to cluster 377 streams which have similar rainfall lines into 7 groups via the genetic algorithm. Then, support vector machines (SVM) were applied to setup the critical rainfall line for debris flows. SVM is a machine learning approach proposed based on statistical learning theory and has been widely used on pattern recognition and regression. This theory raises the generalized ability of learning mechanisms according to the minimum structural risk. Therefore, the advantage of using SVM can obtain results of minimized error rates without many training samples. Finally, the experimental results confirm that SVM method performs well in setting a critical rainfall line for each group of debris-flow streams.
Effect of thermodynamic disequilibrium on critical liquid-vapor flow conditions
Bilicki, Z.; Kestin, J.
1989-01-01
In this lecture we characterize the effect of absence of unconstrained thermodynamic equilibrium and onset of a metastable state on the adiabatic flow of a mixture of liquid and its vapor through a convergent-divergent nozzle. We study steady-state flows and emphasize the relations that are present when the flow is choked. In such cases, there exists a cross-section in which the flow is critical and in which the adiabatic wave of small amplitude is stationary. More precisely, the relaxation process which results from the lack of equilibrium causes the system to be dispersive. In such circumstances, the critical velocity is equal to the frozen speed of sound, a/sub f/ corresponding to /omega/ /yields/ /infinity/. The relaxation process displaces the critical cross-section quite far downstream from the throat and places it in the divergent portion of the channel. We present the topological portrait of solutions in a suitably defined state-velocity space and discuss the potential appearance of normal and dispersed shock waves. In extreme cases, the singular point (usually a saddle) which enables the flow to become supercritical is displaced so far that it is located outside the exit. Then, the flow velocity is everywhere subcritical (w < a/sub f/) even though it may exceed the equilibrium speed of sound (w /approx gt/ a/sub e/) beyond a certain cross-section, and in spite of the presence of a throat. 10 refs., 4 figs.
Modelling Canopy Flows over Complex Terrain
NASA Astrophysics Data System (ADS)
Grant, Eleanor R.; Ross, Andrew N.; Gardiner, Barry A.
2016-06-01
Recent studies of flow over forested hills have been motivated by a number of important applications including understanding CO_2 and other gaseous fluxes over forests in complex terrain, predicting wind damage to trees, and modelling wind energy potential at forested sites. Current modelling studies have focussed almost exclusively on highly idealized, and usually fully forested, hills. Here, we present model results for a site on the Isle of Arran, Scotland with complex terrain and heterogeneous forest canopy. The model uses an explicit representation of the canopy and a 1.5-order turbulence closure for flow within and above the canopy. The validity of the closure scheme is assessed using turbulence data from a field experiment before comparing predictions of the full model with field observations. For near-neutral stability, the results compare well with the observations, showing that such a relatively simple canopy model can accurately reproduce the flow patterns observed over complex terrain and realistic, variable forest cover, while at the same time remaining computationally feasible for real case studies. The model allows closer examination of the flow separation observed over complex forested terrain. Comparisons with model simulations using a roughness length parametrization show significant differences, particularly with respect to flow separation, highlighting the need to explicitly model the forest canopy if detailed predictions of near-surface flow around forests are required.
NASA Astrophysics Data System (ADS)
Afanasyev, A.
2011-12-01
Multiphase flows in porous media with a transition between sub- and supercritical thermodynamic conditions occur in many natural and technological processes (e.g. in deep regions of geothermal reservoirs where temperature reaches critical point of water or in gas-condensate fields where subject to critical conditions retrograde condensation occurs and even in underground carbon dioxide sequestration processes at high formation pressure). Simulation of these processes is complicated due to degeneration of conservation laws under critical conditions and requires non-classical mathematical models and methods. A new mathematical model is proposed for efficient simulation of binary mixture flows in a wide range of pressures and temperatures that includes critical conditions. The distinctive feature of the model lies in the methodology for mixture properties determination. Transport equations and Darcy law are solved together with calculation of the entropy maximum that is reached in thermodynamic equilibrium and determines mixture composition. To define and solve the problem only one function - mixture thermodynamic potential - is required. Such approach allows determination not only single-phase states and two-phase states of liquid-gas type as in classical models but also two-phase states of liquid-liquid type and three-phase states. The proposed mixture model was implemented in MUFITS (Multiphase Filtration Transport Simulator) code for hydrodynamic simulations. As opposed to classical approaches pressure, enthalpy and composition variables together with fully implicit method and cascade procedure are used. The code is capable of unstructured grids, heterogeneous porous media, relative permeability and capillary pressure dependence on temperature and pressure, multiphase diffusion, optional number of sink and sources, etc. There is an additional module for mixture properties specification. The starting point for the simulation is a cubic equation of state that is
Modeling information flow in biological networks
NASA Astrophysics Data System (ADS)
Kim, Yoo-Ah; Przytycki, Jozef H.; Wuchty, Stefan; Przytycka, Teresa M.
2011-06-01
Large-scale molecular interaction networks are being increasingly used to provide a system level view of cellular processes. Modeling communications between nodes in such huge networks as information flows is useful for dissecting dynamical dependences between individual network components. In the information flow model, individual nodes are assumed to communicate with each other by propagating the signals through intermediate nodes in the network. In this paper, we first provide an overview of the state of the art of research in the network analysis based on information flow models. In the second part, we describe our computational method underlying our recent work on discovering dysregulated pathways in glioma. Motivated by applications to inferring information flow from genotype to phenotype in a very large human interaction network, we generalized previous approaches to compute information flows for a large number of instances and also provided a formal proof for the method.
Conveying Lava Flow Hazards Through Interactive Computer Models
NASA Astrophysics Data System (ADS)
Thomas, D.; Edwards, H. K.; Harnish, E. P.
2007-12-01
As part of an Information Sciences senior class project, a software package of an interactive version of the FLOWGO model was developed for the Island of Hawaii. The software is intended for use in an ongoing public outreach and hazards awareness program that educates the public about lava flow hazards on the island. The design parameters for the model allow an unsophisticated user to initiate a lava flow anywhere on the island and allow it to flow down-slope to the shoreline while displaying a timer to show the rate of advance of the flow. The user is also able to modify a range of input parameters including eruption rate, the temperature of the lava at the vent, and crystal fraction present in the lava at the source. The flow trajectories are computed using a 30 m digital elevation model for the island and the rate of advance of the flow is estimated using the average slope angle and the computed viscosity of the lava as it cools in either a channel (high heat loss) or lava tube (low heat loss). Even though the FLOWGO model is not intended to, and cannot, accurately predict the rate of advance of a tube- fed or channel-fed flow, the relative rates of flow advance for steep or flat-lying terrain convey critically important hazard information to the public: communities located on the steeply sloping western flanks of Mauna Loa may have no more than a few hours to evacuate in the face of a threatened flow from Mauna Loa's southwest rift whereas communities on the more gently sloping eastern flanks of Mauna Loa and Kilauea may have weeks to months to prepare for evacuation. Further, the model also can show the effects of loss of critical infrastructure with consequent impacts on access into and out of communities, loss of electrical supply, and communications as a result of lava flow implacement. The interactive model has been well received in an outreach setting and typically generates greater involvement by the participants than has been the case with static maps
Entropy in the Bak-Sneppen Model for Self-Organized Criticality
NASA Astrophysics Data System (ADS)
Yang, Chun-Bin
2003-03-01
The distributions of fitness on the sites of one- and two-dimensional lattices are studied for the nearest-neighbour Bak-Sneppen model on self-organized criticality. The distributions show complicated behaviour showing that the system is far from equilibrium. By introducing the ``energy'' of a site, the entropy flow from the system to its environment is investigated.
A critical comparison of two-equation turbulence models
NASA Technical Reports Server (NTRS)
Lang, N. J.; Shih, T. H.
1991-01-01
Several two-equation models were proposed and tested against benchmark flows by various researchers. For each study, different numerical methods or codes were used to obtain the results which were reported to be an improvement over other models. However, these comparisons may be overshadowed by the different numerical schemes used to obtain the results. With this in mind, several existing two-equation turbulence models, including k-epsilon, k-tau, k-omega, and q-omega models, are implemented into a common flow solver code for near wall turbulent flows. The quality of each model is based on several criteria, including robustness and accuracy of predicting the turbulent quantities.
Modelling boundary layer flow over barnacle-fouled surfaces
NASA Astrophysics Data System (ADS)
Sadique, Jasim; Yang, Xiang; Meneveau, Charles; Mittal, Rajat
2014-11-01
Macro-biofouling is a critical concern for the marine industry. However, there is little data on flow and drag over such surfaces. Accurate modelling of such multi-scale flows remains a big challenge. Such simulations are vital in providing insights into the fundamental flow physics, and they can be used to estimate the timing, need and effectiveness of measures used to counteract bio-fouling. This talk focuses on the use of a sharp-interface immersed boundary method coupled with a wall model and large-eddy simulations to carry out accurate simulations of a turbulent boundary layer flow over macro-fouled surfaces. For the current study, high resolution scans of barnacles were used to create simple geometrical representations. Simulations were then carried out to test how well these simpler geometric models mimic the flow over actual barnacles. Simulations of array of modeled barnacles, with different barnacle densities have also been carried out and we present results on the effect distribution density on the flow physics and drag on the surfaces. This work is funded by ONR Grant N00014-12-1-0582.
Critical review of glass performance modeling
Bourcier, W.L.
1994-07-01
Borosilicate glass is to be used for permanent disposal of high-level nuclear waste in a geologic repository. Mechanistic chemical models are used to predict the rate at which radionuclides will be released from the glass under repository conditions. The most successful and useful of these models link reaction path geochemical modeling programs with a glass dissolution rate law that is consistent with transition state theory. These models have been used to simulate several types of short-term laboratory tests of glass dissolution and to predict the long-term performance of the glass in a repository. Although mechanistically based, the current models are limited by a lack of unambiguous experimental support for some of their assumptions. The most severe problem of this type is the lack of an existing validated mechanism that controls long-term glass dissolution rates. Current models can be improved by performing carefully designed experiments and using the experimental results to validate the rate-controlling mechanisms implicit in the models. These models should be supported with long-term experiments to be used for model validation. The mechanistic basis of the models should be explored by using modern molecular simulations such as molecular orbital and molecular dynamics to investigate both the glass structure and its dissolution process.
Turbulent motion of mass flows. Mathematical modeling
NASA Astrophysics Data System (ADS)
Eglit, Margarita; Yakubenko, Alexander; Yakubenko, Tatiana
2016-04-01
New mathematical models for unsteady turbulent mass flows, e.g., dense snow avalanches and landslides, are presented. Such models are important since most of large scale flows are turbulent. In addition to turbulence, the two other important points are taken into account: the entrainment of the underlying material by the flow and the nonlinear rheology of moving material. The majority of existing models are based on the depth-averaged equations and the turbulent character of the flow is accounted by inclusion of drag proportional to the velocity squared. In this paper full (not depth-averaged) equations are used. It is assumed that basal entrainment takes place if the bed friction equals the shear strength of the underlying layer (Issler D, M. Pastor Peréz. 2011). The turbulent characteristics of the flow are calculated using a three-parameter differential model (Lushchik et al., 1978). The rheological properties of moving material are modeled by one of the three types of equations: 1) Newtonian fluid with high viscosity, 2) power-law fluid and 3) Bingham fluid. Unsteady turbulent flows down long homogeneous slope are considered. The flow dynamical parameters and entrainment rate behavior in time as well as their dependence on properties of moving and underlying materials are studied numerically. REFERENCES M.E. Eglit and A.E. Yakubenko, 2014. Numerical modeling of slope flows entraining bottom material. Cold Reg. Sci. Technol., 108, 139-148 Margarita E. Eglit and Alexander E. Yakubenko, 2016. The effect of bed material entrainment and non-Newtonian rheology on dynamics of turbulent slope flows. Fluid Dynamics, 51(3) Issler D, M. Pastor Peréz. 2011. Interplay of entrainment and rheology in snow avalanches; a numerical study. Annals of Glaciology, 52(58), 143-147 Lushchik, V.G., Paveliev, A.A. , and Yakubenko, A.E., 1978. Three-parameter model of shear turbulence. Fluid Dynamics, 13, (3), 350-362
Critical Differences of Asymmetric Magnetic Reconnection from Standard Models
NASA Astrophysics Data System (ADS)
Nitta, S.; Wada, T.; Fuchida, T.; Kondoh, K.
2016-09-01
We have clarified the structure of asymmetric magnetic reconnection in detail as the result of the spontaneous evolutionary process. The asymmetry is imposed as ratio k of the magnetic field strength in both sides of the initial current sheet (CS) in the isothermal equilibrium. The MHD simulation is carried out by the HLLD code for the long-term temporal evolution with very high spatial resolution. The resultant structure is drastically different from the symmetric case (e.g., the Petschek model) even for slight asymmetry k = 2. (1) The velocity distribution in the reconnection jet clearly shows a two-layered structure, i.e., the high-speed sub-layer in which the flow is almost field aligned and the acceleration sub-layer. (2) Higher beta side (HBS) plasma is caught in a lower beta side plasmoid. This suggests a new plasma mixing process in the reconnection events. (3) A new large strong fast shock in front of the plasmoid forms in the HBS. This can be a new particle acceleration site in the reconnection system. These critical properties that have not been reported in previous works suggest that we contribute to a better and more detailed knowledge of the reconnection of the standard model for the symmetric magnetic reconnection system.
Compressor Flow Control Concepts. 2; UEET Compressor Flow Control Modeling
NASA Technical Reports Server (NTRS)
Chima, Rodrick V.
2001-01-01
Several passive flow control devices have been modeled computationally in the Swift CFD code. The models were applied to the first stage rotor and stator of the baseline UEET compressor in an attempt to improve efficiency and/or stall margin. The devices included suction surface bleed, tip injection, self-aspirated rotors, area-ruled casing, and vortex generators. The models and computed results will be described in the presentation. None of the results have shown significant gains in efficiency; however, casing vortex generators have shown potential improvements in stall margin.
Experimental Conditions for the Minimum Critical Reynolds Number in Pipe Flows
NASA Astrophysics Data System (ADS)
Kanda, Hidesada
2004-11-01
A transition model has been constructed for determining a critical Reynolds number Rc between laminar and turbulent flows in circular pipes. (1) Transition occurs in the entrance region under the conditions of natural disturbances. (2) Entrance shape determines Rc with disturbances; Rc takes a minimum value when the contraction ratio at the inlet is minimum or in the case of a straight pipe. (3) In the entrance boundary layer, there exists a normal wall strength which decreases as Re increases; let its rate of doing work be NW. The velocity profile develops from a uniform distribution at the inlet to a parabolic one; KE is the constant required acceleration power. Thus, the occurrence of the transition is determined by RW and KE: when RW > KE, transition never occurs, and when RW < KE, transition occurs. To prove the model, color-dye experiments were carried out. Consequently, two different types of Rc existed: Rc1 from laminar to turbulent and Rc2 from turbulent to laminar. As the ratio of bellmouth diameter to pipe diameter increases, Rc2 increased. Both Rc1 and Rc2 took a minimum value of approximately 2000 in the case of a straight pipe.
A Model for Critical Games Literacy
ERIC Educational Resources Information Center
Apperley, Tom; Beavis, Catherine
2013-01-01
This article outlines a model for teaching both computer games and videogames in the classroom for teachers. The model illustrates the connections between in-game actions and youth gaming culture. The article explains how the out-of-school knowledge building, creation and collaboration that occurs in gaming and gaming culture has an impact on…
Mathematical Models of Continuous Flow Electrophoresis
NASA Technical Reports Server (NTRS)
Saville, D. A.; Snyder, R. S.
1985-01-01
Development of high resolution continuous flow electrophoresis devices ultimately requires comprehensive understanding of the ways various phenomena and processes facilitate or hinder separation. A comprehensive model of the actual three dimensional flow, temperature and electric fields was developed to provide guidance in the design of electrophoresis chambers for specific tasks and means of interpreting test data on a given chamber. Part of the process of model development includes experimental and theoretical studies of hydrodynamic stability. This is necessary to understand the origin of mixing flows observed with wide gap gravitational effects. To insure that the model accurately reflects the flow field and particle motion requires extensive experimental work. Another part of the investigation is concerned with the behavior of concentrated sample suspensions with regard to sample stream stability particle-particle interactions which might affect separation in an electric field, especially at high field strengths. Mathematical models will be developed and tested to establish the roles of the various interactions.
CFD Modeling for Active Flow Control
NASA Technical Reports Server (NTRS)
Buning, Pieter G.
2001-01-01
This presentation describes current work under UEET Active Flow Control CFD Research Tool Development. The goal of this work is to develop computational tools for inlet active flow control design. This year s objectives were to perform CFD simulations of fully gridded vane vortex generators, micro-vortex genera- tors, and synthetic jets, and to compare flowfield results with wind tunnel tests of simple geometries with flow control devices. Comparisons are shown for a single micro-vortex generator on a flat plate, and for flow over an expansion ramp with sidewall effects. Vortex core location, pressure gradient and oil flow patterns are compared between experiment and computation. This work lays the groundwork for evaluating simplified modeling of arrays of devices, and provides the opportunity to test simple flow control device/sensor/ control loop interaction.
Modeling of Turbulent Flow in Electromagnetically Levitated Metal Droplets
NASA Technical Reports Server (NTRS)
Berry, S.; Hyers, R. W.; Abedian, B.; Racz, L. M.; Rose, M. Franklin (Technical Monitor)
2001-01-01
This article details an effort to improve the understanding and prediction of turbulent flow inside a droplet of molten metal levitated in an electromagnetic field. It is shown that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is in the transitional regime between laminar and turbulent flow. Past research efforts have used laminar, enhanced viscosity, and k-epsilon turbulence models to describe the flow. The method highlighted in our study is the renormalization group (RNG) algorithm. We show that an accurate description of the turbulent eddy viscosity is critical in order to obtain realistic velocity fields, and that the turbulent eddy viscosity cannot be uniform in levitated droplets. The RNG method does not impose isotropic length or time scales on the flow field, thus allowing such nonuniform features to be captured. A number of other materials processing applications exhibit similarly complex flow characteristics, such as highly recirculating, transitional, and free surface flows, for which this modeling approach may prove useful.
Grain size reduction in granular flows of spheres - The effects of critical impact energy
NASA Technical Reports Server (NTRS)
Richman, M. W.; Oyediran, A. A.
1992-01-01
Methods employed to derive recent kinetic theories for rapid noncomminuting granular flows are extended to homogeneous flows in which a fraction of the repeated collisions produce tiny fractures on the particles' peripheries and gradually reduce their effective diameters. The theory consists of balance equations for mass, momentum, and energy, as well as constitutive relations for the presence tensor and collisional rates of mass and energy lost. The work of Richman and Chou (1989) is improved by incorporating into the constitutive theory the critical impact energy below which no mass loss occurs in a binary collision. The theory is applied to granular shear flows and, for fixed shear rates, predicts the time variations of the solid fraction granular temperature, and induced stresses, as well as their extreme sensitivities to small changes in the critical impact energy.
Transient Wellbore Fluid Flow Model
Energy Science and Technology Software Center (ESTSC)
1982-04-06
WELBORE is a code to solve transient, one-dimensional two-phase or single-phase non-isothermal fluid flow in a wellbore. The primary thermodynamic variables used in solving the equations are the pressure and specific energy. An equation of state subroutine provides the density, quality, and temperature. The heat loss out of the wellbore is calculated by solving a radial diffusion equation for the temperature changes outside the bore. The calculation is done at each node point in themore » wellbore.« less
Critical fluctuations in the Bak-Sneppen model
NASA Astrophysics Data System (ADS)
Xiao, S. S.; Yang, C. B.
2012-05-01
The critical fluctuation properties of fitness distribution in one-dimensional Bak-Sneppen model are studied in terms of the normalized factorial moments and erraticity moments. For a fitness window above the critical gap, no anomaly scaling of the moments is found, while for a window below the gap intermittent behaviors are observed. The scaling exponent for the BS model is different from that in two dimensional Ising model for second order phase transition.
Neural network model for extracting optic flow.
Tohyama, Kazuya; Fukushima, Kunihiko
2005-01-01
When we travel in an environment, we have an optic flow on the retina. Neurons in the area MST of macaque monkeys are reported to have a very large receptive field and analyze optic flows on the retina. Many MST-cells respond selectively to rotation, expansion/contraction and planar motion of the optic flow. Many of them show position-invariant responses to optic flow, that is, their responses are maintained during the shift of the center of the optic flow. It has long been suggested mathematically that vector-field calculus is useful for analyzing optic flow field. Biologically, plausible neural network models based on this idea, however, have little been proposed so far. This paper, based on vector-field hypothesis, proposes a neural network model for extracting optic flows. Our model consists of hierarchically connected layers: retina, V1, MT and MST. V1-cells measure local velocity. There are two kinds of MT-cell: one is for extracting absolute velocities, the other for extracting relative velocities with their antagonistic inputs. Collecting signals from MT-cells, MST-cells respond selectively to various types of optic flows. We demonstrate through a computer simulation that this simple network is enough to explain a variety of results of neurophysiological experiments. PMID:16112546
Optimization of solver for gas flow modeling
NASA Astrophysics Data System (ADS)
Savichkin, D.; Dodulad, O.; Kloss, Yu
2014-05-01
The main purpose of the work is optimization of the solver for rarefied gas flow modeling based on the Boltzmann equation. Optimization method is based on SIMD extensions for ×86 processors. Computational code is profiled and manually optimized with SSE instructions. Heat flow, shock waves and Knudsen pump are modeled with optimized solver. Dependencies of computational time from mesh sizes and CPU capabilities are provided.
Holistic Flow Model of Spiritual Wellness
ERIC Educational Resources Information Center
Purdy, Melanie; Dupey, Peggy
2005-01-01
The Holistic Flow Model of Spiritual Wellness is a conceptualization of spiritual health and well-being that has implications for clinical practice and research. The model is unique in its placement of the spirit at the center of Life and in its fluid vision of the spirit. The authors present the model after a discussion of spirituality and the…
Software reliability models for critical applications
Pham, H.; Pham, M.
1991-12-01
This report presents the results of the first phase of the ongoing EG G Idaho, Inc. Software Reliability Research Program. The program is studying the existing software reliability models and proposes a state-of-the-art software reliability model that is relevant to the nuclear reactor control environment. This report consists of three parts: (1) summaries of the literature review of existing software reliability and fault tolerant software reliability models and their related issues, (2) proposed technique for software reliability enhancement, and (3) general discussion and future research. The development of this proposed state-of-the-art software reliability model will be performed in the second place. 407 refs., 4 figs., 2 tabs.
Software reliability models for critical applications
Pham, H.; Pham, M.
1991-12-01
This report presents the results of the first phase of the ongoing EG&G Idaho, Inc. Software Reliability Research Program. The program is studying the existing software reliability models and proposes a state-of-the-art software reliability model that is relevant to the nuclear reactor control environment. This report consists of three parts: (1) summaries of the literature review of existing software reliability and fault tolerant software reliability models and their related issues, (2) proposed technique for software reliability enhancement, and (3) general discussion and future research. The development of this proposed state-of-the-art software reliability model will be performed in the second place. 407 refs., 4 figs., 2 tabs.
A compendium of fracture flow models, 1994
Diodato, D.M.
1994-11-01
The report is designed to be used as a decision-making aid for individuals who need to simulate fluid flow in fractured porous media. Fracture flow codes of varying capability in the public and private domain were identified in a survey of government, academia, and industry. The selection and use of an appropriate code requires conceptualization of the geology, physics, and chemistry (for transport) of the fracture flow problem to be solved. Conceptual models that have been invoked to describe fluid flow in fractured porous media include explicit discrete fracture, dual continuum (porosity and/or permeability), discrete fracture network, multiple interacting continua, multipermeability/multiporosity, and single equivalent continuum. The explicit discrete-fracture model is a ``near-field`` representation, the single equivalent continuum model is a ``far-field`` representation, and the dual-continuum model is intermediate to those end members. Of these, the dual-continuum model is the most widely employed. The concept of multiple interacting continua has been applied in a limited number of examples. Multipermeability/multiporosity provides a unified conceptual model. The ability to accurately describe fracture flow phenomena will continue to improve as a result of advances in fracture flow research and computing technology. This improvement will result in enhanced capability to protect the public environment, safety, and health.
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.
Modeling and Prediction of Hot Deformation Flow Curves
NASA Astrophysics Data System (ADS)
Mirzadeh, Hamed; Cabrera, Jose Maria; Najafizadeh, Abbas
2012-01-01
The modeling of hot flow stress and prediction of flow curves for unseen deformation conditions are important in metal-forming processes because any feasible mathematical simulation needs accurate flow description. In the current work, in an attempt to summarize, generalize, and introduce efficient methods, the dynamic recrystallization (DRX) flow curves of a 17-4 PH martensitic precipitation hardening stainless steel, a medium carbon microalloyed steel, and a 304 H austenitic stainless steel were modeled and predicted using (1) a hyperbolic sine equation with strain dependent constants, (2) a developed constitutive equation in a simple normalized stress-normalized strain form and its modified version, and (3) a feed-forward artificial neural network (ANN). These methods were critically discussed, and the ANN technique was found to be the best for the modeling available flow curves; however, the developed constitutive equation showed slightly better performance than that of ANN and significantly better predicted values than those of the hyperbolic sine equation in prediction of flow curves for unseen deformation conditions.
Effect of rolling motion on critical heat flux for subcooled flow boiling in vertical tube
Hwang, J. S.; Park, I. U.; Park, M. Y.; Park, G. C.
2012-07-01
This paper presents defining characteristics of the critical heat flux (CHF) for the boiling of R-134a in vertical tube operation under rolling motion in marine reactor. It is important to predict CHF of marine reactor having the rolling motion in order to increase the safety of the reactor. Marine Reactor Moving Simulator (MARMS) tests are conducted to measure the critical heat flux using R-134a flowing upward in a uniformly heated vertical tube under rolling motion. MARMS was rotated by motor and mechanical power transmission gear. The CHF tests were performed in a 9.5 mm I.D. test section with heated length of 1 m. Mass fluxes range from 285 to 1300 kg m{sup -2}s{sup -1}, inlet subcooling from 3 to 38 deg. C and outlet pressures from 13 to 24 bar. Amplitudes of rolling range from 15 to 40 degrees and periods from 6 to 12 sec. To convert the test conditions of CHF test using R-134a in water, Katto's fluid-to-fluid modeling was used in present investigation. A CHF correlation is presented which accounts for the effects of pressure, mass flux, inlet subcooling and rolling angle over all conditions tested. Unlike existing transient CHF experiments, CHF ratio of certain mass flux and pressure are different in rolling motion. For the mass fluxes below 500 kg m{sup -2}s{sup -1} at 13, 16 (region of relative low mass flux), CHF ratio was decreased but was increased above that mass flux (region of relative high mass flux). Moreover, CHF tend to enhance in entire mass flux at 24 bar. (authors)
A model for insect tracheolar flow
NASA Astrophysics Data System (ADS)
Staples, Anne; Chatterjee, Krishnashis
2015-11-01
Tracheoles are the terminal ends of the microscale tracheal channels present in most insect respiratory systems that transport air directly to the tissue. From a fluid dynamics perspective, tracheolar flow is notable because it lies at the intersection of several specialized fluid flow regimes. The flow through tracheoles is creeping, microscale gas flow in the rarefied regime. Here, we use lubrication theory to model the flow through a single microscale tracheole and take into account fluid-structure interactions through an imposed periodic wall deformation corresponding to the rhythmic abdominal compression found in insects, and rarefaction effects using slip boundary conditions. We compare the pressure, axial pressure gradient, and axial and radial velocities in the channel, and the volumetric flow rate through the channel for no-slip, low slip, and high slip conditions under two different channel deformation regimes. We find that the presence of slip tends to reduce the flow rate through the model tracheole and hypothesize that one of the mechanical functions of tracheoles is to act as a diffuser to decelerate the flow, enhance mixing, and increase the residency time of freshly oxygenated air at the surface of the tissue. This work was funded by the NSF under grant no. 1437387.
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.
Scaled Experimental Modeling of VHTR Plenum Flows
ICONE 15
2007-04-01
Abstract The Very High Temperature Reactor (VHTR) is the leading candidate for the Next Generation Nuclear Power (NGNP) Project in the U.S. which has the goal of demonstrating the production of emissions free electricity and hydrogen by 2015. Various scaled heated gas and water flow facilities were investigated for modeling VHTR upper and lower plenum flows during the decay heat portion of a pressurized conduction-cooldown scenario and for modeling thermal mixing and stratification (“thermal striping”) in the lower plenum during normal operation. It was concluded, based on phenomena scaling and instrumentation and other practical considerations, that a heated water flow scale model facility is preferable to a heated gas flow facility and to unheated facilities which use fluids with ranges of density to simulate the density effect of heating. For a heated water flow lower plenum model, both the Richardson numbers and Reynolds numbers may be approximately matched for conduction-cooldown natural circulation conditions. Thermal mixing during normal operation may be simulated but at lower, but still fully turbulent, Reynolds numbers than in the prototype. Natural circulation flows in the upper plenum may also be simulated in a separate heated water flow facility that uses the same plumbing as the lower plenum model. However, Reynolds number scaling distortions will occur at matching Richardson numbers due primarily to the necessity of using a reduced number of channels connected to the plenum than in the prototype (which has approximately 11,000 core channels connected to the upper plenum) in an otherwise geometrically scaled model. Experiments conducted in either or both facilities will meet the objectives of providing benchmark data for the validation of codes proposed for NGNP designs and safety studies, as well as providing a better understanding of the complex flow phenomena in the plenums.
Flow field mapping in data rack model
NASA Astrophysics Data System (ADS)
Manoch, L.; Matěcha, J.; Pohan, P.
2013-04-01
The main objective of this study was to map the flow field inside the data rack model, fitted with three 1U server models. The server model is based on the common four-processor 1U server. The main dimensions of the data rack model geometry are taken fully from the real geometry. Only the model was simplified with respect to the greatest possibility in the experimental measurements. The flow field mapping was carried out both experimentally and numerically. PIV (Particle Image Velocimetry) method was used for the experimental flow field mapping, when the flow field has been mapped for defined regions within the 2D/3D data rack model. Ansys CFX and OpenFOAM software were used for the numerical solution. Boundary conditions for numerical model were based on data obtained from experimental measurement of velocity profile at the output of the server mockup. This velocity profile was used as the input boundary condition in the calculation. In order to achieve greater consistency of the numerical model with experimental data, the numerical model was modified with regard to the results of experimental measurements. Results from the experimental and numerical measurements were compared and the areas of disparateness were identified. In further steps the obtained proven numerical model will be utilized for the real geometry of data racks and data.
Reduced order model of draft tube flow
NASA Astrophysics Data System (ADS)
Rudolf, P.; Štefan, D.
2014-03-01
Swirling flow with compact coherent structures is very good candidate for proper orthogonal decomposition (POD), i.e. for decomposition into eigenmodes, which are the cornerstones of the flow field. Present paper focuses on POD of steady flows, which correspond to different operating points of Francis turbine draft tube flow. Set of eigenmodes is built using a limited number of snapshots from computational simulations. Resulting reduced order model (ROM) describes whole operating range of the draft tube. ROM enables to interpolate in between the operating points exploiting the knowledge about significance of particular eigenmodes and thus reconstruct the velocity field in any operating point within the given range. Practical example, which employs axisymmetric simulations of the draft tube flow, illustrates accuracy of ROM in regions without vortex breakdown together with need for higher resolution of the snapshot database close to location of sudden flow changes (e.g. vortex breakdown). ROM based on POD interpolation is very suitable tool for insight into flow physics of the draft tube flows (especially energy transfers in between different operating points), for supply of data for subsequent stability analysis or as an initialization database for advanced flow simulations.
Modeling of Turbulent Free Shear Flows
NASA Technical Reports Server (NTRS)
Yoder, Dennis A.; DeBonis, James R.; Georgiadis, Nicolas J.
2013-01-01
The modeling of turbulent free shear flows is crucial to the simulation of many aerospace applications, yet often receives less attention than the modeling of wall boundary layers. Thus, while turbulence model development in general has proceeded very slowly in the past twenty years, progress for free shear flows has been even more so. This paper highlights some of the fundamental issues in modeling free shear flows for propulsion applications, presents a review of past modeling efforts, and identifies areas where further research is needed. Among the topics discussed are differences between planar and axisymmetric flows, development versus self-similar regions, the effect of compressibility and the evolution of compressibility corrections, the effect of temperature on jets, and the significance of turbulent Prandtl and Schmidt numbers for reacting shear flows. Large eddy simulation greatly reduces the amount of empiricism in the physical modeling, but is sensitive to a number of numerical issues. This paper includes an overview of the importance of numerical scheme, mesh resolution, boundary treatment, sub-grid modeling, and filtering in conducting a successful simulation.
Supersonic boundary-layer flow turbulence modeling
NASA Technical Reports Server (NTRS)
Wang, Chi-Rong
1993-01-01
Baldwin-Lomax and kappa-epsilon turbulence models were modified for use in Navier-Stokes numerical computations of Mach 2.9 supersonic turbulent boundary layer flows along compression ramps. The computational results of Reynolds shear stress profiles were compared with experimental data. The Baldwin-Lomax model was modified to account for the Reynolds shear stress amplification within the flow field. A hybrid kappa-epsilon model with viscous sublayer turbulence treatment was constructed to predict the Reynolds shear stress profiles within the entire flow field. These modified turbulence models were effective for the computations of the surface pressure and the skin friction factor variations along an 8 deg ramp surface. The hybrid kappa-epsilon model could improve the predictions of the Reynolds shear stress profile and the skin friction factor near the corner of a 16 deg ramp.
NASA Technical Reports Server (NTRS)
Papell, S. S.
1972-01-01
Buoyancy effects on the critical heat flux and general data trends for a liquid nitrogen internal flow system were determined by comparison of upflow and downflow data under identical test conditions. The test section had a 1.28 cm diameter flow passage and a 30.5 cm heated length which was subjected to uniform heat fluxes through resistance heating. Test conditions covered a range of pressures from 3.4 to 10.2 atm, inlet velocities from 0.23 to 3.51 m/sec, with the liquid nitrogen temperature at saturated inlet conditions. Data comparisons showed that the critical heat flux for downflow could be up to 36 percent lower than for upflow. A nonmonotonic relationship between the critical heat flux and velocity was determined for upflow but not for downflow. A limiting inlet velocity of 4.12 m/sec was determined to be the minimum velocity required to completely suppress the influence of buoyancy on the critical heat flux for this saturated inlet flow system. A correlation of this limiting fluid velocity is presented that was developed from previously published subcooled liquid nitrogen data and the saturated data of this investigation.
Modeling depth distributions of overland flows
NASA Astrophysics Data System (ADS)
Smith, Mark W.; Cox, Nicholas J.; Bracken, Louise J.
2011-02-01
Hydrological and erosion models use water depth to estimate routing velocity and resultant erosion at each spatial element. Yet the shear stress distribution imposed on the soil surface and any resulting flow detachment and rill incision is controlled by the full probability distribution of depths of overland flow. Terrestrial Laser Scanning (TLS) is used in conjunction with simple field-flume experiments to provide high-resolution measures of overland flow depth-distributions for three semi-arid hillslope transects with differing soil properties. A two-parameter gamma distribution is proposed as the optimum model for depths of both interrill and rill flows. The shape and scale parameters are shown to vary consistently with distance downslope reflecting the morphological signature of runoff processes. The scale parameter is related to the general increase of depths with discharge ( P < 0.0001) as flows gradually concentrate; the shape parameter is more related to the soil surface roughness and potentially provides a control on the rate of depth, but also velocity increase with discharge. Such interactions between surface roughness and overland flows are of crucial importance for flow hydraulics and modeling sediment transport.
Assessment of chemistry models for compressible reacting flows
NASA Astrophysics Data System (ADS)
Lapointe, Simon; Blanquart, Guillaume
2014-11-01
Recent technological advances in propulsion and power devices and renewed interest in the development of next generation supersonic and hypersonic vehicles have increased the need for detailed understanding of turbulence-combustion interactions in compressible reacting flows. In numerical simulations of such flows, accurate modeling of the fuel chemistry is a critical component of capturing the relevant physics. Various chemical models are currently being used in reacting flow simulations. However, the differences between these models and their impacts on the fluid dynamics in the context of compressible flows are not well understood. In the present work, a numerical code is developed to solve the fully coupled compressible conservation equations for reacting flows. The finite volume code is based on the theoretical and numerical framework developed by Oefelein (Prog. Aero. Sci. 42 (2006) 2-37) and employs an all-Mach-number formulation with dual time-stepping and preconditioning. The numerical approach is tested on turbulent premixed flames at high Karlovitz numbers. Different chemical models of varying complexity and computational cost are used and their effects are compared.
Reliable Communication Models in Interdependent Critical Infrastructure Networks
Lee, Sangkeun; Chinthavali, Supriya; Shankar, Mallikarjun
2016-01-01
Modern critical infrastructure networks are becoming increasingly interdependent where the failures in one network may cascade to other dependent networks, causing severe widespread national-scale failures. A number of previous efforts have been made to analyze the resiliency and robustness of interdependent networks based on different models. However, communication network, which plays an important role in today's infrastructures to detect and handle failures, has attracted little attention in the interdependency studies, and no previous models have captured enough practical features in the critical infrastructure networks. In this paper, we study the interdependencies between communication network and other kinds of critical infrastructure networks with an aim to identify vulnerable components and design resilient communication networks. We propose several interdependency models that systematically capture various features and dynamics of failures spreading in critical infrastructure networks. We also discuss several research challenges in building reliable communication solutions to handle failures in these models.
GENERALIZED VISCOPLASTIC MODELING OF DEBRIS FLOW.
Chen, Cheng-lung
1988-01-01
The earliest model developed by R. A. Bagnold was based on the concept of the 'dispersive' pressure generated by grain collisions. Some efforts have recently been made by theoreticians in non-Newtonian fluid mechanics to modify or improve Bagnold's concept or model. A viable rheological model should consist both of a rate-independent part and a rate-dependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i. e. , the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical for general use in debris-flow modeling. In fact, Bagnold's model is found to be only a particular case of the GVF model. analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold's simplified assumption of constant grain concentration.
Critical behaviors near the (tri-)critical end point of QCD within the NJL model
NASA Astrophysics Data System (ADS)
Lu, Ya; Du, Yi-Lun; Cui, Zhu-Fang; Zong, Hong-Shi
2015-10-01
We investigate the dynamical chiral symmetry breaking and its restoration at finite density and temperature within the two-flavor Nambu-Jona-Lasinio model, and mainly focus on the critical behaviors near the critical end point (CEP) and tricritical point (TCP) of quantum chromodynamics. The multi-solution region of the Nambu and Wigner ones is determined in the phase diagram for the massive and massless current quark, respectively. We use the various susceptibilities to locate the CEP/TCP and then extract the critical exponents near them. Our calculations reveal that the various susceptibilities share the same critical behaviors for the physical current quark mass, while they show different features in the chiral limit.
Numerical flow modeling of power plant windboxes
LaRose, J.A.; Hopkins, M.W.
1995-12-31
Numerical flow modeling has become an increasingly important design and analysis tool for improving the air distribution to power plant burners. Uniform air distribution allows the burners to perform as designed to achieve the lowest possible emissions and best fuel burn-out. Modifications can be made internal to the existing windbox to improve the burner-to-burner and burner peripheral air distributions. These modifications can include turning vanes, flow splitters, perforated plate, and burner shrouding. Numerical modeling allows the analysis of design trade-offs between adding flow resistance, fan power, and windbox modification construction cost. Numerical modeling has advantages over physical modeling in that actual geometric scales and air temperatures are used. Advantages over a field data based study include the ability to quickly and cheaply analyze a variety of design options without actually modifying the windbox, and the availability of significantly more data with which to interpret the results. Costs to perform a numerical study are generally one-half to one-third of the cost to perform a physical flow model and can be one-forth of the cost to perform a field study. The continued development of affordable, high speed, large memory workstations and reliable, commercially available computation fluid dynamics (CFD) software allows practical analyses of power plant windboxes. This paper discusses (1) the impact of air distribution on burner performance, (2) the methodology used to perform numerical flow modeling of power plant windboxes, and (3) the results from several windbox analyses including available post-modification observations.
Critical Fluctuations in Cortical Models Near Instability
Aburn, Matthew J.; Holmes, C. A.; Roberts, James A.; Boonstra, Tjeerd W.; Breakspear, Michael
2012-01-01
Computational studies often proceed from the premise that cortical dynamics operate in a linearly stable domain, where fluctuations dissipate quickly and show only short memory. Studies of human electroencephalography (EEG), however, have shown significant autocorrelation at time lags on the scale of minutes, indicating the need to consider regimes where non-linearities influence the dynamics. Statistical properties such as increased autocorrelation length, increased variance, power law scaling, and bistable switching have been suggested as generic indicators of the approach to bifurcation in non-linear dynamical systems. We study temporal fluctuations in a widely-employed computational model (the Jansen–Rit model) of cortical activity, examining the statistical signatures that accompany bifurcations. Approaching supercritical Hopf bifurcations through tuning of the background excitatory input, we find a dramatic increase in the autocorrelation length that depends sensitively on the direction in phase space of the input fluctuations and hence on which neuronal subpopulation is stochastically perturbed. Similar dependence on the input direction is found in the distribution of fluctuation size and duration, which show power law scaling that extends over four orders of magnitude at the Hopf bifurcation. We conjecture that the alignment in phase space between the input noise vector and the center manifold of the Hopf bifurcation is directly linked to these changes. These results are consistent with the possibility of statistical indicators of linear instability being detectable in real EEG time series. However, even in a simple cortical model, we find that these indicators may not necessarily be visible even when bifurcations are present because their expression can depend sensitively on the neuronal pathway of incoming fluctuations. PMID:22952464
Oscillating line source in a shear flow with a free surface: critical layer-like contributions
NASA Astrophysics Data System (ADS)
Ellingsen, Simen Å.; Tyvand, Peder A.
2016-07-01
The linearized water-wave radiation problem for an oscillating submerged line source in an inviscid shear flow with a free surface is investigated analytically at finite, constant depth in the presence of a shear flow varying linearly with depth. The surface velocity is taken to be zero relative to the oscillating source, so that Doppler effects are absent. The radiated wave out from the source is calculated based on Euler's equation of motion with the appropriate boundary and radiation conditions, and differs substantially from the solution obtained by assuming potential flow. To wit, an additional wave is found in the downstream direction in addition to the previously known dispersive wave solutions; this wave is non-dispersive and we show how it is the surface manifestation of a critical layer-like flow generated by the combination of shear and mass flux at the source, passively advected with the flow. As seen from a system moving at the fluid velocity at the source's depth, streamlines form closed curves in a manner similar to Kelvin's cat's eye vortices. A resonant frequency exists at which the critical wave resonates with the downstream propagating wave, resulting in a downstream wave pattern diverging linearly in amplitude away from the source.
Review and selection of unsaturated flow models
Reeves, M.; Baker, N.A.; Duguid, J.O.
1994-04-04
Since the 1960`s, ground-water flow models have been used for analysis of water resources problems. In the 1970`s, emphasis began to shift to analysis of waste management problems. This shift in emphasis was largely brought about by site selection activities for geologic repositories for disposal of high-level radioactive wastes. Model development during the 1970`s and well into the 1980`s focused primarily on saturated ground-water flow because geologic repositories in salt, basalt, granite, shale, and tuff were envisioned to be below the water table. Selection of the unsaturated zone at Yucca Mountain, Nevada, for potential disposal of waste began to shift model development toward unsaturated flow models. Under the US Department of Energy (DOE), the Civilian Radioactive Waste Management System Management and Operating Contractor (CRWMS M&O) has the responsibility to review, evaluate, and document existing computer models; to conduct performance assessments; and to develop performance assessment models, where necessary. This document describes the CRWMS M&O approach to model review and evaluation (Chapter 2), and the requirements for unsaturated flow models which are the bases for selection from among the current models (Chapter 3). Chapter 4 identifies existing models, and their characteristics. Through a detailed examination of characteristics, Chapter 5 presents the selection of models for testing. Chapter 6 discusses the testing and verification of selected models. Chapters 7 and 8 give conclusions and make recommendations, respectively. Chapter 9 records the major references for each of the models reviewed. Appendix A, a collection of technical reviews for each model, contains a more complete list of references. Finally, Appendix B characterizes the problems used for model testing.
Improved modeling techniques for turbomachinery flow fields
Lakshminarayana, B.; Fagan, J.R. Jr.
1995-12-31
This program has the objective of developing an improved methodology for modeling turbomachinery flow fields, including the prediction of losses and efficiency. Specifically, the program addresses the treatment of the mixing stress tensor terms attributed to deterministic flow field mechanisms required in steady-state Computational Fluid Dynamic (CFD) models for turbomachinery flow fields. These mixing stress tensors arise due to spatial and temporal fluctuations (in an absolute frame of reference) caused by rotor-stator interaction due to various blade rows and by blade-to-blade variation of flow properties. This will be accomplished in a cooperative program by Penn State University and the Allison Engine Company. These tasks include the acquisition of previously unavailable experimental data in a high-speed turbomachinery environment, the use of advanced techniques to analyze the data, and the development of a methodology to treat the deterministic component of the mixing stress tenor.
Groundwater Flow Model for Taos, New Mexico
NASA Astrophysics Data System (ADS)
Burck, P. W.; Barroll, P. W.; Core, A. B.; Rappuhn, D.
2003-12-01
The New Mexico Office of the State Engineer - Hydrology Bureau (OSE) has developed a regional groundwater flow model for Taos, New Mexico. The MODFLOW 2000 model will serve as a tool to evaluate alternatives in settlement negotiations in an on-going water rights adjudication. If current settlement negotiations fail, it is conceivable that the model might be used in support of litigation. OSE produced the model in cooperation with technical representatives of the various parties to the adjudication. Regional hydrogeologic data including well records, aquifer test results, stream flow measurements and seepage studies have been shared relatively freely among the parties. A recent deep drilling program conducted in conjunction with the negotiation effort has added substantially to the hydrogeologic data set. Among the hydrologic processes simulated by the model are mountain front recharge; areal recharge from precipitation; evapotranspiration; discharge from springs; river and stream flow; accretions to groundwater from irrigation return flow, seepage from acequias, canals, and ditches, and deep percolation; and pumping by municipal entities and mutual domestic water users associations. The resulting model files are available for all parties to review and evaluate. Comments are assessed and many have resulted in significant improvements to the model. At this stage, however, it is unclear whether adopting this cooperative approach will increase the likelihood of model acceptance by the parties.
Critical heat-flux experiments under low-flow conditions in a vertical annulus. [PWR; BWR; LMFBR
Mishima, K.; Ishii, M.
1982-03-01
An experimental study was performed on critical heat flux (CHF) at low flow conditions for low pressure steam-water upward flow in an annulus. The test section was transparent, therefore, visual observations of dryout as well as various instrumentations were made. The data indicated that a premature CHF occurred due to flow regime transition from churn-turbulent to annular flow. It is shown that the critical heat flux observed in the experiment is essentially similar to a flooding-limited burnout and the critical heat flux can be well reproduced by a nondimensional correlation derived from the previously obtained criterion for flow regime transition. The observed CHF values are much smaller than the standard high quality CHF criteria at low flow, corresponding to the annular flow film dryout. This result is very significant, because the coolability of a heater surface at low flow rates can be drastically reduced by the occurrence of this mode of CHF.
Analytic Model of Reactive Flow
Souers, P C; Vitello, P
2004-11-15
A simple analytic model allows prediction of rate constants and size effect behavior before a hydrocode run if size effect data exists. At infinite radius, it defines not only detonation velocity but also average detonation rate, pressure and energy. This allows the derivation of a generalized radius, which becomes larger as the explosive becomes more non-ideal. The model is applied to near-ideal PBX 9404, in-between ANFO and most non-ideal AN. The power of the pressure declines from 2.3, 1.5 to 0.8 across this set. The power of the burn fraction, F, is 0.8, 0 and 0, so that an F-term is important only for the ideal explosives. The size effect shapes change from concave-down to nearly straight to concave-up. Failure is associated with ideal explosives when the calculated detonation velocity turns in a double-valued way. The effect of the power of the pressure may be simulated by including a pressure cutoff in the detonation rate. The models allows comparison of a wide spectrum of explosives providing that a single detonation rate is feasible.
Analytic Model of Reactive Flow
Souers, P C; Vitello, P
2004-08-02
A simple analytic model allows prediction of rate constants and size effect behavior before a hydrocode run if size effect data exists. At infinite radius, it defines not only detonation velocity but also average detonation rate, pressure and energy. This allows the derivation of a generalized radius, which becomes larger as the explosive becomes more non-ideal. The model is applied to near-ideal PBX 9404, in-between ANFO and most non-ideal AN. The power of the pressure declines from 2.3, 1.5 to 0.8 across this set. The power of the burn fraction, F, is 0.8, 0 and 0, so that an F-term is important only for the ideal explosives. The size effect shapes change from concave-down to nearly straight to concave-up. Failure is associated with ideal explosives when the calculated detonation velocity turns in a double-valued way. The effect of the power of the pressure may be simulated by including a pressure cutoff in the detonation rate. The models allows comparison of a wide spectrum of explosives providing that a single detonation rate is feasible.
A model for transonic plasma flow
Guazzotto, Luca; Hameiri, Eliezer
2014-02-15
A linear, two-dimensional model of a transonic plasma flow in equilibrium is constructed and given an explicit solution in the form of a complex Laplace integral. The solution indicates that the transonic state can be solved as an elliptic boundary value problem, as is done in the numerical code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)]. Moreover, the presence of a hyperbolic region does not necessarily imply the presence of a discontinuity or any other singularity of the solution.
Advanced Numerical Modeling of Turbulent Atmospheric Flows
NASA Astrophysics Data System (ADS)
Kühnlein, Christian; Dörnbrack, Andreas; Gerz, Thomas
The present chapter introduces the method of computational simulation to predict and study turbulent atmospheric flows. This includes a description of the fundamental approach to computational simulation and the practical implementation using the technique of large-eddy simulation. In addition, selected contributions from IPA scientists to computational model development and various examples for applications are given. These examples include homogeneous turbulence, convective boundary layers, heated forest canopy, buoyant thermals, and large-scale flows with baroclinic wave instability.
Unsaturated zone flow modeling for GWTT-95
Ho, C.K.; Altman, S.J.; McKenna, S.A.; Arnold, B.W.
1995-12-31
In accordance with the Nuclear Regulatory Commission regulation regarding groundwater travel times at geologic repositories, various models of unsaturated flow in fractured tuff have been developed and implemented to assess groundwater travel times at the potential repository at Yucca Mountain, Nevada. Kaplan used one-dimensional models to describe the uncertainty and sensitivity of travel times to various processes at Yucca Mountain. Robey and Arnold et al. used a two-dimensional equivalent continuum model (ECM) with inter- and intra-unit heterogeneity in an attempt to assess fast-flow paths through the unsaturated, fractured tuff at Yucca Mountain (GWTT-94). However, significant flow through the fractures in previous models was not simulated due to the characteristics of the ECM, which requires the matrix to be nearly saturated before flow through the fractures is initiated. In the current study (GWTT-95), four two-dimensional cross-sections at Yucca Mountain are simulated using both the ECM and dual-permeability (DK) models. The properties of both the fracture and matrix domains are geostatistically simulated, yielding completely heterogeneous continua. Then, simulations of flow through the four cross-sections are performed using spatially nonuniform infiltration boundary conditions. Steady-state groundwater travel times from the potential repository to the water table are calculated.
Theory and Low-Order Modeling of Unsteady Airfoil Flows
NASA Astrophysics Data System (ADS)
Ramesh, Kiran
is hypothesized, and verified with experimental and computational data, that LEV formation always occurs at the same critical value of LESP irrespective of motion kinematics. Further, the applicability of the LESP criterion in influencing the occurrence of LEV formation is demonstrated. To model the growth and convection of leading-edge vortices, the unsteady thin-airfoil theory is augmented with discrete-vortex shedding from the leading edge. The LESP criterion is used to predict and modulate the shedding of leading-edge vorticity. Comparisons with experiments and CFD for test-cases with different airfoils, Reynolds numbers and motion kinematics, show that the method performs remarkably well in predicting force coefficients and flowfields for unsteady flows. The use of a single empirical parameter - the critical LESP value, allows the determination of onset, growth and termination of leading-edge vortex shedding. In the final part of the research, the discrete-vortex model is extended to flows where the freestream velocity is varying or small in comparison with motion velocity. With this extension, the method is made applicable to a larger set of 2D flows such as perching and hovering maneuvers, gusts, and sinusoidally varying freestream. Abstractions of perching and hovering are designed as test cases and used to validate the low-order model's performance in highly-unsteady, vortex-dominated flows. Alongside development of the low-order methodology, several features of unsteady flows are studied and analyzed with the aid of CFD and experiments. While remaining computationally inexpensive and retaining the essential flow-physics, the method is seen to be successful in prediction of both force coefficients and flow histories.
Preliminary Saturated-Zone Flow Model
1997-06-10
This milestone consists of an updated fully 3D model of ground-water flow within the saturated zone at Yucca Mountain, Nevada. All electronic files pertaining to this deliverable have been transferred via ftp transmission to Steve Bodnar (M and O) and the technical data base. The model was developed using a flow and transport simulator, FEHMN, developed at Los Alamos National Laboratory, and represents a collaborative effort between staff from the US Geological Survey and Los Alamos National Laboratory. The model contained in this deliverable is minimally calibrated and represents work in progress. The flow model developed for this milestone is designed to feed subsequent transport modeling studies at Los Alamos which also use the FEHMN software. In addition, a general-application parameter estimation routine, PEST, was used in conjunction with FEHMN to reduce the difference between observed and simulated values of hydraulic head through the adjustment of model variables. This deliverable in large part consists of the electronic files for Yucca Mountain Site saturated-zone flow model as it existed as of 6/6/97, including the executable version of FEHMN (accession no. MOL.19970610.0204) used to run the code on a Sun Ultrasparc I workstation. It is expected that users of the contents of this deliverable be knowledgeable about the oration of FEHMN.
Numerical modeling of laser thermal propulsion flows
NASA Technical Reports Server (NTRS)
Mccay, T. D.; Thoenes, J.
1984-01-01
An review of the problems associated with modeling laser thermal propulsion flows, a synopsis of the status of such models, and the attributes of a successful model are presented. The continuous gaseous hydrogen laser-supported combustion wave (LSCW) thruster, for which a high-energy laser system (preferably space-based) should exist by the time the propulsion technology is developed, is considered in particular. The model proposed by Raizer (1970) is based on the assumptions of one-dimensional flow at constant pressure with heat conduction as the principal heat transfer mechanism. Consideration is given to subsequent models which account for radiative transfer into the ambient gas; provide a two-dimensional generalization of Raizer's analysis for the subsonic propagation of laser sparks in air; include the effect of forward plasma radiation in a one-dimensional model; and attempt a time-dependent (elliptic) solution of the full Navier-Stokes equations for the flow in a simple axisymmetric thruster. Attention is also given to thruster and nozzle flow models and thermodynamic and transport properties.
Rarefied-flow Shuttle aerodynamics model
NASA Technical Reports Server (NTRS)
Blanchard, Robert C.; Larman, Kevin T.; Moats, Christina D.
1993-01-01
A rarefied-flow shuttle aerodynamic model spanning the hypersonic continuum to the free molecule-flow regime was formulated. The model development has evolved from the High Resolution Accelerometer Package (HiRAP) experiment conducted on the Orbiter since 1983. The complete model is described in detail. The model includes normal and axial hypersonic continuum coefficient equations as functions of angle-of-attack, body flap deflection, and elevon deflection. Normal and axial free molecule flow coefficient equations as a function of angle-of-attack are presented, along with flight derived rarefied-flow transition bridging formulae. Comparisons are made with data from the Operational Aerodynamic Design Data Book (OADDB), applicable wind-tunnel data, and recent flight data from STS-35 and STS-40. The flight-derived model aerodynamic force coefficient ratio is in good agreement with the wind-tunnel data and predicts the flight measured force coefficient ratios on STS-35 and STS-40. The model is not, however, in good agreement with the OADDB. But, the current OADDB does not predict the flight data force coefficient ratios of either STS-35 or STS-40 as accurately as the flight-derived model. Also, the OADDB differs with the wind-tunnel force coefficient ratio data.
Electropolishing the bore of metal capillary tubes: A technique for adjusting the critical flow.
Stoffels, J J; Ells, D R
1979-12-01
A technique has been developed for electropolishing the bore of metal capillary tubes. Although developed specifically for stainless-steel tubes, the technique should be directly applicable to other metals. Tubes with inside diameter as small as 0.20 mm and 110 mm long have been successfully electropolished. The electropolishing technique can be used to increase the critical flow of a capillary tube in a controllable way. PMID:18699437
A void distribution model-flashing flow
Riznic, J.; Ishii, M.; Afgan, N.
1987-01-01
A new model for flashing flow based on wall nucleations is proposed here and the model predictions are compared with some experimental data. In order to calculate the bubble number density, the bubble number transport equation with a distributed source from the wall nucleation sites was used. Thus it was possible to avoid the usual assumption of a constant bubble number density. Comparisons of the model with the data shows that the model based on the nucleation site density correlation appears to be acceptable to describe the vapor generation in the flashing flow. For the limited data examined, the comparisons show rather satisfactory agreement without using a floating parameter to adjust the model. This result indicated that, at least for the experimental conditions considered here, the mechanistic predictions of the flashing phenomenon is possible on the present wall nucleation based model.
Paten, Jeffrey A; Siadat, Seyed Mohammad; Susilo, Monica E; Ismail, Ebraheim N; Stoner, Jayson L; Rothstein, Jonathan P; Ruberti, Jeffrey W
2016-05-24
The type I collagen monomer is one of nature's most exquisite and prevalent structural tools. Its 300 nm triple-helical motifs assemble into tough extracellular fibers that transition seamlessly across tissue boundaries and exceed cell dimensions by up to 4 orders of magnitude. In spite of extensive investigation, no existing model satisfactorily explains how such continuous structures are generated and grown precisely where they are needed (aligned in the path of force) by discrete, microscale cells using materials with nanoscale dimensions. We present a simple fiber drawing experiment, which demonstrates that slightly concentrated type I collagen monomers can be "flow-crystallized" to form highly oriented, continuous, hierarchical fibers at cell-achievable strain rates (<1 s(-1)) and physiologically relevant concentrations (∼50 μM). We also show that application of tension following the drawing process maintains the structural integrity of the fibers. While mechanical tension has been shown to be a critical factor driving collagen fibril formation during tissue morphogenesis in developing animals, the precise role of force in the process of building tissue is not well understood. Our data directly couple mechanical tension, specifically the extensional strain rate, to collagen fibril assembly. We further derive a "growth equation" which predicts that application of extensional strains, either globally by developing muscles or locally by fibroblasts, can rapidly drive the fusion of already formed short fibrils to produce long-range, continuous fibers. The results provide a pathway to scalable connective tissue manufacturing and support a mechano-biological model of collagen fibril deposition and growth in vivo. PMID:27070851
NASA Astrophysics Data System (ADS)
Sawicki, Jerzy; Paczkowski, Tomasz
2015-05-01
The paper presents the results of experimental studies of electrochemical machining process oriented on occurring in the treatment critical states caused by electrolyte flow hydrodynamic conditions in the gap between electrodes. Material forming in electrochemical machining is carried out by anodic dissolution. In general in ECM process, the essence of the treatment is that the workpiece is the anode and the tool is the cathode. The space between the anode and cathode is filled by electrolyte. The current flow between the electrodes causes anodic dissolution process, resulting in the removal of material from the anode. Choosing in the process of electrochemical machining, respectively: anode and cathode material, electrolyte and processing parameters, such conditions can be created that enable a high process efficiency and smoothness of the surface. Inappropriate selection of machining parameters can cause the emergence of critical states in the ECM, which are mainly related to the flow of the electrolyte in the gap between electrodes. This work is an attempt to assess the occurring critical states in ECM on the example of machining of curved surfaces with any sort of outline and curved rotating surfaces.
Ribosome flow model with positive feedback
Margaliot, Michael; Tuller, Tamir
2013-01-01
Eukaryotic mRNAs usually form a circular structure; thus, ribosomes that terminatae translation at the 3′ end can diffuse with increased probability to the 5′ end of the transcript, initiating another cycle of translation. This phenomenon describes ribosomal flow with positive feedback—an increase in the flow of ribosomes terminating translating the open reading frame increases the ribosomal initiation rate. The aim of this paper is to model and rigorously analyse translation with feedback. We suggest a modified version of the ribosome flow model, called the ribosome flow model with input and output. In this model, the input is the initiation rate and the output is the translation rate. We analyse this model after closing the loop with a positive linear feedback. We show that the closed-loop system admits a unique globally asymptotically stable equilibrium point. From a biophysical point of view, this means that there exists a unique steady state of ribosome distributions along the mRNA, and thus a unique steady-state translation rate. The solution from any initial distribution will converge to this steady state. The steady-state distribution demonstrates a decrease in ribosome density along the coding sequence. For the case of constant elongation rates, we obtain expressions relating the model parameters to the equilibrium point. These results may perhaps be used to re-engineer the biological system in order to obtain a desired translation rate. PMID:23720534
Modeling magnetically insulated devices using flow impedance
Mendel, C.W. Jr.; Rosenthal, S.E. )
1995-04-01
In modern pulsed power systems the electric field stresses at metal surfaces in vacuum transmission lines are so high that negative surfaces are space-charge-limited electron emitters. These electrons do not cause unacceptable losses because magnetic fields due to system currents result in net motion parallel to the electrodes. It has been known for several years that a parameter known as flow impedance is useful for describing these flows. Flow impedance is a measure of the separation between the anode and the mean position of the electron cloud, and it will be shown in this paper that in many situations flow impedance depends upon the geometry of the transmission line upstream of the point of interest. It can be remarkably independent of other considerations such as line currents and voltage. For this reason flow impedance is a valuable design parameter. Models of impedance transitions and voltage adders using flow impedance will be developed. Results of these models will be compared to two-dimensional, time-dependent, particle-in-cell simulations.
Steady flow in abdominal aortic aneurysm models.
Budwig, R; Elger, D; Hooper, H; Slippy, J
1993-11-01
Steady flow in abdominal aortic aneurysm models has been examined for four aneurysm sizes over Reynolds numbers from 500 to 2600. The Reynolds number is based on entrance tube diameter, and the inlet condition is fully developed flow. Experimental and numerical methods have been used to determine: (i) the overall features of the flow, (ii) the stresses on the aneurysm walls in laminar flow, and (iii) the onset and characteristics of turbulent flow. The laminar flow field is characterized by a jet of fluid (passing directly through the aneurysm) surrounded by a recirculating vortex. The wall shear stress magnitude in the recirculation zone is about ten times less than in the entrance tube. Both wall shear stress and wall normal stress profiles exhibit large magnitude peaks near the reattachment point at the distal end of the aneurysm. The onset of turbulence in the model is intermittent for 2000 < Re < 2500. The results demonstrate that a slug of turbulence in the entrance tube grows much more rapidly in the aneurysm than in a corresponding length of uniform cross section pipe. When turbulence is present in the aneurysm the recirculation zone breaks down and the wall shear stress returns to a magnitude comparable to that in the entrance tube. PMID:8309237
A Substance Flow Model for Global Phosphorus
NASA Astrophysics Data System (ADS)
Vaccari, D. A.
2015-12-01
A system-based substance flow model (SFM) for phosphorus is developed based on the global phosphorus substance flow analysis (SFA) of Cordell et al (2009). The model is based strictly on mass balance considerations. It predicts the sensitivity of phosphorus consumption to various interventions intended to conserve reserves, as well as interactions among these efforts, allowing a comparison of their impacts on phosphorus demand. The interventions include control of phosphorus losses from soil erosion, food production and food waste, or phosphorus recycling such as from animal manure or human waste.
Enhanced critical heat flux by capillary driven liquid flow on the well-designed surface
NASA Astrophysics Data System (ADS)
Kim, Dong Eok; Park, Su Cheong; Yu, Dong In; Kim, Moo Hwan; Ahn, Ho Seon
2015-07-01
Based on the unique design of the surface morphology, we investigated the effects of gravity and capillary pressure on Critical heat flux (CHF). The micro-structured surfaces for pool boiling tests were comprised with both the rectangular cavity and microchannel structures. The microcavity structures could intrinsically block the liquid flow by capillary pressure effect, and the capillary flow into the boiling surface was one-dimensionally induced only through the microchannel region. Thus, we could clearly establish the relationship between the CHF and capillary wicking flow. The driving potentials for the liquid inflow can be classified into the hydrostatic head by gravitational force, and the capillary pressure induced by the interactions of vapor bubbles, liquid film, and surface solid structures. Through the analysis of the experimental data and visualization of vapor bubble behaviors, we present that the liquid supplement to maintain the nucleate boiling regime in pool boiling condition is governed by the gravitational pressure head and capillary pressure effect.
Discharge coefficient correlations for circular-arc venturi flowmeters at critical /sonic/ flow
NASA Technical Reports Server (NTRS)
Arnberg, B. T.; Britton, C. L.; Seidl, W. F.
1973-01-01
Experimental data are analyzed to support theoretical predictions for discharge coefficients in circular-arc venturi flow meters operating in the critical sonic flow regime at throat Reynolds numbers above 150 thousand. The data tend to verify the predicted 0.25% decrease in the discharge coefficient during transition from a laminar to turbulent boundary layer. Four different test gases and three flow measurement facilities were used in the experiments with 17 venturis with throat sizes from 0.15 to 1.37 in. and Beta ratios ranging from 0.014 to 0.25. Recommendations are given as to how the effectiveness of future studies in the field could be improved.
Critical heat flux of subcooled flow boiling with water for high heat flux application
NASA Astrophysics Data System (ADS)
Inasaka, Fujio; Nariai, Hideki
1993-11-01
Subcooled flow boiling in water is thought to be advantageous in removing high heat load of more than 10 MW/m2. Characteristics of the critical heat flux (CHF), which determines the upper limit of heat removal, are very important for the design of cooling systems. In this paper, studies on subcooled flow boiling CHF, which have been conducted by the authors, are reported. Experiments were conducted using direct current heating of stainless steel tube. For uniform heating conditions, CHF increment in small diameter tubes (1 - 3 mm inside diameter) and the CHF characteristics in tubes with internal twisted tapes were investigated, and also the existing CHF correlations for ordinary tubes (more than 3 mm inside diameter) were evaluated. For peripherally non-uniform heating conditions using the tube, whose wall thickness was partly reduced, the CHF for swirl flow was higher than the CHF under uniform heating conditions with an increase of the non-uniformity factor.
Modelling fluid flow in a reciprocating compressor
NASA Astrophysics Data System (ADS)
Tuhovcak, Jan; Hejčík, Jiří; Jícha, Miroslav
2015-05-01
Efficiency of reciprocating compressor is strongly dependent on the valves characteristics, which affects the flow through the suction and discharge line. Understanding the phenomenon inside the compressor is necessary step in development process. Commercial CFD tools offer wide capabilities to simulate the flow inside the reciprocating compressor, however they are too complicated in terms of computational time and mesh creation. Several parameters describing compressor could be therefore examined without the CFD analysis, such is valve characteristic, flow through the cycle and heat transfer. The aim of this paper is to show a numerical tool for reciprocating compressor based on the energy balance through the cycle, which provides valve characteristics, flow through the cycle and heat losses from the cylinder. Spring-damping-mass model was used for the valve description. Boundary conditions were extracted from the performance test of 4-cylinder semihermetic compressor and numerical tool validation was performed with indicated p-V diagram comparison.
Modeling of Multi-Scale Channeling Phenomena in Porous Flow
NASA Astrophysics Data System (ADS)
Räss, Ludovic; Omlin, Samuel; Yarushina, Viktoriya; Simon, Nina; Podladchikov, Yuri
2015-04-01
Predictive modeling of fluid percolation through tight porous rocks is critical to evaluate environmental risks associated with waste storage and reservoir operations. To understand the evolution of two-phase mixtures of fluid and solid it is insufficient to only combine single-phase fluid flow methods and solid mechanics. A proper coupling of these two different multi-scales physical processes is required to describe the complex evolution of permeability and porosity in space and in time. We conduct numerical modeling experiments in geometrically simple but physically complex systems of stressed rocks containing self-focusing porous flow. Our model is physically and thermodynamically consistent and describes the formation and evolution of fluid pathways. The model consists of a system of coupled equations describing poro-elasto-viscous deformation and flow. Nonlinearity of the solid rheology is also taken into account. We have developed a numerical application based on an iterative finite difference scheme that runs on mutli-GPUs cluster in parallel. In order to validate these models, we consider the largest CO2 sequestration project in operation at the Sleipner field in the Norwegian North Sea. Attempts to match the observations at Sleipner using conventional reservoir simulations fail to capture first order observations, such as the seemingly effortless vertical flow of CO2 through low permeability shale layers and the formation of focused flow channels or chimneys. Conducted high-resolution three-dimensional numerical simulations predict the formation of dynamically evolving high porosity and permeability pathways as a natural outcome of porous flow nonlinearly coupled with rock deformation, which may trigger leakage through low permeability barriers.
Bootstrapping Critical Ising Model on Three Dimensional Real Projective Space.
Nakayama, Yu
2016-04-01
Given conformal data on a flat Euclidean space, we use crosscap conformal bootstrap equations to numerically solve the Lee-Yang model as well as the critical Ising model on a three dimensional real projective space. We check the rapid convergence of our bootstrap program in two dimensions from the exact solutions available. Based on the comparison, we estimate that our systematic error on the numerically solved one-point functions of the critical Ising model on a three dimensional real projective space is less than 1%. Our method opens up a novel way to solve conformal field theories on nontrivial geometries. PMID:27104697
Bootstrapping Critical Ising Model on Three Dimensional Real Projective Space
NASA Astrophysics Data System (ADS)
Nakayama, Yu
2016-04-01
Given conformal data on a flat Euclidean space, we use crosscap conformal bootstrap equations to numerically solve the Lee-Yang model as well as the critical Ising model on a three dimensional real projective space. We check the rapid convergence of our bootstrap program in two dimensions from the exact solutions available. Based on the comparison, we estimate that our systematic error on the numerically solved one-point functions of the critical Ising model on a three dimensional real projective space is less than 1%. Our method opens up a novel way to solve conformal field theories on nontrivial geometries.
Modelling of critical phenomena for ignition of metal particles
NASA Astrophysics Data System (ADS)
Shchepakina, E.; Sobolev, V.
2008-11-01
The singularly perturbed system of differential equations describing metal particle ignition is analyzed. Two general types of combustion reactions are recognized, namely subcritical and supercritical. It is shown that there exists a so-called critical regime, which separates domains of subcritical and supercritical regimes. The critical regime is modelled by the special phase trajectory of the system, which includes the repulsive part with slow variation of the temperature of a metal particle and the thickness of the oxide film. The asymptotic formulae for the calculation of critical values of the modified Semenov number are obtained in the cases of parabolic and cubic laws of oxidation kinetics.
SATURATED ZONE FLOW AND TRANSPORT MODEL ABSTRACTION
B.W. ARNOLD
2004-10-27
The purpose of the saturated zone (SZ) flow and transport model abstraction task is to provide radionuclide-transport simulation results for use in the total system performance assessment (TSPA) for license application (LA) calculations. This task includes assessment of uncertainty in parameters that pertain to both groundwater flow and radionuclide transport in the models used for this purpose. This model report documents the following: (1) The SZ transport abstraction model, which consists of a set of radionuclide breakthrough curves at the accessible environment for use in the TSPA-LA simulations of radionuclide releases into the biosphere. These radionuclide breakthrough curves contain information on radionuclide-transport times through the SZ. (2) The SZ one-dimensional (I-D) transport model, which is incorporated in the TSPA-LA model to simulate the transport, decay, and ingrowth of radionuclide decay chains in the SZ. (3) The analysis of uncertainty in groundwater-flow and radionuclide-transport input parameters for the SZ transport abstraction model and the SZ 1-D transport model. (4) The analysis of the background concentration of alpha-emitting species in the groundwater of the SZ.