Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fast Fluid Dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Topography and tectonics of mid-oceanic ridge axes

Numerical fluid dynamic models of mid-oceanic ridge axes were constructed using distributions of material properties constrained by seismic studies and thermal calculations. The calculations indicate that spreading is passive except for forces caused by density differences due to thermal expansion and partial melt. Except for geometric differences due to temperature distribution, one set of mechanical properties can adequately explain central rifts of slow ridges and central peaks of fast ridges. Viscous head loss in the upwelling material dominates at low spreading rates where material ascends through a narrow conduit. Thermal expansion and partial melting dominate at high spreading rates where a wide low viscosity crustal magma chamber is present. The preferred rheology is 10 to the 20th poise for the upwelling lithosphere; less than 5 x 10 to the 17th for the crustal magma chamber and axial intrusion zone at fast ridges, and a yield stress of 200 bars for the lithosphere. The calculation correctly predicts the existence of central peaks at 'hot-spot' ridges, where seismic evidence indicates a large magma chamber.

Sleep, N. H.↗

Ionization effects in three-dimensional solar granulation simulations

These numerical studies show that ionization influences both the transport and dynamical properties of compressible convection near the surface of the Sun. About two-thirds of the enthalpy transported by convective motions in the region of partial hydrogen ionization is carried as latent heat. The role of fast downflow plumes in total convective transport is substantially elevated by this contribution. Instability of the thermal boundary layer is strongly enhanced by temperature sensitive variations in the radiative properties of the fluid, and this provides a mechanism for plume initiation and cell fragmentation in the surface layers. As the plumes descend, temperature fluctuations and associated buoyancy forces are maintained because of the increased specific heat of the partially ionized material. This can result is supersonic vertical flows. At greater depths, ionization effects diminish, and the plumes are decelerated by significant entrainment of surrounding fluid.

Rast, Mark P.↗

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗

CFD Analysis for Assessing the Effect of Wind on the Thermal Control of the Mars Science Laboratory Curiosity Rover

The challenging range of landing sites for which the Mars Science Laboratory Rover was designed, requires a rover thermal management system that is capable of keeping temperatures controlled across a wide variety of environmental conditions. On the Martian surface where temperatures can be as cold as -123 C and as warm as 38 C, the rover relies upon a Mechanically Pumped Fluid Loop (MPFL) Rover Heat Rejection System (RHRS) and external radiators to maintain the temperature of sensitive electronics and science instruments within a -40 C to 50 C range. The RHRS harnesses some of the waste heat generated from the rover power source, known as the Multi Mission Radioisotope Thermoelectric Generator (MMRTG), for use as survival heat for the rover during cold conditions. The MMRTG produces 110 W of electrical power while generating waste heat equivalent to approximately 2000 W. Heat exchanger plates (hot plates) positioned close to the MMRTG pick up this survival heat from it by radiative heat transfer. Winds on Mars can be as fast as 15 m/s for extended periods. They can lead to significant heat loss from the MMRTG and the hot plates due to convective heat pick up from these surfaces. Estimation of this convective heat loss cannot be accurately and adequately achieved by simple textbook based calculations because of the very complicated flow fields around these surfaces, which are a function of wind direction and speed. Accurate calculations necessitated the employment of sophisticated Computational Fluid Dynamics (CFD) computer codes. This paper describes the methodology and results of these CFD calculations. Additionally, these results are compared to simple textbook based calculations that served as benchmarks and sanity checks for them. And finally, the overall RHRS system performance predictions will be shared to show how these results affected the overall rover thermal performance.

wind↗

Hybrid nanocapsules for in situ TEM imaging of gas evolution reactions in confined liquids

Liquid cell transmission electron microscopy (TEM) enables the direct observation of dynamic physical and chemical processes in liquids at the nanoscale. Quantitative investigations into reactions with fast kinetics and/or multiple reagents will benefit from further advances in liquid cell design that facilitate rapid in situ mixing and precise control over reagent volumes and concentrations. This work reports the development of inorganic–organic nanocapsules for high-resolution TEM imaging of nanoscale reactions in liquids with well-defined zeptoliter volumes. These hybrid nanocapsules, with 48 nm average diameter, consist of a thin layer of gold coating a lipid vesicle. As a model reaction, the nucleation, growth, and diffusion of nanobubbles generated by the radiolysis of water is investigated inside the nanocapsules. When the nanobubbles are sufficiently small (10–25 nm diameter), they are mobile in the nanocapsules, but their movement deviates from Brownian motion, which may result from geometric confinement by the nanocapsules. Gases and fluids can be transported between two nanocapsules when they fuse, demonstrating in situ mixing without using complex microfluidic schemes. The ability to synthesize nanocapsules with controlled sizes and to monitor dynamics simultaneously inside multiple nanocapsules provides opportunities to investigate nanoscale processes such as single nanoparticle synthesis in confined volumes and biological processes such as biomineralization and membrane dynamics.

36 MATERIALS SCIENCE↗

Remote Visualization and Remote Collaboration On Computational Fluid Dynamics

A new technology has been developed for remote visualization that provides remote, 3D, high resolution, dynamic, interactive viewing of scientific data (such as fluid dynamics simulations or measurements). Based on this technology, some World Wide Web sites on the Internet are providing fluid dynamics data for educational or testing purposes. This technology is also being used for remote collaboration in joint university, industry, and NASA projects in computational fluid dynamics and wind tunnel testing. Previously, remote visualization of dynamic data was done using video format (transmitting pixel information) such as video conferencing or MPEG movies on the Internet. The concept for this new technology is to send the raw data (e.g., grids, vectors, and scalars) along with viewing scripts over the Internet and have the pixels generated by a visualization tool running on the viewer's local workstation. The visualization tool that is currently used is FAST (Flow Analysis Software Toolkit).

Watson, Val↗

Interactive Visualization of Computational Fluid Dynamics using Mosaic

The Web provides new Methods for accessing Information world-wide, but the current text-and-pictures approach neither utilizes all the Web's possibilities not provides for its limitations. While the inclusion of pictures and animations in a paper communicates more effectively than text alone, It Is essentially an extension of the concept of "publication." Also, as use of the Web increases putting images and animations online will quickly load even the "Information Superhighway." We need to find forms of communication that take advantage of the special nature of the Web. This paper presents one approach: the use of the Internet and the Mosaic interface for data sharing and collaborative analysis. We will describe (and In the presentation, demonstrate) our approach: using FAST (Flow Analysis Software Toolkit), a scientific visualization package, as a data viewer and interactive tool called from MOSAIC. Our intent is to stimulate the development of other tools that utilize the unique nature of electronic communication.

Clucas, Jean↗

Optimization of Scrap Melting Using an Electric Arc in Steel Manufacturing

Steel industry is crucial to the national economy and security. Around 67% of crude steel in the U.S is produced in electric arc furnaces (EAF), which is energy intensive. Around 140 EAFs operate in the U.S., consuming about 8.6x10 7 MMBtu/year of electricity. One of major challenges for EAFs includes maximizing the efficiency of the electrical energy provided in the form of electric arcs to melt various scrap mixes. To address this issue, a computational fluid dynamics (CFD) methodology is chosen to analyze scrap melting using the electric arc. Due to complex furnace phenomena and the wide variety of potential scenarios, high performance computing (HPC) is essential to yield comprehensive and detailed CFD analyses and systematic parametric studies for optimized EAF operation. The objectives are to 1) simulate scrap melting using electric arc, 2) evaluate electrode/arc position for optimum scrap melting and 3) establish reduced order model for CFD data-base for fast model calculation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Comparison of Computational Results with a Low-g, Nitrogen Slosh and Boiling Experiment

This paper compares a fluid/thermal simulation, in Fluent, with a low-g, nitrogen slosh and boiling experiment. In 2010, the French Space Agency, CNES, performed cryogenic nitrogen experiments in a low-g aircraft campaign. From one parabolic flight, a low-g interval was simulated that focuses on low-g motion of nitrogen liquid and vapor with significant condensation, evaporation, and boiling. The computational results are compared with high-speed video, pressure data, heat transfer, and temperature data from sensors on the axis of the cylindrically shaped tank. These experimental and computational results compare favorably. The initial temperature stratification is in good agreement, and the two-phase fluid motion is qualitatively captured. Temperature data is matched except that the temperature sensors are unable to capture fast temperature transients when the sensors move from wet to dry (liquid to vapor) operation. Pressure evolution is approximately captured, but condensation and evaporation rate modeling and prediction need further theoretical analysis.

Cryogenic Fluid Storage↗

Coupling coarse mesh CFD with fine mesh CFD for modeling for Molten Salt Reactors in the Virtual Test Bed VTB

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims at developing a simulation tool kit to accelerate the development and deployment of nuclear power technologies. NEAMS multiphysics tools have been designed to provide numerical simulation support for the design and licensing of GEN IV reactors. Pronghorn is NEAMS’s coarse-mesh computational fluid dynamics (CFD) tool, which is designed to run 3D core transients in GEN IV reactors at a reduced computational cost. To increase their accuracy, coarse-mesh CFD simulations require calibrated closure coefficients. One way of computing these coefficients is via the Nek5000, NEAMS’s high-fidelity CFD tool. This article discusses our current research lines in informing Pronghorn closure coefficients via Nek5000 to enable multiphysics simulations of the core cavity of the molten-salt fast reactor. We present an application in which Pronghorn mixing length turbulent viscosity has been calibrated from Nek5000 simulations. The resulting Pronghorn thermal-hydraulics model is then coupled to Griffin, the NEAMS neutron transport solver, to solve for the steady-state configuration of the molten-salt fast reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A drift kinetic model for the expander region of a magnetic mirror

Here, we present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ~5T e /e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Internal Mixing Studied for GE/ARL Ejector Nozzle

To achieve jet noise reduction goals for the High Speed Civil Transport aircraft, researchers have been investigating the mixer-ejector nozzle concept. For this concept, a primary nozzle with multiple chutes is surrounded by an ejector. The ejector mixes low-momentum ambient air with the hot engine exhaust to reduce the jet velocity and, hence, the jet noise. It is desirable to mix the two streams as fast as possible in order to minimize the length and weight of the ejector. An earlier model of the mixer-ejector nozzle was tested extensively in the Aerodynamic Research Laboratory (ARL) of GE Aircraft Engines at Cincinnati, Ohio. While testing was continuing with later generations of the nozzle, the earlier model was brought to the NASA Lewis Research Center for relatively fundamental measurements. Goals of the Lewis study were to obtain details of the flow field to aid computational fluid dynamics (CFD) efforts and obtain a better understanding of the flow mechanisms, as well as to experiment with mixing enhancement devices, such as tabs. The measurements were made in an open jet facility for cold (unheated) flow without a surrounding coflowing stream.

Zaman, Khairul↗

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗

Automated bidirectional coupling of multiscale models of aerosol dosimetry: Validation with subject-specific deposition data

Assessing the toxicity of airborne particulate matter or the efficacy of inhaled drug depends upon accurate estimates of deposited fraction of inhaled materials. In silico approaches can provide important insights into site- or airway-specific deposition of inhaled aerosols in the respiratory system. In this study, we improved on our recently developed 3D/1D model that simulate aerosol transport and deposition in the whole lung over multiple breath cycles. A subject-specific multiscale lung model of a healthy male subject using computational fluid-particle dynamics (CFPD) in a 3D model of the oral cavity through the large bronchial airways entering each lobe was bidirectionally coupled with a recently improved Multiple Path Particle Dosimetry (MPPD) model to predict aerosol deposition over the entire respiratory tract over multiple breaths for four conditions matching experimental aerosol exposures in the same subject from which the model was developed. These include two particle sizes (1 and 2.9 μm) and two subject-specific breathing rates of ~300 ml/s (slow breathing) and ~750 ml/s (fast breathing) at a target tidal volume of 1 L. In silico predictions of retained fraction were 0.31 and 0.29 for 1 μm and 0.66 and 0.62 for 2.9 μm during slow and fast breathing, respectively, and compared well with experimental data (1 μm: 0.31 ± 0.01 (slow) and 0.27 ± 0.01 (fast), 2.9 μm: 0.63 ± 0.03 (slow) and 0.68 ± 0.02 (fast)). Furthermore, these results provide a great deal of confidence in the validity and reliability of our approach.

60 APPLIED LIFE SCIENCES↗

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗

A numerical model for dynamic wave rotor analysis

A numerical model has been developed which can predict the dynamic (and steady state) performance of a wave rotor, given the geometry and time dependent boundary conditions. The one-dimensional, perfect gas, CFD based code tracks the gasdynamics in each of the wave rotor passages as they rotate past the various ducts. The model can operate both on and off-design, allowing dynamic behavior to be studied throughout the operating range of the wave rotor. The model accounts for several major loss mechanisms including finite passage opening time, fluid friction, heat transfer to and from the passage walls, and leakage to and from the passage ends. In addition, it can calculate the amount of work transferred to and from the fluid when the flow in the ducts is not aligned with the passages such as occurs in off-design operation. Since it is one-dimensional, the model runs reasonably fast on a typical workstation. This paper will describe the model and present the results of some transient calculations for a conceptual four port wave rotor designed as a topping cycle for a small gas turbine engine.

Paxson, D. E.↗

Supersonic shock wave/vortex interaction

Although shock wave/vortex interaction is a basic and important fluid dynamics problem, very little research has been conducted on this topic. Therefore, a detailed experimental study of the interaction between a supersonic streamwise turbulent vortex and a shock wave was carried out at the Penn State Gas Dynamics Laboratory. A vortex is produced by replaceable swirl vanes located upstream of the throat of various converging-diverging nozzles. The supersonic vortex is then injected into either a coflowing supersonic stream or ambient air. The structure of the isolated vortex is investigated in a supersonic wind tunnel using miniature, fast-response, five-hole and total temperature probes and in a free jet using laser Doppler velocimetry. The cases tested have unit Reynolds numbers in excess of 25 million per meter, axial Mach numbers ranging from 2.5 to 4.0, and peak tangential Mach numbers from 0 (i.e., a pure jet) to about 0.7. The results show that the typical supersonic wake-like vortex consists of a non-isentropic, rotational core, where the reduced circulation distribution is self similar, and an outer isentropic, irrotational region. The vortex core is also a region of significant turbulent fluctuations. Radial profiles of turbulent kinetic energy and axial-tangential Reynolds stress are presented. The interactions between the vortex and both oblique and normal shock waves are investigated using nonintrusive optical diagnostics (i.e. schlieren, planar laser scattering, and laser Doppler velocimetry). Of the various types, two Mach 2.5 overexpanded-nozzle Mach disc interactions are examined in detail. Below a certain vortex strength, a 'weak' interaction exists in which the normal shock is perturbed locally into an unsteady 'bubble' shock near the vortex axis, but vortex breakdown (i.e., a stagnation point) does not occur. For stronger vortices, a random unsteady 'strong' interaction results that causes vortex breakdown. The vortex core reforms downstream of the rear stagnation point, and the reduced circulation distribution once again becomes self-similar in this region. A-new model of this interaction is proposed. Finally, a curve defining the approximate limits of supersonic vortex breakdown is presented.

Settles, G. S.↗

Liquid–Liquid Microfluidic Flows for Ultrafast 5-Hydroxymethyl Furfural Extraction

We study the effects of liquid–liquid flow patterns on the extraction of 5-hydroxymethyl furfural from the aqueous phase into an organic phase (methyl isobutyl ketone, ethyl acetate, and 2-pentanol) at ultrashort residence times of 0.2–20 s in a coiled capillary microchannel. Slug, droplet, slug-droplet, dispersed, parallel, annular, and irregular two-phase flow patterns are generated. We evaluate the extraction efficiency and the volumetric mass transfer coefficient via experiments and computational fluid dynamics (CFD). The CFD predictions are in quantitative agreement with the experimental data. The convective and diffusive contributions to the mass transfer of different flow patterns are analyzed. For all solvent pairs, the extraction efficiency varies irregularly with increasing flow rate due to changes in the mass transfer mechanism caused by a transition between flow regimes. Overall, we demonstrate >90% extraction efficiencies with high volumetric mass transfer coefficients ranging from 0.006 to 2.17 s–1. Furthermore, a merit index is defined to assess the trade-off between the fast mass transfer and the pressure drop in biphasic microchannels. In general, the mass transfer rate increases significantly at high flow rates at the expense of a high-pressure drop. Furthermore, the ethyl acetate/water irregular flow pattern is optimal due to a large throughput, extraction efficiency, and mass transfer coefficient, with a minimal pressure drop.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗