Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex fluid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

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.↗

Surface Systems and Interface Standardization

Space exploration on planetary surfaces will require the use of various surface systems which will likely need to interface with one another. These systems must communicate and share data, as well as distribute power and transfer fluids for sustainable surface operations. Such systems can range from landers, surface habitats, mobility systems, cargo, and In-Situ Resource Utilizations (ISRUs). Previous and current programs have developed space interoperability standards that help aid in reducing potential risks of interface integration. One example of this is the International Deep Space Interoperability Standards (IDSIS) which focuses on deep space exploration. A goal of the Artemis program is to develop a sustained human presence on the lunar surface that would be a training ground for future Mars exploration. Therefore, it will be beneficial to identify interface standards between the surface assets which will help mitigate risk and reduce complexity in the harsh environments of space exploration. A key contribution to sustainability is commonality between hardware and software interfaces. Generic interfaces for data, power, and fluids will reduce risk, promote interoperability, and define standard interfaces across surface exploration projects and programs. Standardized interfaces would be advantageous for improving efficiency and reducing overall complexity, which are critical considerations for future space exploration. Furthermore, it will provide cost reductions to the Artemis programs over their life cycle (in Operations & Maintenance (O&M) and Logistics). Unique proprietary interfaces if considered or allowed would not only increase complexity but also add cost to the programs. Of course, understanding what will work and not work effectively in these unique environments such as the lunar surface is important. These unique environments require much needed intelligent design, prototyping, comprehensive testing, and field experience, utilizing consensus on common interface solutions.

Jaime Gomez↗

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↗

Chemical reacting flows

Future aerospace propulsion concepts involve the combustion of liquid or gaseous fuels in a highly turbulent internal airstream. Accurate predictive computer codes which can simulate the fluid mechanics, chemistry, and turbulence-combustion interaction of these chemical reacting flows will be a new tool that is needed in the design of these future propulsion concepts. Experimental and code development research is being performed at LeRC to better understand chemical reacting flows with the long-term goal of establishing these reliable computer codes. Our approach to understand chemical reacting flows is to look at separate, more simple parts of this complex phenomenon as well as to study the full turbulent reacting flow process. As a result, we are engaged in research on the fluid mechanics associated with chemical reacting flows. We are also studying the chemistry of fuel-air combustion. Finally, we are investigating the phenomenon of turbulence-combustion interaction. Research, both experimental and analytical, is highlighted in each of these three major areas.

Mularz, Edward J.↗

Mathematical description of complex chemical kinetics and application to CFD modeling codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗

Mathematical Description of Complex Chemical Kinetics and Application to CFD Modeling Codes

A major effort in combustion research at the present time is devoted to the theoretical modeling of practical combustion systems. These include turbojet and ramjet air-breathing engines as well as ground-based gas-turbine power generating systems. The ability to use computational modeling extensively in designing these products not only saves time and money, but also helps designers meet the quite rigorous environmental standards that have been imposed on all combustion devices. The goal is to combine the very complex solution of the Navier-Stokes flow equations with realistic turbulence and heat-release models into a single computer code. Such a computational fluid-dynamic (CFD) code simulates the coupling of fluid mechanics with the chemistry of combustion to describe the practical devices. This paper will focus on the task of developing a simplified chemical model which can predict realistic heat-release rates as well as species composition profiles, and is also computationally rapid. We first discuss the mathematical techniques used to describe a complex, multistep fuel oxidation chemical reaction and develop a detailed mechanism for the process. We then show how this mechanism may be reduced and simplified to give an approximate model which adequately predicts heat release rates and a limited number of species composition profiles, but is computationally much faster than the original one. Only such a model can be incorporated into a CFD code without adding significantly to long computation times. Finally, we present some of the recent advances in the development of these simplified chemical mechanisms.

Bittker, D. A.↗

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↗

A spatially adaptive high-order meshless method for fluid–structure interactions

We present a scheme implementing an a posteriori refinement strategy in the context of a high-order meshless method for problems involving point singularities and fluid–solid interfaces. The generalized moving least squares (GMLS) discretization used in this work has been previously demonstrated to provide high-order compatible discretization of the Stokes and Darcy problems, offering a high-fidelity simulation tool for problems with moving boundaries. The meshless nature of the discretization is particularly attractive for adaptive h -refinement, especially when resolving the near-field aspects of variables and point singularities governing lubrication effects in fluid–structure interactions. We demonstrate that the resulting spatially adaptive GMLS method is able to achieve optimal convergence in the presence of singularities for both the div-grad and Stokes problems. Furthermore, we present a series of simulations for flows of colloid suspensions, in which the refinement strategy efficiently achieved highly accurate solutions, particularly for colloids with complex geometries.

42 ENGINEERING↗

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↗

An algorithm for computing the 2D structure of fast rotating stars

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditioner of the Jacobian matrix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Permeability Variations Associated With Fault Reactivation in a Claystone Formation Investigated by Field Experiments and Numerical Simulations

We studied the relation between rupture and changes in permeability within a fault zone intersecting the Opalinus Clay formation at 300 m depth in the Mont Terri Underground Research Laboratory (Switzerland). A series of water injection experiments were performed in a borehole straddle interval set within the damage zone of the main fault. A three-component displacement sensor allowed an estimation of the displacement of a minor fault plane reactivated during a succession of step rate pressure tests. The experiment reveals that the fault hydromechanical (HM) behavior is different from one test to the other with varying pressure levels needed to trigger rupture and different slip behavior under similar pressure conditions. Numerical simulations were performed to better understand the reason for such different behavior and to investigate the relation between rupture nucleation, permeability change, pressure diffusion, and rupture propagation. Here, our main findings are as follows: (i) a rate frictional law and a rate-and-state permeability law can reproduce the first test, but it appears that the rate constitutive parameters must be pressure dependent to reproduce the complex HM behavior observed during the successive injection tests; (ii) almost similar ruptures can create or destroy the fluid diffusion pathways; (iii) a too high or too low diffusivity created by the main rupture prevents secondary rupture events from occurring whereas “intermediate” diffusivity favors the nucleation of a secondary rupture associated with the fluid diffusion. However, because rupture may in certain cases destroy permeability, this succession of ruptures may not necessarily create a continuous hydraulic pathway.

58 GEOSCIENCES↗

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↗

Numerical Simulation of Slag Eye Formation and Slag Entrapment in a Bottom-Blown Argon-Stirred Ladle

A transient mathematical model is developed for simulating the bubble-steel-slag-top gas four-phase flow in a bottom-blown argon-stirred ladle with a 70-ton capacity. The Lagrangian discrete phase model (DPM) is used for describing the moving behavior of bubbles in the steel and slag. To observe the formation process of slag eye, the volume of fluid (VOF) model is used to track the interfaces between three incompressible phases: metal/slag, metal/gas, and slag/gas. The complex multiphase turbulent flow induced by bubble-liquid interactions is solved by a large eddy simulation (LES) model. Slag eye area and slag droplet dispersion are investigated under different gas flow rates. The results show that the movement of bubbles, formation and collapse of slag eye, volatility of steel/slag interface and behavior of slag entrapment can be properly predicted in the current model. When the gas flow rate is 300 L/min, the circulation driven by the bubble plume will stir the entire ladle adequately and form a slag eye of the right size. At the same time, it will not cause strong erosion to the ladle wall, and the fluctuation of the interface is of adequate intensity, which will be helpful for improving the desulfurization efficiency; the slag entrapment behavior can also be decreased. Interestingly, with the motion of liquid steel circulation, the collision and coalescence of dispersed slag droplets occur during the floating process in the vicinity of the wall.

36 MATERIALS SCIENCE↗

A fully 3-D numerical model to predict flood wave propagation and assess efficiency of flood protection measures

Discussed are the development and validation of a 3-D, non-hydrostatic, Reynolds-Averaged Navier–Stokes (RANS) model using the Volume of Fluid (VOF) approach to simulate flooding events in natural environments. A series of validation simulations in simple and complex geometries are performed to test the different modules of the flow solver. The paper also discusses the differences between simulation results obtained using the 3-D model and those obtained using a standard, 2-D hydrostatic, depth-averaged model. The predictions by the two models show fairly large differences in regions where 3-D effects are expected to be significant (e.g., in regions containing highly-curved river reaches and near hydraulic structures). Then, the validated 3-D model is employed to investigate the efficiency of using flood-protection walls to reduce adverse effects of flooding in the Iowa River near Iowa City, Iowa, USA. Finally, the propagation of a flood wave in a reach of the Iowa River containing a bridge is simulated and the effects of the flow becoming pressurized beneath the bridge deck are discussed. For such applications, in which the regime can change locally between open channel and pressurized flow, the use of fully 3-D, non-hydrostatic models may be the only viable option to accurately predict flow hydrodynamics and its effects on bed shear stresses and sediment entrainment potential, especially around locations where the flow becomes pressurized.

Water Resources↗

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↗

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↗

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↗