Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid interfaces”

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

Mesoscale Polymer Surfactants: Photolithographic Production and Localization at Droplet Interfaces

Stabilization of fluid droplets, classically as oil-in-water or water-in-oil emulsions, is typically conducted using molecular surfactants or small particulates that localize at oil-water interfaces. In this paper, we describe a method whereby thin polymer films are converted photolithographically to ribbon-like mesoscale objects, which in turn adsorb to fluid inter-faces where they extend as appendages, or arms, from the droplet surface. These “mesoscale polymer surfactants”, or MPSs, were prepared from thin polymer films containing reactive functional moieties, including coumarin for photo-crosslinking, triphenylsulfonium for photoacid generation, and tert-butyl ester for solubility switching. The resultant MPSs, prepared initially on Si substrates, were released into water to reveal an exquisite shape sensitivity (forming straight, bent, or helical structures) and affinity for droplet interfaces based on their preparation conditions and the properties of the surrounding liquid. Notably, the lithographic techniques employed were amenable to differentiating the wettability of MPS segments, affording access to diblock-like MPSs which adhered to dispersed droplets via their hydrophobic segments, allowing their hydrophilic segments to extend into the continuous phase.

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

New insights into Mn 2+ and Mg 2+ inhibition of calcite growth

Impurity ion and isotope partitioning into carbonate minerals provide a window into the molecular processes occurring at the fluid-mineral interface during crystal growth. Here, we employ calcium isotope fractionation together with process-based modeling to elucidate the mechanisms by which two divalent cations with starkly contrasting compatibility, magnesium and manganese, inhibit calcite growth and incorporate into the mineral lattice. Calcite growth inhibition by Mg 2+ is log-linear and K Mg is on the order of 0.02–0.03 throughout the range of {Mg 2+ }/{Ca 2+ } studied here (0.01–2.6). Mn 2+ exhibits much stronger log-linear growth rate inhibition at low Mn 2+ concentrations (fluid {Mn 2+ }/{Ca 2+ } = 0.001–0.02). Mn 2+ is readily incorporated into the calcite lattice to form a calcite-rhodochrosite solid solution, with large partition coefficients (K Mn 4.6–15.6) inversely correlated to growth rate. For both Mn 2+ and Mg 2+ , calcium isotope fractionation is found to be invariant with {Me 2+ }/{Ca 2+ } despite more than an order of magnitude decline in growth rate. This invariant Δ 44/40 Ca suggests that the presence of Mn 2 + or Mg 2+ does not significantly change the relative rates of Ca 2+ attachment and detachment at kink sites during growth, indicative of a dominantly kink blocking inhibition mechanism. Because the partitioning behavior dictates that Mn 2+ must attach to the surface significantly faster than Ca 2+ , attachment of Mn 2+ is likely to be as a non-monomer species such as an ion pair or possibly a larger polynuclear cluster. We propose that calcite growth rate inhibition by Mn is determined by the kinetics of carbonate attachment at Mn-occupied kink sites, potentially due to slow re-orientation kinetics of carbonate ions that have formed an inner-sphere complex with Mn 2+ at the surface but must reorient to incorporate into the lattice. We demonstrate that patterns in Mg 2+ partitioning and inhibition behavior are broadly consistent with growth inhibition driven by slow Mg 2+ -aquo complex dehydration relative to Ca 2+ but argue that this mechanism likely represents one endmember scenario, seen in Mg-calcite growth at low supersaturations and net precipitation rates. During growth at faster net precipitation rates, some portion of Mg 2+ is likely incorporated as a partially hydrated or otherwise complexed species, but calcite growth remains significantly inhibited by the kinetics of CO 3 2- attachment at Mg 2+ kink sites. These findings suggest a hybrid classical/nonclassical growth mechanism whereby Ca 2+ incorporates largely as a free ion at kink sites while Mn 2+ and some portion of Mg 2+ are incorporated via non-monomer attachment. This pattern may be generalizable; trace constituent cations with aquo-complex desolvation rates significantly slower than the mineral growth rate preferentially incorporate as a non-monomer species during otherwise classical crystal growth.

58 GEOSCIENCES↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
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 OSTI2025

A simplified CFD approach for modeling mass transport in catalytic open-cell foams

A simplified macroscopic CFD approach is presented to model mass transport including chemical reactions in washcoated open-cell foams. The foam is treated as a porous medium. Species conversion during chemical reactions is modeled using appropriate source terms based on reaction rate expressions and modified to account for the mass transport resistances occurring at the fluid-washcoat interfaces and within the washcoat layers. As example, the catalytic CO oxidation over platinum is studied. The simulation results show good agreement with experimental data from literature. A parametric study on washcoat parameters, such as thickness, tortuosity, porosity, and size, is carried out. Increasing the washcoat thickness from 5 to 100 µm or decreasing the tortuosity to porosity ratio from 5 to 20 decreases the CO conversion by 10 %. The proposed model is found to be reliable and has the advantage of lower computational cost, making it a suitable tool for foam-based catalytic reactor design.

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

Validation of SPH code Spheral to model interacting solid bodies in a supersonic flow

Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.

Airburst↗
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2025

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗
DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗
DOE OSTI2025

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Foreword: Special Section on Multiphysics Aspects of Power Electronics Packaging - Power Die, Power Module, and Converter Level - Part 1

Power electronics are increasingly being used to condition electricity for a wide array of applications, such as transportation (on land, air, and water), data centers, radio frequency, directed energy, wind, solar, and grid-tied applications. To increase power density, performance, efficiency, and reliability and reduce cost, innovations and developments are needed in the multiphysics packaging of power electronics at a die, module, and converter level. This includes fundamental research and development related to emerging high-voltage, high-temperature, and high-switching-frequency power electronics; packaging materials; thermal materials and interfaces; fluid-based thermal management technologies; reliability; condition monitoring; and prognostics. To address these important aspects, this Special Section on Multiphysics Aspects of Power Electronics Packaging includes several articles to be published in two parts. More details are given below on the articles included in the first part.

condition monitoring↗
DOE OSTI2024

Foreword: Special Section on Multiphysics Aspects of Power Electronics Packaging—Power Die, Power Module, and Converter Level: Part 2

Power electronics are increasingly being used to condition electricity for a wide array of applications, such as transportation (on land, air, and water), data centers, radio frequency, directed energy, wind, solar, and grid-tied applications. Here, to increase power density, performance, efficiency, and reliability-as well as to reduce cost-innovations and developments are needed in the multiphysics packaging of power electronics at a die, module, and converter level. This includes fundamental R&D related to emerging high-voltage, high-temperature, and high-switching-frequency power electronics, packaging materials, thermal materials and interfaces, fluid-based thermal management technologies, reliability, condition monitoring, and prognostics. Latest developments in this area are published as a Special Section on Multiphysics Aspects of Power Electronics Packaging. The first part was published in the May 2024 issue of the IEEE Transactions on Components, Packaging and Manufacturing Technology (Volume 14, Issue 5). The second part of that Special Section is being published in this issue. A brief summary of the papers included in the second part are given below.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

SalineMoose

SalineMoose is a MOOSE wrapper for the open-source Saline code, which is an interface for fluid property databsae files. This code is meant to provide access to Saline, which can easily be integrated into MOOSE-based codes that need access to the fluid thermophysical properties that Saline can provide. Wrapping Saline in SalineMoose allows MOOSE codes to use the native MOOSE build system for linking Saline to their MOOSE application.

Salko, Robert↗
DOE OSTI2022

Effect of Packing Nonuniformity at the Fiber Bundle–Case Interface on Performance of Hollow Fiber Membrane Gas Separation Modules

High-fidelity simulations of momentum and mass transfer within a hollow fiber gas separation membrane module are here reported. The simulations capture the potential detrimental effects of poor fiber packing at the bundle–case interface on fluid distribution and performance. Results are presented for both circular and planar fiber bundles. The length over which bundle–case gaps affects flow is determined. The length increases dramatically with increasing fiber packing fraction. As the packing fraction approaches 0.6, the impact extends over the entire bundle diameter for small modules (<1000 fibers). The results clearly demonstrate the detrimental effect of poor packing along the case and can be used to develop module manufacturing guidelines. To reduce computational costs, an equivalent planar bundle module approximation is developed. The approximate simulations agree well with results from full 3-D simulations and can reduce computational costs without sacrificing fidelity.

42 ENGINEERING↗
DOE OSTI2023

Aeroelastic Analysis Using Deforming Cartesian Grids

Ongoing work in air-vehicle design illustrates the potential of advanced concepts to provide significant improvements in efficiency; but with their incorporation of lightweight flexible structures, such configurations may require active control systems to ensure reliability and safety. However, many contemporary analysis methods are inefficient for aeroelastic analysis and design of such configurations. This paper describes the development of a new approach that automates the geometry setup, mesh generation, and assembly of fluid–structural coupling interfaces to enable efficient aeroelastic and aeroservoelastic analysis of advanced concepts. The core elements for this approach are a cut-cell Cartesian grid-based computational fluid dynamics solver, a nonlinear beam element structural model, a conservative fluid–structural interface treatment, and the formulation and implementation of a new deforming grid capability within the cut-cell Cartesian grid solver. In this paper, emphasis is on this latter component with detailed description given of the mesh motion strategy, evaluation of fluxes and structural loads at the surface, and computation of geometrical properties such as cell volume, directed face areas, centroids, and motion-induced fluxes for deforming Cartesian grids required to advance the flow states. Aeroelastic simulations exercising the capability show favorable agreement with data and predictions in the literature for subsonic and supersonic applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Interface models for three-dimensional Rayleigh–Taylor instability

We derive interface models for three-dimensional Rayleigh–Taylor instability (RTI), making use of a novel asymptotic expansion in the non-locality of the fluid flow. These interface models are derived for the purpose of studying universal features associated with RTI such as the Froude number in single-mode RTI, the predicted quadratic growth of the interface amplitude under multi-mode random perturbations, the optimal (viscous) mixing rates induced by the RTI and the self-similarity of horizontally averaged density profiles and the remarkable stabilization of the mixing layer growth rate which arises for the three-fluid two-interface heavy–light–heavy configuration, in which the addition of a third fluid bulk slows the growth of the mixing layer to a linear rate. Our interface models can capture the formation of small-scale structures induced by severe interface roll-up, reproduce experimental data in a number of different regimes and study the effects of multiple interface interactions even as the interface separation distance becomes exceedingly small. Compared with traditional numerical schemes used to study such phenomena, our models provide a computational speed-up of at least two orders of magnitude.

Mechanics↗