AFDM: An Advanced Fluid-Dynamics Model
Bohl, W.R.; Parker, F.R. ); Wilhelm, D. . Inst. fuer Neutronenphysik und Reaktortechnik); Berthier, J. ); Goutagny, L. . Inst. de Protection et de Surete Nucleaire); Ninokata,
1990-09-01
AFDM, or the Advanced Fluid-Dynamics Model, is a computer code that investigates new approaches simulating the multiphase-flow fluid-dynamics aspects of severe accidents in fast reactors. The AFDM formalism starts with differential equations similar to those in the SIMMER-II code. These equations are modified to treat three velocity fields and supplemented with a variety of new models. The AFDM code has 12 topologies describing what material contacts are possible depending on the presence or absence of a given material in a computational cell, on the dominant liquid, and on the continuous phase. Single-phase, bubbly, churn-turbulent, cellular, and dispersed flow regimes are permitted for the pool situations modeled. Virtual mass terms are included for vapor in liquid-continuous flow. Interfacial areas between the continuous and discontinuous phases are convected to allow some tracking of phenomenological histories. Interfacial areas are also modified by models of nucleation, dynamic forces, turbulence, flashing, coalescence, and mass transfer. Heat transfer is generally treated using engineering correlations. Liquid-vapor phase transitions are handled with the nonequilibrium, heat-transfer-limited model, whereas melting and freezing processes are based on equilibrium considerations. Convection is treated using a fractional-step method of time integration, including a semi-implicit pressure iteration. A higher-order differencing option is provided to control numerical diffusion. The Los Alamos SESAME equation-of-state has been implemented using densities and temperatures as the independent variables. AFDM programming has vectorized all computational loops consistent with the objective of producing an exportable code. 24 refs., 4 figs.
AFDM: An advanced fluid-dynamics model. Volume 6: EOS-AFDM interface
Henneges, G.; Kleinheins, S.
1994-01-01
This volume of the Advanced Fluid-Dynamics Model (AFDM) documents the modeling of the equation of state (EOS) in the code. The authors present an overview of the basic concepts underlying the thermodynamics modeling and resulting EOS, which is a set of relations between the thermodynamic properties of materials. The AFDM code allows for multiphase-multimaterial systems, which they explore in three phase models: two-material solid, two-material liquid, and three-material vapor. They describe and compare two ways of specifying the EOS of materials: (1) as simplified analytic expressions, or (2) as tables that precisely describe the properties of materials and their interactions for mechanical equilibrium. Either of the two EOS models implemented in AFDM can be selected by specifying the option when preprocessing the source code for compilation. Last, the authors determine thermophysical properties such as surface tension, thermal conductivities, and viscosities in the model for the intracell exchanges of AFDM. Specific notations, routines, EOS data, plots, test results, and corrections to the code are available in the appendices.
NASA Astrophysics Data System (ADS)
Gulamali, M. Y.; Saunders, J. H.; Jackson, M. D.; Pain, C. C.
2009-04-01
We present results from a new computational multi-fluid dynamics code, designed to model the transport of heat, mass and chemical species during flow of single or multiple immiscible fluid phases through porous media, including gravitational effects and compressibility. The model also captures the electrical phenomena which may arise through electrokinetic, electrochemical and electrothermal coupling. Building on the advanced computational technology of the Imperial College Ocean Model, this new development leads the way towards a complex multiphase code using arbitrary unstructured and adaptive meshes, and domains decomposed to run in parallel over a cluster of workstations or a dedicated parallel computer. These facilities will allow efficient and accurate modelling of multiphase flows which capture large- and small-scale transport phenomena, while preserving the important geology and/or surface topology to make the results physically meaningful and realistic. Applications include modelling of contaminant transport in aquifers, multiphase flow during hydrocarbon production, migration of carbon dioxide during sequestration, and evaluation of the design and safety of nuclear reactors. Simulations of the streaming potential resulting from multiphase flow in laboratory- and field-scale models demonstrate that streaming potential signals originate at fluid fronts, and at geologic boundaries where fluid saturation changes. This suggests that downhole measurements of streaming potential may be used to inform production strategies in oil and gas reservoirs. As water encroaches on an oil production well, the streaming-potential signal associated with the water front encompasses the well even when the front is up to 100 m away, so the potential measured at the well starts to change significantly relative to a distant reference electrode. Variations in the geometry of the encroaching water front could be characterized using an array of electrodes positioned along the well
Advanced Multigrid Solvers for Fluid Dynamics
NASA Technical Reports Server (NTRS)
Brandt, Achi
1999-01-01
The main objective of this project has been to support the development of multigrid techniques in computational fluid dynamics that can achieve "textbook multigrid efficiency" (TME), which is several orders of magnitude faster than current industrial CFD solvers. Toward that goal we have assembled a detailed table which lists every foreseen kind of computational difficulty for achieving it, together with the possible ways for resolving the difficulty, their current state of development, and references. We have developed several codes to test and demonstrate, in the framework of simple model problems, several approaches for overcoming the most important of the listed difficulties that had not been resolved before. In particular, TME has been demonstrated for incompressible flows on one hand, and for near-sonic flows on the other hand. General approaches were advanced for the relaxation of stagnation points and boundary conditions under various situations. Also, new algebraic multigrid techniques were formed for treating unstructured grid formulations. More details on all these are given below.
AIR INGRESS ANALYSIS: COMPUTATIONAL FLUID DYNAMIC MODELS
Chang H. Oh; Eung S. Kim; Richard Schultz; Hans Gougar; David Petti; Hyung S. Kang
2010-08-01
The Idaho National Laboratory (INL), under the auspices of the U.S. Department of Energy, is performing research and development that focuses on key phenomena important during potential scenarios that may occur in very high temperature reactors (VHTRs). Phenomena Identification and Ranking Studies to date have ranked an air ingress event, following on the heels of a VHTR depressurization, as important with regard to core safety. Consequently, the development of advanced air ingress-related models and verification and validation data are a very high priority. Following a loss of coolant and system depressurization incident, air will enter the core of the High Temperature Gas Cooled Reactor through the break, possibly causing oxidation of the in-the core and reflector graphite structure. Simple core and plant models indicate that, under certain circumstances, the oxidation may proceed at an elevated rate with additional heat generated from the oxidation reaction itself. Under postulated conditions of fluid flow and temperature, excessive degradation of the lower plenum graphite can lead to a loss of structural support. Excessive oxidation of core graphite can also lead to the release of fission products into the confinement, which could be detrimental to a reactor safety. Computational fluid dynamic model developed in this study will improve our understanding of this phenomenon. This paper presents two-dimensional and three-dimensional CFD results for the quantitative assessment of the air ingress phenomena. A portion of results of the density-driven stratified flow in the inlet pipe will be compared with results of the experimental results.
Fluid-dynamical model for antisurfactants
NASA Astrophysics Data System (ADS)
Conn, Justin J. A.; Duffy, Brian R.; Pritchard, David; Wilson, Stephen K.; Halling, Peter J.; Sefiane, Khellil
2016-04-01
We construct a fluid-dynamical model for the flow of a solution with a free surface at which surface tension acts. This model can describe both classical surfactants, which decrease the surface tension of the solution relative to that of the pure solvent, and antisurfactants (such as many salts when added to water, and small amounts of water when added to alcohol) which increase it. We demonstrate the utility of the model by considering the linear stability of an infinitely deep layer of initially quiescent fluid. In particular, we predict the occurrence of an instability driven by surface-tension gradients, which occurs for antisurfactant, but not for surfactant, solutions.
Direct modeling for computational fluid dynamics
NASA Astrophysics Data System (ADS)
Xu, Kun
2015-06-01
All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes (NS) equations. The current computational fluid dynamics (CFD) focuses on the numerical solution of partial differential equations (PDEs), and its aim is to get the accurate solution of these governing equations. Under such a CFD practice, it is hard to develop a unified scheme that covers flow physics from kinetic to hydrodynamic scales continuously because there is no such governing equation which could make a smooth transition from the Boltzmann to the NS modeling. The study of fluid dynamics needs to go beyond the traditional numerical partial differential equations. The emerging engineering applications, such as air-vehicle design for near-space flight and flow and heat transfer in micro-devices, do require further expansion of the concept of gas dynamics to a larger domain of physical reality, rather than the traditional distinguishable governing equations. At the current stage, the non-equilibrium flow physics has not yet been well explored or clearly understood due to the lack of appropriate tools. Unfortunately, under the current numerical PDE approach, it is hard to develop such a meaningful tool due to the absence of valid PDEs. In order to construct multiscale and multiphysics simulation methods similar to the modeling process of constructing the Boltzmann or the NS governing equations, the development of a numerical algorithm should be based on the first principle of physical modeling. In this paper, instead of following the traditional numerical PDE path, we introduce direct modeling as a principle for CFD algorithm development. Since all computations are conducted in a discretized space with limited cell resolution, the flow physics to be modeled has to be done in the mesh size and time step scales. Here, the CFD is more or less a direct
Computational fluid dynamics modelling in cardiovascular medicine
Morris, Paul D; Narracott, Andrew; von Tengg-Kobligk, Hendrik; Silva Soto, Daniel Alejandro; Hsiao, Sarah; Lungu, Angela; Evans, Paul; Bressloff, Neil W; Lawford, Patricia V; Hose, D Rodney; Gunn, Julian P
2016-01-01
This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards ‘digital patient’ or ‘virtual physiological human’ representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges. PMID:26512019
Computational fluid dynamics modelling in cardiovascular medicine.
Morris, Paul D; Narracott, Andrew; von Tengg-Kobligk, Hendrik; Silva Soto, Daniel Alejandro; Hsiao, Sarah; Lungu, Angela; Evans, Paul; Bressloff, Neil W; Lawford, Patricia V; Hose, D Rodney; Gunn, Julian P
2016-01-01
This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards 'digital patient' or 'virtual physiological human' representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges. PMID:26512019
COMPUTATIONAL FLUID DYNAMICS MODELING ANALYSIS OF COMBUSTORS
Mathur, M.P.; Freeman, Mark; Gera, Dinesh
2001-11-06
In the current fiscal year FY01, several CFD simulations were conducted to investigate the effects of moisture in biomass/coal, particle injection locations, and flow parameters on carbon burnout and NO{sub x} inside a 150 MW GEEZER industrial boiler. Various simulations were designed to predict the suitability of biomass cofiring in coal combustors, and to explore the possibility of using biomass as a reburning fuel to reduce NO{sub x}. Some additional CFD simulations were also conducted on CERF combustor to examine the combustion characteristics of pulverized coal in enriched O{sub 2}/CO{sub 2} environments. Most of the CFD models available in the literature treat particles to be point masses with uniform temperature inside the particles. This isothermal condition may not be suitable for larger biomass particles. To this end, a stand alone program was developed from the first principles to account for heat conduction from the surface of the particle to its center. It is envisaged that the recently developed non-isothermal stand alone module will be integrated with the Fluent solver during next fiscal year to accurately predict the carbon burnout from larger biomass particles. Anisotropy in heat transfer in radial and axial will be explored using different conductivities in radial and axial directions. The above models will be validated/tested on various fullscale industrial boilers. The current NO{sub x} modules will be modified to account for local CH, CH{sub 2}, and CH{sub 3} radicals chemistry, currently it is based on global chemistry. It may also be worth exploring the effect of enriched O{sub 2}/CO{sub 2} environment on carbon burnout and NO{sub x} concentration. The research objective of this study is to develop a 3-Dimensional Combustor Model for Biomass Co-firing and reburning applications using the Fluent Computational Fluid Dynamics Code.
Tentner, Adrian; Lo, Simon; Ioilev, Andrey; Melnikov, Vladimir; Samigulin, Maskhud; Ustinenko, Vasily; Kozlov, Valentin
2006-07-01
A new code, CFD-BWR, is being developed for the simulation of two-phase flow phenomena inside a BWR fuel bundle. These phenomena include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. CFD-BWR is a specialized module built on the foundation of the commercial CFD code STAR-CD which provides general two-phase flow modeling capabilities. New models describing the inter-phase mass, momentum, and energy transfer phenomena specific for BWRs have been developed and implemented in the CFD-BWR module. A set of experiments focused on two-phase flow and phase-change phenomena has been identified for the validation of the CFD-BWR code and results of two experiment analyses focused on the radial void distribution are presented. The close agreement between the computed results, the measured data and the correlation results provides confidence in the accuracy of the models. (authors)
Modeling Tools Predict Flow in Fluid Dynamics
NASA Technical Reports Server (NTRS)
2010-01-01
"Because rocket engines operate under extreme temperature and pressure, they present a unique challenge to designers who must test and simulate the technology. To this end, CRAFT Tech Inc., of Pipersville, Pennsylvania, won Small Business Innovation Research (SBIR) contracts from Marshall Space Flight Center to develop software to simulate cryogenic fluid flows and related phenomena. CRAFT Tech enhanced its CRUNCH CFD (computational fluid dynamics) software to simulate phenomena in various liquid propulsion components and systems. Today, both government and industry clients in the aerospace, utilities, and petrochemical industries use the software for analyzing existing systems as well as designing new ones."
Computational fluid dynamic modelling of cavitation
NASA Technical Reports Server (NTRS)
Deshpande, Manish; Feng, Jinzhang; Merkle, Charles L.
1993-01-01
Models in sheet cavitation in cryogenic fluids are developed for use in Euler and Navier-Stokes codes. The models are based upon earlier potential-flow models but enable the cavity inception point, length, and shape to be determined as part of the computation. In the present paper, numerical solutions are compared with experimental measurements for both pressure distribution and cavity length. Comparisons between models are also presented. The CFD model provides a relatively simple modification to an existing code to enable cavitation performance predictions to be included. The analysis also has the added ability of incorporating thermodynamic effects of cryogenic fluids into the analysis. Extensions of the current two-dimensional steady state analysis to three-dimensions and/or time-dependent flows are, in principle, straightforward although geometrical issues become more complicated. Linearized models, however offer promise of providing effective cavitation modeling in three-dimensions. This analysis presents good potential for improved understanding of many phenomena associated with cavity flows.
Computational fluid dynamics modeling for emergency preparedness and response
Lee, R.L.; Albritton, J.R.; Ermak, D.L.; Kim, J.
1995-02-01
Computational fluid dynamics (CFD) has (CFD) has played an increasing in the improvement of atmospheric dispersion modeling. This is because many dispersion models are now driven by meteorological fields generated from CFD models or, in numerical weather prediction`s terminology, prognostic models. Whereas most dispersion models typically involve one or a few scalar, uncoupled equations, the prognostic equations are a set of highly-couple equations whose solution requires a significant level of computational power. Recent advances in computer hardware and software have enabled modestly-priced, high performance, workstations to exhibit the equivalent computation power of some mainframes. Thus desktop-class machines that were limited to performing dispersion calculations driven by diagnostic wind fields may now be used to calculate complex flows using prognostic CFD models. The Release and Advisory Capability (ARAC) program at Lawrence Livermore National Laboratory (LLNL) has, for the past several years, taken advantage of the improvements in hardware technology to develop a national emergency response capability based on executing diagnostic models on workstations. Diagnostic models that provide wind fields are, in general, simple to implement, robust and require minimal time for execution. Because these models typically contain little physics beyond mass-conservation, their performance is extremely sensitive to the quantity and quality of input meteorological data and, in spite of their utility, can be applied with confidence to only modestly complex flows. We are now embarking on a development program to incorporate prognostic models to generate, in real-time, the meteorological fields for the dispersion models. In contrast to diagnostic models, prognostic models are physically-based and are capable of incorporating many physical processes to treat highly complex flow scenarios.
Particle hopping vs. fluid-dynamical models for traffic flow
Nagel, K.
1995-12-31
Although particle hopping models have been introduced into traffic science in the 19509, their systematic use has only started recently. Two reasons for this are, that they are advantageous on modem computers, and that recent theoretical developments allow analytical understanding of their properties and therefore more confidence for their use. In principle, particle hopping models fit between microscopic models for driving and fluiddynamical models for traffic flow. In this sense, they also help closing the conceptual gap between these two. This paper shows connections between particle hopping models and traffic flow theory. It shows that the hydrodynamical limits of certain particle hopping models correspond to the Lighthill-Whitham theory for traffic flow, and that only slightly more complex particle hopping models produce already the correct traffic jam dynamics, consistent with recent fluid-dynamical models for traffic flow. By doing so, this paper establishes that, on the macroscopic level, particle hopping models are at least as good as fluid-dynamical models. Yet, particle hopping models have at least two advantages over fluid-dynamical models: they straightforwardly allow microscopic simulations, and they include stochasticity.
Computational fluid dynamics modeling for emergency preparedness & response
Lee, R.L.; Albritton, J.R.; Ermak, D.L.; Kim, J.
1995-07-01
Computational fluid dynamics (CFD) has played an increasing role in the improvement of atmospheric dispersion modeling. This is because many dispersion models are now driven by meteorological fields generated from CFD models or, in numerical weather prediction`s terminology, prognostic models. Whereas most dispersion models typically involve one or a few scalar, uncoupled equations, the prognostic equations are a set of highly-coupled, nonlinear equations whose solution requires a significant level of computational power. Until recently, such computer power could be found only in CRAY-class supercomputers. Recent advances in computer hardware and software have enabled modestly-priced, high performance, workstations to exhibit the equivalent computation power of some mainframes. Thus desktop-class machines that were limited to performing dispersion calculations driven by diagnostic wind fields may now be used to calculate complex flows using prognostic CFD models. The Atmospheric Release and Advisory Capability (ARAC) program at Lawrence Livermore National Laboratory (LLNL) has, for the past several years, taken advantage of the improvements in hardware technology to develop a national emergency response capability based on executing diagnostic models on workstations. Diagnostic models that provide wind fields are, in general, simple to implement, robust and require minimal time for execution. Such models have been the cornerstones of the ARAC operational system for the past ten years. Kamada (1992) provides a review of diagnostic models and their applications to dispersion problems. However, because these models typically contain little physics beyond mass-conservation, their performance is extremely sensitive to the quantity and quality of input meteorological data and, in spite of their utility, can be applied with confidence to only modestly complex flows.
Modeling quantum fluid dynamics at nonzero temperatures
Berloff, Natalia G.; Brachet, Marc; Proukakis, Nick P.
2014-01-01
The detailed understanding of the intricate dynamics of quantum fluids, in particular in the rapidly growing subfield of quantum turbulence which elucidates the evolution of a vortex tangle in a superfluid, requires an in-depth understanding of the role of finite temperature in such systems. The Landau two-fluid model is the most successful hydrodynamical theory of superfluid helium, but by the nature of the scale separations it cannot give an adequate description of the processes involving vortex dynamics and interactions. In our contribution we introduce a framework based on a nonlinear classical-field equation that is mathematically identical to the Landau model and provides a mechanism for severing and coalescence of vortex lines, so that the questions related to the behavior of quantized vortices can be addressed self-consistently. The correct equation of state as well as nonlocality of interactions that leads to the existence of the roton minimum can also be introduced in such description. We review and apply the ideas developed for finite-temperature description of weakly interacting Bose gases as possible extensions and numerical refinements of the proposed method. We apply this method to elucidate the behavior of the vortices during expansion and contraction following the change in applied pressure. We show that at low temperatures, during the contraction of the vortex core as the negative pressure grows back to positive values, the vortex line density grows through a mechanism of vortex multiplication. This mechanism is suppressed at high temperatures. PMID:24704874
Green Algae as Model Organisms for Biological Fluid Dynamics
NASA Astrophysics Data System (ADS)
Goldstein, Raymond E.
2015-01-01
In the past decade, the volvocine green algae, spanning from the unicellular Chlamydomonas to multicellular Volvox, have emerged as model organisms for a number of problems in biological fluid dynamics. These include flagellar propulsion, nutrient uptake by swimming organisms, hydrodynamic interactions mediated by walls, collective dynamics and transport within suspensions of microswimmers, the mechanism of phototaxis, and the stochastic dynamics of flagellar synchronization. Green algae are well suited to the study of such problems because of their range of sizes (from 10 μm to several millimeters), their geometric regularity, the ease with which they can be cultured, and the availability of many mutants that allow for connections between molecular details and organism-level behavior. This review summarizes these recent developments and highlights promising future directions in the study of biological fluid dynamics, especially in the context of evolutionary biology, that can take advantage of these remarkable organisms.
Green Algae as Model Organisms for Biological Fluid Dynamics*
Goldstein, Raymond E.
2015-01-01
In the past decade the volvocine green algae, spanning from the unicellular Chlamydomonas to multicellular Volvox, have emerged as model organisms for a number of problems in biological fluid dynamics. These include flagellar propulsion, nutrient uptake by swimming organisms, hydrodynamic interactions mediated by walls, collective dynamics and transport within suspensions of microswimmers, the mechanism of phototaxis, and the stochastic dynamics of flagellar synchronization. Green algae are well suited to the study of such problems because of their range of sizes (from 10 μm to several millimetres), their geometric regularity, the ease with which they can be cultured and the availability of many mutants that allow for connections between molecular details and organism-level behavior. This review summarizes these recent developments and highlights promising future directions in the study of biological fluid dynamics, especially in the context of evolutionary biology, that can take advantage of these remarkable organisms. PMID:26594068
AIR INGRESS ANALYSIS: PART 2 – COMPUTATIONAL FLUID DYNAMIC MODELS
Chang H. Oh; Eung S. Kim; Richard Schultz; Hans Gougar; David Petti; Hyung S. Kang
2011-01-01
The Idaho National Laboratory (INL), under the auspices of the U.S. Department of Energy, is performing research and development that focuses on key phenomena important during potential scenarios that may occur in very high temperature reactors (VHTRs). Phenomena Identification and Ranking Studies to date have ranked an air ingress event, following on the heels of a VHTR depressurization, as important with regard to core safety. Consequently, the development of advanced air ingress-related models and verification and validation data are a very high priority. Following a loss of coolant and system depressurization incident, air will enter the core of the High Temperature Gas Cooled Reactor through the break, possibly causing oxidation of the in-the core and reflector graphite structure. Simple core and plant models indicate that, under certain circumstances, the oxidation may proceed at an elevated rate with additional heat generated from the oxidation reaction itself. Under postulated conditions of fluid flow and temperature, excessive degradation of the lower plenum graphite can lead to a loss of structural support. Excessive oxidation of core graphite can also lead to the release of fission products into the confinement, which could be detrimental to a reactor safety. Computational fluid dynamic model developed in this study will improve our understanding of this phenomenon. This paper presents two-dimensional and three-dimensional CFD results for the quantitative assessment of the air ingress phenomena. A portion of results of the density-driven stratified flow in the inlet pipe will be compared with results of the experimental results.
Advances in computational fluid dynamics solvers for modern computing environments
NASA Astrophysics Data System (ADS)
Hertenstein, Daniel; Humphrey, John R.; Paolini, Aaron L.; Kelmelis, Eric J.
2013-05-01
EM Photonics has been investigating the application of massively multicore processors to a key problem area: Computational Fluid Dynamics (CFD). While the capabilities of CFD solvers have continually increased and improved to support features such as moving bodies and adjoint-based mesh adaptation, the software architecture has often lagged behind. This has led to poor scaling as core counts reach the tens of thousands. In the modern High Performance Computing (HPC) world, clusters with hundreds of thousands of cores are becoming the standard. In addition, accelerator devices such as NVIDIA GPUs and Intel Xeon Phi are being installed in many new systems. It is important for CFD solvers to take advantage of the new hardware as the computations involved are well suited for the massively multicore architecture. In our work, we demonstrate that new features in NVIDIA GPUs are able to empower existing CFD solvers by example using AVUS, a CFD solver developed by the Air Force Research Labratory (AFRL) and the Volcanic Ash Advisory Center (VAAC). The effort has resulted in increased performance and scalability without sacrificing accuracy. There are many well-known codes in the CFD space that can benefit from this work, such as FUN3D, OVERFLOW, and TetrUSS. Such codes are widely used in the commercial, government, and defense sectors.
Supersonic diode pumped alkali lasers: Computational fluid dynamics modeling
NASA Astrophysics Data System (ADS)
Rosenwaks, Salman; Yacoby, Eyal; Waichman, Karol; Sadot, Oren; Barmashenko, Boris D.
2015-10-01
We report on recent progress on our three-dimensional computational fluid dynamics (3D CFD) modeling of supersonic diode pumped alkali lasers (DPALs), taking into account fluid dynamics and kinetic processes in the lasing medium. For a supersonic Cs DPAL with laser section geometry and resonator parameters similar to those of the 1-kW flowing-gas subsonic Cs DPAL [A.V. Bogachev et al., Quantum Electron. 42, 95 (2012)] the maximum achievable output power, ~ 7 kW, is 25% higher than that achievable in the subsonic case. Comparison between semi-analytical and 3D CFD models for Cs shows that the latter predicts much higher maximum achievable output power than the former. Optimization of the laser parameters using 3D CFD modeling shows that very high power and optical-to-optical efficiency, 35 kW and 82%, respectively, can be achieved in a Cs supersonic device pumped by a collimated cylindrical (0.5 cm diameter) beam. Application of end- or transverse-pumping by collimated rectangular (large cross section ~ 2 - 4 cm2) beam makes it possible to obtain even higher output power, > 250 kW, for ~ 350 kW pumping power. The main processes limiting the power of Cs supersonic DPAL are saturation of the D2 transition and large ~ 40% losses of alkali atoms due to ionization, whereas the influence of gas heating is negligibly small. For supersonic K DPAL both gas heating and ionization effects are shown to be unimportant and the maximum achievable power, ~ 40 kW and 350 kW, for pumping by ~ 100 kW cylindrical and ~ 700 kW rectangular beam, respectively, are higher than those achievable in the Cs supersonic laser. The power achieved in the supersonic K DPAL is two times higher than for the subsonic version with the same resonator and K density at the gas inlet, the maximum optical-to-optical efficiency being 82%.
Greg Weirs; Hyung Lee
2011-09-01
V&V and UQ are the primary means to assess the accuracy and reliability of M&S and, hence, to establish confidence in M&S. Though other industries are establishing standards and requirements for the performance of V&V and UQ, at present, the nuclear industry has not established such standards or requirements. However, the nuclear industry is beginning to recognize that such standards are needed and that the resources needed to support V&V and UQ will be very significant. In fact, no single organization has sufficient resources or expertise required to organize, conduct and maintain a comprehensive V&V and UQ program. What is needed is a systematic and standardized approach to establish and provide V&V and UQ resources at a national or even international level, with a consortium of partners from government, academia and industry. Specifically, what is needed is a structured and cost-effective knowledge base that collects, evaluates and stores verification and validation data, and shows how it can be used to perform V&V and UQ, leveraging collaboration and sharing of resources to support existing engineering and licensing procedures as well as science-based V&V and UQ processes. The Nuclear Energy Knowledge base for Advanced Modeling and Simulation (NE-KAMS) is being developed at the Idaho National Laboratory in conjunction with Bettis Laboratory, Sandia National Laboratories, Argonne National Laboratory, Utah State University and others with the objective of establishing a comprehensive and web-accessible knowledge base to provide V&V and UQ resources for M&S for nuclear reactor design, analysis and licensing. The knowledge base will serve as an important resource for technical exchange and collaboration that will enable credible and reliable computational models and simulations for application to nuclear power. NE-KAMS will serve as a valuable resource for the nuclear industry, academia, the national laboratories, the U.S. Nuclear Regulatory Commission (NRC) and
Benchmarking computational fluid dynamics models for lava flow simulation
NASA Astrophysics Data System (ADS)
Dietterich, Hannah; Lev, Einat; Chen, Jiangzhi
2016-04-01
Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, and COMSOL. Using the new benchmark scenarios defined in Cordonnier et al. (Geol Soc SP, 2015) as a guide, we model viscous, cooling, and solidifying flows over horizontal and sloping surfaces, topographic obstacles, and digital elevation models of natural topography. We compare model results to analytical theory, analogue and molten basalt experiments, and measurements from natural lava flows. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We can apply these models to reconstruct past lava flows in Hawai'i and Saudi Arabia using parameters assembled from morphology, textural analysis, and eruption observations as natural test cases. Our study highlights the strengths and weaknesses of each code, including accuracy and computational costs, and provides insights regarding code selection.
Modelling reacting localized air pollution using Computational Fluid Dynamics (CFD)
NASA Astrophysics Data System (ADS)
Karim, A. A.; Nolan, P. F.
2011-02-01
A Computational Fluid Dynamics (CFD) approach is used to model reacting NO 2 dispersion of vehicle pollutants released from a dual carriageway in Maidstone, UK. The simulations are carried out over the course of one full day during January, 2008. The developed CFD model utilizes a modified k- ɛ turbulence model and Arrhenius reaction kinetics with source terms for the reactions which include a photo-stationary set with peroxy radicals. An approach is taken whereby the reactions are solved specific to the rush hour period corresponding to the availability of certain hydrocarbons released from the vehicles. The results of the simulation are compared with field measurements taken at the site which is made up of several, different sized buildings on varying terrain in Maidstone UK. The predictions and field measurements are considered over a 12 h period with averaged hourly results. It was found that the reactive pollutant approach greatly improves the predictions as compared to the experiments. Furthermore the effect of peroxy radicals during rush hour periods is found to be a major disturbance to the photo-stationary set and its inclusion improved the predictions further.
Computational fluid dynamics framework for aerodynamic model assessment
NASA Astrophysics Data System (ADS)
Vallespin, D.; Badcock, K. J.; Da Ronch, A.; White, M. D.; Perfect, P.; Ghoreyshi, M.
2012-07-01
This paper reviews the work carried out at the University of Liverpool to assess the use of CFD methods for aircraft flight dynamics applications. Three test cases are discussed in the paper, namely, the Standard Dynamic Model, the Ranger 2000 jet trainer and the Stability and Control Unmanned Combat Air Vehicle. For each of these, a tabular aerodynamic model based on CFD predictions is generated along with validation against wind tunnel experiments and flight test measurements. The main purpose of the paper is to assess the validity of the tables of aerodynamic data for the force and moment prediction of realistic aircraft manoeuvres. This is done by generating a manoeuvre based on the tables of aerodynamic data, and then replaying the motion through a time-accurate computational fluid dynamics calculation. The resulting forces and moments from these simulations were compared with predictions from the tables. As the latter are based on a set of steady-state predictions, the comparisons showed perfect agreement for slow manoeuvres. As manoeuvres became more aggressive some disagreement was seen, particularly during periods of large rates of change in attitudes. Finally, the Ranger 2000 model was used on a flight simulator.
Computational Fluid Dynamics Model for Saltstone Vault 4 Vapor Sapce
Lee, Si Young
2005-06-27
Computational fluid dynamics (CFD) methods have been used to estimate the flow patterns for vapor space inside the Saltstone Vault No.4 under different operating scenarios. The purpose of this work is to examine the gas motions inside the vapor space under the current vault configurations. A CFD model took three-dimensional transient momentum-energy coupled approach for the vapor space domain of the vault. The modeling calculations were based on prototypic vault geometry and expected normal operating conditions as defined by Waste Solidification Engineering. The modeling analysis was focused on the air flow patterns near the ventilated corner zones of the vapor space inside the Saltstone vault. The turbulence behavior and natural convection mechanism used in the present model were benchmarked against the literature information and theoretical results. The verified model was applied to the Saltstone vault geometry for the transient assessment of the air flow patterns inside the vapor space of the vault region using the boundary conditions as provided by the customer. The present model considered two cases for the estimations of the flow patterns within the vapor space. One is the reference baseline case. The other is for the negative temperature gradient between the roof inner and top grout surface temperatures intended for the potential bounding condition. The flow patterns of the vapor space calculated by the CFD model demonstrate that the ambient air comes into the vapor space of the vault through the lower-end ventilation hole, and it gets heated up by the Benard-cell type circulation before leaving the vault via the higher-end ventilation hole. The calculated results are consistent with the literature information.
Computational fluid dynamic modeling of fluidized-bed polymerization reactors
Rokkam, Ram
2012-01-01
Polyethylene is one of the most widely used plastics, and over 60 million tons are produced worldwide every year. Polyethylene is obtained by the catalytic polymerization of ethylene in gas and liquid phase reactors. The gas phase processes are more advantageous, and use fluidized-bed reactors for production of polyethylene. Since they operate so close to the melting point of the polymer, agglomeration is an operational concern in all slurry and gas polymerization processes. Electrostatics and hot spot formation are the main factors that contribute to agglomeration in gas-phase processes. Electrostatic charges in gas phase polymerization fluidized bed reactors are known to influence the bed hydrodynamics, particle elutriation, bubble size, bubble shape etc. Accumulation of electrostatic charges in the fluidized-bed can lead to operational issues. In this work a first-principles electrostatic model is developed and coupled with a multi-fluid computational fluid dynamic (CFD) model to understand the effect of electrostatics on the dynamics of a fluidized-bed. The multi-fluid CFD model for gas-particle flow is based on the kinetic theory of granular flows closures. The electrostatic model is developed based on a fixed, size-dependent charge for each type of particle (catalyst, polymer, polymer fines) phase. The combined CFD model is first verified using simple test cases, validated with experiments and applied to a pilot-scale polymerization fluidized-bed reactor. The CFD model reproduced qualitative trends in particle segregation and entrainment due to electrostatic charges observed in experiments. For the scale up of fluidized bed reactor, filtered models are developed and implemented on pilot scale reactor.
This paper discusses the status and application of Computational Fluid Dynamics (CFD) models to address challenges for modeling human exposures to air pollutants around urban building microenvironments. There are challenges for more detailed understanding of air pollutant sour...
A FRAMEWORK FOR FINE-SCALE COMPUTATIONAL FLUID DYNAMICS AIR QUALITY MODELING AND ANALYSIS
Fine-scale Computational Fluid Dynamics (CFD) simulation of pollutant concentrations within roadway and building microenvironments is feasible using high performance computing. Unlike currently used regulatory air quality models, fine-scale CFD simulations are able to account rig...
Advanced feedback control of indoor air quality using real-time computational fluid dynamics
Ratnam, E.; Campbell, T.; Bradley, R.
1998-10-01
This paper describes the partial implementation of a novel method of controlling indoor air quality (IAQ) for critical applications. The proposed method uses a numerical modeling technique known as computational fluid dynamics (CFD) for modeling the effect of variable ventilation rates for intelligent and rapid control of air contamination in space. This paper describes how a CFD model is made to run in real time linked to a feedback control loop. The technique was simulated in a graphical programming language. The simulation results indicate that a quasi-transient potential flow CFD model is a viable technique for feedback control of IAQ, and it is currently being implemented in an experimental validation.
Computational fluid dynamics models and congenital heart diseases
Pennati, Giancarlo; Corsini, Chiara; Hsia, Tain-Yen; Migliavacca, Francesco
2013-01-01
Mathematical modeling is a powerful tool to investigate hemodynamics of the circulatory system. With improving imaging techniques and detailed clinical investigations, it is now possible to construct patient-specific models of reconstructive surgeries for the treatment of congenital heart diseases. These models can help clinicians to better understand the hemodynamic behavior of different surgical options for a treated patient. This review outlines recent advances in mathematical modeling in congenital heart diseases, the discoveries and limitations these models present, and future directions that are on the horizon. PMID:24432298
Simulation of Tailrace Hydrodynamics Using Computational Fluid Dynamics Models
Cook, Christopher B.; Richmond, Marshall C.
2001-05-01
This report investigates the feasibility of using computational fluid dynamics (CFD) tools to investigate hydrodynamic flow fields surrounding the tailrace zone below large hydraulic structures. Previous and ongoing studies using CFD tools to simulate gradually varied flow with multiple constituents and forebay/intake hydrodynamics have shown that CFD tools can provide valuable information for hydraulic and biological evaluation of fish passage near hydraulic structures. These studies however are incapable of simulating the rapidly varying flow fields that involving breakup of the free-surface, such as those through and below high flow outfalls and spillways. Although the use of CFD tools for these types of flow are still an active area of research, initial applications discussed in this report show that these tools are capable of simulating the primary features of these highly transient flow fields.
Sittitavornwong, Somsak; Waite, Peter D.; Shih, Alan M.; Cheng, Gary C.; Koomullil, Roy; Ito, Yasushi; Cure, Joel K; Harding, Susan M.; Litaker, Mark
2013-01-01
Purpose Evaluate the soft tissue change of the upper airway after maxillomandibular advancement (MMA) by computational fluid dynamics (CFD). Materials and Methods Eight OSAS patients who required MMA were recruited into this study. All participants had pre- and post-operative computed tomography (CT) and underwent MMA by a single oral and maxillofacial surgeon. Upper airway CT data sets for these 8 participants were created with high-fidelity 3-D numerical models for computational fluid dynamics (CFD). The 3-D models were simulated and analyzed to study how changes in airway anatomy affects pressure effort required for normal breathing. Airway dimensions, skeletal changes, Apnea-Hypopnea Index (AHI), and pressure efforts of pre- and post-operative 3-D models were compared and correlations interpreted. Results After MMA, laminar and turbulent air flow was significantly decreased at every level of the airway. The cross-sectional areas at the soft palate and tongue base were significantly increased. Conclusions This study shows that MMA increases airway dimensions by the increasing the occipital base (Base) - pogonion (Pg) distance. An increase of the Base-Pg distance showed a significant correlation with an AHI improvement and a decreased pressure effort of the upper airway. Decreasing the pressure effort will decrease the breathing workload. This improves the condition of OSAS. PMID:23642544
NASA Astrophysics Data System (ADS)
Emelyanov, V. N.; Karpenko, A. G.; Volkov, K. N.
2015-06-01
Modern graphics processing units (GPU) provide architectures and new programming models that enable to harness their large processing power and to design computational fluid dynamics (CFD) simulations at both high performance and low cost. Possibilities of the use of GPUs for the simulation of internal fluid flows are discussed. The finite volume method is applied to solve three-dimensional (3D) unsteady compressible Euler and Navier-Stokes equations on unstructured meshes. Compute Inified Device Architecture (CUDA) technology is used for programming implementation of parallel computational algorithms. Solution of some fluid dynamics problems on GPUs is presented and approaches to optimization of the CFD code related to the use of different types of memory are discussed. Speedup of solution on GPUs with respect to the solution on central processor unit (CPU) is compared with the use of different meshes and different methods of distribution of input data into blocks. Performance measurements show that numerical schemes developed achieve 20 to 50 speedup on GPU hardware compared to CPU reference implementation. The results obtained provide promising perspective for designing a GPU-based software framework for applications in CFD.
A FRAMEWORK FOR FINE-SCALE COMPUTATIONAL FLUID DYNAMICS AIR QUALITY MODELING AND ANALYSIS
This paper discusses a framework for fine-scale CFD modeling that may be developed to complement the present Community Multi-scale Air Quality (CMAQ) modeling system which itself is a computational fluid dynamics model. A goal of this presentation is to stimulate discussions on w...
REMOVAL OF TANK AND SEWER SEDIMENT BY GATE FLUSHING: COMPUTATIONAL FLUID DYNAMICS MODEL STUDIES
This presentation will discuss the application of a computational fluid dynamics 3D flow model to simulate gate flushing for removing tank/sewer sediments. The physical model of the flushing device was a tank fabricated and installed at the head-end of a hydraulic flume. The fl...
Applicability of fluid-dynamical modeling of nucleus-nucleus collisions at relativistic energies
NASA Astrophysics Data System (ADS)
Hazineh, Dean; Auvinen, Jussi; Nahrgang, Marlene; Bass, Steffen
2015-10-01
At sufficiently high temperatures and densities, similar to the conditions found in the early universe, QCD matter forms a deconfined state called the quark gluon plasma (QGP). This state of matter can be created in collisions of ultra-relativistic heavy-ions, and RHIC data suggests that this QGP behaves similar to an ideal fluid. Viscous relativistic fluid dynamics therefore is one of the preferred theoretical tools to model the time-evolution and properties of the QGP. As the collision energy or the system size is decreased, the range of applicability of viscous fluid dynamics becomes smaller as the length scale of the interaction among the basic constituents is similar to the overall scale of the collision system itself. In order to investigate the validity of fluid-dynamical modeling of proton-nucleus and nucleus-nucleus collisions at LHC and RHIC, we conduct an analysis of the spatial and temporal evolution of the Knudsen number, i.e. the ratio of the microscopic mean free path to the macroscopic length scale of the system. We show results for large and small collision systems, as a function of the specific shear viscosity, and discuss the range of applicability of fluid-dynamical modeling in relativistic proton-nucleus and nucleus-nucleus collisions at different energies.
Modeling the Fluid Dynamics in a Human Stomach to Gain Insight of Food Digestion
Ferrua, MJ; Singh, RP
2010-01-01
During gastric digestion, food is disintegrated by a complex interaction of chemical and mechanical effects. Although the mechanisms of chemical digestion are usually characterized by using in vitro analysis, the difficulty in reproducing the stomach geometry and motility has prevented a good understanding of the local fluid dynamics of gastric contents. The goal of this study was to use computational fluid dynamics (CFD) to develop a 3-D model of the shape and motility pattern of the stomach wall during digestion, and use it to characterize the fluid dynamics of gastric contents of different viscosities. A geometrical model of an averaged-sized human stomach was created, and its motility was characterized by a series of antral-contraction waves of up to 80% relative occlusion. The flow field within the model (predicted using the software Fluent™) strongly depended on the viscosity of gastric contents. By increasing the viscosity, the formation of the 2 flow patterns commonly regarded as the main mechanisms driving digestion (i.e., the retropulsive jet-like motion and eddy structures) was significantly diminished, while a significant increase of the pressure field was predicted. These results were in good agreement with experimental data previously reported in the literature, and suggest that, contrary to the traditional idea of a rapid and complete homogenization of the meal, gastric contents associated with high viscous meals are poorly mixed. This study illustrates the capability of CFD to provide a unique insight into the fluid dynamics of the gastric contents, and points out its potential to develop a fundamental understanding and modeling of the mechanisms involved in the digestion process. Practical Application This study illustrates the capability of computational fluid dynamic techniques to provide a unique insight into the dynamics of the gastric contents, pointing out its potential to develop a fundamental understanding and modeling of the human
Computational Fluid Dynamics Modeling of Nickel Hydrogen Batteries
NASA Technical Reports Server (NTRS)
Cullion, R.; Gu, W. B.; Wang, C. Y.; Timmerman, P.
2000-01-01
An electrochemical Ni-H2 battery model has been expanded to include thermal effects. A thermal energy conservation equation was derived from first principles. An electrochemical and thermal coupled model was created by the addition of this equation to an existing multiphase, electrochemical model. Charging at various rates was investigated and the results validated against experimental data. Reaction currents, pressure changes, temperature profiles, and concentration variations within the cell are predicted numerically and compared with available data and theory.
Zhang, Jiafeng; Zhang, Pei; Fraser, Katharine H.; Griffith, Bartley P.; Wu, Zhongjun J.
2012-01-01
With the recent advances in computer technology, computational fluid dynamics (CFD) has become an important tool to design and improve blood contacting artificial organs, and to study the device-induced blood damage. Commercial CFD software packages are readily available, and multiple CFD models are provided by CFD software developers. However, the best approach of using CFD effectively to characterize fluid flow and to predict blood damage in these medical devices remains debatable. This study aimed to compare these CFD models and provide useful information on the accuracy of each model in modeling blood flow in circulatory assist devices. The laminar and five turbulence models (Spalart-Allmaras, k-ε (k-epsilon), k-ω (k-omega), SST (Menter’s Shear Stress Transport), and Reynolds Stress) were implemented to predict blood flow in a clinically used circulatory assist device, CentriMag® centrifugal blood pump (Thoratec, MA). In parallel, a transparent replica of the CentriMag® pump was constructed and selected views of the flow fields were measured with digital particle image velocimetry (DPIV). CFD results were compared with the DPIV experimental results. Compared with the experiment, all the selected CFD models predicted the flow pattern fairly well except the area of the outlet. However, quantitatively, the laminar model results were the most deviated from the experimental data. On the other hand, k-ε RNG models and Reynolds Stress model are the most accurate. In conclusion, for the circulatory assist devices, turbulence models provide more accurate results than laminar model. Among the selected turbulence models, k-ε and Reynolds Stress Method models are recommended. PMID:23441681
Zhang, Jiafeng; Zhang, Pei; Fraser, Katharine H; Griffith, Bartley P; Wu, Zhongjun J
2013-04-01
With the recent advances in computer technology, computational fluid dynamics (CFDs) has become an important tool to design and improve blood-contacting artificial organs, and to study the device-induced blood damage. Commercial CFD software packages are readily available, and multiple CFD models are provided by CFD software developers. However, the best approach of using CFD effectively to characterize fluid flow and to predict blood damage in these medical devices remains debatable. This study aimed to compare these CFD models and provide useful information on the accuracy of each model in modeling blood flow in circulatory assist devices. The laminar and five turbulence models (Spalart-Allmaras, k-ε (k-epsilon), k-ω (k-omega), SST [Menter's Shear Stress Transport], and Reynolds Stress) were implemented to predict blood flow in a clinically used circulatory assist device, the CentriMag centrifugal blood pump. In parallel, a transparent replica of the CentriMag pump was constructed and selected views of the flow fields were measured with digital particle image velocimetry (DPIV). CFD results were compared with the DPIV experimental results. Compared with the experiment, all the selected CFD models predicted the flow pattern fairly well except the area of the outlet. However, quantitatively, the laminar model results were the most deviated from the experimental data. On the other hand, k-ε renormalization group theory models and Reynolds Stress model are the most accurate. In conclusion, for the circulatory assist devices, turbulence models provide more accurate results than the laminar model. Among the selected turbulence models, k-ε and Reynolds Stress Method models are recommended. PMID:23441681
Phantom model of physiologic intracranial pressure and cerebrospinal fluid dynamics.
Bottan, Simone; Poulikakos, Dimos; Kurtcuoglu, Vartan
2012-06-01
We describe herein a novel life-size phantom model of the intracranial cavity and its validation. The cerebrospinal fluid (CSF) domains including ventricular, cysternal, and subarachnoid spaces were derived via magnetic resonance imaging. Brain mechanical properties and cranio-spinal compliance were set based on published data. Both bulk and pulsatile physiologic CSF flow were modeled. Model validation was carried out by comparisons of flow and pressure measurements in the phantom with published in vivo data of healthy subjects. Physiologic intracranial pressure with 10 mmHg mean and 0.4 mmHg peak pulse amplitude was recorded in the ventricles. Peak CSF flow rates of 0.2 and 2 ml/s were measured in the cerebral aqueduct and subarachnoid space, respectively. The phantom constitutes a first-of-its-kind approach to modeling physiologic intracranial dynamics in vitro. Herein, we describe the phantom design and manufacturing, definition and implementation of its operating parameters, as well as the validation of the modeled dynamics. PMID:22333981
Mathematical modeling of fluid dynamics in pulsatile cardiopulmonary bypass.
Pennati, Giancarlo; Fiore, Gianfranco B; Laganà, Katia; Fumero, Roberto
2004-02-01
The design criteria of an extracorporeal circuit suitable for pulsatile flow are quite different and more entangled than for steady flow. The time and costs of the design process could be reduced if mutual influences between the pulsatile pump and other extracorporeal devices were considered without experimental trial-and-error activities. With this in mind, we have developed a new lumped-parameter mathematical model of the hydraulic behavior of the arterial side of an extracorporeal circuit under pulsatile flow conditions. Generally, components feature a resistant-inertant-compliant behavior and the most relevant nonlinearities are accounted for. Parameter values were derived either by experimental tests or by analytical analysis. The pulsatile pump is modeled as a pure pulsatile flow generator. Model predictions were compared with flow rate and pressure tracings measured during hydraulic tests on two different circuits at various flow rates and pulse frequencies. The normalized root mean square error did not exceed 24% and the model accurately describes the changes that occur in the basic features of the pressure and flow wave propagating from the pulsatile pump to the arterial cannula. PMID:14961960
Computational Fluid Dynamics Modeling of the John Day Dam Tailrace
Rakowski, Cynthia L.; Perkins, William A.; Richmond, Marshall C.; Serkowski, John A.
2010-07-08
US Army Corps of Engineers - Portland District required that a two-dimensional (2D) depth-averaged and a three-dimensional (3D) free-surface numerical models to be developed and validated for the John Day tailrace. These models were used to assess potential impact of a select group of structural and operational alternatives to tailrace flows aimed at improving fish survival at John Day Dam. The 2D model was used for the initial assessment of the alternatives in conjunction with a reduced-scale physical model of the John Day Project. A finer resolution 3D model was used to more accurately model the details of flow in the stilling basin and near-project tailrace hydraulics. Three-dimensional model results were used as input to the Pacific Northwest National Laboratory particle tracking software, and particle paths and times to pass a downstream cross section were used to assess the relative differences in travel times resulting from project operations and structural scenarios for multiple total river flows. Streamlines and neutrally-buoyant particles were seeded in all turbine and spill bays with flows. For a Total River of 250 kcfs running with the Fish Passage Plan spill pattern and a spillwall, the mean residence times for all particles were little changed; however the tails of the distribution were truncated for both spillway and powerhouse release points, and, for the powerhouse releases, reduced the residence time for 75% of the particles to pass a downstream cross section from 45.5 minutes to 41.3 minutes. For a total river of 125 kcfs configured with the operations from the Fish Passage Plan for the temporary spillway weirs and for a proposed spillwall, the neutrally-buoyant particle tracking data showed that the river with a spillwall in place had the overall mean residence time increase; however, the residence time for 75% of the powerhouse-released particles to pass a downstream cross section was reduced from 102.4 min to 89 minutes.
Singular Limits of Voigt Models in Fluid Dynamics
NASA Astrophysics Data System (ADS)
Coti Zelati, Michele; Gal, Ciprian G.
2015-06-01
We investigate the long-term behavior, as a certain regularization parameter vanishes, of the three-dimensional Navier-Stokes-Voigt model of a viscoelastic incompressible fluid. We prove the existence of global and exponential attractors of optimal regularity. We then derive explicit upper bounds for the dimension of these attractors in terms of the three-dimensional Grashof number and the regularization parameter. Finally, we also prove convergence of the (strong) global attractor of the 3D Navier-Stokes-Voigt model to the (weak) global attractor of the 3D Navier-Stokes equation. Our analysis improves and extends recent results obtained by Kalantarov and Titi (Chin Ann Math Ser B 30:697-714, 2009).
Fluid dynamic modeling of nano-thermite reactions
NASA Astrophysics Data System (ADS)
Martirosyan, Karen S.; Zyskin, Maxim; Jenkins, Charles M.; Yuki Horie, Yasuyuki
2014-03-01
This paper presents a direct numerical method based on gas dynamic equations to predict pressure evolution during the discharge of nanoenergetic materials. The direct numerical method provides for modeling reflections of the shock waves from the reactor walls that generates pressure-time fluctuations. The results of gas pressure prediction are consistent with the experimental evidence and estimates based on the self-similar solution. Artificial viscosity provides sufficient smoothing of shock wave discontinuity for the numerical procedure. The direct numerical method is more computationally demanding and flexible than self-similar solution, in particular it allows study of a shock wave in its early stage of reaction and allows the investigation of "slower" reactions, which may produce weaker shock waves. Moreover, numerical results indicate that peak pressure is not very sensitive to initial density and reaction time, providing that all the material reacts well before the shock wave arrives at the end of the reactor.
Development and mathematical modeling of biomagnetic fluid dynamics
NASA Astrophysics Data System (ADS)
Haik, Yousef
In this study the human blood as a biomagnetic fluid and the blood flow under the influence of magnetic field is investigated. Experiments were performed to examine the effect of magnetic field on blood flow. The blood flow in a straight pipe when aligned to the magnetic field lines was reduced by 28% when the field strength was 10 tesla. The sedimentation rate of red blood cells in the stagnant container due to gravity was reduced when the whole blood was subjected to a horizontally imposed magnetic field from a permanent magnet of 1.2 tesla. The physical properties of whole blood under the influence of high magnetic field were also investigated. It was found that the magnetic susceptibility of whole blood is in the order of 10-6 when deoxygenated and 10-7 when oxygenated. In the present study the whole blood is treated as a single fluid. The governing equations for the flow under magnetic field are derived based on the modified Stokes' principle. The model coefficients were obtained by performing experiments on the whole blood. The mathematical model is simplified based on the equilibrium approximation and the experimental observations. Numerical simulations were obtained for the situations such as acceleration and retardation of blood flow, removal of blood thrombus from the vein wall by increasing the shear rate and the behavior of the blood in the aneurysm in the vicinity of a magnetic field. It is shown that the blood flow characteristic can change dramatically when the fluid is subjected to magnetic field. Finally different biomagnetic applications are proposed. Among these proposed applications is a device to separate the red blood cells from the whole blood for possible treatment of the white cells. Another proposed application is the Magnetic MEMS. Magnetic micro delivery systems that utilizes the magnetic force to deliver a micro machine to a specific site and magnetic micro rotary systems that uses the magnetic torque to drive the micromachine are
JACKSON VL
2011-08-31
The primary purpose of the tank mixing and sampling demonstration program is to mitigate the technical risks associated with the ability of the Hanford tank farm delivery and celtification systems to measure and deliver a uniformly mixed high-level waste (HLW) feed to the Waste Treatment and Immobilization Plant (WTP) Uniform feed to the WTP is a requirement of 24590-WTP-ICD-MG-01-019, ICD-19 - Interface Control Document for Waste Feed, although the exact definition of uniform is evolving in this context. Computational Fluid Dynamics (CFD) modeling has been used to assist in evaluating scaleup issues, study operational parameters, and predict mixing performance at full-scale.
Parallel technology for numerical modeling of fluid dynamics problems by high-accuracy algorithms
NASA Astrophysics Data System (ADS)
Gorobets, A. V.
2015-04-01
A parallel computation technology for modeling fluid dynamics problems by finite-volume and finite-difference methods of high accuracy is presented. The development of an algorithm, the design of a software implementation, and the creation of parallel programs for computations on large-scale computing systems are considered. The presented parallel technology is based on a multilevel parallel model combining various types of parallelism: with shared and distributed memory and with multiple and single instruction streams to multiple data flows.
NASA Technical Reports Server (NTRS)
Noor, Ahmed K. (Editor); Dwoyer, Douglas L. (Editor)
1988-01-01
Recent advances in computational structural and fluid dynamics are discussed in reviews and reports. Topics addressed include fluid-structure interaction and aeroelasticity, CFD techniques for reacting flows, micromechanics, stability and eigenproblems, probabilistic methods and chaotic dynamics, and perturbation and spectral methods. Consideration is given to finite-element, finite-volume, and boundary-element methods; adaptive methods; parallel processing machines and applications; and visualization, mesh generation, and AI interfaces.
Good modelling practice in applying computational fluid dynamics for WWTP modelling.
Wicklein, Edward; Batstone, Damien J; Ducoste, Joel; Laurent, Julien; Griborio, Alonso; Wicks, Jim; Saunders, Stephen; Samstag, Randal; Potier, Olivier; Nopens, Ingmar
2016-01-01
Computational fluid dynamics (CFD) modelling in the wastewater treatment (WWT) field is continuing to grow and be used to solve increasingly complex problems. However, the future of CFD models and their value to the wastewater field are a function of their proper application and knowledge of their limits. As has been established for other types of wastewater modelling (i.e. biokinetic models), it is timely to define a good modelling practice (GMP) for wastewater CFD applications. An International Water Association (IWA) working group has been formed to investigate a variety of issues and challenges related to CFD modelling in water and WWT. This paper summarizes the recommendations for GMP of the IWA working group on CFD. The paper provides an overview of GMP and, though it is written for the wastewater application, is based on general CFD procedures. A forthcoming companion paper to provide specific details on modelling of individual wastewater components forms the next step of the working group. PMID:26942517
NASA Astrophysics Data System (ADS)
Avramova, Maria
In the past few decades the need for improved nuclear reactor safety analyses has led to a rapid development of advanced methods for multidimensional thermal-hydraulic analyses. These methods have become progressively more complex in order to account for the many physical phenomena anticipated during steady state and transient Light Water Reactor (LWR) conditions. The advanced thermal-hydraulic subchannel code COBRA-TF (Thurgood, M. J. et al., 1983) is used worldwide for best-estimate evaluations of the nuclear reactor safety margins. In the framework of a joint research project between the Pennsylvania State University (PSU) and AREVA NP GmbH, the theoretical models and numerics of COBRA-TF have been improved. Under the name F-COBRA-TF, the code has been subjected to an extensive verification and validation program and has been applied to variety of LWR steady state and transient simulations. To enable F-COBRA-TF for industrial applications, including safety margins evaluations and design analyses, the code spacer grid models were revised and substantially improved. The state-of-the-art in the modeling of the spacer grid effects on the flow thermal-hydraulic performance in rod bundles employs numerical experiments performed by computational fluid dynamics (CFD) calculations. Because of the involved computational cost, the CFD codes cannot be yet used for full bundle predictions, but their capabilities can be utilized for development of more advanced and sophisticated models for subchannel-level analyses. A subchannel code, equipped with improved physical models, can be then a powerful tool for LWR safety and design evaluations. The unique contributions of this PhD research are seen as development, implementation, and qualification of an innovative spacer grid model by utilizing CFD results within a framework of a subchannel analysis code. Usually, the spacer grid models are mostly related to modeling of the entrainment and deposition phenomena and the heat
Studies of turbulence models in a computational fluid dynamics model of a blood pump.
Song, Xinwei; Wood, Houston G; Day, Steven W; Olsen, Don B
2003-10-01
Computational fluid dynamics (CFD) is used widely in design of rotary blood pumps. The choice of turbulence model is not obvious and plays an important role on the accuracy of CFD predictions. TASCflow (ANSYS Inc., Canonsburg, PA, U.S.A.) has been used to perform CFD simulations of blood flow in a centrifugal left ventricular assist device; a k-epsilon model with near-wall functions was used in the initial numerical calculation. To improve the simulation, local grids with special distribution to ensure the k-omega model were used. Iterations have been performed to optimize the grid distribution and turbulence modeling and to predict flow performance more accurately comparing to experimental data. A comparison of k-omega model and experimental measurements of the flow field obtained by particle image velocimetry shows better agreement than k-epsilon model does, especially in the near-wall regions. PMID:14616539
NASA Astrophysics Data System (ADS)
Tang, H. S.; Qu, K.; Wu, X. G.
2014-09-01
It is now becoming important to develop our capabilities to simulate coastal ocean flows involved with distinct physical phenomena occurring at a vast range of spatial and temporal scales. This paper presents a hybrid modeling system for such simulation. The system consists of a fully three dimensional (3D) fluid dynamics model and a geophysical fluid dynamics model, which couple with each other in two-way and march in time simultaneously. Particularly, in the hybrid system, the solver for incompressible flow on overset meshes (SIFOM) resolves fully 3D small-scale local flow phenomena, while the unstructured grid finite volume coastal ocean model (FVCOM) captures large-scale background flows. The integration of the two models are realized via domain decomposition implemented with an overset grid method. Numerical experiments on performance of the system in resolving flow patterns and solution convergence rate show that the SIFOM-FVCOM system works as intended, and its solutions compare reasonably with data obtained with measurements and other computational approaches. Its unparalleled capabilities to predict multiphysics and multiscale phenomena with high-fidelity are demonstrated by three typical applications that are beyond the reach of other currently existing models. It is anticipated that the SIFOM-FVCOM system will serve as a new platform to study many emerging coastal ocean problems.
NASA Technical Reports Server (NTRS)
Blotzer, Michael J.; Woods, Jody L.
2009-01-01
This viewgraph presentation reviews computational fluid dynamics as a tool for modelling the dispersion of carbon monoxide at the Stennis Space Center's A3 Test Stand. The contents include: 1) Constellation Program; 2) Constellation Launch Vehicles; 3) J2X Engine; 4) A-3 Test Stand; 5) Chemical Steam Generators; 6) Emission Estimates; 7) Located in Existing Test Complex; 8) Computational Fluid Dynamics; 9) Computational Tools; 10) CO Modeling; 11) CO Model results; and 12) Next steps.
Computational Fluid Dynamics (CFD) simulations provide a number of unique opportunities for expanding and improving capabilities for modeling exposures to environmental pollutants. The US Environmental Protection Agency's National Exposure Research Laboratory (NERL) has been c...
Bedrock Channel and Cave Evolution Models Based on Computational Fluid Dynamics
NASA Astrophysics Data System (ADS)
Perne, M.; Covington, M. D.; Cooper, M.
2014-12-01
Models of bedrock channel cross-section evolution typically rely on simple approximations of boundary shear stress to calculate erosion rates across the channel. While such models provide a useful tool for gaining general insight into channel dynamics, they also exhibit a narrower range of behaviors than seen in nature and scale experiments. Recent computational advances enable use of computational fluid dynamics (CFD) to relax many of the assumptions used in these simple models by simulating the full 3D flow field and resulting erosion. We have developed a model of bedrock channel evolution at the reach scale, using CFD, that alternates flow simulation steps with channel evolution steps and evolves the channel in time according to shear stresses calculated from the CFD runs. Caves provide an ideal field setting for studying bedrock channel dynamics, because long records of incision are often preserved in the form of channel widths, meander patterns, and sculpted forms, such as scallops, that indicate flow velocity and direction. However, most existing numerical models of cave formation investigate processes on larger scales, treat conduits as simple shapes, such as cylinders, and deal with the early stages of speleogenesis when sediment transport and erosion mechanisms other than dissolution do not have to be taken into account. Therefore, initial applications of the CFD model focus on the dynamics of cave channels, and particularly on the controls of channel width. While discharge, base level, sediment supply, and the ratio of dissolution to mechanical erosion, are likely to play important roles in determining channel width, we lack a quantitative understanding for the importance of these various factors. Notches in passage walls are thought to result from lateral erosion during periods of increased sediment load when the bed is armored. Modeling is used to check the plausibility of this explanation, and examine whether other mechanisms may also produce notches
Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P.; Russell, MacKenzie R.; Jones, Robert M.; King, Matt; Betterton, Eric A.; Sáez, A. Eduardo
2014-01-01
Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition. PMID:25621085
NASA Astrophysics Data System (ADS)
Stovern, Michael; Felix, Omar; Csavina, Janae; Rine, Kyle P.; Russell, MacKenzie R.; Jones, Robert M.; King, Matt; Betterton, Eric A.; Sáez, A. Eduardo
2014-09-01
Mining operations are potential sources of airborne particulate metal and metalloid contaminants through both direct smelter emissions and wind erosion of mine tailings. The warmer, drier conditions predicted for the Southwestern US by climate models may make contaminated atmospheric dust and aerosols increasingly important, due to potential deleterious effects on human health and ecology. Dust emissions and dispersion of dust and aerosol from the Iron King Mine tailings in Dewey-Humboldt, Arizona, a Superfund site, are currently being investigated through in situ field measurements and computational fluid dynamics modeling. These tailings are heavily contaminated with lead and arsenic. Using a computational fluid dynamics model, we model dust transport from the mine tailings to the surrounding region. The model includes gaseous plume dispersion to simulate the transport of the fine aerosols, while individual particle transport is used to track the trajectories of larger particles and to monitor their deposition locations. In order to improve the accuracy of the dust transport simulations, both regional topographical features and local weather patterns have been incorporated into the model simulations. Results show that local topography and wind velocity profiles are the major factors that control deposition.
Recent advances in the development of implicit schemes for the equations of fluid dynamics
NASA Technical Reports Server (NTRS)
Warming, R. F.; Beam, R. M.
1981-01-01
Innovations and extensions of implicit schemes for equations of fluid dynamics are presented. The notation and theory for linear multistep methods are reviewed, and extensions of work by Beam and Warming (1979) include the implementation of one-leg methods, ADI methods for equations with mixed derivatives, flux vector splitting, the P-dimensional wave equation, and boundary conditions. Numerical experiments indicate that implicit treatment of the boundary conditions is necessary for unconditional stability, and the improvement and implementation of the boundary condition theory should improve the implicit algorithms for gas dynamic equations.
Technical Review of the CENWP Computational Fluid Dynamics Model of the John Day Dam Forebay
Rakowski, Cynthia L.; Serkowski, John A.; Richmond, Marshall C.
2010-12-01
The US Army Corps of Engineers Portland District (CENWP) has developed a computational fluid dynamics (CFD) model of the John Day forebay on the Columbia River to aid in the development and design of alternatives to improve juvenile salmon passage at the John Day Project. At the request of CENWP, Pacific Northwest National Laboratory (PNNL) Hydrology Group has conducted a technical review of CENWP's CFD model run in CFD solver software, STAR-CD. PNNL has extensive experience developing and applying 3D CFD models run in STAR-CD for Columbia River hydroelectric projects. The John Day forebay model developed by CENWP is adequately configured and validated. The model is ready for use simulating forebay hydraulics for structural and operational alternatives. The approach and method are sound, however CENWP has identified some improvements that need to be made for future models and for modifications to this existing model.
Improvement performance of secondary clarifiers by a computational fluid dynamics model
NASA Astrophysics Data System (ADS)
Ghawi, Ali G.; Kriš, J.
2011-12-01
Secondary clarifier is one of the most commonly used unit operations in wastewater treatment plants. It is customarily designed to achieve the separation of solids from biologically treated effluents through the clarification of biological solids and the thickening of sludge. As treatment plants receive increasingly high wastewater flows, conventional sedimentation tanks suffer from overloading problems, which result in poor performance. Modification of inlet baffles through the use of an energy dissipating inlet (EDI) was proposed to enhance the performance in the circular clarifiers at the Al-Dewanyia wastewater treatment plant. A 3-dimensional fully mass conservative clarifier model, based on modern computational fluid dynamics theory, was applied to evaluate the proposed tank modification and to estimate the maximum capacity of the existing and modified clarifiers. A Computational Fluid Dynamics (CFD) model was formulated to describe the tank is performance, and design parameters were obtained based on the experimental results. The study revealed that velocity and (suspended solids) SS is a better parameter than TS (total solids), (Biochemical Oxygen Demand) BOD, (Chemical Oxygen Demand) COD to evaluate the performance of sedimentation tanks and that the removal efficiencies of the suspended solids, biochemical oxygen demand, and chemical oxygen demand were higher in the baffle.
An advanced combustion research facility for validating computational fluid dynamics codes
NASA Astrophysics Data System (ADS)
Bullard, J. B.; Hurley, C. D.; Eccles, N. C.
1991-12-01
The Sector Combustion Rig (SCR), built to obtain experimental data which could be used to verify computational fluid dynamic programs and to investigate the formation and consumption of combustion products through a combustor, is described. This rig was designed to accommodate sectors of full size engine combustion chambers and to test them at real or simulated engine operating conditions. Changes made to improve the operating, measurement, and data handling capabilities of the rig as a result of experience from several years of operations are described together with some of the features which contribute to the uniqueness of the SCR. The SCR gas analysis system and instrumentation are described. Extracts from some results obtained during a recent program of tests on a Rolls-Royce RB211 combustor are given.
Computational Modeling and Analysis of the Fluid Dynamics of Competitive Swimming
NASA Astrophysics Data System (ADS)
Mittal, Rajat
2009-11-01
In order to swim efficiently and/or fast, a swimmer needs to master the subtle cause-and-effect relationship that exists between his/her movements and the surrounding fluid. This is what makes swimming one of the most technical of all sports. For the most part, science has played little if any role in helping swimmers and coaches improve swimming techniques or even to better understand the fluid dynamics of human swimming. Experiments of free swimming humans are extremely difficult to conduct and computational modeling approaches have, in the past, been unable to address this very complex problem. However, the development of a new class of numerical methods, coupled with unique animation and analysis tools is making it possible to analyze swimming strokes in all their complexity. The talk will focus on describing a relatively new numerical method that has been developed to solve flows with highly complex, moving/deforming boundaries. Numerical simulations are used to perform a detailed analysis of the dolphin kick. This stroke has emerged as an important component of competitive swimming in recent years and our analysis has allowed us to extract some useful insights into the fluid dynamics of this stroke. In addition, we also address the continuing debate about the role of lift versus drag in thrust production for human swimming.
Bennett, James S; Feigley, Charles E; Khan, Jamil; Hosni, Mohammad H
2003-01-01
Understanding source behavior is important in controlling exposure to airborne contaminants. Industrial hygienists are often asked to infer emission information from room concentration data. This is not easily done, but models that make simplifying assumptions regarding contaminant transport are frequently used. The errors resulting from these assumptions are not yet well understood. This study compares emission estimates from the single-zone completely mixed (CM-1), two-zone completely mixed (CM-2), and uniform diffusivity (UD) models with the emissions set as boundary conditions in computational fluid dynamic (CFD) simulations of a workplace. The room airflow and concentration fields were computed using Fluent 4. These numerical experiments were factorial combinations of three source locations, five receptor locations, three dilution airflow rates, and two generation rate profiles, constant and time-varying. The aim was to compute plausible concentration fields, not to simulate exactly the processes in a real workroom. Thus, error is defined here as the difference between model and CFD predictions. For the steady-state case the UD model had the lowest error. When the source near-field contained the breathing zone receptor, the CM-2 model was applied. Then, in decreasing agreement with CFD were UD, CM-2, and CM-1. Averaging over all source and receptor locations (CM-2 applied for only one), in decreasing order of agreement with CFD were UD, CM-1, and CM-2. Source and receptor location had large effects on emission estimates using the CM-1 model and some effect using the UD model. A location-specific mixing factor (location factor) derived from steady-state concentration gradients was used to build a more accurate time-dependent emission model, CM-L. Total mass emitted from a time-varying source was modeled most accurately by CM-L, followed by CM-1 and CM-2. PMID:14674808
Computational fluid dynamics modeling of coal gasification in a pressurized spout-fluid bed
Zhongyi Deng; Rui Xiao; Baosheng Jin; He Huang; Laihong Shen; Qilei Song; Qianjun Li
2008-05-15
Computational fluid dynamics (CFD) modeling, which has recently proven to be an effective means of analysis and optimization of energy-conversion processes, has been extended to coal gasification in this paper. A 3D mathematical model has been developed to simulate the coal gasification process in a pressurized spout-fluid bed. This CFD model is composed of gas-solid hydrodynamics, coal pyrolysis, char gasification, and gas phase reaction submodels. The rates of heterogeneous reactions are determined by combining Arrhenius rate and diffusion rate. The homogeneous reactions of gas phase can be treated as secondary reactions. A comparison of the calculated and experimental data shows that most gasification performance parameters can be predicted accurately. This good agreement indicates that CFD modeling can be used for complex fluidized beds coal gasification processes. 37 refs., 7 figs., 5 tabs.
Dr. Chenn Zhou
2008-10-15
Pulverized coal injection (PCI) into the blast furnace (BF) has been recognized as an effective way to decrease the coke and total energy consumption along with minimization of environmental impacts. However, increasing the amount of coal injected into the BF is currently limited by the lack of knowledge of some issues related to the process. It is therefore important to understand the complex physical and chemical phenomena in the PCI process. Due to the difficulty in attaining trus BF measurements, Computational fluid dynamics (CFD) modeling has been identified as a useful technology to provide such knowledge. CFD simulation is powerful for providing detailed information on flow properties and performing parametric studies for process design and optimization. In this project, comprehensive 3-D CFD models have been developed to simulate the PCI process under actual furnace conditions. These models provide raceway size and flow property distributions. The results have provided guidance for optimizing the PCI process.
Epstein, Joshua M.; Pankajakshan, Ramesh; Hammond, Ross A.
2011-01-01
We introduce a novel hybrid of two fields—Computational Fluid Dynamics (CFD) and Agent-Based Modeling (ABM)—as a powerful new technique for urban evacuation planning. CFD is a predominant technique for modeling airborne transport of contaminants, while ABM is a powerful approach for modeling social dynamics in populations of adaptive individuals. The hybrid CFD-ABM method is capable of simulating how large, spatially-distributed populations might respond to a physically realistic contaminant plume. We demonstrate the overall feasibility of CFD-ABM evacuation design, using the case of a hypothetical aerosol release in Los Angeles to explore potential effectiveness of various policy regimes. We conclude by arguing that this new approach can be powerfully applied to arbitrary population centers, offering an unprecedented preparedness and catastrophic event response tool. PMID:21687788
An integrated muscle mechanic-fluid dynamic model of lamprey swimming
NASA Astrophysics Data System (ADS)
Hsu, Chia-Yu; Tytell, Eric; Fauci, Lisa
2009-11-01
In an effort towards a detailed understanding of the generation and control of vertebrate locomotion, including the role of the CPG and its interactions with reflexive feedback, muscle mechanics, and external fluid dynamics, we study a simple vertebrate, the lamprey. Lamprey body undulations are a result of a wave of neural activation that passes from head to tail, causing a wave of muscle activation. These active forces are mediated by passive structural forces. We present recent results from a model that fully couples a viscous, incompressible fluid with nonlinear muscle mechanics. We measure the dependence of the phase lag between activation wave and mechanical wave as a function of model parameters, such as body stiffness and muscle strength. Simulation results are compared to experiments utilizing both real and synthetic lamprey.
Computational Fluid Dynamics Modeling of a Lithium/Thionyl Chloride Battery with Electrolyte Flow
Gu, W.B.; Jungst, Rudolph G.; Nagasubramanian, Ganesan; Wang, C.Y.; Weidner, John.
1999-06-11
A two-dimensional model is developed to simulate discharge of a lithium/thionyl chloride primary battery. The model accounts for not only transport of species and charge, but also the electrode porosity variations and the electrolyte flow induced by the volume reduction caused by electrochemical reactions. Numerical simulations are performed using a finite volume method of computational fluid dynamics. The predicted discharge curves for various temperatures are compared to the experimental data with excellent agreement. Moreover, the simulation results. in conjunction with computer visualization and animation techniques, confirm that cell utilization in the temperature and current range of interest is limited by pore plugging or clogging of the front side of the cathode as a result of LiCl precipitation. The detailed two-dimensional flow simulation also shows that the electrolyte is replenished from the cell header predominantly through the separator into the front of the cathode during most parts of the discharge, especially for higher cell temperatures.
The coupling of fluids, dynamics, and controls on advanced architecture computers
NASA Technical Reports Server (NTRS)
Atwood, Christopher
1995-01-01
This grant provided for the demonstration of coupled controls, body dynamics, and fluids computations in a workstation cluster environment; and an investigation of the impact of peer-peer communication on flow solver performance and robustness. The findings of these investigations were documented in the conference articles.The attached publication, 'Towards Distributed Fluids/Controls Simulations', documents the solution and scaling of the coupled Navier-Stokes, Euler rigid-body dynamics, and state feedback control equations for a two-dimensional canard-wing. The poor scaling shown was due to serialized grid connectivity computation and Ethernet bandwidth limits. The scaling of a peer-to-peer communication flow code on an IBM SP-2 was also shown. The scaling of the code on the switched fabric-linked nodes was good, with a 2.4 percent loss due to communication of intergrid boundary point information. The code performance on 30 worker nodes was 1.7 (mu)s/point/iteration, or a factor of three over a Cray C-90 head. The attached paper, 'Nonlinear Fluid Computations in a Distributed Environment', documents the effect of several computational rate enhancing methods on convergence. For the cases shown, the highest throughput was achieved using boundary updates at each step, with the manager process performing communication tasks only. Constrained domain decomposition of the implicit fluid equations did not degrade the convergence rate or final solution. The scaling of a coupled body/fluid dynamics problem on an Ethernet-linked cluster was also shown.
Fluid dynamics of moving fish in a two-dimensional multiparticle collision dynamics model.
Reid, Daniel A P; Hildenbrandt, H; Padding, J T; Hemelrijk, C K
2012-02-01
The fluid dynamics of animal locomotion, such as that of an undulating fish, are of great interest to both biologists and engineers. However, experimentally studying these fluid dynamics is difficult and time consuming. Model studies can be of great help because of their simpler and more detailed analysis. Their insights may guide empirical work. Particularly the recently introduced multiparticle collision dynamics method may be suitable for the study of moving organisms because it is computationally fast, simple to implement, and has a continuous representation of space. As regards the study of hydrodynamics of moving organisms, the method has only been applied at low Reynolds numbers (below 120) for soft, permeable bodies, and static fishlike shapes. In the present paper we use it to study the hydrodynamics of an undulating fish at Reynolds numbers 1100-1500, after confirming its performance for a moving insect wing at Reynolds number 75. We measure (1) drag, thrust, and lift forces, (2) swimming efficiency and spatial structure of the wake, and (3) distribution of forces along the fish body. We confirm the resemblance between the simulated undulating fish and empirical data. In contrast to theoretical predictions, our model shows that for steadily undulating fish, thrust is produced by the rear 2/3 of the body and that the slip ratio U/V (with U the forward swimming speed and V the rearward speed of the body wave) correlates negatively (instead of positively) with the actual Froude efficiency of swimming. Besides, we show that the common practice of modeling individuals while constraining their sideways acceleration causes them to resemble unconstrained fish with a higher tailbeat frequency. PMID:22463238
Fluid Dynamics of Underwater Flight in Sea Butterflies: Insights from Computational Modeling
NASA Astrophysics Data System (ADS)
Zhou, Zhuoyu; Mittal, Rajat; Yen, Jeannette; Webster, Donald
2014-11-01
Sea butterflies such as Limacine helicina swim by flapping their wing-like parapodia, in a stroke that exhibits a clap-and-fling type kinematics as well as a strong interaction between the parapodia and the body of the animal at the end of downstroke. We used numerical simulations based on videogrammetric data to examine the fluid dynamics and force generation associated with this swimming motion. The unsteady lift-generating mechanism of clap-and-fling results in a sawtooth trajectory with a characteristic ``wobble'' in pitch. We employ coupled flow-body-dynamics simulations to model the free-swimming motion of the organism and explore the efficiency of propulsion as well the factors such as shell weight, that affect its sawtooth swimming trajectory. This work is funded by NSF Grant 1246317 from the Division of Polar Programs.
NASA Astrophysics Data System (ADS)
Ivancic, B.; Riedmann, H.; Frey, M.; Knab, O.; Karl, S.; Hannemann, K.
2016-07-01
The paper summarizes technical results and first highlights of the cooperation between DLR and Airbus Defence and Space (DS) within the work package "CFD Modeling of Combustion Chamber Processes" conducted in the frame of the Propulsion 2020 Project. Within the addressed work package, DLR Göttingen and Airbus DS Ottobrunn have identified several test cases where adequate test data are available and which can be used for proper validation of the computational fluid dynamics (CFD) tools. In this paper, the first test case, the Penn State chamber (RCM1), is discussed. Presenting the simulation results from three different tools, it is shown that the test case can be computed properly with steady-state Reynolds-averaged Navier-Stokes (RANS) approaches. The achieved simulation results reproduce the measured wall heat flux as an important validation parameter very well but also reveal some inconsistencies in the test data which are addressed in this paper.
Richmond, Marshall C.; Rakowski, Cynthia L.; Perkins, William A.; Serkowski, John A.; Ebner, Laurie L.; Schlenker, Stephen J.
2009-07-27
At The Dalles Dam, located between Oregon and Washington on the Columbia River, juvenile salmon passing over the spillway have a survival rate that is below acceptable levels. An important factor affecting survival is the egress route fish take through the immediate tailrace of the dam. Passage through the high-energy spillway and stilling basin environment can leave fish disoriented and vulnerable to predators. Egress conditions can be improved through structural and operational modifications that provide flow paths that move fish more rapidly into the thalweg of the river hence reducing their exposure to predators. We used the results from free-surface computational fluid dynamics (CFD) modeling combined with Lagrangian particle tracking to evaluate the tailrace egress conditions at The Dalles Dam for different alignments of a proposed guidance wall and for different spillway discharge scenarios.
Assessment of Computational Fluid Dynamics (CFD) Models for Shock Boundary-Layer Interaction
NASA Technical Reports Server (NTRS)
DeBonis, James R.; Oberkampf, William L.; Wolf, Richard T.; Orkwis, Paul D.; Turner, Mark G.; Babinsky, Holger
2011-01-01
A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses.
Numeric modeling approximation of the fluid dynamics in an optical fiber trap
NASA Astrophysics Data System (ADS)
Hernández Zavala, J. E.; Cerecedo Nuñez, H. H.; Vigueras Zuñiga, M. O.; Padilla Sosa, P.
2014-09-01
This document presents a first approach to study the behavior of a static fluid radiated by infrared light (980nm,100mW) transmitted by a single-mode optical fiber, for this simulation temperature and radiation pressure are calculated based on the intensity delivered by a laser diode. The Computing Fluid Dynamics (CFD) results were based on a mesh Tet/Hybrid, TGrid for a Silica micro-particle and a mesh Hex/Wedge, Cooper for the beam. The results show that as the particle moves along the axis, temperature and pressure decreases, having the points of mayor temperature and pressure around the axis. The conclusion of this work is that it is possible to simulate the interactions between the beam, the micro-particle and the surrounding medium in terms of temperature, velocity and pressure using the energy and viscous model.
Computational fluid dynamics modelling of left valvular heart diseases during atrial fibrillation
Saglietto, Andrea; Gaita, Fiorenzo; Ridolfi, Luca; Anselmino, Matteo
2016-01-01
Background: Although atrial fibrillation (AF), a common arrhythmia, frequently presents in patients with underlying valvular disease, its hemodynamic contributions are not fully understood. The present work aimed to computationally study how physical conditions imposed by pathologic valvular anatomy act on AF hemodynamics. Methods: We simulated AF with different severity grades of left-sided valvular diseases and compared the cardiovascular effects that they exert during AF, compared to lone AF. The fluid dynamics model used here has been recently validated for lone AF and relies on a lumped parameterization of the four heart chambers, together with the systemic and pulmonary circulation. The AF modelling involves: (i) irregular, uncorrelated and faster heart rate; (ii) atrial contractility dysfunction. Three different grades of severity (mild, moderate, severe) were analyzed for each of the four valvulopathies (AS, aortic stenosis, MS, mitral stenosis, AR, aortic regurgitation, MR, mitral regurgitation), by varying–through the valve opening angle–the valve area. Results: Regurgitation was hemodynamically more relevant than stenosis, as the latter led to inefficient cardiac flow, while the former introduced more drastic fluid dynamics variation. Moreover, mitral valvulopathies were more significant than aortic ones. In case of aortic valve diseases, proper mitral functioning damps out changes at atrial and pulmonary levels. In the case of mitral valvulopathy, the mitral valve lost its regulating capability, thus hemodynamic variations almost equally affected regions upstream and downstream of the valve. In particular, the present study revealed that both mitral and aortic regurgitation strongly affect hemodynamics, followed by mitral stenosis, while aortic stenosis has the least impact among the analyzed valvular diseases. Discussion: The proposed approach can provide new mechanistic insights as to which valvular pathologies merit more aggressive treatment of
Wu, Binxin
2010-12-01
In this paper, 12 turbulence models for single-phase non-newtonian fluid flow in a pipe are evaluated by comparing the frictional pressure drops obtained from computational fluid dynamics (CFD) with those from three friction factor correlations. The turbulence models studied are (1) three high-Reynolds-number k-ε models, (2) six low-Reynolds-number k-ε models, (3) two k-ω models, and (4) the Reynolds stress model. The simulation results indicate that the Chang-Hsieh-Chen version of the low-Reynolds-number k-ε model performs better than the other models in predicting the frictional pressure drops while the standard k-ω model has an acceptable accuracy and a low computing cost. In the model applications, CFD simulation of mixing in a full-scale anaerobic digester with pumped circulation is performed to propose an improvement in the effective mixing standards recommended by the U.S. EPA based on the effect of rheology on the flow fields. Characterization of the velocity gradient is conducted to quantify the growth or breakage of an assumed floc size. Placement of two discharge nozzles in the digester is analyzed to show that spacing two nozzles 180° apart with each one discharging at an angle of 45° off the wall is the most efficient. Moreover, the similarity rules of geometry and mixing energy are checked for scaling up the digester. PMID:21047058
Environmental Fluid Dynamics Code
The Environmental Fluid Dynamics Code (EFDC)is a state-of-the-art hydrodynamic model that can be used to simulate aquatic systems in one, two, and three dimensions. It has evolved over the past two decades to become one of the most widely used and technically defensible hydrodyn...
Description of fluid dynamics and coupled transports in models of a laminar flow diffusion chamber.
Trávníčková, Tereza; Havlica, Jaromír; Ždímal, Vladimír
2013-08-14
The aim of this study is to assess how much the results of nucleation experiments in a laminar flow diffusion chamber (LFDC) are influenced by the complexity of the model of the transport properties. The effects of the type of fluid dynamic model (the steady state compressible Navier-Stokes system for an ideal gas/parabolic profile approximation) and the contributions of the coupled terms describing the Dufour effects and thermodiffusion on the predicted magnitude of the nucleation maxima and its location were investigated. This study was performed on the model of the homogeneous nucleation of an n-butanol-He vapor mixture in a LFDC. The isothermal dependencies of the nucleation rate on supersaturation were determined at three nucleation temperatures: 265 K, 270 K, and 280 K. For this purpose, the experimental LFDC data measured by A. P. Hyvärinen et al. [J. Chem. Phys. 124, 224304 (2006)] were reevaluated using transport models at different levels of complexity. Our results indicate that the type of fluid dynamical model affects both the position of the nucleation maxima in the LFDC and the maximum value of the nucleation rate. On the other hand, the Dufour effects and thermodiffusion perceptibly influence only the value of the maximal nucleation rate. Its position changes only marginally. The dependence of the maximum experimental nucleation rate on the saturation ratio and nucleation temperature was acquired for each case. Based on this dependence, we presented a method for the comparison and evaluation of the uncertainties of simpler models' solutions for the results, where we assumed that the model with Navier-Stokes equations and both coupled effects taken into account was the basis. From this comparison, it follows that an inappropriate choice of mathematical models could lead to relative errors of the order of several hundred percent in the maximum experimental nucleation rate. In the conclusion of this study, we also provide some general recommendations
NASA Astrophysics Data System (ADS)
Dietterich, H. R.; Lev, E.; Chen, J.; Cashman, K. V.; Honor, C.
2015-12-01
Recent eruptions in Hawai'i, Iceland, and Cape Verde highlight the need for improved lava flow models for forecasting and hazard assessment. Existing models used for lava flow simulation range in assumptions, complexity, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess the capabilities of existing models and test the development of new codes, we conduct a benchmarking study of computational fluid dynamics models for lava flows, including VolcFlow, OpenFOAM, Flow3D, and COMSOL. Using new benchmark scenarios defined in Cordonnier et al. (2015) as a guide, we model Newtonian, Herschel-Bulkley and cooling flows over inclined planes, obstacles, and digital elevation models with a wide range of source conditions. Results are compared to analytical theory, analogue and molten basalt experiments, and measurements from natural lava flows. Our study highlights the strengths and weakness of each code, including accuracy and computational costs, and provides insights regarding code selection. We apply the best-fit codes to simulate the lava flows in Harrat Rahat, a predominately mafic volcanic field in Saudi Arabia. Input parameters are assembled from rheology and volume measurements of past flows using geochemistry, crystallinity, and present-day lidar and photogrammetric digital elevation models. With these data, we use our verified models to reconstruct historic and prehistoric events, in order to assess the hazards posed by lava flows for Harrat Rahat.
Experimental Validation of a Computational Fluid Dynamics Model of Copper Electrowinning
NASA Astrophysics Data System (ADS)
Leahy, Martin J.; Philip Schwarz, M.
2010-12-01
The hydrodynamics that occur in the space between the electrode plates in copper electrowinning (EW) are simulated using a computational fluid dynamics model (CFD). The model solves for the phases of gas oxygen bubbles and electrolyte using the Navier-Stokes equations in a CFD framework. An oxygen source is added to the anode, which sets up a recirculation pattern. The gradients in copper near the cathode lead to buoyancy forces, which result in an uplift in the electrolyte close to the cathode. This study investigates the experimental validation of the CFD model using a small/medium-scale real EW system. The predicted fluid velocity profiles are compared with the experimental values, which have been measured along various cross sections of the gap between the anode and the cathode. The results show that the CFD model accurately predicts the velocity profile at several heights in the plate pair. The CFD model prediction of the gas hold-up and the recirculation pattern is compared with visualizations from the experiment. The CFD model prediction is shown to be good across several different operating conditions and geometries, showing that the fundamental underlying equations used in the CFD model transfer to these cases without adjusting the model parameters.
NASA Astrophysics Data System (ADS)
Dean, Cayla; Soloviev, Alexander; Hirons, Amy; Frank, Tamara; Wood, Jon
2016-02-01
Recent studies (Dewar et al., 2006; Wilhelmus and Dabiri, 2014) suggest that diel vertical migrations (DVM) of zooplankton (or other migrating organisms) may have an impact on ocean mixing, though details are not completely clear. Zooplankton that undergo DVM can have an impact on oil transport through the water column, and oil and dispersants can have a negative or even lethal effect on the organisms. Kunze et al. (2006) reported an increase of dissipation rate of turbulent kinetic energy, ε, by four to five orders of magnitude during DVM of zooplankton over background turbulence in Saanich Inlet, British Columbia, Canada. However, the effect was not observed in the same area by Rousseau et al. (2010) and was later reassessed by Kunze (2011). In our work, an 11-month data set obtained in the Straits of Florida with a bottom-mounted acoustic Doppler current profiler revealed strong sound scattering layers undergoing DVM. We used a 3-D non-hydrostatic computational fluid dynamics model with Lagrangian particle injections (a proxy for migrating organisms) via a discrete phase model to simulate the effect of turbulence generation by DVM. We tested a range of organism concentrations from 1000 to 10,000 organisms/m3 based on measurements by Greenlaw (1979) and Mackie and Mills (1983) in Saanich Inlet. At a concentration close to the upper limit, the simulation showed an increase in ε by two to three orders of magnitude during DVM over background turbulence, 10-9 W kg-1. At a concentration of 1000 organisms/m3, almost no turbulence above the background level was produced in the model. These results suggest that the Kunze et al. (2006) observations could have been performed at a larger concentration of migrating zooplankton than those reported by Rousseau et al. (2010). No exact zooplankton concentrations data were provided in either work. The difference between observations and the model can, in part, be explained by the fact that Kunze et al. (2006) measured
Rakowski, Cynthia L.; Richmond, Marshall C.; Serkowski, John A.
2006-12-01
A computational fluid dynamics (CFD) model was used in an investigation into the suppression of a surface vortex that forms and the south-most spilling bay at The Dalles Project. The CFD work complemented work at the prototype and the reduced-scale physical models. The CFD model was based on a model developed for other work in the forebay but had additional resolution added near the spillway. Vortex suppression devices (VSDs) were to placed between pier noses and/or in the bulkhead slot of the spillway bays. The simulations in this study showed that placing VSD structures or a combination of structures to suppress the vortex would still result in near-surface flows to be entrained in a vortex near the downstream spillwall. These results were supported by physical model and prototype studies. However, there was a consensus of the fish biologists at the physical model that the fish would most likely move north and if the fish went under the VSD it would immediately exit the forebay through the tainter gate and not get trapped between VSDs or the VSDs and the tainter gate if the VSDs were deep enough.
Park, Wonhyoung; Song, Yunsun; Park, Kye Jin; Koo, Hae-Won; Yang, Kuhyun
2016-01-01
Purpose Hemodynamic factors are considered to play an important role in initiation and progression of the recurrence after endosaccular coiling of the intracranial aneurysms. We made paired virtual models of completely coiled aneurysms which were subsequently recanalized and compared to identify hemodynamic characteristics related to the recurred aneurysmal sac. Materials and Methods We created paired virtual models of computational fluid dynamics (CFD) in five aneurysms which were initially regarded as having achieved complete occlusion and then recurred during follow-up. Paired virtual models consisted of the CFD model of 3D rotational angiography obtained in the recurred aneurysm and the control model of the initial, parent artery after artificial removal of the coiled and recanalized aneurysm. Using the CFD analysis of the virtual model, we analyzed the hemodynamic characteristics on the neck of each aneurysm before and after its recurrence. Results High wall shear stress (WSS) was identified at the cross-sectionally identified aneurysm neck at which recurrence developed in all cases. A small vortex formation with relatively low velocity in front of the neck was also identified in four cases. The aneurysm recurrence locations corresponded to the location of high WSS and/or small vortex formation. Conclusion Recanalized aneurysms revealed increased WSS and small vortex formation at the cross-sectional neck of the aneurysm. This observation may partially explain the hemodynamic causes of future recanalization after coil embolization. PMID:26958410
NASA Astrophysics Data System (ADS)
Laan, Nick; de Bruin, Karla G.; Slenter, Denise; Wilhelm, Julie; Jermy, Mark; Bonn, Daniel
2015-06-01
Bloodstain Pattern Analysis is a forensic discipline in which, among others, the position of victims can be determined at crime scenes on which blood has been shed. To determine where the blood source was investigators use a straight-line approximation for the trajectory, ignoring effects of gravity and drag and thus overestimating the height of the source. We determined how accurately the location of the origin can be estimated when including gravity and drag into the trajectory reconstruction. We created eight bloodstain patterns at one meter distance from the wall. The origin’s location was determined for each pattern with: the straight-line approximation, our method including gravity, and our method including both gravity and drag. The latter two methods require the volume and impact velocity of each bloodstain, which we are able to determine with a 3D scanner and advanced fluid dynamics, respectively. We conclude that by including gravity and drag in the trajectory calculation, the origin’s location can be determined roughly four times more accurately than with the straight-line approximation. Our study enables investigators to determine if the victim was sitting or standing, or it might be possible to connect wounds on the body to specific patterns, which is important for crime scene reconstruction.
Supersonic Retro-Propulsion Experimental Design for Computational Fluid Dynamics Model Validation
NASA Technical Reports Server (NTRS)
Berry, Scott A.; Laws, Christopher T.; Kleb, W. L.; Rhode, Matthew N.; Spells, Courtney; McCrea, Andrew C.; Truble, Kerry A.; Schauerhamer, Daniel G.; Oberkampf, William L.
2011-01-01
The development of supersonic retro-propulsion, an enabling technology for heavy payload exploration missions to Mars, is the primary focus for the present paper. A new experimental model, intended to provide computational fluid dynamics model validation data, was recently designed for the Langley Research Center Unitary Plan Wind Tunnel Test Section 2. Pre-test computations were instrumental for sizing and refining the model, over the Mach number range of 2.4 to 4.6, such that tunnel blockage and internal flow separation issues would be minimized. A 5-in diameter 70-deg sphere-cone forebody, which accommodates up to four 4:1 area ratio nozzles, followed by a 10-in long cylindrical aftbody was developed for this study based on the computational results. The model was designed to allow for a large number of surface pressure measurements on the forebody and aftbody. Supplemental data included high-speed Schlieren video and internal pressures and temperatures. The run matrix was developed to allow for the quantification of various sources of experimental uncertainty, such as random errors due to run-to-run variations and bias errors due to flow field or model misalignments. Some preliminary results and observations from the test are presented, although detailed analyses of the data and uncertainties are still on going.
Three-Dimensional Computational Fluid Dynamics Modeling of Solid Oxide Electrolysis Cells and Stacks
Grant Hawkes; James O'Brien; Carl Stoots; Stephen Herring
2008-07-01
A three-dimensional computational fluid dynamics (CFD) electrochemical model has been created for detailed analysis of a high-temperature electrolysis stack (solid oxide fuel cells operated as electrolyzers). Inlet and outlet plenum flow distributions are discussed. Maldistribution of plena flow show deviations in per-cell operating conditions due to non-uniformity of species concentrations. Models have also been created to simulate experimental conditions and for code validation. Comparisons between model predictions and experimental results are discussed. Mass, momentum, energy, and species conservation and transport are provided via the core features of the commercial CFD code FLUENT. A solid-oxide fuel cell (SOFC) model adds the electrochemical reactions and loss mechanisms and computation of the electric field throughout the cell. The FLUENT SOFC user-defined subroutine was modified for this work to allow for operation in the electrolysis mode. Model results provide detailed profiles of temperature, Nernst potential, operating potential, activation over-potential, anode-side gas composition, cathode-side gas composition, current density and hydrogen production over a range of stack operating conditions. Variations in flow distribution, and species concentration are discussed. End effects of flow and per-cell voltage are also considered. Predicted mean outlet hydrogen and steam concentrations vary linearly with current density, as expected. Contour plots of local electrolyte temperature, current density, and Nernst potential indicate the effects of heat transfer, reaction cooling/heating, and change in local gas composition.
Jenny, Patrick Torrilhon, Manuel; Heinz, Stefan
2010-02-20
In this paper, a stochastic model is presented to simulate the flow of gases, which are not in thermodynamic equilibrium, like in rarefied or micro situations. For the interaction of a particle with others, statistical moments of the local ensemble have to be evaluated, but unlike in molecular dynamics simulations or DSMC, no collisions between computational particles are considered. In addition, a novel integration technique allows for time steps independent of the stochastic time scale. The stochastic model represents a Fokker-Planck equation in the kinetic description, which can be viewed as an approximation to the Boltzmann equation. This allows for a rigorous investigation of the relation between the new model and classical fluid and kinetic equations. The fluid dynamic equations of Navier-Stokes and Fourier are fully recovered for small relaxation times, while for larger values the new model extents into the kinetic regime. Numerical studies demonstrate that the stochastic model is consistent with Navier-Stokes in that limit, but also that the results become significantly different, if the conditions for equilibrium are invalid. The application to the Knudsen paradox demonstrates the correctness and relevance of this development, and comparisons with existing kinetic equations and standard solution algorithms reveal its advantages. Moreover, results of a test case with geometrically complex boundaries are presented.
Dudley, Peter N; Bonazza, Riccardo; Jones, T Todd; Wyneken, Jeanette; Porter, Warren P
2014-01-01
As global temperatures increase throughout the coming decades, species ranges will shift. New combinations of abiotic conditions will make predicting these range shifts difficult. Biophysical mechanistic niche modeling places bounds on an animal's niche through analyzing the animal's physical interactions with the environment. Biophysical mechanistic niche modeling is flexible enough to accommodate these new combinations of abiotic conditions. However, this approach is difficult to implement for aquatic species because of complex interactions among thrust, metabolic rate and heat transfer. We use contemporary computational fluid dynamic techniques to overcome these difficulties. We model the complex 3D motion of a swimming neonate and juvenile leatherback sea turtle to find power and heat transfer rates during the stroke. We combine the results from these simulations and a numerical model to accurately predict the core temperature of a swimming leatherback. These results are the first steps in developing a highly accurate mechanistic niche model, which can assists paleontologist in understanding biogeographic shifts as well as aid contemporary species managers about potential range shifts over the coming decades. PMID:25354303
Dudley, Peter N.; Bonazza, Riccardo; Jones, T. Todd; Wyneken, Jeanette; Porter, Warren P.
2014-01-01
As global temperatures increase throughout the coming decades, species ranges will shift. New combinations of abiotic conditions will make predicting these range shifts difficult. Biophysical mechanistic niche modeling places bounds on an animal’s niche through analyzing the animal’s physical interactions with the environment. Biophysical mechanistic niche modeling is flexible enough to accommodate these new combinations of abiotic conditions. However, this approach is difficult to implement for aquatic species because of complex interactions among thrust, metabolic rate and heat transfer. We use contemporary computational fluid dynamic techniques to overcome these difficulties. We model the complex 3D motion of a swimming neonate and juvenile leatherback sea turtle to find power and heat transfer rates during the stroke. We combine the results from these simulations and a numerical model to accurately predict the core temperature of a swimming leatherback. These results are the first steps in developing a highly accurate mechanistic niche model, which can assists paleontologist in understanding biogeographic shifts as well as aid contemporary species managers about potential range shifts over the coming decades. PMID:25354303
A new approach for fluid dynamics simulation: The Short-lived Water Cuboid Particle model
NASA Astrophysics Data System (ADS)
Qiao, Changjian; Li, Jiansong; Tian, Zongshun
2016-09-01
There are many researches to simulate the fluid which adopt the traditional particle-based approach and the grid-based approach. However, it needs massive storage in the traditional particle-based approach and it is very complicated to design the grid-based approach with the Navier-Stokes Equations or the Shallow Water Equations (SWEs) because of the difficulty of solving equations. This paper presents a new model called the Short-lived Water Cuboid Particle model. It updates the fluid properties (mass and momentum) recorded in the fixed Cartesian grids by computing the weighted sum of the water cuboid particles with a time step life. Thus it is a two-type-based approach essentially, which not only owns efficient computation and manageable memory like the grid-based approach, but also deals with the discontinuous water surface (wet/dry fronts, boundary conditions, etc.) with high accuracy as well as the particle-based approach. The proposed model has been found capable to simulate the fluid excellently for three laboratory experimental cases and for the field case study of the Malpasset dam-break event occurred in France in 1959. The obtained results show that the model is proved to be an alternative approach to simulate the fluid dynamics with a fair accuracy.
NASA Astrophysics Data System (ADS)
Derakhshani, S. M.; Schott, D. L.; Lodewijks, G.
2013-06-01
Dust emissions can have significant effects on the human health, environment and industry equipment. Understanding the dust generation process helps to select a suitable dust preventing approach and also is useful to evaluate the environmental impact of dust emission. To describe these processes, numerical methods such as Computational Fluid Dynamics (CFD) are widely used, however nowadays particle based methods like Discrete Element Method (DEM) allow researchers to model interaction between particles and fluid flow. In this study, air flow over a stockpile, dust emission, erosion and surface deformation of granular material in the form of stockpile are studied by using DEM and CFD as a coupled method. Two and three dimensional simulations are respectively developed for CFD and DEM methods to minimize CPU time. The standard κ-ɛ turbulence model is used in a fully developed turbulent flow. The continuous gas phase and the discrete particle phase link to each other through gas-particle void fractions and momentum transfer. In addition to stockpile deformation, dust dispersion is studied and finally the accuracy of stockpile deformation results obtained by CFD-DEM modelling will be validated by the agreement with the existing experimental data.
Curtas, Anthony R; Wood, Houston G; Allaire, Paul E; McDaniel, James C; Day, Steven W; Olsen, Don B
2002-01-01
To finalize the design of the next generation of the HeartQuest left ventricular assist device, a suitable impeller had to be designed and tested. The new prototype was based on calculations and test results of previous designs, but required several changes to decrease the size. For most pump designs, this is a simple matter of altering impeller geometry and rotational speed to achieve the desired pressure rise and flow rate. However, this particular pump was limited by housing geometry and the magnetic bearings that support the impeller. Without much freedom in the overall impeller size, the only parameters open to the designers were the blade profiles and the rotating speed. Rather than build several candidates and test them in a rig at enormous cost, computational models of several designs were tested and analyzed. This not only saved money, but also sped up the development time for the project. The computer models were developed in TASCflow, a computational fluid dynamics software package from AEA Technologies. This paper analyzes the data from several of the selected models, paying close attention to pumping performance and general trends from specific design changes. PMID:12296578
The non-equilibrium statistical mechanics of a simple geophysical fluid dynamics model
NASA Astrophysics Data System (ADS)
Verkley, Wim; Severijns, Camiel
2014-05-01
Lorenz [1] has devised a dynamical system that has proved to be very useful as a benchmark system in geophysical fluid dynamics. The system in its simplest form consists of a periodic array of variables that can be associated with an atmospheric field on a latitude circle. The system is driven by a constant forcing, is damped by linear friction and has a simple advection term that causes the model to behave chaotically if the forcing is large enough. Our aim is to predict the statistics of Lorenz' model on the basis of a given average value of its total energy - obtained from a numerical integration - and the assumption of statistical stationarity. Our method is the principle of maximum entropy [2] which in this case reads: the information entropy of the system's probability density function shall be maximal under the constraints of normalization, a given value of the average total energy and statistical stationarity. Statistical stationarity is incorporated approximately by using `stationarity constraints', i.e., by requiring that the average first and possibly higher-order time-derivatives of the energy are zero in the maximization of entropy. The analysis [3] reveals that, if the first stationarity constraint is used, the resulting probability density function rather accurately reproduces the statistics of the individual variables. If the second stationarity constraint is used as well, the correlations between the variables are also reproduced quite adequately. The method can be generalized straightforwardly and holds the promise of a viable non-equilibrium statistical mechanics of the forced-dissipative systems of geophysical fluid dynamics. [1] E.N. Lorenz, 1996: Predictability - A problem partly solved, in Proc. Seminar on Predictability (ECMWF, Reading, Berkshire, UK), Vol. 1, pp. 1-18. [2] E.T. Jaynes, 2003: Probability Theory - The Logic of Science (Cambridge University Press, Cambridge). [3] W.T.M. Verkley and C.A. Severijns, 2014: The maximum entropy
NASA Astrophysics Data System (ADS)
Shen, Yi; Diplas, Panayiotis
2008-01-01
SummaryComplex flow patterns generated by irregular channel topography, such as boulders, submerged large woody debris, riprap and spur dikes, provide unique habitat for many aquatic organisms. Numerical modeling of the flow structures surrounding these obstructions is challenging, yet it represents an important tool for aquatic habitat assessment. In this study, the ability of two- (2-D) and three-dimensional (3-D) computational fluid dynamics models to reproduce these localized complex flow features is examined. The 3-D model is validated with laboratory data obtained from the literature for the case of a flow around a hemisphere under emergent and submerged conditions. The performance of the 2-D and 3-D models is then evaluated by comparing the numerical results with field measurements of flow around several boulders located at a reach of the Smith River, a regulated mountainous stream, obtained at base and peak flows. Close agreement between measured values and the velocity profiles predicted by the two models is obtained outside the wakes behind the hemisphere and boulders. However, the results suggest that in the vicinity of these obstructions the 3-D model is better suited for reproducing the circulation flow behavior at both low and high discharges. Application of the 2-D and 3-D models to meso-scale stream flows of ecological significance is furthermore demonstrated by using a recently developed spatial hydraulic metric to quantify flow complexity surrounding a number of brown trout spawning sites. It is concluded that the 3-D model can provide a much more accurate description of the heterogeneous velocity patterns favored by many aquatic species over a broad range of flows, especially under deep flow conditions when the various obstructions are submerged. Issues pertaining to selection of appropriate models for a variety of flow regimes and potential implication of the 3-D model on the development of better habitat suitability criteria are discussed. The
Rakowski, Cynthia L.; Richmond, Marshall C.; Serkowski, John A.; Johnson, Gary E.
2005-03-10
Computational fluid dynamics (CFD) models were developed to support the siting and design of a behavioral guidance system (BGS) structure in The Dalles Dam (TDA) forebay on the Columbia River. The work was conducted by Pacific Northwest National Laboratory for the U.S. Army Corps of Engineers, Portland District (CENWP). The CFD results were an invaluable tool for the analysis, both from a Regional and Agency perspective (for the fish passage evaluation) and a CENWP perspective (supporting the BGS design and location). The new CFD model (TDA forebay model) included the latest bathymetry (surveyed in 1999) and a detailed representation of the engineered structures (spillway, powerhouse main, fish, and service units). The TDA forebay model was designed and developed in a way that future studies could easily modify or, to a large extent, reuse large portions of the existing mesh. This study resulted in these key findings: (1) The TDA forebay model matched well with field-measured velocity data. (2) The TDA forebay model matched observations made at the 1:80 general physical model of the TDA forebay. (3) During the course of this study, the methodology typically used by CENWP to contour topographic data was shown to be inaccurate when applied to widely-spaced transect data. Contouring methodologies need to be revisited--especially before such things as modifying the bathymetry in the 1:80 general physical model are undertaken. Future alignments can be evaluated with the model staying largely intact. The next round of analysis will need to address fish passage demands and navigation concerns. CFD models can be used to identify the most promising locations and to provide quantified metrics for biological, hydraulic, and navigation criteria. The most promising locations should then be further evaluated in the 1:80 general physical model.
Computational fluid dynamics modeling of bun baking process under different oven load conditions.
Tank, A; Chhanwal, N; Indrani, D; Anandharamakrishnan, C
2014-09-01
A computational fluid dynamics (CFD) model was developed to study the temperature profile of the bun during baking process. Evaporation-condensation mechanism and effect of the latent heat during phase change of water was incorporated in this model to represent actual bun baking process. Simulation results were validated with experimental measurements of bun temperature at two different positions. Baking process is completed within 20 min, after the temperature of crumb become stable at 98 °C. Further, this study was extended to investigate the effect of partially (two baking trays) loaded and fully loaded (eight baking trays) oven on temperature profile of bun. Velocity and temperature profile differs in partially loaded and fully loaded oven. Bun placed in top rack showed rapid baking while bun placed in bottom rack showed slower baking due to uneven temperature distribution in the oven. Hence, placement of bun inside the oven affects temperature of bun and consequently, the quality of the product. PMID:25190860
Wind Turbine Modeling for Computational Fluid Dynamics: December 2010 - December 2012
Tossas, L. A. M.; Leonardi, S.
2013-07-01
With the shortage of fossil fuel and the increasing environmental awareness, wind energy is becoming more and more important. As the market for wind energy grows, wind turbines and wind farms are becoming larger. Current utility-scale turbines extend a significant distance into the atmospheric boundary layer. Therefore, the interaction between the atmospheric boundary layer and the turbines and their wakes needs to be better understood. The turbulent wakes of upstream turbines affect the flow field of the turbines behind them, decreasing power production and increasing mechanical loading. With a better understanding of this type of flow, wind farm developers could plan better-performing, less maintenance-intensive wind farms. Simulating this flow using computational fluid dynamics is one important way to gain a better understanding of wind farm flows. In this study, we compare the performance of actuator disc and actuator line models in producing wind turbine wakes and the wake-turbine interaction between multiple turbines. We also examine parameters that affect the performance of these models, such as grid resolution, the use of a tip-loss correction, and the way in which the turbine force is projected onto the flow field.
NASA Astrophysics Data System (ADS)
Huda, Nazmul; Naser, Jamal; Brooks, Geoffrey; Reuter, Markus A.; Matusewicz, Robert W.
2012-02-01
Slag fuming is a reductive treatment process for molten zinciferous slags for extracting zinc in the form of metal vapor by injecting or adding a reductant source such as pulverized coal or lump coal and natural gas. A computational fluid dynamic (CFD) model was developed to study the zinc slag fuming process from imperial smelting furnace (ISF) slag in a top-submerged lance furnace and to investigate the details of fluid flow, reaction kinetics, and heat transfer in the furnace. The model integrates combustion phenomena and chemical reactions with the heat, mass, and momentum interfacial interaction between the phases present in the system. A commercial CFD package AVL Fire 2009.2 (AVL, Graz, Austria) coupled with a number of user-defined subroutines in FORTRAN programming language were used to develop the model. The model is based on three-dimensional (3-D) Eulerian multiphase flow approach, and it predicts the velocity and temperature field of the molten slag bath, generated turbulence, and vortex and plume shape at the lance tip. The model also predicts the mass fractions of slag and gaseous components inside the furnace. The model predicted that the percent of ZnO in the slag bath decreases linearly with time and is consistent broadly with the experimental data. The zinc fuming rate from the slag bath predicted by the model was validated through macrostep validation process against the experimental study of Waladan et al. The model results predicted that the rate of ZnO reduction is controlled by the mass transfer of ZnO from the bulk slag to slag-gas interface and rate of gas-carbon reaction for the specified simulation time studied. Although the model is based on zinc slag fuming, the basic approach could be expanded or applied for the CFD analysis of analogous systems.
NASA Astrophysics Data System (ADS)
Raffaeli, Leah R.
Techno-economic and systems studies on microalgal growth scenarios to date are abbreviated and missing a number of important variables. By including these variables in a detailed model integrating biology, chemistry, engineering, and financial aspects, a more defined systems analysis is possible. Through optimizing the model productivity based on the resulting net profit, the system analysis results in a more accurate assessment of environmental and economic sustainability of specific algal growth scenarios. Photobioreactor algal growth scenario optimization in the system model has resulted in realistic engineering design requirements based on algal growth requirements and fluid dynamics analysis. Results show feasibility for photobioreactor growth scenarios to be economically sustainable when co-products are included, but definite technological advancements and productivity improvements must be made. The main factors inhibiting a cost effective photobioreactor growth scenario are culture density, temperature, and lighting distribution for solar illuminated photobioreactors, and lighting cost for artificially illuminated photobioreactors. Open pond algal growth scenarios do not show any prospect of economic or environmental sustainability with current technology due to the large amount of surface area required, inefficient water use, and low culture density. All algal growth scenarios are inferior to petro-diesel regarding energy inputs, carbon emissions, and environmental sustainability. No algal growth scenarios analyzed in this study meet the U.S. requirement of biofuel emitting at least 20% less carbon emissions than diesel from crude oil.
Atmospheric and Oceanic Fluid Dynamics
NASA Astrophysics Data System (ADS)
Vallis, Geoffrey K.
2006-11-01
Fluid dynamics is fundamental to our understanding of the atmosphere and oceans. Although many of the same principles of fluid dynamics apply to both the atmosphere and oceans, textbooks tend to concentrate on the atmosphere, the ocean, or the theory of geophysical fluid dynamics (GFD). This textbook provides a comprehensive unified treatment of atmospheric and oceanic fluid dynamics. The book introduces the fundamentals of geophysical fluid dynamics, including rotation and stratification, vorticity and potential vorticity, and scaling and approximations. It discusses baroclinic and barotropic instabilities, wave-mean flow interactions and turbulence, and the general circulation of the atmosphere and ocean. Student problems and exercises are included at the end of each chapter. Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation will be an invaluable graduate textbook on advanced courses in GFD, meteorology, atmospheric science and oceanography, and an excellent review volume for researchers. Additional resources are available at www.cambridge.org/9780521849692. Includes end of chapter review questions to aid understanding Unified and comprehensive treatment of both atmospheric and oceanic fluid dynamics Covers many modern topics and provides up to date knowledge
Computational fluid dynamic applications
Chang, S.-L.; Lottes, S. A.; Zhou, C. Q.
2000-04-03
The rapid advancement of computational capability including speed and memory size has prompted the wide use of computational fluid dynamics (CFD) codes to simulate complex flow systems. CFD simulations are used to study the operating problems encountered in system, to evaluate the impacts of operation/design parameters on the performance of a system, and to investigate novel design concepts. CFD codes are generally developed based on the conservation laws of mass, momentum, and energy that govern the characteristics of a flow. The governing equations are simplified and discretized for a selected computational grid system. Numerical methods are selected to simplify and calculate approximate flow properties. For turbulent, reacting, and multiphase flow systems the complex processes relating to these aspects of the flow, i.e., turbulent diffusion, combustion kinetics, interfacial drag and heat and mass transfer, etc., are described in mathematical models, based on a combination of fundamental physics and empirical data, that are incorporated into the code. CFD simulation has been applied to a large variety of practical and industrial scale flow systems.
Jacob, Rick E.; Lamm, W. J.
2011-11-08
Pulmonary computational fluid dynamics models require 3D images to be acquired over multiple points in the dynamic breathing cycle, with no breath holds or changes in ventilatory mechanics. With small animals, these requirements result in long imaging times ({approx}90 minutes), over which lung mechanics, such as compliance, can gradually change if not carefully monitored and controlled. These changes, caused by derecruitment of parenchymal tissue, are manifested as an upward drift in peak inspiratory pressure or by changes in the pressure waveform and/or lung volume over the course of the experiment. We demonstrate highly repeatable mechanical ventilation in anesthetized rats over a long duration for pulmonary CT imaging throughout the dynamic breathing cycle. We describe significant updates to a basic commercial ventilator that was acquired for these experiments. Key to achieving consistent results was the implementation of periodic deep breaths, or sighs, of extended duration to maintain lung recruitment. In addition, continuous monitoring of breath-to-breath pressure and volume waveforms and long-term trends in peak inspiratory pressure and flow provide diagnostics of changes in breathing mechanics.
Jacob, Richard E.; Lamm, Wayne J.
2011-01-01
Pulmonary computational fluid dynamics models require that three-dimensional images be acquired over multiple points in the dynamic breathing cycle without breath holds or changes in ventilatory mechanics. With small animals, these requirements can result in long imaging times (∼90 minutes), over which lung mechanics, such as compliance, may gradually change if not carefully monitored and controlled. These changes, caused by derecruitment of parenchymal tissue, are manifested as an upward drift in peak inspiratory pressure (PIP) or by changes in the pressure waveform and/or lung volume over the course of the experiment. We demonstrate highly repeatable mechanical ventilation in anesthetized rats over a long duration for dynamic lung x-ray computed tomography (CT) imaging. We describe significant updates to a basic commercial ventilator that was acquired for these experiments. Key to achieving consistent results was the implementation of periodic deep breaths, or sighs, of extended duration to maintain lung recruitment. In addition, continuous monitoring of breath-to-breath pressure and volume waveforms and long-term trends in PIP and flow provide diagnostics of changes in breathing mechanics. PMID:22087338
Computational fluid dynamic control
NASA Technical Reports Server (NTRS)
Hartley, Tom T.; Deabreu-Garcia, Alex
1989-01-01
A general technique is presented for modeling fluid, or gas, dynamic systems specifically for the development of control systems. The numerical methods which are generally used in computational fluid dynamics are borrowed to create either continuous-time or discrete-time models of the particular fluid system. The resulting equations can be either left in a nonlinear form, or easily linearized about an operating point. As there are typically very many states in these systems, the usual linear model reduction methods can be used on them to allow a low-order controller to be designed. A simple example is given which typifies many internal flow control problems. The resulting control is termed computational fluid dynamic control.
Computational fluid dynamics modeling of transport and deposition of pesticides in an aircraft cabin
Isukapalli, Sastry S.; Mazumdar, Sagnik; George, Pradeep; Wei, Binnian; Jones, Byron; Weisel, Clifford P.
2015-01-01
Spraying of pesticides in aircraft cabins is required by some countries as part of a disinsection process to kill insects that pose a public health threat. However, public health concerns remain regarding exposures of cabin crew and passengers to pesticides in aircraft cabins. While large scale field measurements of pesticide residues and air concentrations in aircraft cabins scenarios are expensive and time consuming, Computational Fluid Dynamics (CFD) models provide an effective alternative for characterizing concentration distributions and exposures. This study involved CFD modeling of a twin-aisle 11 row cabin mockup with heated manikins, mimicking a part of a fully occupied Boeing 767 cabin. The model was applied to study the flow and deposition of pesticides under representative scenarios with different spraying patterns (sideways and overhead) and cabin air exchange rates (low and high). Corresponding spraying experiments were conducted in the cabin mockup, and pesticide deposition samples were collected at the manikin’s lap and seat top for a limited set of five seats. The CFD model performed well for scenarios corresponding to high air exchange rates, captured the concentration profiles for middle seats under low air exchange rates, and underestimated the concentrations at window seats under low air exchange rates. Additionally, both the CFD and experimental measurements showed no major variation in deposition characteristics between sideways and overhead spraying. The CFD model can estimate concentration fields and deposition profiles at very high resolutions, which can be used for characterizing the overall variability in air concentrations and surface loadings. Additionally, these model results can also provide a realistic range of surface and air concentrations of pesticides in the cabin that can be used to estimate potential exposures of cabin crew and passengers to these pesticides. PMID:25642134
Computational fluid dynamics modeling of transport and deposition of pesticides in an aircraft cabin
NASA Astrophysics Data System (ADS)
Isukapalli, Sastry S.; Mazumdar, Sagnik; George, Pradeep; Wei, Binnian; Jones, Byron; Weisel, Clifford P.
2013-04-01
Spraying of pesticides in aircraft cabins is required by some countries as part of a disinsection process to kill insects that pose a public health threat. However, public health concerns remain regarding exposures of cabin crew and passengers to pesticides in aircraft cabins. While large scale field measurements of pesticide residues and air concentrations in aircraft cabins scenarios are expensive and time consuming, Computational Fluid Dynamics (CFD) models provide an effective alternative for characterizing concentration distributions and exposures. This study involved CFD modeling of a twin-aisle 11 row cabin mockup with heated manikins, mimicking a part of a fully occupied Boeing 767 cabin. The model was applied to study the flow and deposition of pesticides under representative scenarios with different spraying patterns (sideways and overhead) and cabin air exchange rates (low and high). Corresponding spraying experiments were conducted in the cabin mockup, and pesticide deposition samples were collected at the manikin's lap and seat top for a limited set of five seats. The CFD model performed well for scenarios corresponding to high air exchange rates, captured the concentration profiles for middle seats under low air exchange rates, and underestimated the concentrations at window seats under low air exchange rates. Additionally, both the CFD and experimental measurements showed no major variation in deposition characteristics between sideways and overhead spraying. The CFD model can estimate concentration fields and deposition profiles at very high resolutions, which can be used for characterizing the overall variability in air concentrations and surface loadings. Additionally, these model results can also provide a realistic range of surface and air concentrations of pesticides in the cabin that can be used to estimate potential exposures of cabin crew and passengers to these pesticides.
Computational fluid dynamics modeling of two-phase flow in a BWR fuel assembly. Final CRADA Report.
Tentner, A.; Nuclear Engineering Division
2009-10-13
A direct numerical simulation capability for two-phase flows with heat transfer in complex geometries can considerably reduce the hardware development cycle, facilitate the optimization and reduce the costs of testing of various industrial facilities, such as nuclear power plants, steam generators, steam condensers, liquid cooling systems, heat exchangers, distillers, and boilers. Specifically, the phenomena occurring in a two-phase coolant flow in a BWR (Boiling Water Reactor) fuel assembly include coolant phase changes and multiple flow regimes which directly influence the coolant interaction with fuel assembly and, ultimately, the reactor performance. Traditionally, the best analysis tools for this purpose of two-phase flow phenomena inside the BWR fuel assembly have been the sub-channel codes. However, the resolution of these codes is too coarse for analyzing the detailed intra-assembly flow patterns, such as flow around a spacer element. Advanced CFD (Computational Fluid Dynamics) codes provide a potential for detailed 3D simulations of coolant flow inside a fuel assembly, including flow around a spacer element using more fundamental physical models of flow regimes and phase interactions than sub-channel codes. Such models can extend the code applicability to a wider range of situations, which is highly important for increasing the efficiency and to prevent accidents.
Santoro, Domenico; Raisee, Mehrdad; Moghaddami, Mostafa; Ducoste, Joel; Sasges, Micheal; Liberti, Lorenzo; Notarnicola, Michele
2010-08-15
Advanced Oxidation Processes (AOPs) promoted by ultraviolet light are innovative and potentially cost-effective solutions for treating persistent pollutants recalcitrant to conventional water and wastewater treatment. While several studies have been performed during the past decade to improve the fundamental understanding of the UV-H(2)O(2) AOP and its kinetic modeling, Computational Fluid Dynamics (CFD) has only recently emerged as a powerful tool that allows a deeper understanding of complex photochemical processes in environmental and reactor engineering applications. In this paper, a comprehensive kinetic model of UV-H(2)O(2) AOP was coupled with the Reynolds averaged Navier-Stokes (RANS) equations using CFD to predict the oxidation of tributyl phosphate (TBP) and tri(2-chloroethtyl) phosphate (TCEP) in two different photoreactors: a parallel- and a cross-flow UV device employing a UV lamp emitting primarily 253.7 nm radiation. CFD simulations, obtained for both turbulent and laminar flow regimes and compared with experimental data over a wide range of UV doses, enabled the spatial visualization of hydrogen peroxide and hydroxyl radical distributions in the photoreactor. The annular photoreactor displayed consistently better oxidation performance than the cross-flow system due to the absence of recirculation zones, as confirmed by the hydroxyl radical dose distributions. Notably, such discrepancy was found to be strongly dependent on and directly correlated with the hydroxyl radical rate constant becoming relevant for conditions approaching diffusion-controlled reaction regimes (k(C,OH) > 10(9) M(-1) s(-1)). PMID:20704221
Hariharan, Prasanna; D'Souza, Gavin; Horner, Marc; Malinauskas, Richard A; Myers, Matthew R
2015-09-01
As part of an ongoing effort to develop verification and validation (V&V) standards for using computational fluid dynamics (CFD) in the evaluation of medical devices, we have developed idealized flow-based verification benchmarks to assess the implementation of commonly cited power-law based hemolysis models in CFD. Verification process ensures that all governing equations are solved correctly and the model is free of user and numerical errors. To perform verification for power-law based hemolysis modeling, analytical solutions for the Eulerian power-law blood damage model (which estimates hemolysis index (HI) as a function of shear stress and exposure time) were obtained for Couette and inclined Couette flow models, and for Newtonian and non-Newtonian pipe flow models. Subsequently, CFD simulations of fluid flow and HI were performed using Eulerian and three different Lagrangian-based hemolysis models and compared with the analytical solutions. For all the geometries, the blood damage results from the Eulerian-based CFD simulations matched the Eulerian analytical solutions within ∼1%, which indicates successful implementation of the Eulerian hemolysis model. Agreement between the Lagrangian and Eulerian models depended upon the choice of the hemolysis power-law constants. For the commonly used values of power-law constants (α = 1.9-2.42 and β = 0.65-0.80), in the absence of flow acceleration, most of the Lagrangian models matched the Eulerian results within 5%. In the presence of flow acceleration (inclined Couette flow), moderate differences (∼10%) were observed between the Lagrangian and Eulerian models. This difference increased to greater than 100% as the beta exponent decreased. These simplified flow problems can be used as standard benchmarks for verifying the implementation of blood damage predictive models in commercial and open-source CFD codes. The current study only used power-law model as an illustrative example to emphasize the need
O'Brien, Haley D; Bourke, Jason
2015-12-01
In the mammalian order Artiodactyla, the majority of arterial blood entering the intracranial cavity is supplied by a large arterial meshwork called the carotid rete. This vascular structure functionally replaces the internal carotid artery. Extensive experimentation has demonstrated that the artiodactyl carotid rete drives one of the most effective selective brain cooling mechanisms among terrestrial vertebrates. Less well understood is the impact that the unique morphology of the carotid rete may have on the hemodynamics of blood flow to the cerebrum. It has been hypothesized that, relative to the tubular internal carotid arteries of most other vertebrates, the highly convoluted morphology of the carotid rete may increase resistance to flow during extreme changes in cerebral blood pressure, essentially protecting the brain by acting as a resistor. We test this hypothesis by employing simple and complex physical models to a 3D surface rendering of the carotid rete of the domestic goat, Capra hircus. First, we modeled the potential for increased resistance across the carotid rete using an electrical circuit analog. The extensive branching of the rete equates to a parallel circuit that is bound in series by single tubular arteries, both upstream and downstream. This method calculated a near-zero increase in resistance across the rete. Because basic equations do not incorporate drag, shear-stress, and turbulence, we used computational fluid dynamics to simulate the impact of these computationally intensive factors on resistance. Ultimately, both simple and complex models demonstrated negligible changes in resistance and blood pressure across the arterial meshwork. We further tested the resistive potential of the carotid rete by simulating blood pressures known to occur in giraffes. Based on these models, we found resistance (and blood pressure mitigation as a whole) to be an unlikely function for the artiodactyl carotid rete. PMID:26403501
Walters, D Keith; Luke, William H
2011-01-01
Computational fluid dynamics (CFD) has emerged as a useful tool for the prediction of airflow and particle transport within the human lung airway. Several published studies have demonstrated the use of Eulerian finite-volume CFD simulations coupled with Lagrangian particle tracking methods to determine local and regional particle deposition rates in small subsections of the bronchopulmonary tree. However, the simulation of particle transport and deposition in large-scale models encompassing more than a few generations is less common, due in part to the sheer size and complexity of the human lung airway. Highly resolved, fully coupled flowfield solution and particle tracking in the entire lung, for example, is currently an intractable problem and will remain so for the foreseeable future. This paper adopts a previously reported methodology for simulating large-scale regions of the lung airway (Walters, D. K., and Luke, W. H., 2010, "A Method for Three-Dimensional Navier-Stokes Simulations of Large-Scale Regions of the Human Lung Airway," ASME J. Fluids Eng., 132(5), p. 051101), which was shown to produce results similar to fully resolved geometries using approximate, reduced geometry models. The methodology is extended here to particle transport and deposition simulations. Lagrangian particle tracking simulations are performed in combination with Eulerian simulations of the airflow in an idealized representation of the human lung airway tree. Results using the reduced models are compared with those using the fully resolved models for an eight-generation region of the conducting zone. The agreement between fully resolved and reduced geometry simulations indicates that the new method can provide an accurate alternative for large-scale CFD simulations while potentially reducing the computational cost of these simulations by several orders of magnitude. PMID:21186893
Computational fluid dynamics modeling of laboratory flames and an industrial flare.
Singh, Kanwar Devesh; Gangadharan, Preeti; Chen, Daniel H; Lou, Helen H; Li, Xianchang; Richmond, Peyton
2014-11-01
A computational fluid dynamics (CFD) methodology for simulating the combustion process has been validated with experimental results. Three different types of experimental setups were used to validate the CFD model. These setups include an industrial-scale flare setups and two lab-scale flames. The CFD study also involved three different fuels: C3H6/CH/Air/N2, C2H4/O2/Ar and CH4/Air. In the first setup, flare efficiency data from the Texas Commission on Environmental Quality (TCEQ) 2010 field tests were used to validate the CFD model. In the second setup, a McKenna burner with flat flames was simulated. Temperature and mass fractions of important species were compared with the experimental data. Finally, results of an experimental study done at Sandia National Laboratories to generate a lifted jet flame were used for the purpose of validation. The reduced 50 species mechanism, LU 1.1, the realizable k-epsilon turbulence model, and the EDC turbulence-chemistry interaction model were usedfor this work. Flare efficiency, axial profiles of temperature, and mass fractions of various intermediate species obtained in the simulation were compared with experimental data and a good agreement between the profiles was clearly observed. In particular the simulation match with the TCEQ 2010 flare tests has been significantly improved (within 5% of the data) compared to the results reported by Singh et al. in 2012. Validation of the speciated flat flame data supports the view that flares can be a primary source offormaldehyde emission. PMID:25509554
Computational fluid dynamics research
NASA Technical Reports Server (NTRS)
Chandra, Suresh; Jones, Kenneth; Hassan, Hassan; Mcrae, David Scott
1992-01-01
The focus of research in the computational fluid dynamics (CFD) area is two fold: (1) to develop new approaches for turbulence modeling so that high speed compressible flows can be studied for applications to entry and re-entry flows; and (2) to perform research to improve CFD algorithm accuracy and efficiency for high speed flows. Research activities, faculty and student participation, publications, and financial information are outlined.
Perspectives in Fluid Dynamics
NASA Astrophysics Data System (ADS)
Batchelor, G. K.; Moffatt, H. K.; Worster, M. G.
2002-12-01
With applications ranging from modelling the environment to automotive design and physiology to astrophysics, conventional textbooks cannot hope to give students much information on what topics in fluid dynamics are currently being researched, or how to choose between them. This book rectifies matters. It consists of eleven chapters that introduce and review different branches of the subject for graduate-level courses, or for specialists seeking introductions to other areas. Hb ISBN (2001): 0-521-78061-6
NASA Astrophysics Data System (ADS)
Chitta, Varun
Modeling of complex flows involving the combined effects of flow transition and streamline curvature using two advanced turbulence models, one in the Reynolds-averaged Navier-Stokes (RANS) category and the other in the hybrid RANS-Large eddy simulation (LES) category is considered in this research effort. In the first part of the research, a new scalar eddy-viscosity model (EVM) is proposed, designed to exhibit physically correct responses to flow transition, streamline curvature, and system rotation effects. The four equation model developed herein is a curvature-sensitized version of a commercially available three-equation transition-sensitive model. The physical effects of rotation and curvature (RC) enter the model through the added transport equation, analogous to a transverse turbulent velocity scale. The eddy-viscosity has been redefined such that the proposed model is constrained to reduce to the original transition-sensitive model definition in nonrotating flows or in regions with negligible RC effects. In the second part of the research, the developed four-equation model is combined with a LES technique using a new hybrid modeling framework, dynamic hybrid RANS-LES. The new framework is highly generalized, allowing coupling of any desired LES model with any given RANS model and addresses several deficiencies inherent in most current hybrid models. In the present research effort, the DHRL model comprises of the proposed four-equation model for RANS component and the MILES scheme for LES component. Both the models were implemented into a commercial computational fluid dynamics (CFD) solver and tested on a number of engineering and generic flow problems. Results from both the RANS and hybrid models show successful resolution of the combined effects of transition and curvature with reasonable engineering accuracy, and for only a small increase in computational cost. In addition, results from the hybrid model indicate significant levels of turbulent
Zuo, Wangda; Chen, Qingyan
2011-06-01
To design a healthy indoor environment, it is important to study airborne particle distribution indoors. As an intermediate model between multizone models and computational fluid dynamics (CFD), a fast fluid dynamics (FFD) model can be used to provide temporal and spatial information of particle dispersion in real time. This study evaluated the accuracy of the FFD for predicting transportation of particles with low Stokes number in a duct and in a room with mixed convection. The evaluation was to compare the numerical results calculated by the FFD with the corresponding experimental data and the results obtained by the CFD. The comparison showed that the FFD could capture major pattern of particle dispersion, which is missed in models with well-mixed assumptions. Although the FFD was less accurate than the CFD partially due to its simplification in numeric schemes, it was 53 times faster than the CFD.
Computational Fluid Dynamics Study on the Effects of RATO Timing on the Scale Model Acoustic Test
NASA Technical Reports Server (NTRS)
Nielsen, Tanner; Williams, B.; West, Jeff
2015-01-01
The Scale Model Acoustic Test (SMAT) is a 5% scale test of the Space Launch System (SLS), which is currently being designed at Marshall Space Flight Center (MSFC). The purpose of this test is to characterize and understand a variety of acoustic phenomena that occur during the early portions of lift off, one being the overpressure environment that develops shortly after booster ignition. The SLS lift off configuration consists of four RS-25 liquid thrusters on the core stage, with two solid boosters connected to each side. Past experience with scale model testing at MSFC (in ER42), has shown that there is a delay in the ignition of the Rocket Assisted Take Off (RATO) motor, which is used as the 5% scale analog of the solid boosters, after the signal to ignite is given. This delay can range from 0 to 16.5ms. While this small of a delay maybe insignificant in the case of the full scale SLS, it can significantly alter the data obtained during the SMAT due to the much smaller geometry. The speed of sound of the air and combustion gas constituents is not scaled, and therefore the SMAT pressure waves propagate at approximately the same speed as occurs during full scale. However, the SMAT geometry is much smaller allowing the pressure waves to move down the exhaust duct, through the trench, and impact the vehicle model much faster than occurs at full scale. To better understand the effect of the RATO timing simultaneity on the SMAT IOP test data, a computational fluid dynamics (CFD) analysis was performed using the Loci/CHEM CFD software program. Five different timing offsets, based on RATO ignition delay statistics, were simulated. A variety of results and comparisons will be given, assessing the overall effect of RATO timing simultaneity on the SMAT overpressure environment.
NASA Technical Reports Server (NTRS)
Nielsen, Tanner; West, Jeff
2015-01-01
The Scale Model Acoustic Test (SMAT) is a 5% scale test of the Space Launch System (SLS), which is currently being designed at Marshall Space Flight Center (MSFC). The purpose of this test is to characterize and understand a variety of acoustic phenomena that occur during the early portions of lift off, one being the overpressure environment that develops shortly after booster ignition. The pressure waves that propagate from the mobile launcher (ML) exhaust hole are defined as the ignition overpressure (IOP), while the portion of the pressure waves that exit the duct or trench are the duct overpressure (DOP). Distinguishing the IOP and DOP in scale model test data has been difficult in past experiences and in early SMAT results, due to the effects of scaling the geometry. The speed of sound of the air and combustion gas constituents is not scaled, and therefore the SMAT pressure waves propagate at approximately the same speed as occurs in full scale. However, the SMAT geometry is twenty times smaller, allowing the pressure waves to move down the exhaust hole, through the trench and duct, and impact the vehicle model much faster than occurs at full scale. The DOP waves impact portions of the vehicle at the same time as the IOP waves, making it difficult to distinguish the different waves and fully understand the data. To better understand the SMAT data, a computational fluid dynamics (CFD) analysis was performed with a fictitious geometry that isolates the IOP and DOP. The upper and lower portions of the domain were segregated to accomplish the isolation in such a way that the flow physics were not significantly altered. The Loci/CHEM CFD software program was used to perform this analysis.
State-of-the-art review of computational fluid dynamics modeling for fluid-solids systems
NASA Astrophysics Data System (ADS)
Lyczkowski, R. W.; Bouillard, J. X.; Ding, J.; Chang, S. L.; Burge, S. W.
1994-05-01
As the result of 15 years of research (50 staff years of effort) Argonne National Laboratory (ANL), through its involvement in fluidized-bed combustion, magnetohydrodynamics, and a variety of environmental programs, has produced extensive computational fluid dynamics (CFD) software and models to predict the multiphase hydrodynamic and reactive behavior of fluid-solids motions and interactions in complex fluidized-bed reactors (FBR's) and slurry systems. This has resulted in the FLUFIX, IRF, and SLUFIX computer programs. These programs are based on fluid-solids hydrodynamic models and can predict information important to the designer of atmospheric or pressurized bubbling and circulating FBR, fluid catalytic cracking (FCC) and slurry units to guarantee optimum efficiency with minimum release of pollutants into the environment. This latter issue will become of paramount importance with the enactment of the Clean Air Act Amendment (CAAA) of 1995. Solids motion is also the key to understanding erosion processes. Erosion rates in FBR's and pneumatic and slurry components are computed by ANL's EROSION code to predict the potential metal wastage of FBR walls, intervals, feed distributors, and cyclones. Only the FLUFIX and IRF codes will be reviewed in the paper together with highlights of the validations because of length limitations. It is envisioned that one day, these codes with user-friendly pre- and post-processor software and tailored for massively parallel multiprocessor shared memory computational platforms will be used by industry and researchers to assist in reducing and/or eliminating the environmental and economic barriers which limit full consideration of coal, shale, and biomass as energy sources; to retain energy security; and to remediate waste and ecological problems.
A microfluidic model to study fluid dynamics of mucus plug rupture in small lung airways
Hu, Yingying; Bian, Shiyao; Grotberg, John; Filoche, Marcel; White, Joshua; Takayama, Shuichi; Grotberg, James B.
2015-01-01
Fluid dynamics of mucus plug rupture is important to understand mucus clearance in lung airways and potential effects of mucus plug rupture on epithelial cells at lung airway walls. We established a microfluidic model to study mucus plug rupture in a collapsed airway of the 12th generation. Mucus plugs were simulated using Carbopol 940 (C940) gels at concentrations of 0.15%, 0.2%, 0.25%, and 0.3%, which have non-Newtonian properties close to healthy and diseased lung mucus. The airway was modeled with a polydimethylsiloxane microfluidic channel. Plug motion was driven by pressurized air. Global strain rates and shear stress were defined to quantitatively describe plug deformation and rupture. Results show that a plug needs to overcome yield stress before deformation and rupture. The plug takes relatively long time to yield at the high Bingham number. Plug length shortening is the more significant deformation than shearing at gel concentration higher than 0.15%. Although strain rates increase dramatically at rupture, the transient shear stress drops due to the shear-thinning effect of the C940 gels. Dimensionless time-averaged shear stress, Txy, linearly increases from 3.7 to 5.6 times the Bingham number as the Bingham number varies from 0.018 to 0.1. The dimensionless time-averaged shear rate simply equals to Txy/2. In dimension, shear stress magnitude is about one order lower than the pressure drop, and one order higher than yield stress. Mucus with high yield stress leads to high shear stress, and therefore would be more likely to cause epithelial cell damage. Crackling sounds produced with plug rupture might be more detectable for gels with higher concentration. PMID:26392827
Advances in Fluid Dynamics of Subsurface Flow of Groundwater, Hydrocarbons, and CO2
NASA Astrophysics Data System (ADS)
Weyer, K. U.
2015-12-01
In the past, the chemical methods of contaminant hydrogeology have dominated much of hydrogeological thinking. In their wake, understanding the physics of subsurface fluid flow and its application to practice and science seemingly has played a secondary role and it often has been replaced by numerical modelling only. Building an understanding of the actual physics of subsurface flow beyond numerical modelling, however, is a confusing experience exposing one to conflicting statements from the sides of engineers, hydrogeologists, and, for a decade or more, by the followers of free convection and density-driven flow. Within the physics of subsurface flow a number of questions arise, such as: Is water really incompressible as assumed in engineering hydraulics? How does buoyancy work? Are underground buoyancy forces generally directed vertically upwards or downwards? What is the consequential difference between hydrostatic and hydrodynamic conditions? What are the force fields causing subsurface flow for water, hydrocarbons and CO2? Is fluid flow really driven by pressure gradients as assumed in reservoir engineering? What is the effect of geothermal gradients on subsurface flow? Do convection cells and free convection exist on-shore? How does variable density flow work? Can today's numerical codes adequately determine variable density flow? Does saltwater really sink to the bottom of geologic systems due to its higher density? Aquitards create confining conditions and thereby confine fluid movements to aquifers? Does more water flow in aquifers than aquitards? The presentation will shed light on the maze of conflicting statements issued within engineering hydraulics and groundwater dynamics. It will also present a field case and its numerical modelling of variable density flow at a major industrial landfill site. The presentation will thereby foster the understanding of the correct physics involved and how this physics can be beneficially applied to practical cases
Lee, Juhyun; Moghadam, Mahdi Esmaily; Kung, Ethan; Cao, Hung; Beebe, Tyler; Miller, Yury; Roman, Beth L.; Lien, Ching-Ling; Chi, Neil C.; Marsden, Alison L.; Hsiai, Tzung K.
2013-01-01
Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis. PMID:24009714
Craven, Brent A; Paterson, Eric G; Settles, Gary S; Lawson, Michael J
2009-09-01
The canine nasal cavity contains a complex airway labyrinth, dedicated to respiratory air conditioning, filtering of inspired contaminants, and olfaction. The small and contorted anatomical structure of the nasal turbinates has, to date, precluded a proper study of nasal airflow in the dog. This study describes the development of a high-fidelity computational fluid dynamics (CFD) model of the canine nasal airway from a three-dimensional reconstruction of high-resolution magnetic resonance imaging scans of the canine anatomy. Unstructured hexahedral grids are generated, with large grid sizes ((10-100) x 10(6) computational cells) required to capture the details of the nasal airways. High-fidelity CFD solutions of the nasal airflow for steady inspiration and expiration are computed over a range of physiological airflow rates. A rigorous grid refinement study is performed, which also illustrates a methodology for verification of CFD calculations on complex unstructured grids in tortuous airways. In general, the qualitative characteristics of the computed solutions for the different grid resolutions are fairly well preserved. However, quantitative results such as the overall pressure drop and even the regional distribution of airflow in the nasal cavity are moderately grid dependent. These quantities tend to converge monotonically with grid refinement. Lastly, transient computations of canine sniffing were carried out as part of a time-step study, demonstrating that high temporal accuracy is achievable using small time steps consisting of 160 steps per sniff period. Here we demonstrate that acceptable numerical accuracy (between approximately 1% and 15%) is achievable with practical levels of grid resolution (approximately 100 x 10(6) computational cells). Given the popularity of CFD as a tool for studying flow in the upper airways of humans and animals, based on this work we recommend the necessity of a grid dependence study and quantification of numerical error when
NASA Astrophysics Data System (ADS)
Huang, L.; Gong, S. L.; Gordon, M.; Liggio, J.; Staebler, R.; Stroud, C. A.; Lu, G.; Mihele, C.; Brook, J. R.; Jia, C. Q.
2014-12-01
Many studies have shown that on-road vehicle emissions are the dominant source of ultrafine particles (UFPs; diameter < 100 nm) in urban areas and near-roadway environments. In order to advance our knowledge on the complex interactions and competition among atmospheric dilution, dispersion, and dynamics of UFPs, an aerosol dynamics-computational fluid dynamics (CFD) coupled model is developed and validated against field measurements. A unique approach of applying periodic boundary conditions is proposed to model pollutant dispersion and dynamics in one unified domain from the tailpipe level to the ambient near-road environment. This approach significantly reduces the size of the computational domain, and therefore allows fast simulation of multiple scenarios. The model is validated against measured turbulent kinetic energy (TKE) and horizontal gradient of pollution concentrations perpendicular to a major highway. Through a model sensitivity analysis, the relative importance of individual aerosol dynamical processes on the total particle number concentration (N) and particle number-size distribution (PSD) near a highway is investigated. The results demonstrate that (1) coagulation has a negligible effect on N and particle growth, (2) binary homogeneous nucleation (BHN) of H2SO4-H2O is likely responsible for elevated N closest to the road, and (3) N and particle growth are very sensitive to the condensation of semi-volatile organics (SVOCs), particle dry deposition, and the interaction between these processes. The results also indicate that, without the proper treatment of the atmospheric boundary layer (i.e., its wind profile and turbulence quantities), the nucleation rate would be underestimated by a factor of 5 in the vehicle wake region due to overestimated dilution. Therefore, introducing atmospheric boundary layer (ABL) conditions to activity-based emission models may potentially improve their performance in estimating UFP traffic emissions.
Progress in geophysical fluid dynamics
NASA Astrophysics Data System (ADS)
Robinson, Allan R.
today are powerful enough to allow realistic simulations of turbulent and planetary flows. A school of scientists and philosophers regard such simulations of computational physics as representing the first major advance in scientific methodology in centuries; scientific enterprise is now tripartite, with simulation on a par with theory and experimentation. Data assimilation involves the continual blending of observational data with dynamical model output for the best overall representation of reality. The conceptual model of nature implied is novel. The named discipline of geophysical fluid dynamics is barely three decades old. Scientifically it is an interesting time in the history of human development on earth as aspects of the dynamics of our atmosphere and oceans become solved problems. Geophysical fluid dynamicists are ready to deal with interactive and whole-earth problems, and to continue to expand the horizons of their science via the opportunities provided by space exploration. Progress is occurring in understanding climate and climate change processes which involve dynamical coupling of the oceans and the atmosphere and which cause profound biological and economic effects. Applied geophysical fluid dynamics is essential for the potential success of the International Geosphere-Biosphere Program which seeks to unite earth scientists in the next decade in the pursuit of global change research dedicated to a more habitable planet.
NASA Technical Reports Server (NTRS)
Thorp, Scott A.
1992-01-01
This presentation will discuss the development of a NASA Geometry Exchange Specification for transferring aerodynamic surface geometry between LeRC systems and grid generation software used for computational fluid dynamics research. The proposed specification is based on a subset of the Initial Graphics Exchange Specification (IGES). The presentation will include discussion of how the NASA-IGES standard will accommodate improved computer aided design inspection methods and reverse engineering techniques currently being developed. The presentation is in viewgraph format.
National Ignition Facility computational fluid dynamics modeling and light fixture case studies
Martin, R.; Bernardin, J.; Parietti, L.; Dennison, B.
1998-02-01
This report serves as a guide to the use of computational fluid dynamics (CFD) as a design tool for the National Ignition Facility (NIF) program Title I and Title II design phases at Lawrence Livermore National Laboratory. In particular, this report provides general guidelines on the technical approach to performing and interpreting any and all CFD calculations. In addition, a complete CFD analysis is presented to illustrate these guidelines on a NIF-related thermal problem.
NASA Astrophysics Data System (ADS)
Recent advances in experimental and computational fluid mechanics are discussed in a series of review essays. Topics addressed include transitions to complex flow in thermal convection, optimum hypersonic wings and wave riders, relativistic hydrodynamics, and wind-tunnel wall corrections for unsteady flow (steady wall adaptation and CFD techniques). Consideration is given to axisymmetric laminar interacting boundary layers, differential forms and fluid dynamics, breaking water waves, strong temperature gradients in turbulent wakes, and liquid-crystal 'blue' phases.
NASA Astrophysics Data System (ADS)
Takahashi, Y. O.; Takehiro, S.; Sugiyama, K.; Odaka, M.; Ishiwatari, M.; Sasaki, Y.; Nishizawa, S.; Ishioka, K.; Nakajima, K.; Hayashi, Y.
2012-12-01
) is a collection of various sample programs using ``SPML''. These sample programs provide the basekit for simple numerical experiments of geophysical fluid dynamics. For example, SPMODEL includes 1-dimensional KdV equation model, 2-dimensional barotropic, shallow water, Boussinesq models, 3-dimensional MHD dynamo models in rotating spherical shells. These models are written in the common style in harmony with SPML functions. ``Deepconv'' (Sugiyama et al., 2010) and ``Dcpam'' are a cloud resolving model and a general circulation model for the purpose of applications to the planetary atmospheres, respectively. ``Deepconv'' includes several physical processes appropriate for simulations of Jupiter and Mars atmospheres, while ``Dcpam'' does for simulations of Earth, Mars, and Venus-like atmospheres. ``Rdoc-f95'' is a automatic generator of reference manuals of Fortran90/95 programs, which is an extension of ruby documentation tool kit ``rdoc''. It analyzes dependency of modules, functions, and subroutines in the multiple program source codes. At the same time, it can list up the namelist variables in the programs.
NASA Astrophysics Data System (ADS)
Darby, S. E.; Rinaldi, M.; Rossi Romanelli, L.; Spyropoulos, E.
2003-12-01
River bank erosion often significantly contributes to the catchment sediment yield. Knowledge of the rates & controls on bank erosion events is therefore important in understanding sediment flux. In recent years progress has been made in understanding processes controlling large-scale mass failure (MF) of stream banks, but less attention has been paid to the role that direct fluvial erosion (FE) plays in bank retreat. This is an important omission, not only because FE is a significant process in its own right, but because FE also often triggers mass failure. FE models are typically of the form: E = k(τ - τ c)b where E is the bank erosion rate, τ is the applied fluid shear stress, τ c is the critical stress for entrainment of the bank material, k is an empirically-derived erodibility parameter, and b is an empirically-derived exponent, often assumed to be close to unity. To apply this model, accurate observations of applied fluid stresses, FE rates & bank erodibility are required. Recent developments in bank erosion monitoring technology [e.g. Lawler, 1993], and in the quantification of the bank erodibility parameters k and τ c using jet-testing devices [e.g. Hanson and Simon, 2001; Dapporto, 2001], offer the means of determining FE rates and bank erodibility. Nevertheless, the problem of collecting the high-resolution spatially-distributed data needed to characterise near-bank fluid stresses remains. One possible solution is to use Computational Fluid Dynamics (CFD) models as a substitute for empirical data. CFD simulations potentially offer a means of acquiring near-bank, distributed, boundary shear stress data at very high spatial resolution. In contrast, empirical data sets of comparable spatial extent and resolution are very difficult to obtain, particularly during the large (competent) flows of interest here. The critical question is therefore whether CFD-derived data are sufficiently accurate for this purpose. Herein we evaluate a series of 3-dimensional
Salim, S M; Viswanathan, Shekar; Ray, Madhumita Bhowmick
2006-12-01
Dispersion of airborne contaminants in indoor air was evaluated employing physical measurement, empirical models, and computer simulation methods. Field data collected from a tray of evaporating solvent in the laboratory were compared with computational fluid dynamics (CFD) simulations coupled with evaporation models. The results indicated that mathematical models of evaporation can be coupled with CFD simulations to produce reasonable qualitative predictions of airborne contaminant levels. The airflow pattern within a room is primarily determined by the room layout and the position of the air supply diffusers. Variations in ventilation rate did not alter the airflow pattern, thus generating a characteristic concentration profile of the airborne contaminants. PMID:17050350
Relativistic fluid dynamics. Proceedings.
NASA Astrophysics Data System (ADS)
Anile, A. M.; Choquet-Bruhat, Y.
Contents: 1. Covariant theory of conductivity in ideal fluid or solid media (B. Carter). 2. Hamiltonian techniques for relativistic fluid dynamics and stability theory (D. D. Holm). 3. Covariant fluid mechanics and thermodynamics: an introduction (W. Israel). 4. Relativistic plasmas (H. Weitzner). 5. An improved relativistic warm plasma model (A. M. Anile, S. Pennisi). 6. Relativistic extended thermodynamics II (I. Müller). 7. Relativistic extended thermodynamics: general assumptions and mathematical procedure (T. Ruggeri). 8. Relativistic hydrodynamics and heavy ion reactions (D. Strottman). 9. Some problems in relativistic hydrodynamics (C. G. van Weert).
Saho, Tatsunori; Onishi, Hideo
2015-07-01
In this study, we evaluated hemodynamics using simulated models and determined how cerebral aneurysms develop in simulated and patient-specific models based on medical images. Computational fluid dynamics (CFD) was analyzed by use of OpenFOAM software. Flow velocity, stream line, and wall shear stress (WSS) were evaluated in a simulated model aneurysm with known geometry and in a three-dimensional angiographic model. The ratio of WSS at the aneurysm compared with that at the basilar artery was 1:10 in simulated model aneurysms with a diameter of 10 mm and 1:18 in the angiographic model, indicating similar tendencies. Vortex flow occurred in both model aneurysms, and the WSS decreased in larger model aneurysms. The angiographic model provided accurate CFD information, and the tendencies of simulated and angiographic models were similar. These findings indicate that hemodynamic effects are involved in the development of aneurysms. PMID:25911446
Defraeye, Thijs; Blocken, Bert; Koninckx, Erwin; Hespel, Peter; Carmeliet, Jan
2010-08-26
This study aims at assessing the accuracy of computational fluid dynamics (CFD) for applications in sports aerodynamics, for example for drag predictions of swimmers, cyclists or skiers, by evaluating the applied numerical modelling techniques by means of detailed validation experiments. In this study, a wind-tunnel experiment on a scale model of a cyclist (scale 1:2) is presented. Apart from three-component forces and moments, also high-resolution surface pressure measurements on the scale model's surface, i.e. at 115 locations, are performed to provide detailed information on the flow field. These data are used to compare the performance of different turbulence-modelling techniques, such as steady Reynolds-averaged Navier-Stokes (RANS), with several k-epsilon and k-omega turbulence models, and unsteady large-eddy simulation (LES), and also boundary-layer modelling techniques, namely wall functions and low-Reynolds number modelling (LRNM). The commercial CFD code Fluent 6.3 is used for the simulations. The RANS shear-stress transport (SST) k-omega model shows the best overall performance, followed by the more computationally expensive LES. Furthermore, LRNM is clearly preferred over wall functions to model the boundary layer. This study showed that there are more accurate alternatives for evaluating flow around bluff bodies with CFD than the standard k-epsilon model combined with wall functions, which is often used in CFD studies in sports. PMID:20488446
Computational Fluid Dynamic Modeling of Rocket Based Combined Cycle Engine Flowfields
NASA Technical Reports Server (NTRS)
Daines, Russell L.; Merkle, Charles L.
1994-01-01
Computational Fluid Dynamic techniques are used to study the flowfield of a fixed geometry Rocket Based Combined Cycle engine operating in rocket ejector mode. Heat addition resulting from the combustion of injected fuel causes the subsonic engine flow to choke and go supersonic in the slightly divergent combustor-mixer section. Reacting flow computations are undertaken to predict the characteristics of solutions where the heat addition is determined by the flowfield. Here, adaptive gridding is used to improve resolution in the shear layers. Results show that the sonic speed is reached in the unheated portions of the flow first, while the heated portions become supersonic later. Comparison with results from another code show reasonable agreement. The coupled solutions show that the character of the combustion-based thermal choking phenomenon can be controlled reasonably well such that there is opportunity to optimize the length and expansion ratio of the combustor-mixer.
Measurement-based quantum lattice gas model of fluid dynamics in 2+1 dimensions.
Micci, Michael M; Yepez, Jeffrey
2015-09-01
Presented are quantum simulation results using a measurement-based quantum lattice gas algorithm for Navier-Stokes fluid dynamics in 2+1 dimensions. Numerical prediction of the kinematic viscosity was measured by the decay rate of an initial sinusoidal flow profile. Due to local quantum entanglement in the quantum lattice gas, the minimum kinematic viscosity in the measurement-based quantum lattice gas is lower than achievable in a classical lattice gas. The numerically predicted viscosities precisely match the theoretical predictions obtained with a mean field approximation. Uniform flow profile with double shear layers, on a 16K×8K lattice, leads to the Kelvin-Helmholtz instability, breaking up the shear layer into pairs of counter-rotating vortices that eventually merge via vortex fusion and dissipate because of the nonzero shear viscosity. PMID:26465581
Measurement-based quantum lattice gas model of fluid dynamics in 2+1 dimensions
NASA Astrophysics Data System (ADS)
Micci, Michael M.; Yepez, Jeffrey
2015-09-01
Presented are quantum simulation results using a measurement-based quantum lattice gas algorithm for Navier-Stokes fluid dynamics in 2+1 dimensions. Numerical prediction of the kinematic viscosity was measured by the decay rate of an initial sinusoidal flow profile. Due to local quantum entanglement in the quantum lattice gas, the minimum kinematic viscosity in the measurement-based quantum lattice gas is lower than achievable in a classical lattice gas. The numerically predicted viscosities precisely match the theoretical predictions obtained with a mean field approximation. Uniform flow profile with double shear layers, on a 16 K ×8 K lattice, leads to the Kelvin-Helmholtz instability, breaking up the shear layer into pairs of counter-rotating vortices that eventually merge via vortex fusion and dissipate because of the nonzero shear viscosity.
NASA Technical Reports Server (NTRS)
Fowlis, W. W.
1981-01-01
Systematic scaling or dimensional analysis reveals that certain scales of geophysical fluid flows (such as stellar, ocean, and planetary atmosphere circulations) can be accurately modeled in the laboratory using a procedure which differs from conventional engineering modeling. Rather than building a model to obtain numbers for a specific design problem, the relative effects of the significant forces are systematically varied in an attempt to deepen understanding of the effects of these forces. Topics covered include: (1) modeling a large-scale planetary atmospheric flow in a rotating cylindrical annulus; (2) achieving a radial dielectric body force; (3) spherical geophysical fluid dynamics experiments for Spacelab flights; (4) measuring flow and temperature; and (5) the possible effect of rotational or precessional disturbances on the flow in the rotating spherical containers.
NASA Astrophysics Data System (ADS)
Stockton, Gregory R.
2011-05-01
Over the last 10 years, very large government, military, and commercial computer and data center operators have spent millions of dollars trying to optimally cool data centers as each rack has begun to consume as much as 10 times more power than just a few years ago. In fact, the maximum amount of data computation in a computer center is becoming limited by the amount of available power, space and cooling capacity at some data centers. Tens of millions of dollars and megawatts of power are being annually spent to keep data centers cool. The cooling and air flows dynamically change away from any predicted 3-D computational fluid dynamic modeling during construction and as time goes by, and the efficiency and effectiveness of the actual cooling rapidly departs even farther from predicted models. By using 3-D infrared (IR) thermal mapping and other techniques to calibrate and refine the computational fluid dynamic modeling and make appropriate corrections and repairs, the required power for data centers can be dramatically reduced which reduces costs and also improves reliability.
Laurent, J; Samstag, R W; Ducoste, J M; Griborio, A; Nopens, I; Batstone, D J; Wicks, J D; Saunders, S; Potier, O
2014-01-01
To date, computational fluid dynamics (CFD) models have been primarily used for evaluation of hydraulic problems at wastewater treatment plants (WWTPs). A potentially more powerful use, however, is to simulate integrated physical, chemical and/or biological processes involved in WWTP unit processes on a spatial scale and to use the gathered knowledge to accelerate improvement in plant models for everyday use, that is, design and optimized operation. Evolving improvements in computer speed and memory and improved software for implementing CFD, as well as for integrated processes, has allowed for broader usage of this tool for understanding, troubleshooting, and optimal design of WWTP unit processes. This paper proposes a protocol for an alternative use of CFD in process modelling, as a way to gain insight into complex systems leading to improved modelling approaches used in combination with the IWA activated sludge models and other kinetic models. PMID:25429444
Otto Laporte Lecture: Fluid Dynamics Prize Talk: Simple Models for Turbulent Flows
NASA Astrophysics Data System (ADS)
Pope, Stephen B.
2009-11-01
We focus on the modeling of two turbulent flows: dispersion from a line source in grid turbulence; and, a lifted non-premixed turbulent jet flame. Stochastic Lagrangian models and PDF methods are described, and are shown to model these flows satisfactorily. For the line source, a Lagrangian approach is taken, with the Langevin equation modeling the velocity following a fluid particle, and with a simple relaxation model for the particle temperature. Comparison with experimental data shows that the resulting model describes accurately the dispersion from single and multiple line sources. These simple stochastic Lagrangian models are then applied to the much more challenging case of a lifted non-premixed jet flame. The stochastic Lagrangian models form the basis for a particle/mesh numerical method for solving a modeled transport equation for the Eulerian joint probability density function (PDF) of velocity and composition. The PDF calculations are in excellent agreement with the experimental data, and exhibit the observed extreme sensitivity of the flame to the temperature of the co-flow. The PDF model calculations presented clearly demonstrate that simple models can be very useful, even though aspects of their behavior may be inaccurate or incomplete. The shortcomings of the Langevin equation are examined, and more advanced models (designed to overcome some of these shortcomings) are described. These include models for fluid-particle acceleration, including the effects of intermittency; models accounting for mean shear, which are correct in the rapid- distortion limit; and models designed for use in conjunction with large-eddy simulations (LES).
NASA Technical Reports Server (NTRS)
Connolly, Joseph W.; Friedlander, David; Kopasakis, George
2015-01-01
This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.
NASA Technical Reports Server (NTRS)
Connolly, Joseph W.; Friedlander, David; Kopasakis, George
2014-01-01
This paper covers the development of an integrated nonlinear dynamic simulation for a variable cycle turbofan engine and nozzle that can be integrated with an overall vehicle Aero-Propulso-Servo-Elastic (APSE) model. A previously developed variable cycle turbofan engine model is used for this study and is enhanced here to include variable guide vanes allowing for operation across the supersonic flight regime. The primary focus of this study is to improve the fidelity of the model's thrust response by replacing the simple choked flow equation convergent-divergent nozzle model with a MacCormack method based quasi-1D model. The dynamic response of the nozzle model using the MacCormack method is verified by comparing it against a model of the nozzle using the conservation element/solution element method. A methodology is also presented for the integration of the MacCormack nozzle model with the variable cycle engine.
Sverdlova, N S; Arkali, F; Witzel, U; Perry, S F
2013-10-01
Respiratory evaporative cooling is an important mechanism of temperature control in bird. A computational simulation of the breathing cycle, heat and water loss in anatomical avian trachea/air sac model has not previously been conducted. We report a first attempt to simulate a breathing cycle in a three-dimensional model of avian trachea and air sacs (domestic fowl) using transient computational fluid dynamics. The airflow in the trachea of the model is evoked by changing the volume of the air sacs based on the measured tidal volume and inspiratory/expiratory times for the domestic fowl. We compare flow parameters and heat transfer results with in vivo data and with our previously reported results for a two-dimensional model. The total respiratory heat loss corresponds to about 13-19% of the starvation metabolic rate of domestic fowl. The present study can lend insight into a possible thermoregulatory function in species with long necks and/or a very long trachea, as found in swans and birds of paradise. Assuming the structure of the sauropod dinosaur respiratory system was close to avian, the simulation of the respiratory temperature control (using convective and evaporative cooling) in the extensively experimentally studied domestic fowl may also help in making simulations of respiratory heat control in these extinct animals. PMID:23797184
Poon, Eric K W; Hayat, Umair; Thondapu, Vikas; Ooi, Andrew S H; Ul Haq, Muhammad Asrar; Moore, Stephen; Foin, Nicolas; Tu, Shengxian; Chin, Cheng; Monty, Jason P; Marusic, Ivan; Barlis, Peter
2015-08-01
Percutaneous coronary intervention (PCI) has shown a high success rate in the treatment of coronary artery disease. The decision to perform PCI often relies on the cardiologist's visual interpretation of coronary lesions during angiography. This has inherent limitations, particularly due to the low resolution and two-dimensional nature of angiography. State-of-the-art modalities such as three-dimensional quantitative coronary angiography, optical coherence tomography and invasive fractional flow reserve (FFR) may improve clinicians' understanding of both the anatomical and physiological importance of coronary lesions. While invasive FFR is the gold standard technique for assessment of the haemodynamic significance of coronary lesions, recent studies have explored a surrogate for FFR derived solely from three-dimensional reconstruction of the invasive angiogram, and therefore eliminating need for a pressure wire. Utilizing advanced computational fluid dynamics research, this virtual fractional flow reserve (vFFR) has demonstrated reasonable correlation with invasive measurements and remains an intense area of ongoing study. However, at present, several limitations and computational fluid dynamic assumptions may preclude vFFR from widespread clinical use. This review demonstrates the tight integration of advanced three-dimensional imaging techniques and vFFR in assessing coronary artery disease, reviews the advantages and disadvantages of such techniques and attempts to provide a glimpse of how such advances may benefit future clinical decision-making during PCI. PMID:26247271
Ranganathan, Panneerselvam; Gu, Sai
2016-08-01
The present work concerns with CFD modelling of biomass fast pyrolysis in a fluidised bed reactor. Initially, a study was conducted to understand the hydrodynamics of the fluidised bed reactor by investigating the particle density and size, and gas velocity effect. With the basic understanding of hydrodynamics, the study was further extended to investigate the different kinetic schemes for biomass fast pyrolysis process. The Eulerian-Eulerian approach was used to model the complex multiphase flows in the reactor. The yield of the products from the simulation was compared with the experimental data. A good comparison was obtained between the literature results and CFD simulation. It is also found that CFD prediction with the advanced kinetic scheme is better when compared to other schemes. With the confidence obtained from the CFD models, a parametric study was carried out to study the effect of biomass particle type and size and temperature on the yield of the products. PMID:26927234
Computational Fluid Dynamics Symposium on Aeropropulsion
NASA Technical Reports Server (NTRS)
1991-01-01
Recognizing the considerable advances that have been made in computational fluid dynamics, the Internal Fluid Mechanics Division of NASA Lewis Research Center sponsored this symposium with the objective of providing a forum for exchanging information regarding recent developments in numerical methods, physical and chemical modeling, and applications. This conference publication is a compilation of 4 invited and 34 contributed papers presented in six sessions: algorithms one and two, turbomachinery, turbulence, components application, and combustors. Topics include numerical methods, grid generation, chemically reacting flows, turbulence modeling, inlets, nozzles, and unsteady flows.
Computational fluid dynamics symposium on aeropropulsion
Not Available
1991-01-01
Recognizing the considerable advances that have been made in computational fluid dynamics, the Internal Fluid Mechanics Division of NASA Lewis Research Center sponsored this symposium with the objective of providing a forum for exchanging information regarding recent developments in numerical methods, physical and chemical modeling, and applications. This conference publication is a compilation of 4 invited and 34 contributed papers presented in six sessions: algorithms one and two, turbomachinery, turbulence, components application, and combustors. Topics include numerical methods, grid generation, chemically reacting flows, turbulence modeling, inlets, nozzles, and unsteady flows.
NASA Technical Reports Server (NTRS)
Gonzalez, Marianne; Quinn, Jacqueline; Captain, Janine; Santiago-Bond, Josephine; Starr, Stanley
2015-01-01
The Resource Prospector (RP) mission with the Regolith and Environment Science and Oxygen Lunar Volatile Extraction (RESOLVE) payload aims to show the presence of water in lunar regolith, and establish a proving ground for NASAs mission to Mars. One of the analysis is performed by the Lunar Advanced Volatiles Analysis (LAVA) subsystem, which consists of a fluid network that facilitates the transport of volatile samples to a gas chromatograph and mass spectrometer (GC-MS) instrument. The understanding of fluid dynamics directed from the GC to the MS is important due to the influence of flow rates and pressures that affect the accuracy of and prevent the damage to the overall GC-MS instrument. The micro-scale capillary fluid network within the GC alone has various lengths and inner-diameters; therefore, determination of pressure differentials and flow rates are difficult to model computationally, with additional complexity from the vacuum conditions in space and lack of a lunar atmosphere. A series of tests were performed on an experimental set-up of the system where the inner diameters of the GC transfer line connecting to the MS were varied. The effect on chromatography readings were also studied by applying these lines onto a GC instrument. It was found that a smaller inner diameter transfer line resulted in a lower flow rate, as well as a lower pressure differential across the thermal conductivity detector (TCD) unit of the GC and a negligible pressure drop across the mock-up capillary column. The chromatography was affected with longer retention times and broader peak integrations. It was concluded that a 0.050 mm inner diameter line still proved most suitable for the systems flow rate preferences. In addition, it was evident that this small transfer line portrayed some expense to GC signal characteristics and the wait time for steady-state operation.
Computational fluid dynamics uses in fluid dynamics/aerodynamics education
NASA Technical Reports Server (NTRS)
Holst, Terry L.
1994-01-01
The field of computational fluid dynamics (CFD) has advanced to the point where it can now be used for the purpose of fluid dynamics physics education. Because of the tremendous wealth of information available from numerical simulation, certain fundamental concepts can be efficiently communicated using an interactive graphical interrogation of the appropriate numerical simulation data base. In other situations, a large amount of aerodynamic information can be communicated to the student by interactive use of simple CFD tools on a workstation or even in a personal computer environment. The emphasis in this presentation is to discuss ideas for how this process might be implemented. Specific examples, taken from previous publications, will be used to highlight the presentation.
ERIC Educational Resources Information Center
Bird, R. Byron
1980-01-01
Problems in polymer fluid dynamics are described, including development of constitutive equations, rheometry, kinetic theory, flow visualization, heat transfer studies, flows with phase change, two-phase flow, polymer unit operations, and drag reduction. (JN)
Guyot, Y; Luyten, F P; Schrooten, J; Papantoniou, I; Geris, L
2015-12-01
Bone tissue engineering strategies use flow through perfusion bioreactors to apply mechanical stimuli to cells seeded on porous scaffolds. Cells grow on the scaffold surface but also by bridging the scaffold pores leading a fully filled scaffold following the scaffold's geometric characteristics. Current computational fluid dynamic approaches for tissue engineering bioreactor systems have been mostly carried out for empty scaffolds. The effect of 3D cell growth and extracellular matrix formation (termed in this study as neotissue growth), on its surrounding fluid flow field is a challenge yet to be tackled. In this work a combined approach was followed linking curvature driven cell growth to fluid dynamics modeling. The level-set method (LSM) was employed to capture neotissue growth driven by curvature, while the Stokes and Darcy equations, combined in the Brinkman equation, provided information regarding the distribution of the shear stress profile at the neotissue/medium interface and within the neotissue itself during growth. The neotissue was assumed to be micro-porous allowing flow through its structure while at the same time allowing the simulation of complete scaffold filling without numerical convergence issues. The results show a significant difference in the amplitude of shear stress for cells located within the micro-porous neo-tissue or at the neotissue/medium interface, demonstrating the importance of taking along the neotissue in the calculation of the mechanical stimulation of cells during culture.The presented computational framework is used on different scaffold pore geometries demonstrating its potential to be used a design as tool for scaffold architecture taking into account the growing neotissue. Biotechnol. Bioeng. 2015;112: 2591-2600. © 2015 Wiley Periodicals, Inc. PMID:26059101
Dresser, Thomas J.; Dotson, Curtis L.; Fisher, Richard K.; Graf, Michael J.; Richmond, Marshall C.; Rakowski, Cynthia L.; Carlson, Thomas J.; Mathur, Dilip; Heisey, Paul G.
2007-10-10
This paper, the second part of a 2 part paper, discusses the use of Computational Fluid Dynamics (CFD) to gain further insight into the results of fish release testing conducted to evaluate the modifications made to upgrade Unit 8 at Wanapum Dam. Part 1 discusses the testing procedures and fish passage survival. Grant PUD is working with Voith Siemens Hydro (VSH) and the Pacific Northwest National Laboratory (PNNL) of DOE and Normandeau Associates in this evaluation. VSH has prepared the geometry for the CFD analysis corresponding to the four operating conditions tested with Unit 9, and the 5 operating conditions tested with Unit 8. Both VSH and PNNL have conducting CFD simulations of the turbine intakes, stay vanes, wicket gates, turbine blades and draft tube of the units. Primary objectives of the analyses were: • determine estimates of where the inserted fish passed the turbine components • determine the characteristics of the flow field along the paths calculated for pressure, velocity gradients and acceleration associated with fish sized bodies • determine the velocity gradients at the structures where fish to structure interaction is predicted. • correlate the estimated fish location of passage with observed injuries • correlate the calculated pressure and acceleration with the information recorded with the sensor fish • utilize the results of the analysis to further interpret the results of the testing. This paper discusses the results of the CFD analyses made to assist the interpretation of the fish test results.
Pathapati, Subbu-Srikanth; Sansalone, John J
2011-07-01
Computational fluid dynamics (CFD) is emerging as a model for resolving the fate of particulate matter (PM) by unit operations subject to rainfall-runoff loadings. However, compared to steady flow CFD models, there are greater computational requirements for unsteady hydrodynamics and PM loading models. Therefore this study examines if integrating a stepwise steady flow CFD model can reproduce PM separation by common unit operations loaded by unsteady flow and PM loadings, thereby reducing computational effort. Utilizing monitored unit operation data from unsteady events as a metric, this study compares the two CFD modeling approaches for a hydrodynamic separator (HS), a primary clarifier (PC) tank, and a volumetric clarifying filtration system (VCF). Results indicate that while unsteady CFD models reproduce PM separation of each unit operation, stepwise steady CFD models result in significant deviation for HS and PC models as compared to monitored data; overestimating the physical size requirements of each unit required to reproduce monitored PM separation results. In contrast, the stepwise steady flow approach reproduces PM separation by the VCF, a combined gravitational sedimentation and media filtration unit operation that provides attenuation of turbulent energy and flow velocity. PMID:21644537
Ong, Robert H.; King, Andrew J. C.; Mullins, Benjamin J.; Cooper, Timothy F.; Caley, M. Julian
2012-01-01
We present Computational Fluid Dynamics (CFD) models of the coupled dynamics of water flow, heat transfer and irradiance in and around corals to predict temperatures experienced by corals. These models were validated against controlled laboratory experiments, under constant and transient irradiance, for hemispherical and branching corals. Our CFD models agree very well with experimental studies. A linear relationship between irradiance and coral surface warming was evident in both the simulation and experimental result agreeing with heat transfer theory. However, CFD models for the steady state simulation produced a better fit to the linear relationship than the experimental data, likely due to experimental error in the empirical measurements. The consistency of our modelling results with experimental observations demonstrates the applicability of CFD simulations, such as the models developed here, to coral bleaching studies. A study of the influence of coral skeletal porosity and skeletal bulk density on surface warming was also undertaken, demonstrating boundary layer behaviour, and interstitial flow magnitude and temperature profiles in coral cross sections. Our models compliment recent studies showing systematic changes in these parameters in some coral colonies and have utility in the prediction of coral bleaching. PMID:22701582
Ferroni, Marco; Giusti, Serena; Nascimento, Diana; Silva, Ana; Boschetti, Federica; Ahluwalia, Arti
2016-08-01
The architecture and dynamic physical environment of tissues can be recreated in-vitro by combining 3D porous scaffolds and bioreactors able to apply controlled mechanical stimuli on cells. In such systems, the entity of the stimuli and the distribution of nutrients within the engineered construct depend on the micro-structure of the scaffolds. In this work, we present a new approach for optimizing computational fluid-dynamics (CFD) models for the investigation of fluid-induced forces generated by cyclic squeeze pressure within a porous construct, coupled with oxygen consumption of cardiomyocytes. A 2D axial symmetric macro-scaled model of a squeeze pressure bioreactor chamber was used as starting point for generating time dependent pressure profiles. Subsequently the fluid movement generated by the pressure fields was coupled with a complete 3D micro-scaled model of a porous protein cryogel. Oxygen transport and consumption inside the scaffold was evaluated considering a homogeneous distribution of cardiomyocytes throughout the structure, as confirmed by preliminary cell culture experiments. The results show that a 3D description of the system, coupling a porous geometry and time dependent pressure driven flow with fluid-structure-interaction provides an accurate and meaningful description of the microenvironment in terms of shear stress and oxygen distribution than simple stationary 2D models. PMID:27189671
NASA Astrophysics Data System (ADS)
Ogilvie, Gordon I.
2016-06-01
These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is `frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.
NASA Astrophysics Data System (ADS)
Ogilvie, Gordon I.
2016-06-01
> These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is `frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.
Martin, Kelly J; Picioreanu, Cristian; Nerenberg, Robert
2013-09-01
A two-dimensional, particle-based biofilm model coupled with mass transport and computational fluid dynamics was developed to simulate autotrophic denitrification in a spiral-wound membrane biofilm reactor (MBfR), where hydrogen is supplied via hollow-fiber membrane fabric. The spiral-wound configuration consists of alternating layers of plastic spacer net and membrane fabric that create rows of flow channels, with the top and bottom walls comprised of membranes. The transversal filaments of the spacer partially obstruct the channel flow, producing complex mixing and shear patterns that require multidimensional representation. This study investigated the effect of hydrogen and nitrate concentrations, as well as spacer configuration, on biofilm development and denitrification fluxes. The model results indicate that the cavity spacer filaments, which rest on the bottom membranes, cause uneven biofilm growth. Most biofilm resided on the bottom membranes, only in the wake of the filaments where low shear zones formed. In this way, filament configuration may help achieve a desired biofilm thickness. For the conditions tested in this study, the highest nitrate fluxes were attained by minimizing the filament diameter and maximizing the filament spacing. This lowered the shear stress at the top membranes, allowing for more biofilm growth. For the scenarios studied, biomass limitation at the top membranes hindered performance more significantly than diffusion limitation in the thick biofilms at the bottom membranes. The results also highlighted the importance of two-dimensional modeling to capture uneven biofilm growth on a substratum with geometrical complexity. PMID:23774188
Santoro, Domenico; Crapulli, Ferdinando; Raisee, Mehrdad; Raspa, Giuseppe; Haas, Charles N
2015-06-16
Wastewater disinfection processes are typically designed according to heuristics derived from batch experiments in which the interaction among wastewater quality, reactor hydraulics, and inactivation kinetics is often neglected. In this paper, a computational fluid dynamics (CFD) study was conducted in a nondeterministic (ND) modeling framework to predict the Escherichia coli inactivation by peracetic acid (PAA) in municipal contact tanks fed by secondary settled wastewater effluent. The extent and variability associated with the observed inactivation kinetics were both satisfactorily predicted by the stochastic inactivation model at a 95% confidence level. Moreover, it was found that (a) the process variability induced by reactor hydraulics is negligible when compared to the one caused by inactivation kinetics, (b) the PAA dose required for meeting regulations is dictated equally by the fixed limit of the microbial concentration as well as its probability of occurrence, and (c) neglecting the probability of occurrence during process sizing could lead to an underestimation of the PAA dose required by as much as 100%. Finally, the ND-CFD model was used to generate sizing information in the form of probabilistic disinfection curves relating E. coli inactivation and probability of occurrence with the average PAA dose and PAA residual concentration at the outlet of the contact tank. PMID:25938730
Liu, Huolong; Li, Mingzhong
2014-11-20
In this work a two-compartmental population balance model (TCPBM) was proposed to model a pulsed top-spray fluidized bed granulation. The proposed TCPBM considered the spatially heterogeneous granulation mechanisms of the granule growth by dividing the granulator into two perfectly mixed zones of the wetting compartment and drying compartment, in which the aggregation mechanism was assumed in the wetting compartment and the breakage mechanism was considered in the drying compartment. The sizes of the wetting and drying compartments were constant in the TCPBM, in which 30% of the bed was the wetting compartment and 70% of the bed was the drying compartment. The exchange rate of particles between the wetting and drying compartments was determined by the details of the flow properties and distribution of particles predicted by the computational fluid dynamics (CFD) simulation. The experimental validation has shown that the proposed TCPBM can predict evolution of the granule size and distribution within the granulator under different binder spray operating conditions accurately. PMID:25181553
Modeling Tools for Propulsion Analysis and Computational Fluid Dynamics on the Internet
NASA Technical Reports Server (NTRS)
Muss, J. A.; Johnson, C. W.; Gotchy, M. B.
2000-01-01
The existing RocketWeb(TradeMark) Internet Analysis System (httr)://www.iohnsonrockets.com/rocketweb) provides an integrated set of advanced analysis tools that can be securely accessed over the Internet. Since these tools consist of both batch and interactive analysis codes, the system includes convenient methods for creating input files and evaluating the resulting data. The RocketWeb(TradeMark) system also contains many features that permit data sharing which, when further developed, will facilitate real-time, geographically diverse, collaborative engineering within a designated work group. Adding work group management functionality while simultaneously extending and integrating the system's set of design and analysis tools will create a system providing rigorous, controlled design development, reducing design cycle time and cost.
Kuttler, Andreas; Dimke, Thomas; Kern, Steven; Helmlinger, Gabriel; Stanski, Donald; Finelli, Luca A
2010-12-01
We introduce how biophysical modeling in pharmaceutical research and development, combining physiological observations at the tissue, organ and system level with selected drug physiochemical properties, may contribute to a greater and non-intuitive understanding of drug pharmacokinetics and therapeutic design. Based on rich first-principle knowledge combined with experimental data at both conception and calibration stages, and leveraging our insights on disease processes and drug pharmacology, biophysical modeling may provide a novel and unique opportunity to interactively characterize detailed drug transport, distribution, and subsequent therapeutic effects. This innovative approach is exemplified through a three-dimensional (3D) computational fluid dynamics model of the spinal canal motivated by questions arising during pharmaceutical development of one molecular therapy for spinal cord injury. The model was based on actual geometry reconstructed from magnetic resonance imaging data subsequently transformed in a parametric 3D geometry and a corresponding finite-volume representation. With dynamics controlled by transient Navier-Stokes equations, the model was implemented in a commercial multi-physics software environment established in the automotive and aerospace industries. While predictions were performed in silico, the underlying biophysical models relied on multiple sources of experimental data and knowledge from scientific literature. The results have provided insights into the primary factors that can influence the intrathecal distribution of drug after lumbar administration. This example illustrates how the approach connects the causal chain underlying drug distribution, starting with the technical aspect of drug delivery systems, through physiology-driven drug transport, then eventually linking to tissue penetration, binding, residence, and ultimately clearance. Currently supporting our drug development projects with an improved understanding of systems
NASA Technical Reports Server (NTRS)
Dougherty, Sam; West, Jeff; Droege, Alan; Wilson, Josh; Liever, Peter; Slaby, Matthew
2006-01-01
This paper discusses the Space Shuttle Lift-off CFD model developed for potential Lift-off Debris transport for return-to-flight. The Lift-off portion of the flight is defined as the time starting with tanking of propellants until tower clear, approximately T0+6 seconds, where interactions with the launch pad cease. A CFD model containing the Space Shuttle and launch Pad geometry has been constructed and executed. Simplifications required in the construction of the model are presented and discussed. A body-fitted overset grid of up to 170 million grid points was developed which allowed positioning of the Vehicle relative to the Launch Pad over the first six seconds of Climb-Out. The CFD model works in conjunction with a debris particle transport model and a debris particle impact damage tolerance model. These models have been used to assess the interactions of the Space Shuttle plumes, the wind environment, and their interactions with each other and the Launch Pad and their ultimate effect on potential debris during Lift-off.
Computational fluid dynamics modeling of gas dispersion in multi impeller bioreactor.
Ahmed, Syed Ubaid; Ranganathan, Panneerselvam; Pandey, Ashok; Sivaraman, Savithri
2010-06-01
In the present study, experiments have been carried out to identify various flow regimes in a dual Rushton turbines stirred bioreactor for different gas flow rates and impeller speeds. The hydrodynamic parameters like fractional gas hold-up, power consumption and mixing time have been measured. A two fluid model along with MUSIG model to handle polydispersed gas flow has been implemented to predict the various flow regimes and hydrodynamic parameters in the dual turbines stirred bioreactor. The computational model has been mapped on commercial solver ANSYS CFX. The flow regimes predicted by numerical simulations are validated with the experimental results. The present model has successfully captured the flow regimes as observed during experiments. The measured gross flow characteristics like fractional gas hold-up, and mixing time have been compared with numerical simulations. Also the effect of gas flow rate and impeller speed on gas hold-up and power consumption have been investigated. PMID:20471599
NASA Technical Reports Server (NTRS)
Benton, E. R.
1986-01-01
A spherical harmonic representation of the geomagnetic field and its secular variation for epoch 1980, designated GSFC(9/84), is derived and evaluated. At three epochs (1977.5, 1980.0, 1982.5) this model incorporates conservation of magnetic flux through five selected patches of area on the core/mantle boundary bounded by the zero contours of vertical magnetic field. These fifteen nonlinear constraints are included like data in an iterative least squares parameter estimation procedure that starts with the recently derived unconstrained field model GSFC (12/83). Convergence is approached within three iterations. The constrained model is evaluated by comparing its predictive capability outside the time span of its data, in terms of residuals at magnetic observatories, with that for the unconstrained model.
Computational fluid dynamics modeling of the operation of a flame ionization sensor
Huckaby, E.D.; Chorpening, B.; Thornton, J.
2007-01-01
The sensors and controls research group at the United States Department of Energy (DOE) National Energy Technology Laboratory (NETL) is continuing to develop the Combustion Control and Diagnostics Sensor (CCADS) for gas turbine applications. CCADS measures the electrical conduction of the charged species generated during the combustion process to detect flashback and combustion instabilities, and to monitor equivalence ratio. As part of this effort, combustion models are being developed which include the interaction between the electric field and the transport processes of the charged species. The primary combustion process is computed using a flame wrinkling model developed by Weller et al. (1998). A sub-model for the transport of charged species is attached to this model. The formulation of the charged-species model is similar to that applied by Penderson and Brown (1993) for the simulation of laminar flames. Using the above procedure, numerical simulations are performed and the results are compared with experimental current measurements. Quantitative agreement with experiment was not obtained, however the model does display similar sensitivity to flow and operating conditions as observed in experiments.
Computational Fluid Dynamics Modeling of the Operation of a Flame Ionization Sensor
Huckaby, E.D.; Chorpening, B.T.; Thornton, J.D.
2007-03-01
The sensors and controls research group at the United States Department of Energy (DOE) National Energy Technology Laboratory (NETL) is continuing to develop the Combustion Control and Diagnostics Sensor (CCADS) for gas turbine applications. CCADS uses the electrical conduction of the charged species generated during the combustion process to detect combustion instabilities and monitor equivalence ratio. As part of this effort, combustion models are being developed which include the interaction between the electric field and the transport of charged species. The primary combustion process is computed using a flame wrinkling model (Weller et. al. 1998) which is a component of the OpenFOAM toolkit (Jasak et. al. 2004). A sub-model for the transport of charged species is attached to this model. The formulation of the charged-species model similar that applied by Penderson and Brown (1993) for the simulation of laminar flames. The sub-model consists of an additional flux due to the electric field (drift flux) added to the equations for the charged species concentrations and the solution the electric potential from the resolved charge density. The subgrid interactions between the electric field and charged species transport have been neglected. Using the above procedure, numerical simulations are performed and the results compared with several recent CCADS experiments.
Synchronized Chaos in Geophysical Fluid Dynamics and in the Predictive Modeling of Natural Systems
NASA Astrophysics Data System (ADS)
Duane, Gregory S.
2008-03-01
The ubiquitous phenomenon of synchronization among regular oscillators in Nature has been shown, in the past two decades, to extend to chaotic systems. Despite sensitive dependence on initial conditions, two chaotic systems will commonly fall into synchronized motion along their strange attractors when only some of the many degrees of freedom of one system are coupled to corresponding variables in the other. In geophysical fluid models, synchronization can mediate scale interactions, so that coupling of degrees of freedom that describe medium-scale components of the flow can result in synchronization, or partial synchronization, at all scales. Chaos synchronization has been used to interpret non-local "teleconnection" patterns in the Earth's climate system and to predict new ones. In the realm of practical meteorology, the fact that two PDE systems, conceived as "truth" and "model", respectively, can be made to synchronize when coupled at only a discrete set of points, explains how observations at a discrete set of weather stations can be sufficient for weather prediction by a synchronously coupled model. Minimizing synchronization error leads to general recipes for assimilation of observed data into a running model that systematize the treatment of nonlinearities in the dynamical equations. Equations can generally be added to adapt parameters as well as states as the model is running, so that the model "learns". The synchronization view of predictive modelling extends to any translationally- any PDE with constant coefficients, the general form of physical theories. The reliance on synchronicity as an organizing principle in Nature, alternative to causality, has philosophical roots in the collaboration of Carl Jung and Wolfgang Pauli, on the one hand, and in traditions outside of European science, on the other.
Platfoot, J.H.; Wendel, M.W.; Williams, P.T.
1996-10-01
This report describes the simulation of the dispersion and dilution of dissolved or finely suspended contaminants entering the Clinch river from White Oak Creek. The work is accomplished through the application of a commercial computational fluid dynamics (CFD) solver. This study assumes that contaminants originating in the White Oak Creed watershed, which drains Oak Ridge National Laboratory, will eventually reach the mouth of White Oak Creek and be discharged into the clinch River. The numerical model was developed to support the analysis of the off-site consequences of releases from the ORNL liquid low-level waste system. The system contains storage tanks and transfer lines in Bethel Valley and Melton Valley. Under certain failure modes, liquid low-level waste could be released to White Oak Creek or Melton Branch to White Oak Creek and eventually be discharged to the Clinch River. Since the Clinch River has unrestricted access by the public and water usage from the Clinch River is not controlled by the Department of Energy, such a liquid low-level waste spill would create the possibility of public exposure to the contaminant. This study is limited to the dispersion of the contaminants downstream of the confluence of White Oak Creek.
NASA Technical Reports Server (NTRS)
1989-01-01
An overview of computational fluid dynamics (CFD) activities at the Langley Research Center is given. The role of supercomputers in CFD research, algorithm development, multigrid approaches to computational fluid flows, aerodynamics computer programs, computational grid generation, turbulence research, and studies of rarefied gas flows are among the topics that are briefly surveyed.
Helicon thruster plasma modeling: Two-dimensional fluid-dynamics and propulsive performances
NASA Astrophysics Data System (ADS)
Ahedo, Eduardo; Navarro-Cavallé, Jaume
2013-04-01
An axisymmetric macroscopic model of the magnetized plasma flow inside the helicon thruster chamber is derived, assuming that the power absorbed from the helicon antenna emission is known. Ionization, confinement, subsonic flows, and production efficiency are discussed in terms of design and operation parameters. Analytical solutions and simple scaling laws for ideal plasma conditions are obtained. The chamber model is then matched with a model of the external magnetic nozzle in order to characterize the whole plasma flow and assess thruster performances. Thermal, electric, and magnetic contributions to thrust are evaluated. The energy balance provides the power conversion between ions and electrons in chamber and nozzle, and the power distribution among beam power, ionization losses, and wall losses. Thruster efficiency is assessed, and the main causes of inefficiency are identified. The thermodynamic behavior of the collisionless electron population in the nozzle is acknowledged to be poorly known and crucial for a complete plasma expansion and good thrust efficiency.
NASA Astrophysics Data System (ADS)
Le Floch, Francois; Harris, Wesley L.
2009-11-01
A novel methodology has been developed to address sickle cell disease, based on highly descriptive mathematical models for blood flow in the capillaries. Our investigations focus on the coupling between oxygen delivery and red blood cell dynamics, which is crucial to understanding sickle cell crises and is unique to this blood disease. The main part of our work is an extensive study of blood dynamics through simulations of red cells deforming within the capillary vessels, and relies on the use of a large mathematical system of equations describing oxygen transfer, blood plasma dynamics and red cell membrane mechanics. This model is expected to lead to the development of new research strategies for sickle cell disease. Our simulation model could be used not only to assess current researched remedies, but also to spur innovative research initiatives, based on our study of the physical properties coupled in sickle cell disease.
Helicon thruster plasma modeling: Two-dimensional fluid-dynamics and propulsive performances
Ahedo, Eduardo; Navarro-Cavalle, Jaume
2013-04-15
An axisymmetric macroscopic model of the magnetized plasma flow inside the helicon thruster chamber is derived, assuming that the power absorbed from the helicon antenna emission is known. Ionization, confinement, subsonic flows, and production efficiency are discussed in terms of design and operation parameters. Analytical solutions and simple scaling laws for ideal plasma conditions are obtained. The chamber model is then matched with a model of the external magnetic nozzle in order to characterize the whole plasma flow and assess thruster performances. Thermal, electric, and magnetic contributions to thrust are evaluated. The energy balance provides the power conversion between ions and electrons in chamber and nozzle, and the power distribution among beam power, ionization losses, and wall losses. Thruster efficiency is assessed, and the main causes of inefficiency are identified. The thermodynamic behavior of the collisionless electron population in the nozzle is acknowledged to be poorly known and crucial for a complete plasma expansion and good thrust efficiency.
Fluid dynamic modelling of renal pelvic pressure during endoscopic stone removal
NASA Astrophysics Data System (ADS)
Oratis, Alexandros; Subasic, John; Bird, James; Eisner, Brian
2015-11-01
Endoscopic kidney stone removal procedures are known to increase internal pressure in the renal pelvis, the kidney's urinary collecting system. High renal pelvic pressure incites systemic absorption of irrigation fluid, which can increase the risk of postoperative fever and sepsis or the unwanted absorption of electrolytes. Urologists choose the appropriate surgical procedure based on patient history and kidney stone size. However, no study has been conducted to compare the pressure profiles of each procedure, nor is there a precise sense of how the renal pelvic pressure scales with various operational parameters. Here we develop physical models for the flow rates and renal pelvic pressure for various procedures. We show that the results of our models are consistent with existing urological data on each procedure and that the models can predict pressure profiles where data is unavailable.
Response Surface Modeling of Combined-Cycle Propulsion Components using Computational Fluid Dynamics
NASA Technical Reports Server (NTRS)
Steffen, C. J., Jr.
2002-01-01
Three examples of response surface modeling with CFD are presented for combined cycle propulsion components. The examples include a mixed-compression-inlet during hypersonic flight, a hydrogen-fueled scramjet combustor during hypersonic flight, and a ducted-rocket nozzle during all-rocket flight. Three different experimental strategies were examined, including full factorial, fractionated central-composite, and D-optimal with embedded Plackett-Burman designs. The response variables have been confined to integral data extracted from multidimensional CFD results. Careful attention to uncertainty assessment and modeling bias has been addressed. The importance of automating experimental setup and effectively communicating statistical results are emphasized.
NASA Astrophysics Data System (ADS)
Sandbach, S. D.; Hardy, R. J.; Lane, S. N.; Ashworth, P. J.; Parsons, D. R.
2010-12-01
Our understanding of large rivers is limited due to the difficulties in obtaining field data at these large scales. Data rich results may be obtained using computational fluid dynamic (CFD) models permitting the investigation of detailed flow patterns that would otherwise not be available. However, the application of these models to large rivers is not without its own complications and has yet to be fully developed. This is the result of two limiting factors, our inability; i) to design numerically stable meshes for complex topographies at these spatial resolutions; and; ii) to collect high resolution data appropriate for the boundary conditions of the numerical scheme. Here, we demonstrate a five-term mass-flux scaling algorithm (MFSA) for including bed topography in a very large river, where the discretised form of the mass and momentum equations are modified using a numerical blockage. Converged solutions were obtained using the Reynolds-averaged Navier stokes (RANS) equations modelling turbulence with a κ-ɛ RNG turbulence model. The boundary conditions were supplied from a field investigation of the Rio Paraná upstream of the Paraguay-Paraná confluence. A 38 km long reach was investigated where topographic and velocity data was collected using an acoustic Doppler current profiler (aDcp) and a single beam echo sounder. The model was validated against the aDcp data and in general showed good agreement. The model was then used to explore the impacts of roughness height upon key characteristics of the 3D flow field in large rivers. The results demonstrate the importance of topographic forcing on determining flow structures including the detection of large helical flow structures.
Xu, Chun; Sin, SangHun; McDonough, Joseph M; Udupa, Jayaram K; Guez, Allon; Arens, Raanan; Wootton, David M
2006-01-01
Computational fluid dynamic (CFD) analysis was used to model the effect of airway geometry on internal pressure in the upper airway of three children with obstructive sleep apnea syndrome (OSAS), and three controls. Model geometry was reconstructed from magnetic resonance images obtained during quiet tidal breathing, meshed with an unstructured grid, and solved at normative peak resting flow. The unsteady Reynolds-averaged Navier-Stokes equations were solved with steady flow boundary conditions in inspiration and expiration, using a two-equation low-Reynolds number turbulence model. Model results were validated using an in-vitro scale model, unsteady flow simulation, and reported nasal resistance measurements in children. Pharynx pressure drop strongly correlated to airway area restriction. Inspiratory pressure drop was primarily proportional to the square of flow, consistent with pressure losses due to convective acceleration caused by area restriction. On inspiration, in OSAS pressure drop occurred primarily between the choanae and the region where the adenoids overlap the tonsils (overlap region) due to airway narrowing, rather than in the nasal passages; in controls the majority of pressure drop was in the nasal passages. On expiration, in OSAS the majority of pressure drop occurred between the oropharynx (posterior to the tongue) and overlap region, and local minimum pressure in the overlap region was near atmospheric due to pressure recovery in the anterior nasopharynx. The results suggest that pharyngeal airway shape in children with OSAS significantly affects internal pressure distribution compared to nasal resistance. The model may also help explain regional dynamic airway narrowing during expiration. PMID:16098533
Turolla, Andrea; Santoro, Domenico; de Bruyn, John R; Crapulli, Ferdinando; Antonelli, Manuela
2016-01-01
A computational fluid dynamic (CFD) model was developed to describe the process performance of a semi-batch annular TiO2-UV photoreactor in an Eulerian framework. The model accounted for the optical behaviour of titanium dioxide (TiO2) suspensions, the flow distribution and the oxalic acid degradation in the reactor. The scattering component of the optical model, explicitly included in the CFD simulations using a TiO2-specific scattering phase function integrated in the radiative transfer equation, was calibrated using an optical goniometer by comparing simulated scattering light profiles against irradiance measurements collected for various TiO2 concentrations and UV wavelengths and subsequently solved by the discrete ordinate (DO) radiation model. Several scattering phase functions were tested against the goniometric measurements confirming that the Henyey-Greenstein (HG) equation was the most appropriate angular distribution function at 254 and 355 nm, irrespective of the TiO2 concentration. Using the calibrated HG function, a new approach for quantifying the absolute values of absorption and scattering coefficients in TiO2 suspensions was proposed. It consists of iteratively solving, using the DO model, the radiative transfer equation for various combinations of absorption and scattering coefficients until the error between observed and predicted angular irradiance measurements is minimized. The accuracy of the optical parameters was verified with independent CFD simulations carried out for an annular photoreactor and already available in the literature. Predicted and simulated irradiance and oxalic acid degradation data were found to be in excellent agreement, confirming the considerable potential of the integrated modelling approach presented in this paper for the design, optimization and scale-up of photocatalytic technologies for water and wastewater treatment applications. PMID:26476682
Kinetic-fluid dynamics modeling of I{sub 2} dissociation in supersonic chemical oxygen-iodine lasers
Waichman, K.; Barmashenko, B. D.; Rosenwaks, S.
2009-09-15
The mechanism of I{sub 2} dissociation in supersonic chemical oxygen-iodine lasers (COILs) is studied applying kinetic-fluid dynamics modeling, where pathways involving the excited species I{sub 2}(X {sup 1}SIGMA{sub g}{sup +},10<=v<25), I{sub 2}(X {sup 1}SIGMA{sub g}{sup +},25<=v<=47), I{sub 2}(A{sup '} {sup 3}PI{sub 2u}), I{sub 2}(A {sup 3}PI{sub 1u}), O{sub 2}(X {sup 3}SIGMA{sub g}{sup -},v), O{sub 2}(a {sup 1}DELTA{sub g},v), O{sub 2}(b {sup 1}SIGMA{sub g}{sup +},v), and I({sup 2}P{sub 1/2}) as intermediate reactants are included. The gist of the model is adding the first reactant and reducing the contribution of the second as compared to previous models. These changes, recently suggested by Azyazov, et al. [J. Chem. Phys. 130, 104306 (2009)], significantly improve the agreement with the measurements of the gain in a low pressure supersonic COIL for all I{sub 2} flow rates that have been tested in the experiments. In particular, the lack of agreement for high I{sub 2} flow rates, which was encountered in previous models, has been eliminated in the present model. It is suggested that future modeling of the COIL operation should take into account the proposed contribution of the above mentioned reactants.
Λ polarization in an exact fluid dynamical model for heavy-ion collisions
NASA Astrophysics Data System (ADS)
Xie, Yilong
2016-05-01
Λ polarization is calculated in an exact analytical, rotating model based on parameters from a high resolution (3+1)D Particle-in-Cell Relativistic hydrodynamics calculation. The polarization is attributed to effects from thermal vorticity and for the first time the effects of the radial and axial acceleration are also studied separately.
Ilegbusi, Olusegun; Li, Ziang; Min, Yugang; Meeks, Sanford; Kupelian, Patrick; Santhanam, Anand P
2012-01-01
The aim of this paper is to model the airflow inside lungs during breathing and its fluid-structure interaction with the lung tissues and the lung tumor using subject-specific elastic properties. The fluid-structure interaction technique simultaneously simulates flow within the airway and anisotropic deformation of the lung lobes. The three-dimensional (3D) lung geometry is reconstructed from the end-expiration 3D CT scan datasets of humans with lung cancer. The lung is modeled as a poro-elastic medium with anisotropic elastic property (non-linear Young's modulus) obtained from inverse lung elastography of 4D CT scans for the same patients. The predicted results include the 3D anisotropic lung deformation along with the airflow pattern inside the lungs. The effect is also presented of anisotropic elasticity on both the spatio-temporal volumetric lung displacement and the regional lung hysteresis. PMID:22356987
Shuhaiber, Jeffrey H.; Niehaus, Justin; Gottliebson, William; Abdallah, Shaaban
2013-01-01
OBJECTIVES The theoretical differences in energy losses as well as coronary flow with different band sizes for branch pulmonary arteries (PA) in hypoplastic left heart syndrome (HLHS) remain unknown. Our objective was to develop a computational fluid dynamic model (CFD) to determine the energy losses and pulmonary-to-systemic flow rates. This study was done for three different PA band sizes. METHODS Three-dimensional computer models of the hybrid procedure were constructed using the standard commercial CFD softwares Fluent and Gambit. The computer models were controlled for bilateral PA reduction to 25% (restrictive), 50% (intermediate) and 75% (loose) of the native branch pulmonary artery diameter. Velocity and pressure data were calculated throughout the heart geometry using the finite volume numerical method. Coronary flow was measured simultaneously with each model. Wall shear stress and the ratio of pulmonary-to-systemic volume flow rates were calculated. Computer simulations were compared at fixed points utilizing echocardiographic and catheter-based metric dimensions. RESULTS Restricting the PA band to a 25% diameter demonstrated the greatest energy loss. The 25% banding model produced an energy loss of 16.76% systolic and 24.91% diastolic vs loose banding at 7.36% systolic and 17.90% diastolic. Also, restrictive PA bands had greater coronary flow compared with loose PA bands (50.2 vs 41.9 ml/min). Shear stress ranged from 3.75 Pascals with restrictive PA banding to 2.84 Pascals with loose banding. Intermediate PA banding at 50% diameter achieved a Qp/Qs (closest to 1) at 1.46 systolic and 0.66 diastolic compared with loose or restrictive banding without excess energy loss. CONCLUSIONS CFD provides a unique platform to simulate pressure, shear stress as well as energy losses of the hybrid procedure. PA banding at 50% provided a balanced pulmonary and systemic circulation with adequate coronary flow but without extra energy losses incurred. PMID:23660734
NASA Astrophysics Data System (ADS)
Peladeau-Pigeon, M.; Coolens, C.
2013-09-01
Dynamic contrast-enhanced computed tomography (DCE-CT) is an imaging tool that aids in evaluating functional characteristics of tissue at different stages of disease management: diagnostic, radiation treatment planning, treatment effectiveness, and monitoring. Clinical validation of DCE-derived perfusion parameters remains an outstanding problem to address prior to perfusion imaging becoming a widespread standard as a non-invasive quantitative measurement tool. One approach to this validation process has been the development of quality assurance phantoms in order to facilitate controlled perfusion ex vivo. However, most of these systems fail to establish and accurately replicate physiologically relevant capillary permeability and exchange performance. The current work presents the first step in the development of a prospective suite of physics-based perfusion simulations based on coupled fluid flow and particle transport phenomena with the goal of enhancing the understanding of clinical contrast agent kinetics. Existing knowledge about a controllable, two-compartmental fluid exchange phantom was used to validate the computational fluid dynamics (CFD) simulation model presented herein. The sensitivity of CFD-derived contrast uptake curves to contrast injection parameters, including injection duration and flow rate, were quantified and found to be within 10% accuracy. The CFD model was employed to evaluate two commonly used clinical kinetic algorithms used to derive perfusion parameters: Fick's principle and the modified Tofts model. Neither kinetic model was able to capture the true transport phenomena it aimed to represent but if the overall contrast concentration after injection remained identical, then successive DCE-CT evaluations could be compared and could indeed reflect differences in regional tissue flow. This study sets the groundwork for future explorations in phantom development and pharmaco-kinetic modelling, as well as the development of novel contrast
Coupled Heat Transfer and Fluid Dynamics Modeling of InSb Solidification
NASA Astrophysics Data System (ADS)
Barvinschi, Paul; Barvinschi, Floricica
2011-10-01
A method for the directional solidification of melted InSb in a silica ampoule is presented and solved with COMSOL Multiphysics. The configuration and initial boundary settings of the model resemble those used in a de-wetting vertical Bridgman configuration [1]. A slightly modified version of the method presented by Voller and Prakash [2] is used to account for solidification of the liquid phase, including convection and conduction heat transfer with mushy region phase change. Axial-symmetric numerical simulations of temperature and velocity fields, under normal gravity, are carried out using different thermal conditions.
Fluid dynamic modeling and numerical simulation of low-density hypersonic flow
NASA Technical Reports Server (NTRS)
Cheng, H. K.; Wong, Eric Y.
1988-01-01
The concept of a viscous shock-layer and several related versions of continuum theories/methods are examined for their adequacy as a viable framework to study flow physics and aerothermodynamics of relevance to sustained hypersonic flights. Considering the flat plate at angle of attack, or the wedge, as a generic example for the major aerodynamic component of a hypersonic vehicle, the relative importance of the molecular-transport effects behind the shock (in the form of the 'shock slip') and the wall-slip effects are studied. In the flow regime where the shock-transition-zone thickness remains small compared to the shock radius of curvature, a quasi-one-dimensional shock structure under the Burnett/thirteen-moment approximation, as well as particulate/collisional models, can be consistently developed. The fully viscous version of the shock-layer model is shown to provide the crucial boundary condition downstream the shock in this case. The gas-kinetic basis of the continuum description for the flow behind the bow shock, and certain features affecting the non-equilibrium flow chemistry, are also discussed.
NASA Astrophysics Data System (ADS)
Gu, S.; Eastwick, C. N.; Simmons, K. A.; McCartney, D. G.
2001-09-01
A computational fluid dynamics (CFD) model is developed to predict gas dynamic behavior in a high-velocity oxy-fuel (HVOF) thermal spray gun in which premixed oxygen and propylene are burnt in a 12 mm combustion chamber linked to a parallel-sided nozzle. The CFD analysis is applied to investigate axisymmetric, steady-state, turbulent, compressible, and chemically combusting flow both within the gun and in a free jet region between the gun and the substrate to be coated. The combustion of oxygen and propylene is modeled using a single-step, finite-rate chemistry model that also allows for dissociation of the reaction products. Results are presented to show the effect of (1) fuel-to-oxygen gas ratio and (2) total gas flow rate on the gas dynamic behavior. Along the centerline, the maximum temperature reached is insensitive to the gas ratio but depends on the total flow. However, the value attained (˜2500 K) is significantly lower than the maximum temperature (˜3200 K) of the annular flame in the combustion chamber. By contrast, the centerline gas velocity depends on both total flow and gas ratio, the highest axial gas velocity being attained with the higher flow and most fuel-rich mixture. The gas Mach number increases through the gun and reaches a maximum value of approximately 1.6 around 5 mm downstream from the nozzle exit. The numerical calculations also show that the residual oxygen level is principally dependent on the fuel-to-oxygen ratio and decreases by approximately fivefold as the ratio is varied from 90 to 69% of the stoichiometric requirement. The CFD model is also used to investigate the effect of changes in combustion chamber size and geometry on gas dynamics, and the results are compared with the nominal 12 mm chamber baseline calculations.
Cheng, Cheng; Zhang, Xiaobing
2013-05-01
In conventional models for two-phase reactive flow of interior ballistic, the dynamic collision phenomenon of particles is neglected or empirically simplified. However, the particle collision between particles may play an important role in dilute two-phase flow because the distribution of particles is extremely nonuniform. The collision force may be one of the key factors to influence the particle movement. This paper presents the CFD-DEM approach for simulation of interior ballistic two-phase flow considering the dynamic collision process. The gas phase is treated as a Eulerian continuum and described by a computational fluid dynamic method (CFD). The solid phase is modeled by discrete element method (DEM) using a soft sphere approach for the particle collision dynamic. The model takes into account grain combustion, particle-particle collisions, particle-wall collisions, interphase drag and heat transfer between gas and solid phases. The continuous gas phase equations are discretized in finite volume form and solved by the AUSM+-up scheme with the higher order accurate reconstruction method. Translational and rotational motions of discrete particles are solved by explicit time integrations. The direct mapping contact detection algorithm is used. The multigrid method is applied in the void fraction calculation, the contact detection procedure, and CFD solving procedure. Several verification tests demonstrate the accuracy and reliability of this approach. The simulation of an experimental igniter device in open air shows good agreement between the model and experimental measurements. This paper has implications for improving the ability to capture the complex physics phenomena of two-phase flow during the interior ballistic cycle and to predict dynamic collision phenomena at the individual particle scale. PMID:24891728
Multiplicity fluctuations at the quark-hadron phase transition from a fluid dynamical model
NASA Astrophysics Data System (ADS)
Herold, Christoph; Nahrgang, Marlene; Yan, Yupeng; Kobdaj, Chinorat
2015-04-01
The region of large net-baryon densities in the QCD phase diagram is expected to exhibit a first-order phase transition. Experimentally, its study will be one of the primary objectives for the upcoming FAIR accelerator. We model the transition between quarks and hadrons in a heavy-ion collision using a fluid which is coupled to the explicit dynamics of the chiral order parameter and a dilaton field. This allows us to investigate signals stemming from the nonequilibrium evolution during the expansion of the hot plasma. Special emphasis is put on an event-by-event analysis of baryon number fluctuations which have long since been claimed to be sensitive to a critical point.
Modeling-gas phase reactions in indoor environments using computational fluid dynamics
NASA Astrophysics Data System (ADS)
Sørensen, Dan Nørtoft; Weschler, Charles J.
This CFD modeling study examines the concentrations of two gaseous compounds that react in an indoor setting to produce a hypothetical product. The reactants are ozone and either d-limonene or α-terpinene (which reacts with ozone about 40 times faster than d-limonene). In addition to two different terpenes, the scenarios include two air exchange rates (0.5 and 2.0 h-1). The terpene is introduced as a floor source with an emission pattern similar to a floor-care product. These four scenarios have been set in a fairly large two-dimensional room (13.6×40.6 m) with a supply at the top of the left wall and an exhaust at the bottom of the right wall. The room has been deliberately scaled so that the Reynolds numbers for key flow regimes match those of a room in which the calculated flow field has been validated against measured data. It has been further assumed that ozone interacts with room surfaces while the terpenes do not. The results show that for all four scenarios, under steady-state conditions, there are large concentration gradients within the room for both reactants and product. To some extent this is due to imperfect mixing. However, it also reflects that reactions occur at different rates across the room (because of varying reactant concentrations) and that the time available for reactions to occur varies with the room location (because the "age of the air" varies from point to point). Locally, within the room, the concentrations calculated by the CFD method differ significantly from those calculated by a one-compartment mass-balance model assuming perfect mixing.
Fluid Dynamics of the Heart and its Valves
NASA Astrophysics Data System (ADS)
Peskin, Charles S.
1997-11-01
The fluid dynamics of the heart involve the interaction of blood, a viscous incompressible fluid, with the flexible, elastic, fiber-reinforced heart valve leaflets that are immersed in that fluid. Neither the fluid motion nor the valve leaflet motion are known in advance: both must be computed simultaneously by solving their coupled equations of motion. This can be done by the immersed boundary method(Peskin CS and McQueen DM: A general method for the computer simulation of biological systems interacting with fluids. In: Biological Fluid Dynamics (Ellington CP and Pedley TJ, eds.), The Company of Biologists Limited, Cambridge UK, 1995, pp. 265-276.), which can be extended to incorporate the contractile fiber architecture of the muscular heart walls as well as the valve leaflets and the blood. In this way we arrive at a three-dimensional computer model of the heart(Peskin CS and McQueen DM: Fluid dynamics of the heart and its valves. In: Case Studies in Mathematical Modeling: Ecology, Physiology, and Cell Biology (Othmer HG, Adler FR, Lewis MA, and Dallon JC, eds.), Prentice-Hall, Englewood Cliffs NJ, 1996, pp. 309-337.), which can be used as a test chamber for the design of prosthetic cardiac valves, and also to study the function of the heart in health and in disease. Numerical solutions of the equations of cardiac fluid dynamics obtained by the immersed boundary method will be presented in the form of a video animation of the beating heart.
Implementing fluid dynamics obtained from GeoPET in reactive transport models
NASA Astrophysics Data System (ADS)
Lippmann-Pipke, Johanna; Eichelbaum, Sebastian; Kulenkampff, Johannes
2016-04-01
Flow and transport simulations in geomaterials are commonly conducted on high-resolution tomograms (μCT) of the pore structure or stochastic models that are calibrated with measured integral quantities, like break through curves (BTC). Yet, there existed virtually no method for experimental verification of the simulated velocity distribution results. Positron emission tomography (PET) has unrivaled sensitivity and robustness for non-destructive, quantitative, spatio-temporal measurement of tracer concentrations in body tissue. In the past decade, we empowered PET for its applicability in opaque/geological media - GeoPET (Kulenkampff et al.; Kulenkampff et al., 2008; Zakhnini et al., 2013) and have developed detailed correction schemes to bring the images into sharp focus. Thereby it is the appropriate method for experimental verification and calibration of computer simulations of pore-scale transport by means of the observed propagation of a tracer pulse, c_PET(x,y,z,t). In parallel, we aimed at deriving velocity and porosity distributions directly from our concentration time series of fluid flow processes in geomaterials. This would allow us to directly benefit from lab scale observations and to parameterize respective numerical transport models. For this we have developed a robust spatiotemporal (3D+t) parameter extraction algorithm. Here, we will present its functionality, and demonstrate the use of obtained velocity distributions in finite element simulations of reactive transport processes on drill core scale. Kulenkampff, J., Gruendig, M., Zakhnini, A., Gerasch, R., and Lippmann-Pipke, J.: Process tomography of diffusion with PET for evaluating anisotropy and heterogeneity, Clay Minerals, in press. Kulenkampff, J., Gründig, M., Richter, M., and Enzmann, F.: Evaluation of positron emission tomography for visualisation of migration processes in geomaterials, Physics and Chemistry of the Earth, 33, 937-942, 2008. Zakhnini, A., Kulenkampff, J., Sauerzapf, S
Tentner, A.; Froehle, P.; Wang, C.; Nuclear Engineering Division
2004-01-01
This work has explored the preliminary design of a Computational Fluid Dynamics (CFD) tool for the analysis of transient vehicle underhood thermo-hydrodynamic events using high performance computing platforms. The goal of this tool will be to extend the capabilities of an existing established CFD code, STAR-CD, allowing the car manufacturers to analyze the impact of transient operational events on the underhood thermal management by exploiting the computational efficiency of modern high performance computing systems. In particular, the project has focused on the CFD modeling of the radiator behavior during a specified transient. The 3-D radiator calculations were performed using STAR-CD, which can perform both steady-state and transient calculations, on the cluster computer available at ANL in the Nuclear Engineering Division. Specified transient boundary conditions, based on experimental data provided by Adapco and DaimlerChrysler were used. The possibility of using STAR-CD in a transient mode for the entire period of time analyzed has been compared with other strategies which involve the use of STAR-CD in a steady-state mode at specified time intervals, while transient heat transfer calculations would be performed for the rest of the time. The results of these calculations have been compared with the experimental data provided by Adapco/DaimlerChrysler and recommendations for future development of an optimal strategy for the CFD modeling of transient thermo-hydrodynamic events have been made. The results of this work open the way for the development of a CFD tool for the transient analysis of underhood thermo-hydrodynamic events, which will allow the integrated transient thermal analysis of the entire cooling system, including both the engine block and the radiator, on high performance computing systems.
Froehle, P.; Tentner, A.; Wang, C.
2003-09-05
This work has explored the preliminary design of a Computational Fluid Dynamics (CFD) tool for the analysis of transient vehicle underhood thermo-hydrodynamic events using high performance computing platforms. The goal of this tool will be to extend the capabilities of an existing established CFD code, STAR-CD, allowing the car manufacturers to analyze the impact of transient operational events on the underhood thermal management by exploiting the computational efficiency of modern high performance computing systems. In particular, the project has focused on the CFD modeling of the radiator behavior during a specified transient. The 3-D radiator calculations were performed using STAR-CD, which can perform both steady-state and transient calculations, on the cluster computer available at ANL in the Nuclear Engineering Division. Specified transient boundary conditions, based on experimental data provided by Adapco and DaimlerChrysler were used. The possibility of using STAR-CD in a transient mode for the entire period of time analyzed has been compared with other strategies which involve the use of STAR-CD in a steady-state mode at specified time intervals, while transient heat transfer calculations would be performed for the rest of the time. The results of these calculations have been compared with the experimental data provided by Adapco/DaimlerChrysler and recommendations for future development of an optimal strategy for the CFD modeling of transient thermo-hydrodynamic events have been made. The results of this work open the way for the development of a CFD tool for the transient analysis of underhood thermo-hydrodynamic events, which will allow the integrated transient thermal analysis of the entire cooling system, including both the engine block and the radiator, on high performance computing systems.
NASA Astrophysics Data System (ADS)
Boutchko, Rostyslav; Rayz, Vitaliy L.; Vandehey, Nicholas T.; O'Neil, James P.; Budinger, Thomas F.; Nico, Peter S.; Druhan, Jennifer L.; Saloner, David A.; Gullberg, Grant T.; Moses, William W.
2012-01-01
This paper presents experimental and modeling aspects of applying nuclear emission tomography to study fluid flow in laboratory packed porous media columns of the type frequently used in geophysics, geochemistry and hydrology research. Positron emission tomography (PET) and single photon emission computed tomography (SPECT) are used as non-invasive tools to obtain dynamic 3D images of radioactive tracer concentrations. Dynamic sequences obtained using 18F-FDG PET are used to trace flow through a 5 cm diameter × 20 cm tall sand packed column with and without an impermeable obstacle. In addition, a custom-made rotating column setup placed in a clinical two-headed SPECT camera is used to image 99mTc-DTPA tracer propagation in a through-flowing column (10 cm diameter × 30 cm tall) packed with recovered aquifer sediments. A computational fluid dynamics software package FLUENT is used to model the observed flow dynamics. Tracer distributions obtained in the simulations in the smaller column uniformly packed with sand and in the column with an obstacle are remarkably similar to the reconstructed images in the PET experiments. SPECT results demonstrate strongly non-uniform flow patterns for the larger column slurry-packed with sub-surface sediment and slow upward flow. In the numerical simulation of the SPECT study, two symmetric channels with increased permeability are prescribed along the column walls, which result in the emergence of two well-defined preferential flow paths. Methods and results of this work provide new opportunities in hydrologic and biogeochemical research. The primary target application for developed technologies is non-destructive, non-perturbing, quantitative imaging of flow dynamics within laboratory scale porous media systems.
NASA Astrophysics Data System (ADS)
Chung, T. J.
2002-03-01
Computational fluid dynamics (CFD) techniques are used to study and solve complex fluid flow and heat transfer problems. This comprehensive text ranges from elementary concepts for the beginner to state-of-the-art CFD for the practitioner. It discusses and illustrates the basic principles of finite difference (FD), finite element (FE), and finite volume (FV) methods, with step-by-step hand calculations. Chapters go on to examine structured and unstructured grids, adaptive methods, computing techniques, and parallel processing. Finally, the author describes a variety of practical applications to problems in turbulence, reacting flows and combustion, acoustics, combined mode radiative heat transfer, multiphase flows, electromagnetic fields, and relativistic astrophysical flows. Students and practitioners--particularly in mechanical, aerospace, chemical, and civil engineering--will use this authoritative text to learn about and apply numerical techniques to the solution of fluid dynamics problems.
Park, Stephen; Li, Yebo
2015-05-01
Microalgal growth and systemic productivity is not only affected by environmental conditions such as temperature, irradiance, and nutrient concentrations, but also by physical processes such as fluid flow and particulate sedimentation. Modeling and simulating the system is a cost-effective way to predict the growth behavior under various environmental and physical conditions while determining effective engineering approaches to maximize productivity. Many mathematical models have been proposed to describe microalgal growth, while computational fluid dynamics (CFD) have been used to model the behavior of many fluid systems. Integrating the growth kinetics into a CFD model can help researchers understand the impact of a variety of parameters and determine what measures can be taken to overcome some obstacles in the aquaculture industry--self-shading, biomass sedimentation, and contamination--which prevent the production of high biomass yields. The aim of this study was to integrate physical and environmental effects to predict space- and time-dependent algal growth in industrial scale raceways. A commercial CFD software, ANSYS-Fluent 14.5, was used to solve the proposed models in regards to fluid flow, heat transfer, and nutrient balance. User-defined functions written in C language were used to incorporate the kinetic equations into a three-dimensional standard k-ε turbulence model of an open channel raceway system driven by a single paddlewheel. Simulated results were compared with light intensity, temperature, nutrient concentration, and algal biomass data acquired for 56 day from an industrial scale raceway pond constructed for the growth of Nannochloropsis salina and were observed to be in good agreement with one another. There was up to a 17.6% increase in simulated productivity when the incoming CO2 concentration was increased from 0.0006 to 0.150 g L(-1), while the effect of paddlewheel velocity was not significant. Sensitivity analysis showed that the model
NASA Technical Reports Server (NTRS)
Strutzenberg, L. L.; Dougherty, N. S.; Liever, P. A.; West, J. S.; Smith, S. D.
2007-01-01
This paper details advances being made in the development of Reynolds-Averaged Navier-Stokes numerical simulation tools, models, and methods for the integrated Space Shuttle Vehicle at launch. The conceptual model and modeling approach described includes the development of multiple computational models to appropriately analyze the potential debris transport for critical debris sources at Lift-Off. The conceptual model described herein involves the integration of propulsion analysis for the nozzle/plume flow with the overall 3D vehicle flowfield at Lift-Off. Debris Transport Analyses are being performed using the Shuttle Lift-Off models to assess the risk to the vehicle from Lift-Off debris and appropriately prioritized mitigation of potential debris sources to continue to reduce vehicle risk. These integrated simulations are being used to evaluate plume-induced debris environments where the multi-plume interactions with the launch facility can potentially accelerate debris particles toward the vehicle.
The use of computers for instruction in fluid dynamics
NASA Technical Reports Server (NTRS)
Watson, Val
1987-01-01
Applications for computers which improve instruction in fluid dynamics are examined. Computers can be used to illustrate three-dimensional flow fields and simple fluid dynamics mechanisms, to solve fluid dynamics problems, and for electronic sketching. The usefulness of computer applications is limited by computer speed, memory, and software and the clarity and field of view of the projected display. Proposed advances in personal computers which will address these limitations are discussed. Long range applications for computers in education are considered.
NASA Technical Reports Server (NTRS)
Beutner, Thomas John
1993-01-01
Porous wall wind tunnels have been used for several decades and have proven effective in reducing wall interference effects in both low speed and transonic testing. They allow for testing through Mach 1, reduce blockage effects and reduce shock wave reflections in the test section. Their usefulness in developing computational fluid dynamics (CFD) codes has been limited, however, by the difficulties associated with modelling the effect of a porous wall in CFD codes. Previous approaches to modelling porous wall effects have depended either upon a simplified linear boundary condition, which has proven inadequate, or upon detailed measurements of the normal velocity near the wall, which require extensive wind tunnel time. The current work was initiated in an effort to find a simple, accurate method of modelling a porous wall boundary condition in CFD codes. The development of such a method would allow data from porous wall wind tunnels to be used more readily in validating CFD codes. This would be beneficial when transonic validations are desired, or when large models are used to achieve high Reynolds numbers in testing. A computational and experimental study was undertaken to investigate a new method of modelling solid and porous wall boundary conditions in CFD codes. The method utilized experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution was then imposed as a pressure boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the effect of porosity changes on the wall was investigated. Also, the effectiveness of this method when only sparse experimental measurements were available has been investigated. The current work demonstrated this approach for low speed flows and compared the results with experimental data obtained from a heavily instrumented variable porosity test section. The approach developed was simple, computationally
NASA Astrophysics Data System (ADS)
Hazer, D.; Schmidt, E.; Unterhinninghofen, R.; Richter, G. M.; Dillmann, R.
2009-08-01
Abnormal hemodynamics and biomechanics of blood flow and vessel wall conditions in the arteries may result in severe cardiovascular diseases. Cardiovascular diseases result from complex flow pattern and fatigue of the vessel wall and are prevalent causes leading to high mortality each year. Computational Fluid Dynamics (CFD), Computational Structure Mechanics (CSM) and Fluid Structure Interaction (FSI) have become efficient tools in modeling the individual hemodynamics and biomechanics as well as their interaction in the human arteries. The computations allow non-invasively simulating patient-specific physical parameters of the blood flow and the vessel wall needed for an efficient minimally invasive treatment. The numerical simulations are based on the Finite Element Method (FEM) and require exact and individual mesh models to be provided. In the present study, we developed a numerical tool to automatically generate complex patient-specific Finite Element (FE) mesh models from image-based geometries of healthy and diseased vessels. The mesh generation is optimized based on the integration of mesh control functions for curvature, boundary layers and mesh distribution inside the computational domain. The needed mesh parameters are acquired from a computational grid analysis which ensures mesh-independent and stable simulations. Further, the generated models include appropriate FE sets necessary for the definition of individual boundary conditions, required to solve the system of nonlinear partial differential equations governed by the fluid and solid domains. Based on the results, we have performed computational blood flow and vessel wall simulations in patient-specific aortic models providing a physical insight into the pathological vessel parameters. Automatic mesh generation with individual awareness in terms of geometry and conditions is a prerequisite for performing fast, accurate and realistic FEM-based computations of hemodynamics and biomechanics in the
Lectures on Geophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Samelson, Roger M.
The fluid kaleidoscope of the Earth's ocean and atmosphere churns and sparkles with jets, gyres, eddies, waves, streams, and cyclones. These vast circulations, essential elements of the physical environment that support human life, are given a special character by the Earth's rotation and by their confinement to a shallow surficial layer, thin relative to the solid Earth in roughly the same proportion as an apple skin is to an apple. Geophysical fluid dynamics exploits this special character to develop a unified theoretical approach to the physics of the ocean and atmosphere.With Lectures on Geophysical Fluid Dynamics, Rick Salmon has added an insightful and provocative volume to the handful of authoritative texts currently available on the subject. The book is intended for first-year graduate students, but advanced students and researchers also will find it useful. It is divided into seven chapters, the first four of these adapted from course lectures. The book is well written and presents a fresh and stimulating perspective that complements existing texts. It would serve equally well either as the main text for a core graduate curriculum or as a supplementary resource for students and teachers seeking new approaches to both classical and contemporary problems. A lively set of footnotes contains many references to very recent work. The printing is attractive, the binding is of high quality, and typographical errors are few.
Body fluid dynamics: back to the future.
Bhave, Gautam; Neilson, Eric G
2011-12-01
Pioneering investigations conducted over a half century ago on tonicity, transcapillary fluid exchange, and the distribution of water and solute serve as a foundation for understanding the physiology of body fluid spaces. With passage of time, however, some of these concepts have lost their connectivity to more contemporary information. Here we examine the physical forces determining the compartmentalization of body fluid and its movement across capillary and cell membrane barriers, drawing particular attention to the interstitium operating as a dynamic interface for water and solute distribution rather than as a static reservoir. Newer work now supports an evolving model of body fluid dynamics that integrates exchangeable Na(+) stores and transcapillary dynamics with advances in interstitial matrix biology. PMID:22034644
Body Fluid Dynamics: Back to the Future
Bhave, Gautam; Neilson, Eric G.
2014-01-01
Pioneering investigations conducted over a half century ago on tonicity, transcapillary fluid exchange, and the distribution of water and solute serve as a foundation for understanding the physiology of body fluid spaces. With passage of time, however, some of these concepts have lost their connectivity to more contemporary information. Here we examine the physical forces determining the compartmentalization of body fluid and its movement across capillary and cell membrane barriers, drawing particular attention to the interstitium operating as a dynamic interface for water and solute distribution rather than as a static reservoir. Newer work now supports an evolving model of body fluid dynamics that integrates exchangeable Na+ stores and transcapillary dynamics with advances in interstitial matrix biology. PMID:22034644
Rakowski, Cynthia L.; Serkowski, John A.; Richmond, Marshall C.; Perkins, William A.
2010-12-01
In 2003, an extension of the existing ice and trash sluiceway was added at Bonneville Powerhouse 2 (B2). This extension started at the existing corner collector for the ice and trash sluiceway adjacent to Bonneville Powerhouse 2 and the new sluiceway was extended to the downstream end of Cascade Island. The sluiceway was designed to improve juvenile salmon survival by bypassing turbine passage at B2, and placing these smolt in downstream flowing water minimizing their exposure to fish and avian predators. In this study, a previously developed computational fluid dynamics model was modified and used to characterized tailrace hydraulics and sluiceway egress conditions for low total river flows and low levels of spillway flow. STAR-CD v4.10 was used for seven scenarios of low total river flow and low spill discharges. The simulation results were specifically examined to look at tailrace hydraulics at 5 ft below the tailwater elevation, and streamlines used to compare streamline pathways for streamlines originating in the corner collector outfall and adjacent to the outfall. These streamlines indicated that for all higher spill percentage cases (25% and greater) that streamlines from the corner collector did not approach the shoreline at the downstream end of Bradford Island. For the cases with much larger spill percentages, the streamlines from the corner collector were mid-channel or closer to the Washington shore as they moved downstream. Although at 25% spill at 75 kcfs total river, the total spill volume was sufficient to "cushion" the flow from the corner collector from the Bradford Island shore, areas of recirculation were modeled in the spillway tailrace. However, at the lowest flows and spill percentages, the streamlines from the B2 corner collector pass very close to the Bradford Island shore. In addition, the very flow velocity flows and large areas of recirculation greatly increase potential predator exposure of the spillway passed smolt. If there is
Nichols, B.D.; Mueller, C.; Necker, G.A.; Travis, J.R.; Spore, J.W.; Lam, K.L.; Royl, P.; Redlinger, R.; Wilson, T.L.
1998-10-01
Los Alamos National Laboratory (LANL) and Forschungszentrum Karlsruhe (FzK) are developing GASFLOW, a three-dimensional (3D) fluid dynamics field code as a best-estimate tool to characterize local phenomena within a flow field. Examples of 3D phenomena include circulation patterns; flow stratification; hydrogen distribution mixing and stratification; combustion and flame propagation; effects of noncondensable gas distribution on local condensation and evaporation; and aerosol entrainment, transport, and deposition. An analysis with GASFLOW will result in a prediction of the gas composition and discrete particle distribution in space and time throughout the facility and the resulting pressure and temperature loadings on the walls and internal structures with or without combustion. A major application of GASFLOW is for predicting the transport, mixing, and combustion of hydrogen and other gases in nuclear reactor containments and other facilities. It has been applied to situations involving transporting and distributing combustible gas mixtures. It has been used to study gas dynamic behavior (1) in low-speed, buoyancy-driven flows, as well as sonic flows or diffusion dominated flows; and (2) during chemically reacting flows, including deflagrations. The effects of controlling such mixtures by safety systems can be analyzed. The code version described in this manual is designated GASFLOW 2.1, which combines previous versions of the United States Nuclear Regulatory Commission code HMS (for Hydrogen Mixing Studies) and the Department of Energy and FzK versions of GASFLOW. The code was written in standard Fortran 90. This manual comprises three volumes. Volume I describes the governing physical equations and computational model. Volume II describes how to use the code to set up a model geometry, specify gas species and material properties, define initial and boundary conditions, and specify different outputs, especially graphical displays. Sample problems are included
Leishear, Robert A.; Lee, Si Y.; Poirier, Michael R.; Steeper, Timothy J.; Ervin, Robert C.; Giddings, Billy J.; Stefanko, David B.; Harp, Keith D.; Fowley, Mark D.; Van Pelt, William B.
2012-10-07
Computational fluid dynamics (CFD) is recognized as a powerful engineering tool. That is, CFD has advanced over the years to the point where it can now give us deep insight into the analysis of very complex processes. There is a danger, though, that an engineer can place too much confidence in a simulation. If a user is not careful, it is easy to believe that if you plug in the numbers, the answer comes out, and you are done. This assumption can lead to significant errors. As we discovered in the course of a study on behalf of the Department of Energy's Savannah River Site in South Carolina, CFD models fail to capture some of the large variations inherent in complex processes. These variations, or scatter, in experimental data emerge from physical tests and are inadequately captured or expressed by calculated mean values for a process. This anomaly between experiment and theory can lead to serious errors in engineering analysis and design unless a correction factor, or safety factor, is experimentally validated. For this study, blending times for the mixing of salt solutions in large storage tanks were the process of concern under investigation. This study focused on the blending processes needed to mix salt solutions to ensure homogeneity within waste tanks, where homogeneity is required to control radioactivity levels during subsequent processing. Two of the requirements for this task were to determine the minimum number of submerged, centrifugal pumps required to blend the salt mixtures in a full-scale tank in half a day or less, and to recommend reasonable blending times to achieve nearly homogeneous salt mixtures. A full-scale, low-flow pump with a total discharge flow rate of 500 to 800 gpm was recommended with two opposing 2.27-inch diameter nozzles. To make this recommendation, both experimental and CFD modeling were performed. Lab researchers found that, although CFD provided good estimates of an average blending time, experimental blending times varied
NASA Astrophysics Data System (ADS)
Trout, Joseph; Manson, J. Russell; Rios, Manny; King, David; Decicco, Nicholas
2015-04-01
Wake Vortex Turbulence is the turbulence generated by an aircraft in flight. This turbulence is created by vortices at the tips of the wing that may decay slowly and persist for several minutes after creation. The strength, formation and lifetime of the turbulence and vortices are effected by many things including the weather. Here we present the preliminary results of an investigation of low level wind fields generated by the Weather Research and Forecasting Model and an analysis of historical data. The simulations are used as inputs for the computational fluid dynamics model (OpenFoam) that will be used to investigate the effect of weather on wake turbulence. The initial results of the OpenFoam model are presented elsewhere. Presented here are the initial results from a research grant, ``A Pilot Project to Investigate Wake Vortex Patterns and Weather Patterns at the Atlantic City Airport by the Richard Stockton College of NJ and the FAA''.
Grigioni, M; Daniele, C; Morbiducci, U; D'Avenio, G; Di Benedetto, G; Del Gaudio, C; Barbaro, V
2002-12-01
Vascular access methods, performed by the insertion of cannulae into vessels, may disturb the physiological flow of blood, giving rise to non-physiological pressure variations and shear stresses. To date, the hydrodynamic behaviour of the cannulae has been evaluated comparing their pressure loss-flow rate relationships, as obtained from in vitro experiments using a monodimensional approach; this methodology neither furnish information about the local fluid dynamics nor the established flow field in specific clinical work conditions. Since the shear stress is a critical factor in the design of artificial circulatory devices, more knowledge should be necessary about the local values assumed by the haemodynamic parameters during cannulation. An alternative way to investigate the fluid dynamic as accurately as possible is given by numeric studies. A 3D model of cannula concentrically placed in a rigid wall vessel is presented, with the finite element methodology used to numerically simulate the steady-state flow field in two different venous cannulation case studies, with two cannulae having a central hole and two or four side holes, respectively, with the same boundary conditions. Lower velocity and shear stress peak values have been computed for the model with four side holes upstream of the central hole, in the region of the cannula where the inlet flows meet and towards cannula's outlet, due to the increased flow symmetry and inlet area with respect to the model with two side holes. Starting from the investigation of different cannula designs, numerically assessing the local fluid dynamics, indications can be drawn to support both the design phase and the device optimal clinical use, in order to limit risks of biomechanical origin. Thus the presence of four side holes implied, as a consequence of the greater inlet area and of the increased symmetry, a less disturbed blood flow, together with reduced shear stress values. Furthermore, results show that the numerical
Computational fluid dynamics in oil burner design
Butcher, T.A.
1997-09-01
In Computational Fluid Dynamics, the differential equations which describe flow, heat transfer, and mass transfer are approximately solved using a very laborious numerical procedure. Flows of practical interest to burner designs are always turbulent, adding to the complexity of requiring a turbulence model. This paper presents a model for burner design.
Lawson, M. J.; Li, Y.; Sale, D. C.
2011-10-01
This paper describes the development of a computational fluid dynamics (CFD) methodology to simulate the hydrodynamics of horizontal-axis tidal current turbines. Qualitative measures of the CFD solutions were independent of the grid resolution. Conversely, quantitative comparisons of the results indicated that the use of coarse computational grids results in an under prediction of the hydrodynamic forces on the turbine blade in comparison to the forces predicted using more resolved grids. For the turbine operating conditions considered in this study, the effect of the computational timestep on the CFD solution was found to be minimal, and the results from steady and transient simulations were in good agreement. Additionally, the CFD results were compared to corresponding blade element momentum method calculations and reasonable agreement was shown. Nevertheless, we expect that for other turbine operating conditions, where the flow over the blade is separated, transient simulations will be required.
Mota, J.P.B.; Esteves, I.A.A.C.; Rostam-Abadi, M.
2004-01-01
A computational fluid dynamics (CFD) software package has been coupled with the dynamic process simulator of an adsorption storage tank for methane fuelled vehicles. The two solvers run as independent processes and handle non-overlapping portions of the computational domain. The codes exchange data on the boundary interface of the two domains to ensure continuity of the solution and of its gradient. A software interface was developed to dynamically suspend and activate each process as necessary, and be responsible for data exchange and process synchronization. This hybrid computational tool has been successfully employed to accurately simulate the discharge of a new tank design and evaluate its performance. The case study presented here shows that CFD and process simulation are highly complementary computational tools, and that there are clear benefits to be gained from a close integration of the two. ?? 2004 Elsevier Ltd. All rights reserved.
Computational Fluid Dynamics Library
Energy Science and Technology Software Center (ESTSC)
2005-03-04
CFDLib05 is the Los Alamos Computational Fluid Dynamics LIBrary. This is a collection of hydrocodes using a common data structure and a common numerical method, for problems ranging from single-field, incompressible flow, to multi-species, multi-field, compressible flow. The data structure is multi-block, with a so-called structured grid in each block. The numerical method is a Finite-Volume scheme employing a state vector that is fully cell-centered. This means that the integral form of the conservation lawsmore » is solved on the physical domain that is represented by a mesh of control volumes. The typical control volume is an arbitrary quadrilateral in 2D and an arbitrary hexahedron in 3D. The Finite-Volume scheme is for time-unsteady flow and remains well coupled by means of time and space centered fluxes; if a steady state solution is required, the problem is integrated forward in time until the user is satisfied that the state is stationary.« less
Application of computational fluid dynamics methods to improve thermal hydraulic code analysis
NASA Astrophysics Data System (ADS)
Sentell, Dennis Shannon, Jr.
A computational fluid dynamics code is used to model the primary natural circulation loop of a proposed small modular reactor for comparison to experimental data and best-estimate thermal-hydraulic code results. Recent advances in computational fluid dynamics code modeling capabilities make them attractive alternatives to the current conservative approach of coupled best-estimate thermal hydraulic codes and uncertainty evaluations. The results from a computational fluid dynamics analysis are benchmarked against the experimental test results of a 1:3 length, 1:254 volume, full pressure and full temperature scale small modular reactor during steady-state power operations and during a depressurization transient. A comparative evaluation of the experimental data, the thermal hydraulic code results and the computational fluid dynamics code results provides an opportunity to validate the best-estimate thermal hydraulic code's treatment of a natural circulation loop and provide insights into expanded use of the computational fluid dynamics code in future designs and operations. Additionally, a sensitivity analysis is conducted to determine those physical phenomena most impactful on operations of the proposed reactor's natural circulation loop. The combination of the comparative evaluation and sensitivity analysis provides the resources for increased confidence in model developments for natural circulation loops and provides for reliability improvements of the thermal hydraulic code.
Chang, Yang; Zhao, Xiao-zhuo; Wang, Cheng; Ning, Fang-gang; Zhang, Guo-an
2015-01-01
Inhalation injury is an important cause of death after thermal burns. This study was designed to simulate the velocity and temperature distribution of inhalation thermal injury in the upper airway in humans using computational fluid dynamics. Cervical computed tomography images of three Chinese adults were imported to Mimics software to produce three-dimensional models. After grids were established and boundary conditions were defined, the simulation time was set at 1 minute and the gas temperature was set to 80 to 320°C using ANSYS software (ANSYS, Canonsburg, PA) to simulate the velocity and temperature distribution of inhalation thermal injury. Cross-sections were cut at 2-mm intervals, and maximum airway temperature and velocity were recorded for each cross-section. The maximum velocity peaked in the lower part of the nasal cavity and then decreased with air flow. The velocities in the epiglottis and glottis were higher than those in the surrounding areas. Further, the maximum airway temperature decreased from the nasal cavity to the trachea. Computational fluid dynamics technology can be used to simulate the velocity and temperature distribution of inhaled heated air. PMID:25412055
NASA Astrophysics Data System (ADS)
Septiani, Eka Lutfi; Widiyastuti, W.; Winardi, Sugeng; Machmudah, Siti; Nurtono, Tantular; Kusdianto
2016-02-01
Flame assisted spray dryer are widely uses for large-scale production of nanoparticles because of it ability. Numerical approach is needed to predict combustion and particles production in scale up and optimization process due to difficulty in experimental observation and relatively high cost. Computational Fluid Dynamics (CFD) can provide the momentum, energy and mass transfer, so that CFD more efficient than experiment due to time and cost. Here, two turbulence models, k-ɛ and Large Eddy Simulation were compared and applied in flame assisted spray dryer system. The energy sources for particle drying was obtained from combustion between LPG as fuel and air as oxidizer and carrier gas that modelled by non-premixed combustion in simulation. Silica particles was used to particle modelling from sol silica solution precursor. From the several comparison result, i.e. flame contour, temperature distribution and particle size distribution, Large Eddy Simulation turbulence model can provide the closest data to the experimental result.
Hanna, S R; Brown, M J; Camelli, F E; Chan, S T; Coirier, W J; Hansen, O R; Huber, A H; Kim, S; Reynolds, R M
2006-03-06
Computational Fluid Dynamics (CFD) model simulations of urban boundary layers have improved so that they are useful in many types of flow and dispersion analyses. The study described here is intended to assist in planning emergency response activities related to releases of chemical or biological agents into the atmosphere in large cities such as New York City. Five CFD models (CFD-Urban, FLACS, FEM3MP, FEFLO-Urban, and Fluent-Urban) have been applied by five independent groups to the same 3-D building data and geographic domain in Manhattan, using approximately the same wind input conditions. Wind flow observations are available from the Madison Square Garden March 2005 (MSG05) field experiment. It is seen from the many side-by-side comparison plots that the CFD models simulations of near-surface wind fields generally agree with each other and with field observations, within typical atmospheric uncertainties of a factor of two. The qualitative results shown here suggest, for example, that transport of a release at street level in a large city could reach a few blocks in the upwind and crosswind directions. There are still key differences seen among the models for certain parts of the domain. Further quantitative examinations of differences among the models and the observations are necessary to understand causal relationships.
Lawson, Mi. J.; Li, Y.; Sale, D. C.
2011-01-01
This paper describes the development of a computational fluid dynamics (CFD) methodology to simulate the hydrodynamics of horizontal-axis tidal current turbines (HATTs). First, an HATT blade was designed using the blade element momentum method in conjunction with a genetic optimization algorithm. Several unstructured computational grids were generated using this blade geometry and steady CFD simulations were used to perform a grid resolution study. Transient simulations were then performed to determine the effect of time-dependent flow phenomena and the size of the computational timestep on the numerical solution. Qualitative measures of the CFD solutions were independent of the grid resolution. Conversely, quantitative comparisons of the results indicated that the use of coarse computational grids results in an under prediction of the hydrodynamic forces on the turbine blade in comparison to the forces predicted using more resolved grids. For the turbine operating conditions considered in this study, the effect of the computational timestep on the CFD solution was found to be minimal, and the results from steady and transient simulations were in good agreement. Additionally, the CFD results were compared to corresponding blade element momentum method calculations and reasonable agreement was shown. Nevertheless, we expect that for other turbine operating conditions, where the flow over the blade is separated, transient simulations will be required.
Liu, Rui; Sun, Wei; Liu, Chun-Zhao
2011-01-01
A two-dimensional axisymmetric computational fluid dynamics (CFD) model based on a porous media model and a discrete population balance model was established to investigate the hydrodynamics and mass transfer behavior in an airlift bioreactor for hairy root culture.During the hairy root culture of Echinacea purpurea, liquid and gas velocity, gas holdup, mass transfer rate, as well as oxygen concentration distribution in the airlift bioreactor were simulated by this CFD model. Simulative results indicated that liquid flow and turbulence played a dominant role in oxygen mass transfer in the growth domain of the hairy root culture. The dissolved oxygen concentration in the hairy root clump increased from the bottom to the top of the bioreactor cultured with the hairy roots, which was verified by the experimental detection of dissolved oxygen concentration in the hairy root clump. This methodology provided insight understanding on the complex system of hairy root culture and will help to eventually guide the bioreactor design and process intensification of large-scale hairy root culture. PMID:22238770
Haro, Alexander J.; Dudley, Robert W.; Chelminski, Michael
2012-01-01
A two-dimensional computational fluid dynamics-habitat suitability (CFD–HSI) model was developed to identify potential zones of shallow depth and high water velocity that may present passage challenges for five anadromous fish species in the Penobscot River, Maine, upstream from two existing dams and as a result of the proposed future removal of the dams. Potential depth-challenge zones were predicted for larger species at the lowest flow modeled in the dam-removal scenario. Increasing flows under both scenarios increased the number and size of potential velocity-challenge zones, especially for smaller species. This application of the two-dimensional CFD–HSI model demonstrated its capabilities to estimate the potential effects of flow and hydraulic alteration on the passage of migratory fish.
Fluid Dynamics with Free Surfaces
Energy Science and Technology Software Center (ESTSC)
1992-02-01
RIPPLE is a two-dimensional, transient, free surface incompressible fluid dynamics program. It allows multiple free surfaces with surface tension and wall adhesion forces and has a partial cell treatment which allows curved boundaries and interior obstacles.
Fluid dynamics of heart development.
Santhanakrishnan, Arvind; Miller, Laura A
2011-09-01
The morphology, muscle mechanics, fluid dynamics, conduction properties, and molecular biology of the developing embryonic heart have received much attention in recent years due to the importance of both fluid and elastic forces in shaping the heart as well as the striking relationship between the heart's evolution and development. Although few studies have directly addressed the connection between fluid dynamics and heart development, a number of studies suggest that fluids may play a key role in morphogenic signaling. For example, fluid shear stress may trigger biochemical cascades within the endothelial cells of the developing heart that regulate chamber and valve morphogenesis. Myocardial activity generates forces on the intracardiac blood, creating pressure gradients across the cardiac wall. These pressures may also serve as epigenetic signals. In this article, the fluid dynamics of the early stages of heart development is reviewed. The relevant work in cardiac morphology, muscle mechanics, regulatory networks, and electrophysiology is also reviewed in the context of intracardial fluid dynamics. PMID:21327946
The Variety of Fluid Dynamics.
ERIC Educational Resources Information Center
Barnes, Francis; And Others
1980-01-01
Discusses three research topics which are concerned with eminently practical problems and deal at the same time with fundamental fluid dynamical problems. These research topics come from the general areas of chemical and biological engineering, geophysics, and pure mathematics. (HM)
Rigano, S; Ferrazzi, E; Boito, S; Pennati, G; Padoan, A; Galan, H
2010-01-01
The primary aim of this pilot study was to study uterine artery (UtA) blood flow volume in uneventful human pregnancies delivered at term, at mid and late gestation by means of 3D and bi-dimensional ultrasound imaging with angio-Doppler combined with fluid-dynamic modeling. Secondary aims were to correlate flow volume to placental site and to UtA Pulsatility Index (PI). Women with singleton, low-risk pregnancies were examined at mid and late gestation. The structure and course of the uterine artery (UtA) was studied in each patient by means of 3D-angio-Doppler and included vessel diameter D, blood flow velocity and PI (measured along the UtA). Fetal weight estimation and placental insertion site were assessed by ultrasound. A robust fluid-dynamic modeling was applied to calculate absolute flow and flow per unit fetal weight. Mean UtA diameter and blood flow velocity increased significantly (p < 0.0001) from mid-gestation to late gestation from 2.6 mm and 67.5 cm/s, to 3.0 mm and 85.3 cm/s, respectively, yielding an increasing absolute flow troughout gestation. h coefficient, derived by fluid-dynamic modeling to calculate mean velocity, increased significantly from 0.52 at mid-gestation to 0.57 at late gestation. UtA blood flow volume ml/min/kg-fetal weight was significantly higher at mid-gestation than at late gestation (535 ml/min/kg vs 193 ml/min/kg; p < 0.0001). In cases with strictly lateral placentas the ipsilateral UtA accommodates at mid and late gestation 63% and 67% of the total UtA flow. In central placentas UtA flow was evenly distributed between the two vessels. An inverse correlation was observed between PI and blood flow volume ml/min/kg (Pearson's coefficient r = -0.54). Our work confirms the technological and methodological limitations in the measurement of uterine artery blood flow. However, Doppler measurements supported by three-dimensional angio imaging of the uterine vessel, high resolution imaging and diameter measurement, and a robust
Fundamentals of Geophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
McWilliams, James C.
2006-07-01
Earth's atmosphere and oceans exhibit complex patterns of fluid motion over a vast range of space and time scales. These patterns combine to establish the climate in response to solar radiation that is inhomogeneously absorbed by the materials comprising air, water, and land. Spontaneous, energetic variability arises from instabilities in the planetary-scale circulations, appearing in many different forms such as waves, jets, vortices, boundary layers, and turbulence. Geophysical fluid dynamics (GFD) is the science of all these types of fluid motion. This textbook is a concise and accessible introduction to GFD for intermediate to advanced students of the physics, chemistry, and/or biology of Earth's fluid environment. The book was developed from the author's many years of teaching a first-year graduate course at the University of California, Los Angeles. Readers are expected to be familiar with physics and mathematics at the level of general dynamics (mechanics) and partial differential equations. Covers the essential GFD required for atmospheric science and oceanography courses Mathematically rigorous, concise coverage of basic theory and applications to both oceans and atmospheres Author is a world expert; this book is based on the course he has taught for many years Exercises are included, with solutions available to instructors from solutions@cambridge.org
Computational fluid dynamics - The coming revolution
NASA Technical Reports Server (NTRS)
Graves, R. A., Jr.
1982-01-01
The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.
Tröbs, Monique; Achenbach, Stephan; Röther, Jens; Redel, Thomas; Scheuering, Michael; Winneberger, David; Klingenbeck, Klaus; Itu, Lucian; Passerini, Tiziano; Kamen, Ali; Sharma, Puneet; Comaniciu, Dorin; Schlundt, Christian
2016-01-01
Invasive fractional flow reserve (FFRinvasive), although gold standard to identify hemodynamically relevant coronary stenoses, is time consuming and potentially associated with complications. We developed and evaluated a new approach to determine lesion-specific FFR on the basis of coronary anatomy as visualized by invasive coronary angiography (FFRangio): 100 coronary lesions (50% to 90% diameter stenosis) in 73 patients (48 men, 25 women; mean age 67 ± 9 years) were studied. On the basis of coronary angiograms acquired at rest from 2 views at angulations at least 30° apart, a PC-based computational fluid dynamics modeling software used personalized boundary conditions determined from 3-dimensional reconstructed angiography, heart rate, and blood pressure to derive FFRangio. The results were compared with FFRinvasive. Interobserver variability was determined in a subset of 25 narrowings. Twenty-nine of 100 coronary lesions were hemodynamically significant (FFRinvasive ≤ 0.80). FFRangio identified these with an accuracy of 90%, sensitivity of 79%, specificity of 94%, positive predictive value of 85%, and negative predictive value of 92%. The area under the receiver operating characteristic curve was 0.93. Correlation between FFRinvasive (mean: 0.84 ± 0.11) and FFRangio (mean: 0.85 ± 0.12) was r = 0.85. Interobserver variability of FFRangio was low, with a correlation of r = 0.88. In conclusion, estimation of coronary FFR with PC-based computational fluid dynamics modeling on the basis of lesion morphology as determined by invasive angiography is possible with high diagnostic accuracy compared to invasive measurements. PMID:26596195
Fluid Dynamics in Sucker Rod Pumps
Cutler, R.P.; Mansure, A.J.
1999-01-14
Sucker rod pumps are installed in approximately 90% of all oil wells in the U.S. Although they have been widely used for decades, there are many issues regarding the fluid dynamics of the pump that have not been fully investigated. A project was conducted at Sandia National Laboratories to develop unimproved understanding of the fluid dynamics inside a sucker rod pump. A mathematical flow model was developed to predict pressures in any pump component or an entire pump under single-phase fluid and pumping conditions. Laboratory flow tests were conducted on instrumented individual pump components and on a complete pump to verify and refine the model. The mathematical model was then converted to a Visual Basic program to allow easy input of fluid, geometry and pump parameters and to generate output plots. Examples of issues affecting pump performance investigated with the model include the effects of viscosity, surface roughness, valve design details, plunger and valve pressure differentials, and pumping rate.
NASA Astrophysics Data System (ADS)
Lou, Wentao; Zhu, Miaoyong
2014-10-01
A computation fluid dynamics-simultaneous reaction model (CFD-SRM) coupled model has been proposed to describe the desulfurization behavior in a gas-stirred ladle. For the desulfurization thermodynamics, different models were investigated to determine sulfide capacity and oxygen activity. For the desulfurization kinetic, the effect of bubbly plume flow, as well as oxygen absorption and oxidation reactions in slag eyes are considered. The thermodynamic and kinetic modification coefficients are proposed to fit the measured data, respectively. Finally, the effects of slag basicity and gas flow rate on the desulfurization efficiency are investigated. The results show that as the interfacial reactions (Al2O3)-(FeO)-(SiO2)-(MnO)-[S]-[O] simultaneous kinetic equilibrium is adopted to determine the oxygen activity, and the Young's model with the modification coefficient R th of 1.5 is adopted to determine slag sulfide capacity, the predicted sulfur distribution ratio LS agrees well with the measured data. With an increase of the gas blowing time, the predicted desulfurization rate gradually decreased, and when the modification parameter R k is 0.8, the predicted sulfur content changing with time in ladle agrees well with the measured data. If the oxygen absorption and oxidation reactions in slag eyes are not considered in this model, then the sulfur removal rate in the ladle would be overestimated, and this trend would become more obvious with an increase of the gas flow rate and decrease of the slag layer height. With the slag basicity increasing, the total desulfurization ratio increases; however, the total desulfurization ratio changes weakly as the slag basicity exceeds 7. With the increase of the gas flow rate, the desulfurization ratio first increases and then decreases. When the gas flow rate is 200 NL/min, the desulfurization ratio reaches a maximum value in an 80-ton gas-stirred ladle.
Fluid Dynamics of Human Phonation and Speech
NASA Astrophysics Data System (ADS)
Mittal, Rajat; Erath, Byron D.; Plesniak, Michael W.
2013-01-01
This article presents a review of the fluid dynamics, flow-structure interactions, and acoustics associated with human phonation and speech. Our voice is produced through the process of phonation in the larynx, and an improved understanding of the underlying physics of this process is essential to advancing the treatment of voice disorders. Insights into the physics of phonation and speech can also contribute to improved vocal training and the development of new speech compression and synthesis schemes. This article introduces the key biomechanical features of the laryngeal physiology, reviews the basic principles of voice production, and summarizes the progress made over the past half-century in understanding the flow physics of phonation and speech. Laryngeal pathologies, which significantly enhance the complexity of phonatory dynamics, are discussed. After a thorough examination of the state of the art in computational modeling and experimental investigations of phonatory biomechanics, we present a synopsis of the pacing issues in this arena and an outlook for research in this fascinating subject.
NASA Technical Reports Server (NTRS)
Eisner, M. (Editor)
1974-01-01
The possible utilization of the zero gravity resource for studies in a variety of fluid dynamics and fluid-dynamic related problems was investigated. A group of experiments are discussed and described in detail; these include experiments in the areas of geophysical fluid models, fluid dynamics, mass transfer processes, electrokinetic separation of large particles, and biophysical and physiological areas.
NASA Astrophysics Data System (ADS)
Schellart, W. P.
2011-06-01
Analogue models of lithospheric deformation and fluid dynamic models of mantle flow mostly use some kind of syrup such as honey or glucose syrup to simulate the low-viscosity sub-lithospheric mantle. This paper describes detailed rheological tests and density measurements of three brands of glucose syrup and three brands of honey. Additional tests have been done for one brand of glucose syrup that was diluted with water to various degrees (2%, 5% and 10% by weight). The rheological tests have been done to test the effect of shear strain, shear rate and temperature on the dynamic viscosity of the syrup. The results show that the viscosity of all glucose syrups and honeys is independent of shear strain (i.e. no strain hardening or softening). The viscosity of the glucose syrups is independent of shear rate (γ˙), i.e. linear-viscous or Newtonian, in the range γ˙=10-4-10 s with stress exponents that are almost identical to one ( n = 0.995-1.004). All the honeys show a very weak, but consistent, decrease in viscosity with increasing shear rate of 7-14% from 10 -3 to 10 0 s -1 and have stress exponents more distinct from one ( n = 1.007-1.026). All syrups have a viscosity that is strongly dependent on temperature in the range 0-50 °C, where viscosity decreases with increasing temperature. Such decrease can be fitted with exponential and Arrhenius functions, with the latter giving the best results. Furthermore, the viscosity of glucose syrup decreases approximately exponentially with increasing water content. Oscillation tests indicate that the rheology of all the syrups is entirely dominated by viscous behaviour and not by elastic behaviour at frequencies of 10 -3-10 2 Hz. Finally, the density investigations show that the density of glucose syrup and honey decreases approximately linearly with increasing temperature in the range 10-70 °C, with coefficients of thermal volumetric expansion at 20 °C of 3.89-3.95 × 10 -4 °C -1 and 4.57-4.81 × 10 -4 °C -1 for
ADDRESSING ENVIRONMENTAL ENGINEERING CHALLENGES WITH COMPUTATIONAL FLUID DYNAMICS
This paper discusses the status and application of Computational Fluid Dynamics )CFD) models to address environmental engineering challenges for more detailed understanding of air pollutant source emissions, atmospheric dispersion and resulting human exposure. CFD simulations ...
Pennati, Giancarlo; Corsini, Chiara; Cosentino, Daria; Hsia, Tain-Yen; Luisi, Vincenzo S.; Dubini, Gabriele; Migliavacca, Francesco
2011-01-01
Cavopulmonary connections are surgical procedures used to treat a variety of complex congenital cardiac defects. Virtual pre-operative planning based on in silico patient-specific modelling might become a powerful tool in the surgical decision-making process. For this purpose, three-dimensional models can be easily developed from medical imaging data to investigate individual haemodynamics. However, the definition of patient-specific boundary conditions is still a crucial issue. The present study describes an approach to evaluate the vascular impedance of the right and left lungs on the basis of pre-operative clinical data and numerical simulations. Computational fluid dynamics techniques are applied to a patient with a bidirectional cavopulmonary anastomosis, who later underwent a total cavopulmonary connection (TCPC). Multi-scale models describing the surgical region and the lungs are adopted, while the flow rates measured in the venae cavae are used at the model inlets. Pre-operative and post-operative conditions are investigated; namely, TCPC haemodynamics, which are predicted using patient-specific pre-operative boundary conditions, indicates that the pre-operative balanced lung resistances are not compatible with the TCPC measured flows, suggesting that the pulmonary vascular impedances changed individually after the surgery. These modifications might be the consequence of adaptation to the altered pulmonary blood flows. PMID:22670201
NASA Technical Reports Server (NTRS)
Spirka, T. A.; Myers, J. G.; Setser, R. M.; Halliburton, S. S.; White, R. D.; Chatzimavroudis, G. P.
2005-01-01
A priority of NASA is to identify and study possible risks to astronauts health during prolonged space missions [l]. The goal is to develop a procedure for a preflight evaluation of the cardiovascular system of an astronaut and to forecast how it will be affected during the mission. To predict these changes, a computational cardiovascular model must be constructed. Although physiology data can be used to make a general model, a more desirable subject-specific model requires anatomical, functional, and flow data from the specific astronaut. MRI has the unique advantage of providing images with all of the above information, including three-directional velocity data which can be used as boundary conditions in a computational fluid dynamics (CFD) program [2,3]. MRI-based CFD is very promising for reproduction of the flow patterns of a specific subject and prediction of changes in the absence of gravity. The aim of this study was to test the feasibility of this approach by reconstructing the geometry of MRI-scanned arterial models and reproducing the MRI-measured velocities using CFD simulations on these geometries.
Hoi, Yiemeng; Ionita, Ciprian N.; Tranquebar, Rekha V.; Hoffmann, Kenneth R.; Woodward, Scott, H.; Taulbee, Dale B.; Meng, Hui; Rudin, Stephen
2011-01-01
An asymmetric stent with low porosity patch across the intracranial aneurysm neck and high porosity elsewhere is designed to modify the flow to result in thrombogenesis and occlusion of the aneurysm and yet to reduce the possibility of also occluding adjacent perforator vessels. The purposes of this study are to evaluate the flow field induced by an asymmetric stent using both numerical and digital subtraction angiography (DSA) methods and to quantify the flow dynamics of an asymmetric stent in an in vivo aneurysm model. We created a vein-pouch aneurysm model on the canine carotid artery. An asymmetric stent was implanted at the aneurysm, with 25% porosity across the aneurysm neck and 80% porosity elsewhere. The aneurysm geometry, before and after stent implantation, was acquired using cone beam CT and reconstructed for computational fluid dynamics (CFD) analysis. Both steady-state and pulsatile flow conditions using the measured waveforms from the aneurysm model were studied. To reduce computational costs, we modeled the asymmetric stent effect by specifying a pressure drop over the layer across the aneurysm orifice where the low porosity patch was located. From the CFD results, we found the asymmetric stent reduced the inflow into the aneurysm by 51%, and appeared to create a stasis-like environment which favors thrombus formation. The DSA sequences also showed substantial flow reduction into the aneurysm. Asymmetric stents may be a viable image guided intervention for treating intracranial aneurysms with desired flow modification features. PMID:21666881
Computational Fluid Dynamics Technology for Hypersonic Applications
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.
2003-01-01
Several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented from code validation and code benchmarking efforts to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified. Highlights of diverse efforts to address these challenges are then discussed. One such effort to re-engineer and synthesize the existing analysis capability in LAURA, VULCAN, and FUN3D will provide context for these discussions. The critical (and evolving) role of agile software engineering practice in the capability enhancement process is also noted.
NASA Astrophysics Data System (ADS)
Santini, Maurizio
2015-11-01
X-ray computed tomography (CT) is a well-known technique nowadays, since its first practical application by Sir. G. Hounsfield (Nobel price for medicine 1979) has continually benefited from optimising improvements, especially in medical applications. Indeed, also application of CT in various engineering research fields provides fundamental informations on a wide range of applications, considering that the technique is not destructive, allowing 3D visualization without perturbation of the analysed material. Nowadays, it is technologically possible to design and realize an equipment that achieve a micrometric resolution and even improve the sensibility in revealing differences in materials having very radiotransparency, allowing i.e. to distinguish between different fluids (with different density) or states of matter (like with two-phase flows). At the University of Bergamo, a prototype of an X-ray microCT system was developed since 2008, so being fully operative from 2012, with specific customizations for investigations in thermal-fluid dynamics and multiphase flow researches. A technical session held at the UIT International Conference in L'Aquila (Italy), at which this paper is referring, has presented some microCT fundamentals, to allow the audience to gain basics to follow the “fil-rouge” that links all the instrumentation developments, till the recent applications. Hereinafter are reported some applications currently developed at Bergamo University at the X-ray computed micro-tomography laboratory.
Viitala, Tapani; Liang, Huamin; Gupta, Mayur; Zwinger, Thomas; Yliperttula, Marjo; Bunker, Alex
2012-07-15
We have used computational fluid dynamics modeling (CFD) to synchronize the flow conditions in the flow channels of two complementary surface-sensitive characterization techniques: surface plasmon resonance (SPR) and quartz crystal microbalance (QCM). Since the footprint of the flow channels of the two devices is specified by their function, the flow behavior can only be varied either by altering the height of the flow channel, or altering the volumetric rate of flow (flow rate) through the channel. The relevant quantity that must be calibrated is the shear strain on the measurement surface (center and bottom) of the flow channel. Our CFD modeling shows that the flow behavior is in the Stokes flow regime. We were thus able to generate a scaling expression with parameters for flow rate and flow channel height for each of the two devices: f(QCM)=2.64f(SPR)(h(QCM)/h(SPR)(2), where f(QCM) and f(SPR) are the flow rates in the SPR and QCM flow channels, respectively, and h(QCM)/h(SPR) is the ratio of the heights of the two channels. We demonstrate the success of our calibration procedure through the combined use of commercially available SPR and QCM flow channel devices on both a biomolecular interaction system of surface immobilized biotin and streptavidin and a targeted drug delivery model system of biotinylated liposomes interacting with a streptavidin functionalized surface. PMID:22579516
Meier, U; Kiefer, M
2001-04-01
The internationally accepted methods of calculating cerebrospinal fluid dynamics proceed from the assumption of a pressure-independent resistance to CSF outflow. Our new model focusses on the pressure-dependency of this resistance. In it, we monitor the entire pressure course over time, p(t) during and after infusion. A comparison of the pressure rise, On(p), during infusion, and the decrease, Off(p), to the same pressure level, permits the creation of all the formulas for C(p) and R(p). The simultaneous measurement of resistance and compliance during a single intervention allows us to minimize patient exertion. In contrast to the classical methods, it is not necessary for the ICP to reach a plateau. Our mathematical model differs from the static examination model by describing a pressure-dependent slope of the function for the resistance. This has been demonstrated in a study using H-Tx rats. In this way, we are able to take the non-linearity of the CSF resorption into consideration. PMID:11388039
Fluid Dynamic Verification Experiments on STS-70
NASA Technical Reports Server (NTRS)
Kleis, Stanley J.
1996-01-01
Fluid dynamic experiments were flown on STS-70 as phase two of the engineering evaluation of the first bioreactor Engineering Development Unit (EDU#1). The phase one experiments were comparative cell cultures in identical units on earth and onboard STS-70. In phase two, two types of fluid dynamic experiments were performed. Qualitative comparisons of the basic flow patterns were evaluated with the use of 'dye' streaklines formed from alternate injections of either a mild acid or base solution into the external flow loop that was then perfused into the vessel. The presence of Bromothymol Blue in the fluid then caused color changes from yellow to blue or vice versa, indicating the basic fluid motions. This reversible change could be repeated as desired. In the absence of significant density differences in the fluid, the flow patterns in space should be the same as on earth. Video tape records of the flow patterns for a wide range of operating conditions were obtained. The second type of fluid dynamic experiment was the quantitative evaluation of the trajectories of solid beads of various densities and sizes. The beads were introduced into the vessel and the paths recorded on video tape, with the vessel operated at various rotation rates and flow perfusion rates. Because of space limitations, the video camera was placed as close as possible to the vessel, resulting in significant optical distortion. This report describes the analysis methods to obtain comparisons between the in-flight fluid dynamics and numerical models of the flow field. The methods include optical corrections to the video images and calculation of the bead trajectories for given operating conditions and initial bead locations.
Computational fluid dynamics in cardiovascular disease.
Lee, Byoung-Kwon
2011-08-01
Computational fluid dynamics (CFD) is a mechanical engineering field for analyzing fluid flow, heat transfer, and associated phenomena, using computer-based simulation. CFD is a widely adopted methodology for solving complex problems in many modern engineering fields. The merit of CFD is developing new and improved devices and system designs, and optimization is conducted on existing equipment through computational simulations, resulting in enhanced efficiency and lower operating costs. However, in the biomedical field, CFD is still emerging. The main reason why CFD in the biomedical field has lagged behind is the tremendous complexity of human body fluid behavior. Recently, CFD biomedical research is more accessible, because high performance hardware and software are easily available with advances in computer science. All CFD processes contain three main components to provide useful information, such as pre-processing, solving mathematical equations, and post-processing. Initial accurate geometric modeling and boundary conditions are essential to achieve adequate results. Medical imaging, such as ultrasound imaging, computed tomography, and magnetic resonance imaging can be used for modeling, and Doppler ultrasound, pressure wire, and non-invasive pressure measurements are used for flow velocity and pressure as a boundary condition. Many simulations and clinical results have been used to study congenital heart disease, heart failure, ventricle function, aortic disease, and carotid and intra-cranial cerebrovascular diseases. With decreasing hardware costs and rapid computing times, researchers and medical scientists may increasingly use this reliable CFD tool to deliver accurate results. A realistic, multidisciplinary approach is essential to accomplish these tasks. Indefinite collaborations between mechanical engineers and clinical and medical scientists are essential. CFD may be an important methodology to understand the pathophysiology of the development and
Neidlin, Michael; Steinseifer, Ulrich; Kaufmann, Tim A S
2014-06-01
Neurological complication often occurs during cardiopulmonary bypass (CPB). One of the main causes is hypoperfusion of the cerebral tissue affected by the position of the cannula tip and diminished cerebral autoregulation (CA). Recently, a lumped parameter approach could describe the baroreflex, one of the main mechanisms of cerebral autoregulation, in a computational fluid dynamics (CFD) study of CPB. However, the cerebral blood flow (CBF) was overestimated and the physiological meaning of the variables and their impact on the model was unknown. In this study, we use a 0-D control circuit representation of the Baroreflex mechanism, to assess the parameters with respect to their physiological meaning and their influence on CBF. Afterwards the parameters are transferred to 3D-CFD and the static and dynamic behavior of cerebral autoregulation is investigated. The parameters of the baroreflex mechanism can reproduce normotensive, hypertensive and impaired autoregulation behavior. Further on, the proposed model can mimic the effects of anesthetic agents and other factors controlling dynamic CA. The CFD simulations deliver similar results of static and dynamic CBF as the 0-D control circuit. This study shows the feasibility of a multiscale 0-D/3-D approach to include patient-specific cerebral autoregulation into CFD studies. PMID:24746017
Lumb, A B; Burns, A D; Figueroa Rosette, J A; Gradzik, K B; Ingham, D B; Pourkashanian, M
2015-05-01
We have used computational fluid dynamic modelling to study the effects of tracheal tube size and position on regional gas flow in the large airways. Using a three-dimensional mathematical model, we simulated flow with and without a tracheal tube, replicating both physiological and artificial breathing. Ventilation through a tracheal tube increased proportional flow to the left lung from 39.5% with no tube to 43.1-47.2%, depending on tube position. Ventilation mode and tube distance from the carina had no effect on flow. Lateral displacement and deflection of the tube increased ventilation to the ipsilateral lung; for example, when deflected 10° to the left of centre, flow to the left lung increased from 43.8 to 53.7%. Because of the small diameter of a tracheal tube relative to the trachea, gas exits a tube at high velocity such that regional ventilation may be affected by changes in the position and angle of the tube. PMID:25581493
Liu, Rui; Sun, Wei; Liu, Chun-Zhao
2011-01-01
Recently, cichoric acid production from hairy roots of Echinacea purpurea was significantly improved by ultrasound stimulation in an airlift bioreactor. In this article, the possible mechanism on ultrasound-intensified hairy root culture of E. purpurea in the bioreactor was elucidated with the help of computational fluid dynamics (CFD) simulation, membrane permeability detection, dissolved oxygen concentration detection, confocal laser-scanning microscopy (LSM) observation, and phenylalanine ammonium lyase (PAL) activity analysis. The CFD model developed in Part I was used to simulate the hydrodynamics and oxygen mass transfer in hairy root bioreactor culture stimulated by ultrasound. A dynamic mesh model combined with a changing Schmidt number method was used for the simulation of the ultrasound field. Simulation results and experimental data illustrated that ultrasound intensified oxygen mass transfer in the hairy root clump, which subsequently stimulated root growth and cichoric acid biosynthesis. Ultrasound increased the hairy root membrane permeability, and a high root membrane permeability of 0.359 h(-1) was observed at the bottom region in the bioreactor. LSM observation showed that the change in the membrane permeability recovered to normal in the further culture after ultrasound stimulation. PAL activity in the hairy roots was stimulated by ultrasound increase and was correlated well to cichoric acid accumulation in the hairy roots of E. purpurea. PMID:21850671
Stellar Astrophysical Fluid Dynamics
NASA Astrophysics Data System (ADS)
Thompson, Michael J.; Christensen-Dalsgaard, Jørgen
2003-05-01
Preface; 1. A selective overview Jørgen Christensen-Dalsgaard and Michael J. Thompson; Part I. Stellar Convection and Oscillations: 2. On the diversity of stellar pulsations Wojciech A. Dziembowski; 3. Acoustic radiation and mode excitation by turbulent convection Günter Houdek; 4. Understanding roAp stars Margarida S. Cunha; 5. Waves in the magnetised solar atmosphere Colin S. Rosenthal; Part II. Stellar Rotation and Magnetic Fields: 6. Stellar rotation: a historical survey Leon Mestel; 7. The oscillations of rapidly rotating stars Michel Rieutord; 8. Solar tachocline dynamics: eddy viscosity, anti-friction, or something in between? Michael E. McIntyre; 9. Dynamics of the solar tachocline Pascale Garaud; 10. Dynamo processes: the interaction of turbulence and magnetic fields Michael Proctor; 11. Dynamos in planets Chris Jones; Part III. Physics and Structure of Stellar Interiors: 12. Solar constraints on the equation of state Werner Däppen; 13. 3He transport and the solar neutrino problem Chris Jordinson; 14. Mixing in stellar radiation zones Jean-Paul Zahn; 15. Element settling and rotation-induced mixing in slowly rotating stars Sylvie Vauclair; Part IV. Helio- and Asteroseismology: 16. Solar structure and the neutrino problem Hiromoto Shibahashi; 17. Helioseismic data analysis Jesper Schou; 18. Seismology of solar rotation Takashi Sekii; 19. Telechronohelioseismology Alexander Kosovichev; Part V. Large-Scale Numerical Experiments: 20. Bridges between helioseismology and models of convection zone dynamics Juri Toomre; 21. Numerical simulations of the solar convection zone Julian R. Elliott; 22. Modelling solar and stellar magnetoconvection Nigel Weiss; 23. Nonlinear magnetoconvection in the presence of a strong oblique field Keith Julien, Edgar Knobloch and Steven M. Tobias; 24. Simulations of astrophysical fluids Marcus Brüggen; Part VI. Dynamics: 25. A magic electromagnetic field Donald Lynden-Bell; 26. Continuum equations for stellar dynamics Edward A