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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗
DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
DOE OSTI2014

A penalization technique to model plasma facing components in a tokamak with temperature variations

To properly address turbulent transport in the edge plasma region of a tokamak, it is mandatory to describe the particle and heat outflow on wall components, using an accurate representation of the wall geometry. This is challenging for many plasma transport codes, which use a structured mesh with one coordinate aligned with magnetic surfaces. We propose here a penalization technique that allows modeling of particle and heat transport using such structured mesh, while also accounting for geometrically complex plasma-facing components. Solid obstacles are considered as particle and momentum sinks whereas ionic and electronic temperature gradients are imposed on both sides of the obstacles along the magnetic field direction using delta functions (Dirac). Solutions exhibit plasma velocities (M=1) and temperatures fluxes at the plasma–wall boundaries that match with boundary conditions usually implemented in fluid codes. Grid convergence and error estimates are found to be in agreement with theoretical results obtained for neutral fluid conservation equations. The capability of the penalization technique is illustrated by introducing the non-collisional plasma region expected by the kinetic theory in the immediate vicinity of the interface, that is impossible when considering fluid boundary conditions. Axisymmetric numerical simulations show the efficiency of the method to investigate the large-scale transport at the plasma edge including the separatrix and in realistic complex geometries while keeping a simple structured grid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computational Fluid Dynamics Using the Adaptive Wavelet-Collocation Method

Advancements to the adaptive wavelet-collocation method over the last decade have opened up a number of new possible areas for active research. Volume penalization techniques allow complex immersed boundary conditions to be used with high efficiency for both internal and external flows. Anisotropic methods make it possible to use body-fitted meshes while still taking advantage of the dynamic adaptability properties wavelet-based methods provide. The parallelization of the approach has made it possible to perform large high-resolution simulations of detonation initiation and fluid instabilities to uncover new physical insights that would otherwise be difficult to discover. Other developments include space-time adaptive methods and nonreflecting boundary conditions. This article summarizes the work performed using the adaptive wavelet-collocation method developed by Vasilyev and coworkers over the past decade.

42 ENGINEERING↗
DOE OSTI2020

Smart Proxy Modeling Application of Artificial Intelligence & Machine Learning in Computational Fluid Dynamics

Smart proxy technology leverages the art of artificial intelligence and machine learning in order to build accurate and very fast proxy models for highly complex numerical simulation models. In this project, smart proxy technology is used to replicate the results of a highly complex single phase CFD simulation with reasonable degree of accuracy while reducing the computational cost associated with such CFD simulations. The CFD model under study simulates the combustion of natural gas under various conditions such as varying natural gas composition and flow rate, inlet air flow rate and temperature, and outlet pressure in a High-Pressure Combustion Facility (B6 Combustor) with more than 4 million simulation cells. Only eight CFD simulation runs were used to create a smart proxy model that replicates the detail distribution of Pressure, Temperature, Nitrogen, Oxygen and Carbon-dioxide concentration in the CFD simulation model in seconds with less than 10% error.

03 NATURAL GAS↗
DOE OSTI2026

Additional Data Focused on Phase 1 Geared Toward Computational Fluid Dynamics (CFD) validation

A new validation campaign was developed within the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) to better understand the complex interactions between components of a floating wind system (e.g., columns, pontoons, etc.) in a more generic sense, and to validate modeling approaches for a variety of floating designs. A set of cylinders based on the Offshore Code Comparison Collaboration Continuation (OC4)-DeepCwind semidesign were tested under wave loading, both individually and in different combinations (including the attachment of heave plates of different sizes). The cylinders were held rigid, and pressure sensors were used to measure the distributed loading on the structures, while load cells measured the total integrated hydrodynamic loads.

17 WIND ENERGY↗
NASA NTRS1989

1989 Computational Fluid Dynamics Highlights

This video presents highlights of 1989's CFD graphics, which show shuttle flight problems, F-18 flows, artificial heart, and rotorstrator with more complex blades.

Source record↗
NASA NTRS1975

Forced-flow once-through boilers

A compilation and review of NASA-sponsored research on boilers for use in spacecraft electrical power generation systems is presented. Emphasis is on the heat-transfer and fluid-flow problems. In addition to space applications, much of the boiler technology is applicable to terrestrial and marine uses such as vehicular power, electrical power generation, vapor generation, and heating and cooling. Related research areas are discussed such as condensation, cavitation, line and boiler dynamics, the SNAP-8 project (Mercury-Rankine cycle), and conventional terrestrial boilers (either supercritical or gravity-assisted liquid-vapor separation types). The research effort was directed at developing the technology for once-through compact boilers with high heat fluxes to generate dry vapor stably, without utilizing gravity for phase separations. A background section that discusses, tutorially, the complex aspects of the boiling process is presented. Discussions of tests on alkali metals are interspersed with those on water and other fluids on a phenomenological basis.

Stone, J. R.↗
NASA NTRS1989

Volumetric visualization of 3D data

In recent years, there has been a rapid growth in the ability to obtain detailed data on large complex structures in three dimensions. This development occurred first in the medical field, with CAT (computer aided tomography) scans and now magnetic resonance imaging, and in seismological exploration. With the advances in supercomputing and computational fluid dynamics, and in experimental techniques in fluid dynamics, there is now the ability to produce similar large data fields representing 3D structures and phenomena in these disciplines. These developments have produced a situation in which currently there is access to data which is too complex to be understood using the tools available for data reduction and presentation. Researchers in these areas are becoming limited by their ability to visualize and comprehend the 3D systems they are measuring and simulating.

Russell, Gregory↗
DOE OSTI2014

Anomalous diffusion of field lines and charged particles in Arnold-Beltrami-Childress force-free magnetic fields

The cosmic magnetic fields in regions of low plasma pressure and large currents, such as in interstellar space and gaseous nebulae, are force-free in the sense that the Lorentz force vanishes. The three-dimensional Arnold-Beltrami-Childress (ABC) field is an example of a force-free, helical magnetic field. In fluid dynamics, ABC flows are steady state solutions of the Euler equation. The ABC magnetic field lines exhibit a complex and varied structure that is a mix of regular and chaotic trajectories in phase space. The characteristic features of field line trajectories are illustrated through the phase space distribution of finite-distance and asymptotic-distance Lyapunov exponents. In regions of chaotic trajectories, an ensemble-averaged variance of the distance between field lines reveals anomalous diffusion—in fact, superdiffusion—of the field lines. The motion of charged particles in the force-free ABC magnetic fields is different from the flow of passive scalars in ABC flows. The particles do not necessarily follow the field lines and display a variety of dynamical behavior depending on their energy, and their initial pitch-angle. There is an overlap, in space, of the regions in which the field lines and the particle orbits are chaotic. The time evolution of an ensemble of particles, in such regions, can be divided into three categories. For short times, the motion of the particles is essentially ballistic; the ensemble-averaged, mean square displacement is approximately proportional to t{sup 2}, where t is the time of evolution. The intermediate time region is defined by a decay of the velocity autocorrelation function—this being a measure of the time after which the collective dynamics is independent of the initial conditions. For longer times, the particles undergo superdiffusion—the mean square displacement is proportional to t{sup α}, where α > 1, and is weakly dependent on the energy of the particles. These super-diffusive characteristics, both of magnetic field lines and of particles moving in these fields, strongly suggest that theories of transport in three-dimensional chaotic magnetic fields need a shift from the usual paradigm of quasilinear diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Comprehensive Evaluation of Fast-Response, Reynolds-Averaged Navier–Stokes, and Large-Eddy Simulation Methods Against High-Spatial-Resolution Wind-Tunnel Data in Step-Down Street Canyons

Three computational fluid dynamics (CFD) methods with different levels of flow-physics modelling are comprehensively evaluated against high-spatial-resolution wind-tunnel velocity data from step-down street canyons (i.e., a short building downwind of a tall building). The first method is a semi-empirical fast-response approach using the Quick Urban Industrial Complex (QUIC-URB) model. The second method solves the Reynolds-averaged Navier–Stokes (RANS) equations, and the third one utilizes a fully-coupled fluid-structure interaction large-eddy simulation (LES) model with a grid-turbulence inflow generator. Unlike typical point-by-point evaluation comparisons, here the entire two-dimensional wind-tunnel dataset is used to evaluate the dynamics of dominant flow topological features in the street canyon. Each CFD method is scrutinized for several geometric configurations by varying the downwind-to-upwind building-height ratio (H d /H u ) and street canyon-width to building-width aspect ratio (S / W) for inflow winds perpendicular to the upwind building front face. Disparities between the numerical results and experimental data are quantified in terms of their ability to capture flow topological features for different geometric configurations. Ultimately, all three methods qualitatively predict the primary flow topological features, including a saddle point and a primary vortex. But, the secondary flow topological features, namely an in-canyon separation point and secondary vortices, are only well represented by the LES method despite its failure for taller downwind building cases. Misrepresentation of flow-regime transitions, exaggeration of the coherence of recirculation zones and wake fields, and overestimation of downwards vertical velocity into the canyon are the main defects in QUIC-URB, RANS and LES results, respectively. All three methods underestimate the updrafts and, surprisingly, QUIC-URB outperforms RANS for the streamwise velocity component, while RANS is superior to QUIC-URB for the vertical velocity component in the street canyon.

58 GEOSCIENCES↗
DOE OSTI2023

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Fluid–Solid Reaction Kinetics for Solids of Nonbasic Geometries and Determination of the Appropriate Shape Factors

The rate analysis involving a solid of nonbasic geometry in general requires a complex numerical solution. A closed-form expression of the reaction rate for such a solid would be very convenient and greatly simplify the analysis. A method for this purpose was developed in this work for solids with shapes of finite cylinders and finite rectangular parallelepipeds reacting under the rate control by chemical reaction or mass transfer, which may be extended to other irregular shapes. In this method, an appropriate shape factor is first obtained under the conditions of a topochemical reaction, which is then applied to the unified rate expressions previously obtained for solids of three basic shapes. Good agreements between the rates predicted with the appropriate shape factors and the exact or numerical solutions were obtained for both reaction-controlled and pore-diffusion-controlled reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Altering wettability to recover more oil from tight formations

We describe here a method for chemically modifying fracturing fluids and overflushes to chemically increase oil recovery from tight formations. Oil wetting of tight formations is usually controlled by adhesion to illite, kerogen, or both; adhesion to carbonate minerals may also play a role. Oil-illite adhesion is sensitive to salinity, dissolved divalent cation content, and pH. We measure oil-rock adhesion with middle Bakken formation oil and core to verify a surface complexation model of reservoir wettability. The agreement between the model and experiments suggests that wettability trends in tight formations can be quantitatively predicted and that fracturing fluid and overflush compositions can be individually tailored to increase oil recovery.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Andrieslombaardite, RhSbS, a new platinum-group mineral from the platiniferous Onverwacht Pipe, Republic of South Africa

A hundred years after the discovery of the Merensky Reef in 1924, it is appropriate to present the new mineral andrieslombaardite in honour of Andries Frederik Lombaard who was instrumental in its discovery. Andrieslombaardite, RhSbS, was first described as an unknown mineral from placer deposits associated with the Tulameen Alaskan-Uralian type complex, British Colombia, Canada (Raicevic and Cabri, 1976) but has since been reported from several other deposits including the platiniferous Driekop, Mooihoek, and Onverwacht pipes in the eastern Bushveld Complex, South Africa. The mineral and the name were approved by the Commission on New Minerals Nomenclature and Classification (CNMNC) of the International Mineralogical Association (IMA no. 2022-076) based on data in the co-type samples from Onverwacht and a co-type sample from the Yubdo stream, Birbir River, Ethiopia. Andrieslombaardite in the Onverwacht sample is a single 8 x 20 μm grain attached to laurite in a matrix of altered silicate and Fe-oxyhydroxide minerals. In the Yubdo samples, there are many grains of pale brownish gray andrieslombaardite up to 25 x 55 μm in size, included in Pt-Fe alloys, some associated with erlichmanite, and others attached to bornite and chalcopyrite. The reflectance values (R%) measured in air and in oil at the COM wavelengths are 48.3 and 33.0 (470 nm), 49.3 and 34.0 (546 nm), 51.0 and 35.9 (589 nm), and 51.8 and 36.7 (650 nm). The colour values x, y, Y, λd, and Pe in air are 0.317, 0.322, 50.3, 580, and 3.2, and in oil are 0.319, 0.324, 35.6, 579, and 4.5. The composition of andrieslombaardite is ideally RhSbS, but it contains variable amounts of Fe, Pt, Pd, and Ir that may substitute for Rh. The mineral is cubic with unit-cell dimensions of a = 6.0278(4) Å, V = 219.01(6) Å 3 and Z = 4. It was synthesised at 400 and 550°C using stoichiometric elemental amounts. It is a member of the cobaltite group. The mineralisation of the intrusive dunite pipes was probably introduced at high temperatures, under magmatic conditions. The primary assemblages were to a certain degree overprinted and redistributed by low-temperature hydrothermal fluids. The Pt-Fe alloys from Yubdo containing PGM inclusions such as andrieslombaardite in the Yubdo-Alaskan-type complex were formed at some post-magmatic stage owing to PGE remobilisation during hydrothermal or metamorphic episodes.

58 GEOSCIENCES↗
DOE OSTI2021

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗
DOE OSTI2016

MATLAB/Simulink Framework for Modeling Complex Coolant Flow Configurations of Advanced Automotive Thermal Management Systems

The National Renewable Energy Laboratory's (NREL's) CoolSim MATLAB/Simulink modeling framework was extended by including a newly developed coolant loop solution method aimed at reducing the simulation effort for arbitrarily complex thermal management systems. The new approach does not require the user to identify specific coolant loops and their flow. The user only needs to connect the fluid network elements in a manner consistent with the desired schematic. Using the new solution method, a model of NREL's advanced combined coolant loop system for electric vehicles was created that reflected the test system architecture. This system was built using components provided by the MAHLE Group and included both air conditioning and heat pump modes. Validation with test bench data and verification with the previous solution method were performed for 10 operating points spanning a range of ambient temperatures between -2 degrees C and 43 degrees C. The largest root mean square difference between pressure, temperature, energy and mass flow rate data and simulation results was less than 7%.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2016

All Nuclear Hybrid Superheater Design Method

The 'All Nuclear Hybrid' (ANH) concept combines Light Water (LWR) and Generation IV (Gen IV) reactor technologies to produce superheated steam. LWR technology characteristically produces steam temperatures of ∼300 C. Process steam demands such as refineries and hydrogen production can require steam temperatures that exceed 600 C. Therefore, an additional heat source must be added to provide the heat to increase or 'boost' the temperature of the steam generated from by the LWR. A recent study demonstrated that application of the ANH concept to a refinery optimized the use of nuclear technologies and maximized carbon dioxide emission reductions. The objective of the present work is to size the heat exchangers (Superheaters) required to apply the ANH concept to a Westinghouse AP1000 LWR to produce 675 C superheated steam. Two GenIV reactor concepts are considered to provide the additional heat source and high primary coolant temperatures required to produce the superheated steam: a High temperature Gas Reactor (HTGR) using helium as primary coolant, and a Molten Salt reactor. The details of the GenIV reactor designs are not the focus of this study; however, the primary coolant thermal conditions from these concepts are utilized directly in the superheater sizing calculations. Also, a tube and shell superheater design is assumed for the superheater due to its extensive use in the nuclear and other industries. The ANH concept generates superheated steam by first utilizing LWR technology to accomplish the phase change from liquid water to steam within the steam generator. This is a well-established technology with several decades of successful safe operation. Next, the concept utilizes GenIV reactor technology to superheat the steam. Since both the working fluids in the superheater are single phase, the heat transfer and pressure drop calculations required to size the superheater can be performed without complex two-phase computer analyses (i.e., RELAP5). This study presents an approach for the sizing the ANH superheater using an iterative LMTD heat exchanger design method. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗