SEARCH · Search NASA

Results for “complex fluid”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,477 records · Page 82

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Monte Carlo Methodology Serves Up a Software Success

Widely used for the modeling of gas flows through the computation of the motion and collisions of representative molecules, the Direct Simulation Monte Carlo method has become the gold standard for producing research and engineering predictions in the field of rarefied gas dynamics. Direct Simulation Monte Carlo was first introduced in the early 1960s by Dr. Graeme Bird, a professor at the University of Sydney, Australia. It has since proved to be a valuable tool to the aerospace and defense industries in providing design and operational support data, as well as flight data analysis. In 2002, NASA brought to the forefront a software product that maintains the same basic physics formulation of Dr. Bird's method, but provides effective modeling of complex, three-dimensional, real vehicle simulations and parallel processing capabilities to handle additional computational requirements, especially in areas where computational fluid dynamics (CFD) is not applicable. NASA's Direct Simulation Monte Carlo Analysis Code (DAC) software package is now considered the Agency s premier high-fidelity simulation tool for predicting vehicle aerodynamics and aerothermodynamic environments in rarified, or low-density, gas flows.

Source record↗

NASA NTRS2005

AirShow 1.0 CFD Software Users' Guide

AirShow is visualization post-processing software for Computational Fluid Dynamics (CFD). Upon reading binary PLOT3D grid and solution files into AirShow, the engineer can quickly see how hundreds of complex 3-D structured blocks are arranged and numbered. Additionally, chosen grid planes can be displayed and colored according to various aerodynamic flow quantities such as Mach number and pressure. The user may interactively rotate and translate the graphical objects using the mouse. The software source code was written in cross-platform Java, C++, and OpenGL, and runs on Unix, Linux, and Windows. The graphical user interface (GUI) was written using Java Swing. Java also provides multiple synchronized threads. The Java Native Interface (JNI) provides a bridge between the Java code and the C++ code where the PLOT3D files are read, the OpenGL graphics are rendered, and numerical calculations are performed. AirShow is easy to learn and simple to use. The source code is available for free from the NASA Technology Transfer and Partnership Office.

Mohler, Stanley R., Jr.↗

DOE OSTI2010

Geochemical engineering design tools for uranium in situ recovery : the HYDROGEOCHEM codes.

Geochemical Engineering Design (GED) is based on applications of the principles and various computer models that describe the biogeochemistry and physics of removal of contaminants from water by adsorption, precipitation and filtration. It can be used to optimize or evaluate the efficiency of all phases of in situ recovery (ISR). The primary tools of GED are reactive transport models; this talk describes the potential application of the HYDROGEOCHEM family of codes to ISR. The codes can describe a complete suite of equilibrium or kinetic aqueous complexation, adsorption-desorption, precipitation-dissolution, redox, and acid-base reactions in variably saturated media with density-dependent fluid flow. Applications to ISR are illustrated with simulations of (1) the effectiveness of a reactive barrier to prevent off-site uranium migration and (2) evaluation of the effect of sorption hysteresis on natural attenuation. In the first example, it can be seen that the apparent effectiveness of the barrier depends on monitoring location and that it changes over time. This is due to changes in pH, saturation of sorption sites, as well as the geometry of the flow field. The second simulation shows how sorption hysteresis leads to observable attenuation of a uranium contamination plume. Different sorption mechanisms including fast (or reversible), slow, and irreversible sorption were simulated. The migration of the dissolved and total uranium plumes for the different cases are compared and the simulations show that when 50-100% of the sites have slow desorption rates, the center of mass of the dissolved uranium plume begins to move upstream. This would correspond to the case in which the plume boundaries begin to shrink as required for demonstration of natural attenuation.

58 GEOSCIENCES↗

DOE OSTI2015

Code Comparison Study Fosters Confidence in the Numerical Simulation of Enhanced Geothermal Systems

Numerical simulation has become a standard analytical tool for scientists and engineers to evaluate the potential and performance of enhanced geothermal systems. A variety of numerical simulators developed by industry, universities, and national laboratories are currently available and being applied to better understand enhanced geothermal systems at the field scale. To yield credible predictions and be of value to site operators, numerical simulators must be able to accurately represent the complex coupled processes induced by producing geothermal systems, such as fracture aperture changes due to thermal stimulation, fracture shear displacement with fluid injection, rate of thermal depletion of reservoir rocks, and permeability alteration with mineral precipitation or dissolution. A suite of numerical simulators was exercised on a series of test problems that considered coupled thermal, hydraulic, geomechanical, and geochemical (THMC) processes. Problems were selected and designed to isolate selected coupled processes, to be executed on workstation class computers, and have simple but illustrative metrics for result comparisons. This paper summarizes the initial suite of seven benchmark problems, describes the code comparison activities, provides example results for problems and documents the capabilities of currently available numerical simulation codes to represent coupled processes that occur during the production of geothermal resources. Code comparisons described in this paper use the ISO (International Organization for Standardization) standard ISO-13538 for proficiency testing of numerical simulators. This approach was adopted for a recent code comparison study within the radiation transfer-modeling field of atmospheric sciences, which was focused on canopy reflectance models. This standard specifies statistical methods for analyzing laboratory data from proficiency testing schemes to demonstrate that the measurement results do not exhibit evidence of an unacceptable level of bias.

benchmark problems↗

DOE OSTI2013

Modeling Solar Wind Flow with the Multi-Scale Fluid-Kinetic Simulation Suite

Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS) is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. We have enhanced the code with additional physical treatments for the transport of turbulence and acceleration of pickup ions in the interplanetary space and at the termination shock. In this article, we present the results of our numerical simulation of the solar wind (SW) interaction with the local interstellar medium (LISM) in different time-dependent and stationary formulations. Numerical results are compared with the Ulysses, Voyager, and OMNI observations. Finally, the SW boundary conditions are derived from in-situ spacecraft measurements and remote observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Declaration of the Chemical Composition of Waste-products for the German Repository Konrad Using Substance Vectors - 20500

The waste acceptance criteria of the German repository Konrad require beside a radiological declaration a compositional declaration to be able to avoid negative impact on the groundwater by a presumed leakage into the groundwater after complete corrosion of the containers in the future. The declaration of the composition has to specify compounds and substances not element analyses, because the chemical behavior and toxicity especially of organic compounds varies strongly The relevant compounds have to be described and quantified using defined substance codes provided by the operator. For simple waste types like pure concrete or steel components, the definition of the composition is rather simple. The composition of more complex waste streams is generally declared using producer specific material vectors, which can be derived for defined waste types (e.g. grouted fluids, ash from combustion, super-compacted mixed waste etc.). Additional substance vectors exist for containers, drums and grout, which together with the waste composition yield the total composition of the waste packages. The substance vectors for homogeneous waste streams like cementitious demolition waste, solidified fluids and ash are mostly based on analyses provided by the waste owner. In some cases, analyses from non-radioactive analogues or data sheets can also be used to demonstrate the composition. For each qualified treatment of radioactive waste streams the material composition is described and balanced on the basis of specific material vectors, which are composed of material codes from the material and cask list provided by the Federal company for radioactive waste disposal. For mixed waste streams like super-compacted or grouted mixed waste or materials simply packed in containers an analyses cannot be used. In these cases, the compositional information from the bookkeeping can be used for the entire waste stream to calculate a mean composition including a compositional range (min / max composition). Only waste packages with an approved material description can be delivered to Konrad. Therefore, the material declaration of all waste streams and the resulting waste products is crucial for all waste owners in Germany to be able to have enough waste packages ready, when the repository is going to open in 2027. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2001

Adjustable Autonomy Testbed

The Adjustable Autonomy Testbed (AAT) is a simulation-based testbed located in the Intelligent Systems Laboratory in the Automation, Robotics and Simulation Division at NASA Johnson Space Center. The purpose of the testbed is to support evaluation and validation of prototypes of adjustable autonomous agent software for control and fault management for complex systems. The AA T project has developed prototype adjustable autonomous agent software and human interfaces for cooperative fault management. This software builds on current autonomous agent technology by altering the architecture, components and interfaces for effective teamwork between autonomous systems and human experts. Autonomous agents include a planner, flexible executive, low level control and deductive model-based fault isolation. Adjustable autonomy is intended to increase the flexibility and effectiveness of fault management with an autonomous system. The test domain for this work is control of advanced life support systems for habitats for planetary exploration. The CONFIG hybrid discrete event simulation environment provides flexible and dynamically reconfigurable models of the behavior of components and fluids in the life support systems. Both discrete event and continuous (discrete time) simulation are supported, and flows and pressures are computed globally. This provides fast dynamic simulations of interacting hardware systems in closed loops that can be reconfigured during operations scenarios, producing complex cascading effects of operations and failures. Current object-oriented model libraries support modeling of fluid systems, and models have been developed of physico-chemical and biological subsystems for processing advanced life support gases. In FY01, water recovery system models will be developed.

Malin, Jane T.↗

NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗

DOE OSTI2022

Latent Space Simulation for Carbon Capture Design Optimization

The CO2 capture efficiency in solvent-based carbon capture systems (CCSs) critically depends on the gas-solvent interfacial area (IA), making maximization of IA a foundational challenge in CCS design. While the IA associated with a particular CCS design can be estimated via a computational fluid dynamics (CFD) simulation, using CFD to derive the IAs associated with numerous CCS designs is prohibitively costly. However, previous works such as Deep Fluids (Kim et al., 2019) show that large simulation speed-ups and low error are achievable by replacing CFD simulators with neural-network (NN) surrogates that mimic the CFD simulation process. This raises the possibility for a fast, accurate replacement for a CFD simulator, and thus computationally feasible IA-based CCS-design optimization. As such, here, we explore whether an existing NN-surrogate approach (and variants we develop) can successfully be applied to our complex carbon-capture CFD simulations, with the ultimate goal of obtaining a fast and accurate simulator for our CCS-design application. Our experiments build on the Deep Fluids approach and find that resulting surrogates can produce large speed ups (4000x) while maintaining IA relative errors as low as 4% on unseen CCS configurations (interpolating between configurations seen during training). Thus, despite less faithfulness to the underlying physics of the problem, NN surrogates may be a promising tool for our CCS design optimization problem. Notably, though, the Deep Fluids approach has limitations for our application (such as model non-transferability to CCS-packing changes), which we discuss. We conclude with potential directions for future work that may, like innovations we introduced here (e.g., transformer-based dynamics prediction), improve performance on our complicated dataset of CCS CFD simulations.

AI Surrogate Simulation, Carbon Capture, design op↗

DOE OSTI2023

Reducing the Flow Maldistribution in Heat Exchangers through a Novel Polymer Manifold: Numerical Evaluation

The maldistribution of working fluid is one of the issues in heat exchangers that causes a reduction in performance of not only the heat exchanger but also the entire HVAC system. One of the methods to reduce such maldistribution is to improve manifold design to evenly distribute the flow. In the present work, an advanced maldistribution reduction manifold, which was based on a preliminary maldistribution reduction manifold, was designed to further improve the flow distribution in the heat exchanger. In the design, spiral baffles are used to create vortices in the tubes to regulate the flow in each tube. The design also keeps the tubes away from the manifold inlet to avoid direct flow from the inlet. Due to the complexity, the design of the advanced maldistribution reduction manifold is for AM only, which cannot be fabricated by traditional manufacturing. To evaluate the design, a computational fluid dynamic model is developed to study flow distribution in heat exchanger manifolds. The simulation results reveal that the relative standard deviation of the tubes in the advanced maldistribution reduction design is half of the preliminary maldistribution reduction design and about 1/20 of the reference design.

42 ENGINEERING↗

DOE OSTI2025

In situ Raman investigation of Dy complexation in Cl-bearing aqueous solutions at 20–300 °C

Raman spectroscopy provides a versatile tool for in situ characterization of aqueous rare earth elements (REE) speciation at the molecular level. Complexation of REE with ligands such as Cl – and OH – is of particular interest for understanding the mobility of REE in NaCl-bearing hydrothermal fluids responsible for enriching REE to economic levels in nature. Raman spectroscopic studies of REE speciation in Cl-bearing aqueous fluids are primarily conducted at ambient temperature, whereas natural systems indicate temperatures of >100–600 °C. In this study, the speciation of Dy in acidic chloride-bearing hydrothermal solutions was investigated using confocal Raman spectroscopy with a new capillary Raman heating stage at 20–300 °C. Background solutions (pure water, NaCl-solutions) and solutions with 0.14–1.8 mol kg –1 dissolved DyCl 3 were sealed in quartz capillary cells. Comparison of the spectra for Dy chloride solutions with those for background solutions and the spectra for reference Dy-bearing solids was used to identify Raman bands specific to Dy–O and Dy–Cl bonds. The Raman band for the Dy–O stretching mode of hydrated Dy 3+ aqua ions was measured at 365–384 cm –1 and a Raman band for the Dy–Cl stretching modes of Dy chloride complexes was measured near 240 cm –1 . The Dy–O band decreases systematically with temperature, whereas the Dy–Cl band systematically increases, indicating a systematic increase in the stability of Dy chloride complexes with temperature. Here, this study provides the framework for expanding the use of in situ Raman spectroscopy to investigate the speciation of REE in aqueous solutions to hydrothermal conditions.

58 GEOSCIENCES↗

NASA NTRS1995

Domain decomposition for aerodynamic and aeroacoustic analyses, and optimization

The overarching theme was the domain decomposition, which intended to improve the numerical solution technique for the partial differential equations at hand; in the present study, those that governed either the fluid flow, or the aeroacoustic wave propagation, or the sensitivity analysis for a gradient-based optimization. The role of the domain decomposition extended beyond the original impetus of discretizing geometrical complex regions or writing modular software for distributed-hardware computers. It induced function-space decompositions and operator decompositions that offered the valuable property of near independence of operator evaluation tasks. The objectives have gravitated about the extensions and implementations of either the previously developed or concurrently being developed methodologies: (1) aerodynamic sensitivity analysis with domain decomposition (SADD); (2) computational aeroacoustics of cavities; and (3) dynamic, multibody computational fluid dynamics using unstructured meshes.

Baysal, Oktay↗

DOE OSTI2014

Discrete models of fluids: spatial averaging, closure and model reduction

We consider semidiscrete ODE models of single-phase fluids and two-fluid mixtures. In the presence of multiple fine-scale heterogeneities, the size of these ODE systems can be very large. Spatial averaging is then a useful tool for reducing computational complexity of the problem. The averages satisfy exact balance equations of mass, momentum, and energy. These equations do not form a satisfactory continuum model because evaluation of stress and heat flux requires solving the underlying ODEs. To produce continuum equations that can be simulated without resolving microscale dynamics, we recently proposed a closure method based on the use of regularized deconvolution. Here we continue the investigation of deconvolution closure with the long term objective of developing consistent computational upscaling for multiphase particle methods. The structure of the fine-scale particle solvers is reminiscent of molecular dynamics. For this reason we use nonlinear averaging introduced for atomistic systems by Noll, Hardy, and Murdoch-Bedeaux. We also consider a simpler linear averaging originally developed in large eddy simulation of turbulence. We present several simple but representative examples of spatially averaged ODEs, where the closure error can be analyzed. Based on this analysis we suggest a general strategy for reducing the relative error of approximate closure. For problems with periodic highly oscillatory material parameters we propose a spectral boosting technique that augments the standard deconvolution and helps to correctly account for dispersion effects. We also conduct several numerical experiments, one of which is a complete mesoscale simulation of a stratified two-fluid flow in a channel. In this simulation, the operation count per coarse time step scales sublinearly with the number of particles.

Multiscale modeling, discrete models, fluid flow↗

DOE OSTI2017

Pressure evolution equation for the particulate phase in inhomogeneous compressible disperse multiphase flows

An analytical expression describing the unsteady pressure evolution of the dispersed phase driven by variations in the carrier phase is presented. In this article, the term “dispersed phase” represents rigid particles, droplets, or bubbles. Letting both the dispersed and continuous phases be inhomogeneous, unsteady, and compressible, the developed pressure equation describes the particle response and its eventual equilibration with that of the carrier fluid. The study involves impingement of a plane traveling wave of a given frequency and subsequent volume-averaged particle pressure calculation due to a single wave. The ambient or continuous fluid's pressure and density-weighted normal velocity are identified as the source terms governing the particle pressure. Analogous to the generalized Faxén theorem, which is applicable to the particle equation of motion, the pressure expression is also written in terms of the surface average of time-varying incoming flow properties. The surface average allows the current formulation to be generalized for any complex incident flow, including situations where the particle size is comparable to that of the incoming flow. Further, the particle pressure is also found to depend on the dispersed-to-continuous fluid density ratio and speed of sound ratio in addition to dynamic viscosities of both fluids. The model is applied to predict the unsteady pressure variation inside an aluminum particle subjected to normal shock waves. The results are compared against numerical simulations and found to be in good agreement. Furthermore, it is shown that, although the analysis is conducted in the limit of negligible flow Reynolds and Mach numbers, it can be used to compute the density and volume of the dispersed phase to reasonable accuracy. Finally, analogous to the pressure evolution expression, an equation describing the time-dependent particle radius is deduced and is shown to reduce to the Rayleigh-Plesset equation in the linear limit.

42 ENGINEERING↗

NASA NTRS1982

Computational fluid dynamics - The coming revolution

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.

Graves, R. A., Jr.↗

DOE OSTI2020

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗

DOE OSTI2012

Flow Accelerated Erosion-Corrosion (FAC) considerations for secondary side piping in the AP1000{sup R} nuclear power plant design

The issue of Flow Accelerated Erosion-Corrosion (FAC) in power plant piping is a known phenomenon that has resulted in material replacements and plant accidents in operating power plants. Therefore, it is important for FAC resistance to be considered in the design of new nuclear power plants. This paper describes the design considerations related to FAC that were used to develop a safe and robust AP1000{sup R} plant secondary side piping design. The primary FAC influencing factors include: - Fluid Temperature - Pipe Geometry/layout - Fluid Chemistry - Fluid Velocity - Pipe Material Composition - Moisture Content (in steam lines) Due to the unknowns related to the relative impact of the influencing factors and the complexities of the interactions between these factors, it is difficult to accurately predict the expected wear rate in a given piping segment in a new plant. This paper provides: - a description of FAC and the factors that influence the FAC degradation rate, - an assessment of the level of FAC resistance of AP1000{sup R} secondary side system piping, - an explanation of options to increase FAC resistance and associated benefits/cost, - discussion of development of a tool for predicting FAC degradation rate in new nuclear power plants. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2015

On the rogue waves propagation in non-Maxwellian complex space plasmas

The implications of the non-Maxwellian electron distributions (nonthermal/or suprathermal/or nonextensive distributions) are examined on the dust-ion acoustic (DIA) rogue/freak waves in a dusty warm plasma. Using a reductive perturbation technique, the basic set of fluid equations is reduced to a nonlinear Schrödinger equation. The latter is used to study the nonlinear evolution of modulationally unstable DIA wavepackets and to describe the rogue waves (RWs) propagation. Rogue waves are large-amplitude short-lived wave groups, routinely observed in space plasmas. The possible region for the rogue waves to exist is defined precisely for typical parameters of space plasmas. It is shown that the RWs strengthen for decreasing plasma nonthermality and increasing superthermality. For nonextensive electrons, the RWs amplitude exhibits a bit more complex behavior, depending on the entropic index q. Moreover, our numerical results reveal that the RWs exist with all values of the ion-to-electron temperature ratio σ for nonthermal and superthermal distributions and there is no limitation for the freak waves to propagate in both two distributions in the present plasma system. But, for nonextensive electron distribution, the bright- and dark-type waves can propagate in this case, which means that there is a limitation for the existence of freak waves. Our systematic investigation should be useful in understanding the properties of DIA solitary waves that may occur in non-Maxwellian space plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗