Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Fluid Dynamics Uncertainty Quantification Challenge Problem: XFOIL vs. MFOIL

Uncertainty quantification (UQ) has become more critical in aerospace engineering due to the growing dependence on computational tools for design optimization and performance analyses of aerospace vehicles. Even though the significance of UQ in assessing the credibility of computational analyses is well recognized, its costs and complexity impede its integration into standard practices, particularly in computational fluid dynamics (CFD) and other fluid analyses. This paper presents a UQ study for low-fidelity computational aerodynamics analyses with XFOIL and mfoil (i.e., the MATLAB version of XFOIL with several implementation modifications); these tools are utilized widely in both research and education. The main contributions of this paper are as follows: 1) improved precision in quantifying the uncertainty of the baseline Monte Carlo results used to benchmark surrogate modeling techniques for UQ, 2) quantification of the effect of the implementation differences between XFOIL and mfoil on solution quantities of interest (QoIs), such as lift and pitching moment coefficients, and 3) development of an open-source UQ library for use with XFOIL and mfoil, which has educational values and helps promote UQ for fluid analyses with aerospace applications. Results and discussions revolve around cases 1-4 of the challenge problem posed by the AIAA Fluid Dynamics Technical Committee’s Uncertainty Quantification Discussion Group (UQDG). In case 3, this work employs CFDverify, an open-source solution verification software, to quantify the discretization error and evaluate the extrapolated QoIs based on the grid convergence index (GCI). This UQ study differentiates itself from previous studies in the rigor of handling baseline Monte Carlo uncertainty and in including mfoil, which is a more accessible alternative to XFOIL. Finally, despite the growing computing power, low-fidelity computational tools remain valuable, such as for aerodynamic shape optimization at Mach numbers below 0.65 and low-to-mid Reynolds numbers.

Lay, Aidan S [University of Tennessee, Knoxville (↗
NASA NTRS1987

Velocimetry with refractive index matching for complex flow configurations, phase 1

The feasibility of obtaining detailed velocity field measurements in large Reynolds number flow of the Space Shuttle Main Engine (SSME) main injector bowl was demonstrated using laser velocimetry and the developed refractive-index-matching technique. An experimental system to provide appropriate flow rates and temperature control of refractive-index-matching fluid was designed and tested. Test results are presented to establish the feasibility of obtaining accurate velocity measurements that map the entire field including the flow through the LOX post bundles: sample mean velocity, turbulence intensity, and spectral results are presented. The results indicate that a suitable fluid and control system is feasible for the representation of complex rocket-engine configurations and that measurements of velocity characteristics can be obtained without the optical access restrictions normally associated with laser velocimetry. The refractive-index-matching technique considered needs to be further developed and extended to represent other rocket-engine flows where current methods either cannot measure with adequate accuracy or they fail.

Thompson, B. E.↗
DOE OSTI2018

An exploration of pressure dynamics using differential equations defined on a fractal geometry

Pressure tests have been used to understand the behavior of naturally fractured reservoirs (NFRs) by establishing parameters such as permeability, connectivity of fractures, system compressibility, fracture network parameters, and others associated with the physical properties of the reservoir. The complex fracture network of a NFR produces anomalous fluid flow, which is shown by the pressure behavior at the well. A plausible model to explain this behavior is to assume that it arises from a fractal fracture system. However, to estimate the parameters that define this fractal fracture system, fluid flow models based on the modification of Darcy’s law are used. On the other hand, in this work, a fractal structure with known fractal dimension: Sierpinski gasket (SG) is proposed. The objective of this work is to solve numerically the diffusion equation on SG, and establish a comparison with flow models used in the oil industry to investigate the fractality assumption. This study is focused on the qualitative behavior of the solution of the diffusion equation defined on a SG supported by the Kigami’s theory of differential equations on fractals.

42 ENGINEERING↗
DOE OSTI2022

Chemical and Reactive Transport Processes Associated with Hydraulic Fracturing of Unconventional Oil/Gas Shales

Hydraulic fracturing of unconventional oil/gas shales has changed the energy landscape of the U.S.. Recovery of hydrocarbons from tight, hydraulically fractured shales is a highly inefficient process, with estimated recoveries of < 25% for natural gas and < 5% for oil. This review focuses on the complex chemical interactions of additives in hydraulic fracturing fluid (HFF) with minerals and organic matter in oil/gas shales. These interactions are intended to increase hydrocarbon recovery by increasing porosities and permeabilities of tight shales. However, fluid-shale interactions result in the dissolution of shale minerals and the release and transport of chemical components. They also result in mineral precipitation in the shale matrix, which can reduce permeability, porosity, and hydrocarbon recovery. Competition between mineral dissolution and mineral precipitation processes influences the amounts of oil and gas recovered. We review the temporal/spatial origins and distribution of unconventional oil/gas shales from mudstones and shales, followed by discussion of their global and U.S. distributions and compositional differences from different U.S. sedimentary basins. We discuss the major types of chemical additives in HFF with their intended purposes, including drilling muds. Fracture distribution, porosity, permeability, and the identity and molecular-level speciation of minerals and organic matter in oil/gas shales throughout the hydraulic fracturing process are discussed. Also discussed are analysis methods used in characterizing oil/gas shales before and after hydraulic fracturing, including permeametry and porosimetry measurements, X-ray diffraction/Rietveld refinement, X-ray computed tomography, scanning/transmission electron microscopy, and laboratory- and synchrotron-based imaging/spectroscopic methods. Reactive transport and spatial scaling are discussed in some detail in order to relate fundamental molecular-scale processes to fluid transport. Furthermore, our review concludes with a discussion of potential environmental impacts of hydraulic fracturing and important knowledge gaps that must be bridged to achieve improved mechanistic understanding of fluid transport in oil/gas shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2016

Derivation of the Hall and extended magnetohydrodynamics brackets

There are several plasma models intermediate in complexity between ideal magnetohydrodynamics (MHD) and two-fluid theory, with Hall and Extended MHD being two important examples. In this paper, we investigate several aspects of these theories, with the ultimate goal of deriving the noncanonical Poisson brackets used in their Hamiltonian formulations. We present fully Lagrangian actions for each, as opposed to the fully Eulerian, or mixed Eulerian-Lagrangian, actions that have appeared previously. As an important step in this process, we exhibit each theory's two advected fluxes (in analogy to ideal MHD's advected magnetic flux), discovering also that with the correct choice of gauge they have corresponding Lie-dragged potentials resembling the electromagnetic vector potential, and associated conserved helicities. Finally, using the Euler-Lagrange map, we show how to derive the noncanonical Eulerian brackets from canonical Lagrangian ones.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Breakup of a liquid rivulet falling over an inclined plate: Identification of a critical Weber number

Here, we have numerically investigated the breakup of a rivulet falling over a smooth inclined plate using the volume of fluid method. The breakup the rivulet is a complex phenomenon that is dictated by many factors, such as solvent properties, contact angle, inertia, plate inclination, etc. An extensive simulation campaign was conducted wherein these factors were systematically investigated. Regimes for a stable rivulet and an unstable rivulet that leads to the breakup and formation of a droplet are examined in terms of a critical value of the Weber number (We cr ) that delineates these regimes. The effect of plate inclination on the breakup of the rivulet shows that the critical Weber number decreases with increased inclination angle (θ) owing to higher liquid velocity. However, the effect is negligible beyond θ > 60°. The impact of solvent properties is characterized using the Kapitza number (Ka). Variation of We cr with Ka shows two trends depending on the Ka value of the solvent. Solvents with lower Ka values, corresponding to high viscosities and/or low surface tensions, show smaller values of the critical Weber number and the variation is linear. While solvents with higher Ka values exhibit higher values of the We cr and the variation in We cr is steep. This behavior is more pronounced with increasing contact angle. Higher contact angles promote rivulet breakup so that inertia must be higher to the breakup. As a result, a phenomenological scaling for a critical Weber number with the Kapitza number and contact angle is presented that can offer insight into rivulet breakup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

An experimental study of the existence regions and non-linear interactions of drift wave and Kelvin–Helmholtz instabilities in a linear magnetized plasma

Experimental observations of the intrinsic excitation and non-linear interactions of drift wave (DW) and Kelvin–Helmholtz (KH) instabilities in a linear magnetized plasma column are presented. The experiments are carried out in the inverse mirror plasma experimental device (IMPED)—a cylindrical, magnetized, linear plasma machine designed to study low-frequency waves and instabilities in plasma. A novel feature of IMPED is the ability to control plasma profiles, such as the density n(r)⁠, electron temperature T e (r)⁠, and plasma potential V p (r) by varying the ratio Rm of the magnetic field in the main chamber to that in the source chamber. At high values of Rm, higher-density gradient scale length promotes the drift wave (DW) instability while lower Rm value results in a higher radial electric field, inducing a sheared poloidal flow that enhances the dominance of the Kelvin–Helmholtz (KH) mode. The background and fluctuating plasma parameters are characterized using various configurations of multiple in situ electric probes at different spatial locations to quantify the local gradients that excite the low-frequency primary instabilities. Statistical, spectral, and bispectral analysis of the density and potential signals help identify these modes in terms of wave number, frequency, phase, and amplitude and also delineate the nature of their non-linear interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

ROM-Based Surrogate Systems Modeling of EBR-II

We report the System Analysis Module (SAM), developed and maintained by Argonne National Laboratory, is designed to provide whole-plant transient safety analysis capabilities for a number of advanced non-light water reactors, including sodium-cooled fast reactor (SFR), lead-cooled fast reactor (LFR), and molten salt reactor (MSR)/fluoride-salt-cooled high-temperature reactor (FHR) designs. SAM is primarily constructed as a systems-level analysis tool, with the potential to incorporate reduced order models from three-dimensional computational fluid dynamics (CFD) simulations to improve characterization of complex, multidimensional physics. It is recognized that the computational expense associated with CFD can be intractable for various engineering analyses, such as uncertainty quantification, inference, and design optimization. This paper explores the reducibility of a SAM model using recent advances in randomized linear algebra techniques, which attempt to find recurring patterns in the various realizations generated by a model after randomly perturbing all its input parameters. The reduction is described in terms of fewer degrees of freedom (DOFs), referred to as the active DOFs, for the model variables such as input model parameters and model responses. The results indicate that there is significant room for additional reduction that may be leveraged for additional computational gains when employing SAM for engineering-intensive analyses that require repeated model executions. Different from physics-based reduction approaches, the proposed approach allows one to estimate upper bounds on the reduction errors, which are rigorously developed in this work. Finally, different methods for surrogate model construction, such as regression and neural network-based training, are employed to correlate the input and output active DOFs, which are related back to the original variables using matrix-based linear transformations.

42 ENGINEERING↗
DOE OSTI2016

Plasma dark matter direct detection

Dark matter in spiral galaxies like the Milky Way may take the form of a dark plasma. Hidden sector dark matter charged under an unbroken U(1)' gauge interaction provides a simple and well defined particle physics model realising this possibility. The assumed U(1)' neutrality of the Universe then implies (at least) two oppositely charged dark matter components with self-interactions mediated via a massless 'dark photon' (the U(1)' gauge boson). In addition to nuclear recoils such dark matter can give rise to keV electron recoils in direct detection experiments. In this context, the detailed physical properties of the dark matter plasma interacting with the Earth is required. This is a complex system, which is here modelled as a fluid governed by the magnetohydrodynamic equations. These equations are numerically solved for some illustrative examples, and implications for direct detection experiments discussed. In particular, the analysis presented here leaves open the intriguing possibility that the DAMA annual modulation signal is due primarily to electron recoils (or even a combination of electron recoils and nuclear recoils). The importance of diurnal modulation (in addition to annual modulation) as a means of probing this kind of dark matter is also emphasised.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Technical Report on NETL's Non Newtonian Multiphase Slurry Workshop: A path forward to understanding non-Newtonian multiphase slurry flows

The Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) sponsored a workshop on non-Newtonian multiphase slurry at NETL’s Morgantown campus August 19 and 20, 2013. The objective of this special two-day meeting of 20-30 invited experts from industry, National Labs and academia was to identify and address technical issues associated with handling non-Newtonian multiphase slurries across various facilities managed by DOE. Particular emphasis during this workshop was placed on applications managed by the Office of Environmental Management (EM). The workshop was preceded by two webinars wherein personnel from ORP and NETL provided background information on the Hanford WTP project and discussed the critical design challenges facing this project. In non-Newtonian fluids, viscosity is not constant and exhibits a complex dependence on applied shear stress or deformation. Many applications under EM’s tank farm mission involve non-Newtonian slurries that are multiphase in nature; tank farm storage and handling, slurry transport, and mixing all involve multiphase flow dynamics, which require an improved understanding of the mechanisms responsible for rheological changes in non-Newtonian multiphase slurries (NNMS). To discuss the issues in predicting the behavior of NNMS, the workshop focused on two topic areas: (1) State-of-the-art in non-Newtonian Multiphase Slurry Flow, and (2) Scaling up with Confidence and Ensuring Safe and Reliable Long-Term Operation.

42 ENGINEERING↗
DOE OSTI2022

PDC Modifications for Analysis of Gas-Cooled Reactors with Direct Helium Brayton Cycle

Capabilities of the Plant Dynamics Code (PDC) have been extended to allow steady-state and transient simulation of graphite-cooled reactors with direct helium Brayton cycle. On the cycle side, the most significant code modification is the addition of helium properties, in the format required by the code’s equations. Since the code was already formulated to handle more complex real gas properties, adding helium as working fluid that behaves like ideal gas was fairly straightforward. A reactor module was added to PDC to simulate a reactor cooled by the working fluid of the Brayton cycle. Two options are supported: channel type, typical for graphite gas-cooled reactors, and pin type, typical for light-water and liquid metal-cooled reactors. The reactor module is an extension of the electrical heater model and simulates heat deposition in the fuel and transfer of this heat from the fuel to the coolant through the matrix and tube materials. The new reactor module becomes the third option in PDC for modeling heat addition to the cycle, besides previously modeled heat addition heat exchanger and electrical heater. In addition to those changes, other minor code modifications and improvements were introduced during the work of expanding PDC to modeling of gas-cooled reactors. These modifications are summarized in the last chapter of this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Rapid CFD Using Machine Learning Algorithms (CRADA Final Report)

This is a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (“LLNL”) and Guardian Glass, LLC ("Guardian Glass") to develop a fast-running emulator of the reactive Computational Fluid Dynamics (“CFD”) simulations needed to understand the complex reactions and flows in the glass melting, fining, and forming subprocesses. This CRADA project is sponsored under the High-Performance Computing for Manufacturing (“HPC4Mfg”) Program of the Department of Energy’s Advanced Manufacturing Office (“AMO”) within the Energy Efficiency and Renewable Energy (“EERE”) Office.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SULI Report - Development of a Molten Salt Circulation Loop for in-situ Spectroscopy

This project supports the development of real time optical monitoring capabilities for molten salt reactor (MSR) environments by designing, testing, and refining a molten salt circulation loop suitable for combined laser induced breakdown spectroscopy (LIBS) and ultraviolet visible (UV Vis) absorption measurements. Online spectroscopic monitoring is increasingly important for nuclear safeguards, corrosion tracking, and material accountancy, yet MSR process fluids present substantial challenges due to their chemical complexity and hazards such as high temperatures and radiation. To address these needs, this work focuses on Phase I, the development of a room temperature aqueous circulation loop that serves as a surrogate platform for evaluating flow behavior, optical access, and component performance prior to high temperature salt operation. Initial testing identified several practical issues—including leaks, obstructions, and two-phase flow through the absorption cell—that were systematically resolved through hardware replacement, flow path redesign, and venturi pressure optimization. Relocating the flow cell upstream of the primary venturi enabled periods of stable single-phase flow, demonstrating the feasibility of integrating optical diagnostics into a circulation system. The results of Phase I provide essential design insight for Phase II, which will incorporate furnace compatible materials and LiCl KCl eutectic salt. Completion of the molten salt system will deliver a reusable testbed for evaluating multimodal spectroscopic techniques, advancing nondestructive, real time monitoring tools for future MSR and nuclear fuel cycle applications.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗
NASA NTRS1982

A basic code for the prediction of transient three-dimensional turbulent flowfields

A primitive pressure-velocity finite difference code has been developed to predict transient three-dimensional turbulent flow. The code is a simplified yet effective prediction procedure for use by persons with little experience in computational fluid dynamics, and into which user-oriented complexities can easily be added. The method is based on the transient two-dimensional Los Alamos SOLA prediction technique for laminar flows. Turbulence is simulated by means of a two-equation turbulence model; species diffusion and buoyancy are also included. Two applications of the code are presented to local destratification near the release structure of a reservoir and to the deflection of a jet entering normally into a uniform cross-flow. Predicted results exhibit good agreement with experimental data, showing that a useful characterization of fully three-dimensional flows is now available.

Busnaina, A. A.↗
NASA NTRS1987

Unsteady three-dimensional Navier-Stokes simulations of turbine rotor-stator interaction

Fluid flows within turbomachinery tend to be extremely complex. Understanding such flows is crucial to efforts to improve current turbomachinery designs, and the computational approach can be used to great advantage in this regard. This study presents a finite-difference, unsteady, thin-layer Navier-Stokes approach to calculating the flow within an axial turbine stage. The relative motion between the stator and rotor airfoils is made possible with the use of patched grids that move relative to each other. The calculation includes end-wall and tip-leakage effects. Results in the form of time-averaged surface pressures, pressure amplitudes (corresponding to the pressure fluctuation in time), near-surface velocity vectors, and pressure contours in the passage areas are presented. The numerical results are compared with experimental data wherever possible, and the agreement between the two is found to be good.

Rai, Man Mohan↗
NASA NTRS1988

Three-dimensional Navier-Stokes simulations of turbine rotor-stator interaction

Fluid flows within turbomachinery tend to be extremely complex in nature. Understanding such flows is crucial to improving current designs of turbomachinery. The computational approach can be used to great advantage in understanding flows in turbomachinery. A finite difference, unsteady, thin layer, Navier-Stokes approach to calculating the flow within an axial turbine stage is presented. The relative motion between the stator and rotor airfoils is made possible with the use of patched grids that move relative to each other. The calculation includes endwall and tip leakage effects. An introduction to the rotor-stator problem and sample results in the form of time averaged surface pressures are presented. The numerical data are compared with experimental data and the agreement between the two is found to be good.

Rai, Man Mohan↗
NASA NTRS1988

Adaptive remeshing for transient problems with moving bodies

The combination of adaptive remeshing techniques, flow solvers for transient problems with moving grids, and integrators for rigid body motion is presented. The resulting scheme allows the economic simulation of fully coupled fluid-rigid body interaction problems of arbitrary geometric complexity. Several results are given to demonstrate the capabilities developed.

Lohner, Rainald↗