Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluids”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effects of Carbonate Minerals on Shale-Hydraulic Fracturing Fluid Interactions in the Marcellus Shale

Natural gas extracted from tight shale formations, such as the Marcellus Shale, represents a significant and developing front in energy exploration. By fracturing these formations using pressurized fracturing fluid, previously unobtainable hydrocarbon reserves may be tapped. While pursuing this resource, hydraulic fracturing operations leave chemically complex fluids in the shale formation for at least two weeks. This provides a substantial opportunity for the hydraulic fracturing fluid (HFF) to react with the shale formation at reservoir temperature and pressure. In this study, we investigated the effects of the carbonates on shale-HFF reactions with a focus on the Marcellus Shale. We performed autoclave experiments at high temperature and pressure reservoir conditions using a carbonate-rich and a decarbonated or carbonate-free version of the same shale sample. We observed that carbonate minerals buffer the pH of the solution, which in turn prevents clay dissolution. Carbonate and bicarbonate ions also scavenge reactive oxidizing species (ROS), which prevents oxidation of shale organic matter and volatile organic compounds (VOCs). Carbonate-free samples also show higher pyrite dissolution compared to the carbonate-rich sample due to chelation reactions. This study demonstrates how carbonate minerals (keeping all other variables constant) affect shale-HFF reactions that can potentially impact porosity, microfracture integrity, and the release of heavy metals and volatile organic contaminants in the produced water.

58 GEOSCIENCES↗
DOE OSTI2022

Rheo-Structural Spectroscopy: Fingerprinting the In Situ Response of Fluids to Arbitrary Flow Fields

The objectives of this project were to develop new sample environments, measurement methodologies and associated modeling tools for characterizing the structural response to arbitrarily complex processing flows using small angle scattering, and to apply these new tools for understanding the fundamental physics governing the structuring of anisotropic particulate and polymeric materials under flow histories and conditions relevant to industrial processing flows. The research resulted in the development and implementation of a new sample environment, the fluidic four roll mill (FFoRM), for in situ small angle neutron and X-ray scattering (SANS/SAXS) measurements. These measurements are capable of generating large data sets that “fingerprint” how a complex fluid responds to a wide range of flow histories involving time variations in deformation type and rate. New modeling tools were developed to extract detailed microstructural information from such data sets, including orientation distribution functions and interparticle correlation functions, as well as reduced-order parametric descriptors of these high-dimensional functions that can be used to readily map, visualize and interpret a fluid’s structural response to its flow history. These new tools were applied to a range of model materials involving elongated particle suspensions in order to provide new insights into the physics of how flow couples with orientational and structural order in complex flows, particularly under non-dilute conditions for which no accurate theories currently exist. Using these investigations, we elucidated a number of new insights into the fundamental phenomena driving such process-structure-property relationships. These findings provide guidance for the further development of rheological models, and ultimately can inform the rational and model-based design of flow processes to achieve optimized orientational ordering that is key to the properties and function of a wide range of energy-relevant materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fluid Mechanics Challenges in Direct-Ink-Writing Additive Manufacturing

Direct-ink writing (DIW) has rapidly become a versatile 3D fabrication method due to its ability to deposit a wide range of complex fluids into customizable 3D geometries. This review highlights key fundamental fluid mechanics and soft matter challenges across the different stages of the DIW printing process. The rheology of fluids and suspensions governs the flow behavior through narrow nozzles, posing questions about extrudability, confined flow dynamics, and clogging mechanisms. Downstream, the formation and deposition of extruded filaments involve extensional flows and potential instabilities, while postdeposition dynamics introduces complexities related to yield stress and structural stability. These stages are inherently interdependent, as optimizing material composition without considering filament stability risks compromising the final structure. As DIW applications expand through advanced ink formulations, developing fundamental fluid mechanics frameworks is essential to replace trial-and-error approaches with predictive design methodologies to enable more precise control and improved reliability of the printing process.

3D printing↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

Recycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO 2 ) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ 34 S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H 2 O-NaCl solutions to calibrate a high-pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO 2 from FMQ-2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO 2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur-rich and high δ 34 S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34 S-rich subducted lithologies.

58 GEOSCIENCES↗
DOE OSTI2021

Topology optimization of additively manufactured fluidic components free of internal support structures

This report integrates projection-based approaches for implementing overhang constraints with fluid topology optimization to design additively manufactured fluidic components that do not require internal support structures. Internal support structures are challenging, potentially impossible to remove in complex fluid networks, yet would degrade fluidic performance if left within the part. This presents a major challenge in coupling topology optimization with additive manufacturing in the design of fluid components. The proposed approach leverages past work in overhang constraints for solid mechanics, including projection formulations and adjoint sensitivity analysis, and considers incompressible Navier–Stokes equations for laminar fluid flow. The approach is demonstrated on 2D and 3D fluid topology optimization problems to minimize pressure drop and constrain mass flow rates in pipes and manifolds. Resulting designs are crisp, logical, satisfy performance constraints, and satisfy overhang constraints, eliminating the need for internal support structures.

42 ENGINEERING↗
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2025

The structure of liquid carbon elucidated by in situ X-ray diffraction

Carbon has a central role in biology and organic chemistry, and its solid allotropes provide the basis of much of our modern technology. However, the liquid form of carbon remains nearly uncharted, and the structure of liquid carbon and most of its physical properties are essentially unknown. But liquid carbon is relevant for modelling planetary interiors and the atmospheres of white dwarfs, as an intermediate state for the synthesis of advanced carbon materials, inertial confinement fusion implosions, hypervelocity impact events on carbon materials and our general understanding of structured fluids at extreme conditions. Here we present a precise structure measurement of liquid carbon at pressures of around 1 million atmospheres obtained by in situ X-ray diffraction at an X-ray free-electron laser. Our results show a complex fluid with transient bonding and approximately four nearest neighbours on average, in agreement with quantum molecular dynamics simulations. The obtained data substantiate the understanding of the liquid state of one of the most abundant elements in the universe and can test models of the melting line. The demonstrated experimental abilities open the path to performing similar studies of the structure of liquids composed of light elements at extreme conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Prototype equation of state for phase transition of confined fluids based on the generalized van der Waals partition function

A simple self-consistent prototype equation of state (EOS) based on the Generalized van der Waals (vdW) partition function has been demonstrated to describe the phase transition of simple fluids in nanopores with uniform size. Different from those commonly presented in the literature, the new EOS does not need an auxiliary equation that is conventionally applied to provide the capillary pressure derived from surface tension. Lastly, the encouraging performance of the EOS calls for further extension to applications with more complex fluids and porous media.

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

A Graphics Processing Unit–Based, Industrial Grade Compositional Reservoir Simulator

Summary Recently, graphics processing units (GPUs) have been demonstrated to provide a significant performance benefit for black-oil reservoir simulation, as well as flash calculations that serve an important role in compositional simulation. A comprehensive approach to compositional simulation based on GPUs has yet to emerge, and the question remains as to whether the benefits observed in black-oil simulation persist with a more complex fluid description. We present a positive answer to this question through the extension of a commercial GPU-based black-oil simulator to include a compositional description based on standard cubic equations of state (EOSs). We describe the motivations for the selected nonlinear formulation, including the choice of primary variables and iteration scheme, and support for both fully implicit methods (FIMs) and adaptive implicit methods (AIMs). We then present performance results on an example sector model and simplified synthetic case designed to allow a detailed examination of runtime and memory scaling with respect to the number of hydrocarbon components and model size, as well as the number of processors. We finally show results from two complex asset models (synthetic and real) and examine performance scaling with respect to GPU generation, demonstrating that performance correlates strongly with GPU memory bandwidth. NOTE: This paper is also published as part of the 2021 SPE Reservoir Simulation Conference Special Issue.

Engineering↗
DOE OSTI2025

Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets

Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Hydrodynamically Controlled Self‐Organization in Mixtures of Active and Passive Colloids

Abstract Active particles are known to exhibit collective behavior and induce structure in a variety of soft‐matter systems. However, many naturally occurring complex fluids are mixtures of active and passive components. The authors examine how activity induces organization in such multi‐component systems. Mixtures of passive colloids and colloidal micromotors are investigated and it is observed that even a small fraction of active particles induces reorganization of the passive components in an intriguing series of phenomena. Experimental observations are combined with large‐scale simulations that explicitly resolve the near‐ and far‐field effects of the hydrodynamic flow and simultaneously accurately treat the fluid–colloid interfaces. It is demonstrated that neither conventional molecular dynamics simulations nor the reduction of hydrodynamic effects to phoretic attractions can explain the observed phenomena, which originate from the flow field that is generated by the active colloids and subsequently modified by the aggregating passive units. These findings not only offer insight into the organization of biological or synthetic active–passive mixtures, but also open avenues to controlling the behavior of passive building blocks by means of small amounts of active particles.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Three-dimensional realizations of flood flow in large-scale rivers using the neural fuzzy-based machine-learning algorithms

Machine learning methods have been extensively used to study the dynamics of complex fluid flows. One such algorithm, known as adaptive neural fuzzy inference system (ANFIS), can generate data-driven predictions for flow fields, but has not been applied to natural geophysical flows in large-scale rivers. Herein, we demonstrate the potential of ANFIS to produce three-dimensional (3D) realizations of the instantaneous flood flow field in several large-scale, virtual meandering rivers. The 3D dynamics of flood flow in large-scale rivers were obtained using large-eddy simulation (LES). The LES results, i.e., the 3D velocity components, were employed to train the learnable coefficients of an ANFIS. Further, the trained ANFIS, along with a few time-steps of LES results (precursor data) were then used to produce 3D realizations of flood flow fields in large-scale rivers with geometries other than the one the ANFIS was trained with. We also used the trained ANFIS to generate 3D realizations of river flow at a discharge other than that the ANFIS was trained with. The flow field results obtained from ANFIS were validated using separate LES runs to assess the accuracy of the 3D instantaneous realizations of the machine learning algorithm. An error analysis was conducted to quantify the discrepancies among the ANFIS and LES results for various flood flow predictions in large-scale rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Combined selection of the dynamic model and modeling error in nonlinear aeroelastic systems using Bayesian Inference

Here, we report a Bayesian framework for concurrent selection of physics-based models and (modeling) error models. We investigate the use of colored noise to capture the mismatch between the predictions of calibrated models and observational data that cannot be explained by measurement error alone within the context of Bayesian estimation for stochastic ordinary differential equations. Proposed models are characterized by the average data-fit, a measure of how well a model fits the measurements, and the model complexity measured using the Kullback–Leibler divergence. The use of a more complex error models increases the average data-fit but also increases the complexity of the combined model, possibly over-fitting the data. Bayesian model selection is used to find the optimal physical model as well as the optimal error model. The optimal model is defined using the evidence, where the average data-fit is balanced by the complexity of the model. The effect of colored noise process is illustrated using a nonlinear aeroelastic oscillator representing a rigid NACA0012 airfoil undergoing limit cycle oscillations due to complex fluid–structure interactions. Several quasi-steady and unsteady aerodynamic models are proposed with colored noise or white noise for the model error. The use of colored noise improves the predictive capabilities of simpler models.

42 ENGINEERING↗
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2021

Flow over an espresso cup: inferring 3-D velocity and pressure fields from tomographic background oriented Schlieren via physics-informed neural networks

Tomographic background oriented Schlieren (Tomo-BOS) imaging measures density or temperature fields in three dimensions using multiple camera BOS projections, and is particularly useful for instantaneous flow visualizations of complex fluid dynamics problems. We propose a new method based on physics-informed neural networks (PINNs) to infer the full continuous three-dimensional (3-D) velocity and pressure fields from snapshots of 3-D temperature fields obtained by Tomo-BOS imaging. The PINNs seamlessly integrate the underlying physics of the observed fluid flow and the visualization data, hence enabling the inference of latent quantities using limited experimental data. In this hidden fluid mechanics paradigm, we train the neural network by minimizing a loss function composed of a data mismatch term and residual terms associated with the coupled Navier–Stokes and heat transfer equations. We first quantify the accuracy of the proposed method based on a two-dimensional synthetic data set for buoyancy-driven flow, and subsequently apply it to the Tomo-BOS data set, where we are able to infer the instantaneous velocity and pressure fields of the flow over an espresso cup based only on the temperature field provided by the Tomo-BOS imaging. Moreover, we conduct an independent PIV experiment to validate the PINN inference for the unsteady velocity field at a centre plane. To explain the observed flow physics, we also perform systematic PINN simulations at different Reynolds and Richardson numbers and quantify the variations in velocity and pressure fields. Furthermore, the results in this paper indicate that the proposed deep learning technique can become a promising direction in experimental fluid mechanics.

97 MATHEMATICS AND COMPUTING↗