Search NASA⌕ Search

SEARCH · Search NASA

Results for “Volume fractions”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Soot Formation and Ignition Characteristics of Ethanol/Gasoline Blends in a Rapid Compression Machine

With the ever-increasing demand for sustainable energy, alcohol fuels have garnered interest for use in heavy duty engines. The significant infrastructure for ethanol production and blending of ethanol with gasoline make these fuels/fuel blends desirable candidates. However, development of heavy duty engine technology that is capable of burning alcohol fuels while retaining the advantages of traditional diesel combustion requires an improved understanding of the soot formation for these fuels under conditions relevant to mixing-controlled combustion. This work uses an extinction diagnostic to study the sooting tendency of ethanol and gasoline/ethanol blends ranging from E10 to E98 during ignition in a homogeneous environment. Experiments were conducted in a rapid compression machine (RCM) for compressed conditions of 20 ± 1 bar and an approximately constant temperature (± 10K) which was unique for each fuel. For a given soot volume fraction, a linear relationship was observed between ethanol content and the equivalence ratio in which that soot volume fraction was formed. Accounting for the oxygenated nature of ethanol, E85 and E98 fuels produced similar amounts of soot at a given Φ ox, , suggesting other factors outside of fuel oxygen content, such as fuel morphology, impact soot formation. Ignition delay data is reported for compressed pressures of 20 ± 1 bar and compressed temperatures ranging from 633 – 670 K for E10 and 771 – 789 K for E98. Varying pressures for E10 and E98 at conditions producing similar soot volume fractions demonstrated a linear dependence of soot formation on pressure, regardless of if the pressure considered was at top dead center or peak combustion pressure. Furthermore, the data gleaned from this work will be used to select soot models and chemical kinetic mechanisms for RCM simulations to ultimately model heavy duty engine technology with the studied fuels.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Unstable Invasion of Sedimenting Granular Suspensions

Here we investigate the development of mobility inversion and fingering when a granular suspension is injected radially between horizontal parallel plates of a cell filled with a miscible fluid. While the suspension spreads uniformly when the suspension and the displaced fluid densities are exactly matched, even a small density difference is found to result in a dense granular front which develops fingers with angular spacing that increase with granular volume fraction and decrease with injection rate. We show that the timescale over which the instability develops is given by the volume fraction dependent settling timescale of the grains in the cell. We then show that the mobility inversion and the nonequilibrium Korteweg surface tension due to granular volume fraction gradients determine the number of fingers at the onset of the instability in these miscible suspensions.

02 PETROLEUM↗
DOE OSTI2020

A Proposed Common Model of Multi-phase Strength and Equation of State for a Tri-laboratory Collaboration (Working Draft 1.1)

The proposed simple common model for multiphase strength and EoS (CMMP) is meant to be sufficiently simple that each of the collaborating labs can share in a common starting point. Another objective is to start with relatively simple assumptions, which will not necessarily capture details of the physical processes, and incrementally add complexity in order to identify the minimal-needed technical detail. Through this co-evolution of model and experiment, we will better learn the importance of various theoretical approximations and where to invest future resources in experiment and model development. This simple framework is based on pressure and temperature equilibrium of all co-existing phases combined with deviatoric stress averaging through a volume fraction weighted ow stress and a volume fraction weighted shear modulus. Implementations of the framework based on equilibrium phase fractions (i.e. instantaneous kinetic rate) and for finite rate transition kinetics are proposed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rheology of bi-disperse dense fiber suspensions

While significant progress has been made in the modeling and simulation of uniform fiber suspensions, no existing model has been validated for industrially-relevant concentrated suspensions containing fibers of multiple aspect ratios. In the present work, we investigate bi-disperse suspensions with two fiber populations in varying aspect ratios in a steady shear flow using direct numerical simulations. Moreover, we measure the suspension viscosity by creating a controlled length bidispersity for nylon fibers suspended in a Newtonian fluid. Here the results showed good agreement between the experimentally measured and numerically predicted viscosity for bi-disperse suspensions. The ratio between the aspect ratio of large to small fibers (size ratio) and the volume fraction of large fibers (composition) in bi-disperse systems strongly affected the rheological behavior of the suspension. The increment of relative viscosity associated with size ratio and composition can be explained by the decrease in the maximum flowable limit or jamming volume fraction. Moreover, the relative viscosity of bi-disperse suspensions collapses, when plotted against the reduced volume fraction, demonstrating the controlling influence of the jamming fraction in bi-disperse fiber suspensions.

42 ENGINEERING↗
DOE OSTI2025

Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods

The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.

Vargas, Arturo [Lawrence Livermore National Labora↗
DOE OSTI2023

Economic competitiveness of pultruded fiber composites for wind turbine applications

Pultrusion manufacturing of fiber reinforced polymers has been shown to yield some of the highest mechanical properties for unidirectional composites, having a high degree of fiber alignment with consistent performance. Pultrusions offer a low-cost manufacturing approach for producing unidirectional composites with a constant cross-section and are used in many applications, including spar caps of wind turbine blades. However, as an intermediate processing step for wind blades, the additional cost of manufacturing pultrusions must be accompanied by sufficient increases in mechanical performance and system benefits. Wind turbine blades are manufactured using vacuum-assisted resin transfer molding with infused unidirectional fiberglass or carbon pultrusions for the spar cap. Infused fiberglass composites are among the most cost-effective structural materials available and replacing this material in the cost-driven wind industry has proven challenging, where infused fiberglass spar caps are still the predominant material system in use. To evaluate alternative material systems in a pultruded composite form, it is necessary to understand the costs for this additional manufacturing step which are shown to add 33%–55% on top of the material costs. A pultrusion cost model has been developed and used to quantify cost sensitivities to various processing parameters. The mechanical performance for pultruded composites is improved versus resin-infusion manufacturing with a 17% increase in design strength at a constant fiber volume fraction, but also enables higher achievable fiber volume fractions. The cost-specific mechanical performance is compared as a function of processing parameters for pultruded composites to identify the opportunities for alternative material and manufacturing approaches for wind turbine spar caps. Finally, four materials are compared in a representative wind turbine blade model to assess the performance of pultruded carbon fiber systems and pultruded fiberglass relative to infused fiberglass, where the pultruded systems produce lower weight blades with various cost distinctions.

42 ENGINEERING↗
DOE OSTI2021

Influence of O 2 / N 2 oxidizer blends on soot formation and radiative heat flux in PMMA-air 2D slab burner for understanding hybrid rocket combustion

The focus of this study is to establish a correlation between the flame to fuel incident radiative heat flux (RHF) in 1D steady-state diffusion flames to non-intrusively measured RHF in 2D slab burner experiments for a range of O 2 /N 2 oxidizer blends. Slab burner experiments are conducted for a range of increasing O 2 concentrations from 17% to 100% in the oxidizer stream. Non-intrusive two-color pyrometry is used to determine flame temperature, soot volume fraction, and flame to fuel surface RHF. Measurements and modeling confirm soot volume fraction attains a local maximum for an O 2 concentration of 45%. Low-strain diffusion flame modeling shows the reason for the local maximum is the relative location of the soot formation/nucleation fronts relative to oxidation which decreases with increasing O 2 concentration. Overall good agreement between modeling and measurements of radiative heat flux to the fuel is observed for the fuel and range of O 2 concentrations considered. Furthermore, these results demonstrate the viability of using flamelet modeling approaches for describing hybrid rocket motors.

2D slab burner↗
DOE OSTI2024

Exploring 2D X-ray diffraction phase fraction analysis with convolutional neural networks: Insights from kinematic-diffraction simulations

Abstract Deep-learning models are effective for analyzing the complex information in 2D X-ray diffraction (XRD) patterns. Accurately collecting parameters of the material sample is crucial during model training, significantly impacting model performance. In this study, we employ a kinematic-diffraction simulator to generate simulated 2D XRD patterns for Ti–6Al–4V alloy, allowing precise control of sample parameters. These simulated patterns are used to train convolutional neural networks, predicting $$\upbeta$$ β -phase volume fractions. The training data set consists exclusively of 2D XRD patterns with pure $$\upalpha$$ α - or pure $$\upbeta$$ β -phase, while the testing set incorporates patterns with intermediate phase volume fraction. In particular, we investigate how the architectures of the model influence prediction reliability and computational performance. Experimental results reveal that, with appropriate training, the convolutional neural network accurately detects intermediate phase volume fractions even trained with only pure-phase patterns, achieving a mean square error accuracy of $$9.4 \times 10^{-4}$$ 9.4 × 10 - 4 . Graphical abstract

Yue, Weiqi↗
DOE OSTI2020

Effect of copper content on the tensile elongation of Al–Cu–Mn–Zr alloys: Experiments and finite element simulations

Microstructures of cast aluminum alloys used in automotive engine applications often consist of intermetallic particles that can impact the tensile elongation of these alloys. Here in this paper, we investigate the effect of intermetallic grain boundary particles on tensile elongation by fabricating a series of alloys with Cu content varying between 6.0 - 9.0 wt% in cast Al–Cu–Mn–Zr (ACMZ) type compositions. The tensile elongation of as-aged ACMZ alloys decreases monotonically with increase in Cu content. While the microstructure within the grains and yield stress of the alloy remains invariant with Cu content, the decrease in tensile elongation correlates well with increase in the size and volume fraction of grain boundary particles. Microstructural observations are combined with finite element simulations to explain the trend in tensile elongation with changing Cu content. Crack initiation is found to occur by brittle fracture of the grain boundary particles. Increase in particle size promotes crack initiation by reduction in size dependent particle fracture strength. Lower inter-particle spacing at higher particle volume fraction further facilitates crack initiation by increasing stress within the particles caused by the interaction between stress fields of neighboring particles. Increase in particle volume fraction also accelerates crack propagation through the formation of macro shear zones in the microstructure. The increase in Cu content of cast ACMZ alloys, therefore, decreases tensile elongation by promoting both crack initiation and crack propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Supersaturating Water with CO 2 via Nanobubble Formation

The subsurface storage of CO 2 is limited by the solubility of CO 2 in brine. One proposed solution to increase the volume fraction of CO 2 in a reservoir is by injecting the CO 2 at high pressure in nanobubble form. These nanobubbles, diameter <50 nm, would be small enough that they would only undergo Brownian motion, preventing them rising under buoyancy and escaping the formation. There has also been some experimental work that suggests the present of these nanobubbles significantly (~2x) increase the solubility of gaseous CO 2 in the aqueous phase. NETL and ExxonMobil will test the viability of this strategy to increase the volume fraction of CO 2 in a subsurface reservoir, determine the stability of these nanobubble systems, and understand the ultimate limits of increasing the CO 2 volume fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Competing Magnetic Interactions and Field-Induced Metamagnetic Transition in Highly Crystalline Phase-Tunable Iron Oxide Nanorods

The inherent existence of multi phases in iron oxide nanostructures highlights the significance of them being investigated deliberately to understand and possibly control the phases. Here, the effects of annealing at 250 °C with a variable duration on the bulk magnetic and structural properties of high aspect ratio biphase iron oxide nanorods with ferrimagnetic Fe3O4 and antiferromagnetic α-Fe2O3 are explored. Increasing annealing time under a free flow of oxygen enhanced the α-Fe2O3 volume fraction and improved the crystallinity of the Fe3O4 phase, identified in changes in the magnetization as a function of annealing time. A critical annealing time of approximately 3 h maximized the presence of both phases, as observed via an enhancement in the magnetization and an interfacial pinning effect. This is attributed to disordered spins separating the magnetically distinct phases which tend to align with the application of a magnetic field at high temperatures. The increased antiferromagnetic phase can be distinguished due to the field-induced metamagnetic transitions observed in structures annealed for more than 3 h and was especially prominent in the 9 h annealed sample. Our controlled study in determining the changes in volume fractions with annealing time will enable precise control over phase tunability in iron oxide nanorods, allowing custom-made phase volume fractions in different applications ranging from spintronics to biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In situ synchrotron diffraction and modeling of non-equilibrium solidification of a MnFeCoNiCu alloy

The solidification mechanism and segregation behavior of laser-melted Mn 35 Fe 5 Co 20 Ni 20 Cu 20 was firstly investigated via in situ synchrotron x-ray diffraction at millisecond temporal resolution. The transient composition evolution of the random solid solution during sequential solidification of dendritic and interdendritic regions complicates the analysis of synchrotron diffraction data via any single conventional tool, such as Rietveld refinement. Therefore, a novel approach combining a hard-sphere approximation model, thermodynamic simulation, thermal expansion measurement and microstructural characterization was developed to assist in a fundamental understanding of the evolution of local composition, lattice parameter, and dendrite volume fraction corresponding to the diffraction data. This methodology yields self-consistent results across different methods. Via this approach, four distinct stages were identified, including: (I) FCC dendrite solidification, (II) solidification of FCC interdendritic region, (III) solid-state interdiffusion and (IV) final cooling with marginal diffusion. It was found out that in Stage I, Cu and Mn were rejected into liquid as Mn 35 Fe 5 Co 20 Ni 20 Cu 20 solidified dendritically. During Stage II, the lattice parameter disparity between dendrite and interdendritic region escalated as Cu and Mn continued segregating into the interdendritic region. After complete solidification, during Stage III, the lattice parameter disparity gradually decreases, demonstrating a degree of composition homogenization. The volume fraction of dendrites slightly grew from 58.3 to 65.5%, based on the evolving composition profile across a dendrite/interdendritic interface in diffusion calculations. Postmortem metallography further confirmed that dendrites have a volume fraction of 64.7 ± 5.3% in the final microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗