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

Destabilization of binary mixing layer in supercritical conditions

Compressible mixing layer instabilities are of importance to a wide range of environmental and industrial applications. Past studies have focused on either ideal-gas or real-fluid thermodynamic regimes of single-species mixing layers. However, mixing layers of binary mixtures at supercritical conditions, commonly encountered in fuel injection systems, introduce additional complexities due to the added compositional degree-of-freedom. Moreover, the effect of strong variations in thermodynamic response functions across the Widom line on the binary mixing layer stability remains poorly understood. Thus, the objective of this study is to examine the coupling between the hydrodynamic instability and the real-fluid thermodynamics across the Widom line and its effects on the overall binary mixing layer dynamics. To this end, we develop a linear stability analysis of the full binary-species compressible transport equations coupled with the PC-SAFT equation of state. Analysis shows the existence of a novel instability mechanism that arises from juxtapositioning of the Widom-line transition and the hydrodynamic inflection point. This novel thermodynamically induced instability mechanism has the net effect of destabilizing the binary mixing layer at lowering supercritical conditions towards the critical pressure point. This is in contrast to previous stability analyses of supercritical single-species mixing layers, where increasing pressure destabilizes the flow due to its effect on reducing the density stratification. The discovered thermodynamically induced instability mechanism of binary mixing flows highlights the need for an extension of classical instability criteria to incorporate the effect of strong variations in the thermodynamic response functions across the Widom line on mixing layer instability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The nonlinear rheology of complex yield stress foods

Many foods have a yield stress that allows them to retain a desired shape at rest, but transition into a viscous fluid when being served or consumed. The determination of the yield stress of the food dictates how the foods are formed and packaged, how they are served, and how they are perceived when being eaten. Oscillatory shearing provides an ideal test protocol to map the rheology across a range of time and flow strength scales. We couple oscillatory shearing and an iterative recovery procedure to show that the yielding process is a continuous transition for two common yield stress foods. We show that unrecoverable processes from oscillatory tests are equivalent to the steady shear flow behavior. We show that this yielding behavior can be well approximated by a recently published model that treats yield stress materials as continuous viscoelastic fluids with a rate-dependent relaxation time and has parameters that can be obtained from the linear viscoelastic oscillatory frequency sweep and the steady shear flow curve.

Griebler, James J.↗
DOE OSTI2025

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗
DOE OSTI2023

Particle resuspension: Challenges and perspectives for future models

Using what has become a celebrated catchphrase, Philip W. Anderson once wrote that “more is different” (Science, Vol. 177, Issue 4047, pp. 393–396, 1972). First formulated in the context of condensed matter, this statement carries far beyond the sole limits of solid-state physics. It emphasizes that collective behavior can be more than the mere sum of what happens for elementary constituents or the mere collation of the evolution of each degree of freedom. Said otherwise, complex phenomena can arise out of the interplay between multiple sub-phenomena each of which can be relatively simple. The process of particle resuspension, in which discrete particles adhering on a surface are pulled off and carried away by a fluid flow, is another example involving a web of phenomena pertaining to fluid mechanics, particle dynamics and interface chemistry whose cross-effects create an intricate topic. The purpose of this review is to analyze the physics at play in particle resuspension in order to bring insights into the rich complexity of this common but challenging concern. Following the more-is-different vision, this is performed by starting from a range of practical observations and experimental data. We then work our way through the investigation of the key mechanisms which play a role in the overall process. In turn, these mechanisms reveal an array of fundamental interactions, such as particle–fluid, particle–particle and particle–surface, whose combined effects create the tapestry of current applications. At the core of this analysis are descriptions of these physical phenomena and the different ways through which they are intertwined to build up various models used to provide quantitative assessment of particle resuspension. The physics of particle resuspension implies to hold together processes occurring at extremely different space and time scales and models are key in providing a single vehicle to lead us through such multiscale journeys. This raises questions on what makes up a model and one objective of the present work is to clarify the essence of a modeling approach. In spite of its ubiquitous nature, particle resuspension is still at the early stages of developments. Many extensions need to be worked out and revisiting the art of modeling is not a moot point. The need to consider more complex objects than small and spherical particles and, moreover, to come up with unified descriptions of mono- and multilayer resuspension put the emphasis on solid model foundations if we are to go beyond current limits. Further, this is very much modeling in the making and new ideas are proposed to stimulate interest into this everyday but challenging issue in physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Structural Heterogeneity and Hydrodynamics of an Intrinsically Disordered Protein Condensate

Biology demonstrates precise control over the free-energy landscape through the selective partitioning of biomacromolecules into membraneless organelles, enabling essential functions such as biochemical transformations, signaling cascades, and mechanical reinforcement. Although the function of these condensates depends on their underlying structure and hydrodynamics, molecular-scale information on these systems remains sparse. Here, in this study, neutron scattering is used to probe the organization and dynamics of the intrinsically disordered N-terminal domain of Galectin-3, an extracellular lectin responsible for facilitating liquid–liquid phase separation on the cellular surface, in both dilute and condensed phases. Dilute solutions contain isolated protein chains in equilibrium with mesoscopic clusters, whereas the condensed phase adopts a bicontinuous, microemulsion-like morphology. The dilute phase behavior is quantitatively described by coarse-grained polymer models from soft-matter physics, demonstrating their predictive power for complex biological proteins. At elevated concentrations, the proteins self-assemble akin to block copolymers, microphase separating through the aggregation of hydrophobic domains along the protein contour. The resulting condensate remains fluid-like despite a 25-fold increase in concentration; its internal hydrodynamics slow by only a factor of 3 relative to dilute protein chains. These results provide a molecular-level framework for how disordered proteins achieve both the structural complexity and dynamic fluidity of biomolecular condensates.

Carrick, Brian R. [Massachusetts Inst. of Technolo↗
DOE OSTI2021

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗
DOE OSTI2025

An experimental study of synthetic Hydroxybastnäsite-(La) solubility and speciation in carbonate bearing aqueous solutions at 175–250 °C

The transport and enrichment of rare earth element (REE) ore bodies are dependent on the stability of aqueous metal ligand complexes and the solubility of REE bearing minerals. REE ores are commonly associated with igneous systems having aqueous fluids with high carbonate concentrations and REE solubilities have been shown to be dependent on temperature and associate anion aqueous ligands present in solution. Furthermore, this work presents solubility experiments of hydroxybastnäsite-(La) at elevated temperatures in aqueous solutions of varying carbonate concentrations. At lower temperatures, hydroxybastnäsite-(La) solubility is controlled by neutral mono-carbonate LaCO 3 OH° but at higher temperatures and activities of carbonate species, charged di-carbonate La(CO 3 ) 2 - increases and predominates. This divergence, and the difference in solubility products of other hydroxybastnäsite-(REE) phases, provides a potential mechanism for REE fractionation in carbonate dominated aqueous solutions. To illustrate one such mechanism the solubility data of hydroxybastnäsite-(La) is compared with previously reported data of hydroxybastnäsite-(Nd) at elevated temperatures.

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

Episodic hydrocarbon charge in tight Mississippian reservoirs of Central Oklahoma, USA: Insights from oil inclusion geochemistry

The unforeseen discovery of substantial hydrocarbons in ultra-tight Mississippian reservoirs in Oklahoma testifies to the complexity of the hydrocarbon charge across the flanks of the Nemaha Uplift. Based on crude oil and fluid-inclusion chemistry, the hydrocarbons at the western side of the Nemaha Uplift consist of highly-mature, thermally-cracked and less-mature uncracked hydrocarbons (0.50–1.4% calculated vitrinite reflectance, R c ). The oils are less mature on the eastern flank of the structure (0.5–0.7% R c ). Based on distributions of biomarkers and higher diamondoid (composed of three adamantane cages or more) distributions, the hydrocarbons entrapped within reservoir fluid-inclusions are sourced from the Woodford Shale and also from Mississippian carbonates. To our knowledge, this is the first work to describe oil correlations using diamondoids inside fluid inclusions. Petroleum compositions in the ultra-tight Mississippian reservoirs are highly variable, which highlights the need for a thorough assessment of mechanisms controlling charge. By integrating oil and fluid inclusion geochemistry and fluid-inclusion microthermometry with basin modeling, we show that the primary controls are: (1) proximity to the source-rock kitchen(s), (2) episodic charging from Permian through Cretaceous time, and (3) mixing.

58 GEOSCIENCES↗
DOE OSTI2021

Modeling Snowflake Divertors in MAST-U Tokamak

We report that in a snowflake (SF) divertor, two magnetic field nulls are placed close to each other, creating four strike points (SPs) compared to two in a standard X-point divertor. In preparation for MAST-U experiments, magnetic configurations with the standard and SF divertors with various locations and separation distances of the nulls were modeled using the two-dimensional multi-fluid code UEDGE with a full plasma transport model featuring charge-state-resolved sputtered carbon impurities. The complex interplay of the plasma transport and magnetic configurations was comprehensively studied using a simple model for the theoretically predicted fast plasma mixing driven by the 'churning' mode instability in the two-null SF region. The modeling results show that (1) all SF-plus configurations and SF-minus configuration with closely located nulls produce the same plasma parameters and heat fluxes at the same SPs; (2) SF divertors approach the outer and inner SP detachment conditions at lower upstream density w.r.t. the standard divertor; (3) heat flux profiles at primary SPs are substantially broadened and peak values are reduced in SF configurations w.r.t. SN divertors; this broadening becomes more pronounced with the fast plasma mixing increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Physically rigorous reduced-order flow models of fractured subsurface environments without explosive computational cost

Fractured media models comprise discontinuities of multiple lengths (e.g. fracture lengths and apertures, wellbore area) that fall into the relatively insignificant length scales spanning millimeter-scale fractures to centimeter-scale wellbores in comparison to the extensions of the field of interest, and challenge the conventional discretization methods imposing highly-fine meshing and formidably large numerical cost. By utilizing the recent developments in the finite element analysis of electromagnetics that allow to represent material properties on a hierarchical geometry, this project develops computational capabilities to model fluid flow, heat conduction, transport and induced polarization in large-scale geologic environments that possess geometrically-complex fractures and man-made infrastructures without explosive computational cost. The computational efficiency and robustness of this multi-physics modeling tool are demonstrated by considering various highly-realistic complex geologic environments that are common in many energy and national security related engineering problems.

42 ENGINEERING↗
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 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↗
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 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↗