Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A generalized wind turbine cross section as a reduced-order model to gain insights in blade aeroelastic challenges

In this work, we present an approach to study the aeroelastic stability of a wind turbine by focusing on the dynamics of a blade cross section. We present a methodology to obtain a reduced-order model of the blade dynamics in the form of generalized cross-sectional quantities that approximates the aerodynamic and structural properties of the full blade. The motivation for the work is to gain a physical understanding of the influence of aerodynamic models such as dynamic wake and dynamic stall on the frequency and damping of the structure using a reduced-order model with low computational cost. The model may be coupled to two-dimensional computational fluid dynamics softwares or engineering unsteady airfoil aerodynamics models accounting for dynamic wake and dynamic stall. In the latter case, we can obtain monolithic state-space forms of the aeroelastic system of equations, which simplifies the determination of the modal parameters and therefore the study of stability. The work investigates wind turbines in operation or at standstill, where vortex-induced vibrations and stall-induced vibrations, respectively, might be an issue. The implementation is made available as part of the open-source Python package WELIB and as part of the open-source unsteady aerodynamic driver of OpenFAST.

17 WIND ENERGY↗
DOE OSTI2024

Workflow for Developing and Operating Subsurface Hydrogen Storage Facilities in Porous Reservoirs

Long-duration (seasonal) storage of natural gas (NG), which primarily consists of methane (CH 4 ), has been practiced for more than a hundred years at underground gas storage (UGS) facilities that use depleted hydrocarbon reservoirs, saline aquifers, and salt caverns. To enable hydrogen (H 2 ) to be used as a long-duration, energy-storage medium, similar facilities are envisioned for underground H 2 storage (UHS) of either H 2 or H 2 /NG mixtures. Experience with UGS can be used to guide recommended practices for developing and operating UHS facilities in porous reservoirs. The most important factors (formation/fluid properties and engineering choices) that influence the performance of UHS reservoirs have been identified and quantified in previous studies. These factors and choices influence phenomena that determine the sweep efficiency of the stored working gas. These phenomena include viscous fingering, hysteretic capillary trapping, and gravity override of the working gas, as well as the upconing of nonproductive fluid that determine the sweep efficiency of the stored working gas. This report describes initial recommended-practices and a project-development workflow for UHS facilities that utilize porous reservoirs, based on the current state-of-knowledge about H 2 behavior in the subsurface. The workflow sequentially addresses all aspects of UHS project development, including the identification of H 2 sources and users, site ranking and down-selection, geologic and reservoir-engineering characterization, reservoir design, testing, risk management, commissioning, operations, and monitoring for a UHS facility. The goal is to enable UHS facilities to be developed in an efficient and timely manner, while carefully managing project risks. This workflow is similar to that which has been developed for UGS facilities (see Figure 1 of API, 2022), with the addition of tasks and subtasks specific to H 2 and UHS. The project-development workflow is broken down into three major stages: (1) define the H 2 use case; (2) rank, down-select, and characterize potential, candidate UHS sites; and (3) reservoir design, integrity testing, risk assessment, commissioning, operations, and monitoring for selected UHS sites. Each major stage is further broken down into tasks and subtasks, which are described at a high level. This report also provides more detailed descriptions of all tasks and subtasks that involve reservoir analysis and testing.

08 HYDROGEN↗
DOE OSTI2025

Beam Shaping and Oscillation for Metals Additive Manufacturing

Additive Manufacturing (AM) is pivotal for next-generation industrial manufacturing but is hampered by process control challenges leading to unwanted defects and sub-optimal material properties. Our research builds on prior work that utilized beam oscillation techniques to mitigate these issues (Wu et al., 2023). We focus on investigating the influence of beam oscillation frequency on melt pool dynamics, utilizing in-situ synchrotron material characterization from the ESRF (Chen et al., 2026). Preliminary findings indicated that beam oscillation results in unique fluid flow patterns and enhanced melting stability. These findings suggested the possibility of optimized AM process control, including minimized defects and tailored thermal histories. We proposed further research to explore the effects of beam oscillation frequencies comprehensively, using the in-situ characterization capabilities at the ESRF. The goal was to quantify this relationship, offering a pathway to more precise AM process control and optimization. Our research is expected to offer new scientific insights into the fluid mechanics and thermal behaviour in AM, contributing to its broader applicability and commercial viability.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗
DOE OSTI2020

Transient flow modeling in fractured media using graphs

In this work, we describe a method to simulate transient fluid flows in fractured media using an approach based on graph theory. Our approach builds on past work where the graph-based approach was successfully used to simulate steady-state fluid flows in fractured media. We find a mean computational speedup of the order of 1400 from an ensemble of a 100 discrete fracture networks in contrast to the O(10 4 ) speedup that was obtained for steady-state flows earlier. However, the transient flows considered here involve an additional degree of complexity that was not present in the steady-state flows considered previously with a graph-based approach, that of time marching and solution of the flow equations within a time-stepping scheme. We verify our method with an analytical test case and demonstrate its use on a practical problem related to fluid flows in hydraulically fractured reservoirs. By enabling the study of transient flows, we create an opportunity for a wide set of possibilities where a steady-state approximation is not sufficient, such as the example motivated by hydraulic fracturing that we present here. This work validates the concept that graphs are able to reliably capture the topological properties of the fracture network and serve as effective surrogates in an uncertainty-quantification framework.

58 GEOSCIENCES↗
DOE OSTI2020

Modeling and simulation of gas separations with spiral–wound membranes

Models for gas separations with spiral–wound membranes are developed and found to exhibit good agreement with experiments performed on N 2 /O 2 mixtures. The two–dimensional (2D) model can be accurately approximated by a one–dimensional (1D) surrogate model when the spacer widths are chosen to make the channel pressure drops small. Subsequently, the separation of propane/propylene mixtures from the recycle purge stream of a polypropylene reactor is investigated. Assuming ideal gas is found to lead to significant overestimations in membrane stage cuts (sometimes more than 10%), an extent comparable to that associated with extrapolating constant olefin permeance from a low–pressure condition. While olefin permeance can change significantly with pressure, using a constant–permeance formulation can result in a small (< 2.5%) underprediction in stage cut if the value for the permeance is taken from the feed condition. Finally, membrane properties and costs necessary for a viable separation process are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Mechanically-, Thermally, and Chemically-Robust High-Temperature Ceramic Composites

In an effort to lower the levelized cost of electricity generated by Concentrated Solar Power (CSP) plants, the heat-to-electricity conversion efficiency may be increased by operating such plants with higher turbine inlet temperatures using closed Brayton cycles with high-pressure, supercritical carbon dioxide (sCO 2 ) as the working fluid. For example, if the turbine inlet temperature is raised to ≥750°C, in combination with sCO2 closed Brayton cycles, the relative heat-to-electricity efficiency may be raised by more than 20% (compared to subcritical steam Rankine cycles operating with turbine inlet temperatures ≤550°C). The associated reduction in the cost of dispatchable electricity from such CSP plants (utilizing high-temperature thermal energy storage) would be an important breakthrough on the path towards direct competition with fossil-fuel-based plants. A key barrier to achieve such high-temperature efficient operation has been the limited thermomechanical performance of metal alloys used in compact, primary printed circuit- type heat exchangers (PCHEXs) for heat transfer to high-pressure sCO 2 . The maximum allowed stresses for the use of conventional stainless steels and nickel-based superalloys at high sCO 2 pressures (≥20 MPa) decline rapidly at temperatures ≥550°C. This project has been focused on demonstrating the desired high-temperature properties and attractive manufacturing characteristics of mechanically-, thermally-, and chemically-robust ceramic/metal composites (cermets), in order to allow for the use of such composites in primary PCHEXs for heat transfer to sCO 2 at ≥750°C and ≥20 MPa in CSP plants. Prior work with high-temperature, co-continuous carbide/refractory metal cermets, such as ZrC/W composites, has already demonstrated: i) the attractive combination of high-temperature properties exhibited by such materials (e.g., high values of stiffness, failure strength, and thermal conductivity at 800°C), and ii) cost-effective processes for manufacturing PCHEX plates comprised of such cermets with tailorable channel patterns. However, ZrC/W cermets were not found to be inherently corrosion resistant in sCO 2 at 750°C and 20 MPa. While the application of a Cu coating to ZrC/W surfaces, along with a modest addition (50 ppm) of CO to CO 2 , rendered such composites resistant to corrosion for 1000 h in such sCO 2 -based fluids at 750°C, such additional steps would increase the design complexity of such cermet-based components. An alternative approach (the focus of this project) is to develop mechanically-robust, cost- effective cermets that are inherently resistant to high-temperature oxidation in air and CO 2 . The overall goals of this project have been to demonstrate: i) the oxidation resistance of at least one cermet at 750°C in air and in CO 2 (i.e., with a projected annual corrosion ≤30 μm), ii) the mechanically robust nature of at least one cermet at 750°C (i.e., with an average failure strength ≥200 MPa), and iii) that the selected cermet can be manufactured via low-cost forming techniques.

14 SOLAR ENERGY↗
DOE OSTI2024

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗
DOE OSTI2021

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Measuring Fault Zone and Host Rock Hydraulic Properties Using Tidal Responses

Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.

58 GEOSCIENCES↗
DOE OSTI2025

Material extrusion additive manufacturing of wood pulp-reinforced epoxy composites

Direct ink writing (DIW) is an extrusion-based form of 3D-printing that has gained popularity over the last decade. DIW uses thixotropic fluid extrusion to form a particular shape. In order to form stable structures, the rheology of the paste is important to allow for extrusion from the syringe, stability of the growing print, and prevention of unwanted seeping flow during jog moves. In this work, we use wood pulp as a bio-based filler that can provide shear thinning properties to the ink, which produces a stable ink for DIW processing. Additionally, the filler imparts improved mechanical and thermal performance compared to neat resin. The wood pulp provided the shear thinning behavior necessary for DIW printing, and pulp loadings greater than 6 wt%, provided sufficient yield stress so that the composite could self-support during printing. Nanoclay was utilized to further improve ink rheology and appearance to enable larger scale printing. Overall, this work showed successful DIW of an epoxy resin with a sustainable filler improving its stiffness and thermal properties and provides an avenue for further development of bio-based inks for DIW towards various applications.

Lamm, Meghan E. [Oak Ridge National Laboratory (OR↗
DOE OSTI2021

Learning intermolecular forces at liquid–vapor interfaces

By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2023

Hydrophobic, Thermal-shock Resistant Lightweight Cement Composites in Hybrid Thermal Energy-Storage Geothermal Well System

Under the DOE EERE Geothermal Technology Office (GTO)-funded project in FY 2019-2022, Brookhaven National Laboratory (BNL) aimed to develop hydrophobic thermal-shock resistant lightweight well cement composites suitable for temporarily storing a curtailed hot working fluid incurred from excess geothermal-, solar-, and wind-generated energies in underground reservoirs thermal energy storage sites. This storage energy systems are designed to provide energy during the periods between energy-harvesting and -requiring periods. One key factor to the success of this project was hydrophobicity of cement composites. The hydrophobic composites not only upgrade thermal insulating properties of water-saturated composites for minimizing heat loss in upper well surface regions, but also provide improved mechanical properties related to thermal-shock resistance and the protection of carbon steel (CS) casing against corrosion.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗
DOE OSTI2020

Deposit geometry and oxygen concentration spatial variations due to composition change in printed functionally graded components

Functionally graded materials (FGMs) with site specific chemical composition are commonly manufactured by directed energy deposition (DED). Although previous work fabricated an FGM with a compositional variation between a ferritic and austenitic alloy, difficulties arose due to variations in deposit shape with composition change. This problem also occurs for FGMs in literature; however, unlike other cases, the thermophysical properties of these two alloys were similar throughout the build. Here, we investigate the role of chemical composition and surface-active elements on deposit geometry during the manufacture of FGMs by laser DED. Single-track experiments for the relevant FGM compositions are analyzed with results from a well-tested, three-dimensional, transient numerical heat transfer and fluid flow model and thermodynamic calculations. Experiments showed deposit shape varied as a function of composition for constant laser power and scan speed. Thermodynamic analysis indicated that the oxygen solubility in the fusion zone varied significantly for each composition used for the FGM. Numerical modeling revealed that the change in fluid flow caused by Marangoni convection due to dissolved oxygen in the fusion zone were mainly responsible for the changes in deposit shape observed in experiments. Finally, because oxygen can be introduced into the fusion zone through the feedstock as well as the surrounding atmosphere, these findings elucidate a previously unconsidered aspect of process control during DED fabrication of FGMs.

36 MATERIALS SCIENCE↗