Search NASA⌕ Search

SEARCH · Search NASA

Results for “Very viscous flows”

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

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

In-situ magnetic alignment model for additive manufacturing of anisotropic bonded magnets

In this work, we report a mathematical framework which predicts the degree of alignment (DoA) in an in-situ aligned additively manufactured 3D printed bonded magnets. A multiphysics model is developed which couples the harmonious interactions of magnetic particles in a viscous flowing polymer under the presence of an external magnetic field. The hydrodynamic fluid-particle interaction is paired with the magnetophoretic force to predict the particle trajectory and distribution during extrusion through a 3D printer nozzle. Succeeding the force balance, a magnetohydrodynamic torque equilibrium analysis is performed to predict the net-orientation of the magnetic particles as a function of the applied field strength, viscous forces, and particle-to-particle interactions (P2P). Experimental validation of the DoA predictions is performed using 65 vol% Nd-Fe-B+Sm-Fe-N in Nylon12 (DoA exp = 0.620 and DoA theory = 0.686) and 15 vol% Sm-Co in PLA (DoA exp = 0.830 and DoA theory = 0.863). A parametric analysis is performed to analyze the effect of operating and design parameters like alignment field strength, magnetic loading fraction, extrusion load, and particle size. The model predicts a competing behavior between particle-fluid and particle-particle interactions under the presence of an applied field. Lastly, the model provides a framework to efficiently predict the DoA in tandem with a functionalized-magnetic 3D printer and allows the user to adjust the operating parameters according to the desired DoA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unraveling the Highly Plastic Behavior of ALD-Aluminum Oxide Encapsulations by Small-Scale Tensile Testing

We present a study directly measuring the electron-beam-induced plasticity of amorphous Al 2 O 3 coatings. Core–shell nanostructures are employed as small-scale model systems for two-dimensional coatings made by atomic layer deposition (ALD). Copper nanowires (NWs) are used as substrates for ALD deposition, representing a model system for interconnects commonly found in integrated circuits. Experiments are performed in situ in a transmission electron microscope (TEM) and further analyzed with electron energy loss spectroscopy (EELS). Our in situ TEM tensile experiments reveal the highly plastic behavior of the ALD shell, which withstands a maximum strain of 188%. Comparable samples under beam-off conditions show a brittle fracture, which underlines the effect of electron irradiation. The electron-beam-activated bond switching within the amorphous network enables compensation of the applied tensile strain, leading to viscous flow. By incorporating an intermediate nanocrystalline layer within the Al 2 O 3 shell, the plasticity is suppressed and brittle fracture occurs. This work directly demonstrates the tuning of mechanical properties in amorphous ALD structures through electron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rigidity‐Driven Structural Isomers in the NaCl–Ga 2 S 3 System: Implications for Energy Storage

Alternative energy sources require the search for innovative materials with promising functionalities. Systems with unusual chemical properties represent an insufficiently explored domain, concealing unexpected features. Using diffraction and Raman spectroscopy over a wide temperature range, supported by first‐principles simulations, a rare phenomenon is unveiled: phase‐dependent chemical interactions between binary components in the NaCl–Ga 2 S 3 system. In this unique occurrence, previously intact binary crystalline species transform upon melting into mixed liquid structural isomers, forming bonds with new partners. The chemical combinatorics appears to be fully reversible for stable crystals and liquids. Despite this, rapidly frozen glasses out of thermodynamic equilibrium remain in a metastable isomeric state, offering remarkable properties, particularly a high room‐temperature Na + conductivity, comparable to the best sodium halide superionic conductors and therefore encouraging for sodium solid‐state batteries and energy applications. A rigidity paradigm is responsible for the observed phenomenon, as the extremely constrained Ga 2 S 3 crystal lattice does not survive viscous flow, breaking up at a short‐range level. The removal of rigidity constraints and dense packing leads to a significant increase in empty space, which is the origin of high sodium diffusivity. Broadly, the rigidity‐driven structural isomerism opens up an inspiring path to the discovery of atypical materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Nonaxisymmetric Shapes of Biological Membranes from Locally Induced Curvature

In various biological processes such as endocytosis and caveolae formation, the cell membrane is locally deformed into curved morphologies. Previous models to study membrane morphologies resulting from locally induced curvature often only consider the possibility of axisymmetric shapes-an indeed unphysical constraint. Past studies predict that the cell membrane buds at low resting tensions and stalls at a flat pit at high resting tensions. In this work, we lift the restriction to axisymmetry to study all possible membrane morphologies. Only if the resting tension of the membrane is low, we reproduce axisymmetric membrane morphologies. When the resting tension is moderate to high, we show that 1) axisymmetric membrane pits are unstable and 2) nonaxisymmetric ridge-shaped structures are energetically favorable. Furthermore, we find the interplay between intramembrane viscous flow and the rate of induced curvature affects the membrane's ability to transition into nonaxisymmetric ridges and axisymmetric buds. In particular, we show that axisymmetric buds are favored when the induced curvature is rapidly increased, whereas nonaxisymmetric ridges are favored when the curvature is slowly increased. Overall, our results hold relevant implications for biological processes such as endocytosis and physical phenomena like phase separation in lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Coupled confined phase behavior and transport of methane in slit nanopores

The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.

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

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Densification of non-radioactive porous siliceous particles loaded with cesium potassium copper hexacyanoferrate by spark plasma sintering

In the development of a conditioning process for cesium-bearing sorbents from the decontamination of aqueous waste streams, the treatment of a granular material composed of (Cs,K) 2 CuFe(CN) 6 loaded onto mesoporous silica particles has been investigated in an additive-free two-step conditioning process involving (i) calcination at 1000 °C to remove volatile species and (ii) densification by spark plasma sintering. A 92.0 % dense pellet was obtained after 10 min sintering at 850 °C under 100 MPa. Density and microstructural measurements indicate that densification is first driven by the viscous flow of a vitreous phase containing Cs, formed during calcination. Additionally, the density of the sample then increases further during sintering by grain boundary diffusion. However, a limited but non negligible part of cesium is lost during densification.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Viscosity-temperature relation based on the evolution of medium-range structures of silica

The viscosity of glass is the most important technological property to glass manufactures and various applications. Practically, finding an accurate equation to express the glass viscosity behavior in the entire glass transition temperature range is tremendously challenge because it spans more than ten orders of magnitude. After a brief review of existing empirical viscosity equations, this work focuses on the correlating silica viscosity behavior with the evolution of glass medium-range structure, based on the recently proposed nanoflake model. From this new model, a new equation is constructed, which correctly describes the Arrhenius-type behavior of silica viscosity η above the melting temperature Tm, and non-Arrhenius-type behavior from Tm to a critical temperature Tc. At temperature lower than Tc, the equation predicts Arrhenius-type behavior again for the viscous flow with increased activation energy. The new equation agrees with experimental data in the entire temperature spanned from extremely high to extremely low. Furthermore, the application of the new equation is shown to extend to sodium silicate glasses as well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi‐Layered Systems for Permanent Geologic Storage of CO 2 at the Gigatonne Scale

Abstract The effectiveness of Carbon Capture and Storage (CCS) as an imperative decarbonization technology relies on the sealing capacity of a fine‐grained caprock to permanently store CO 2 deep underground. Uncertainties in assessing the caprock sealing capacity increase with the spatial and temporal scales and may delay CCS deployment at the gigatonne scale. We have developed a computationally efficient transport model to capture the dynamics of basin‐wide upward CO 2 migration in a multi‐layered setting over geological time scales. We find that massive capillary breakthrough and viscous flow of CO 2 , even through pervasively fractured caprocks, are unlikely to occur and compromise the storage security. Potential leakage from the injection reservoir is hampered by repetitive layering of overlying caprocks. This finding agrees with geologic intuition and should be understandable by the public, contributing to the development of climate policies around this technology with increased confidence that CO 2 will be indefinitely contained in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

Unravel unusual hardening behavior of a Pd–Ni–P metallic glass in its supercooled liquid region

Monolithic metallic glasses often exhibit work softening induced by high atom mobility in the supercooled liquid region in tension. In this work, we report an unusual viscous flow of the Pd 42.5 Ni 42.5 P 15 metallic glass in its supercooled liquid region, which is characterized by a remarkable hardening behavior with decent plasticity during tension. To unravel the causes of this unusual hardening, we conducted structural and thermodynamic analyses by employing several experiments, including dynamic mechanical analysis, differential scanning calorimetry, high-resolution transmission electron microscopy, and in situ heating synchrotron high-energy x-ray diffraction. It was found that the unusual hardening behavior was attributed to structural ordering during deformation, rather than nanocrystallization or structural relaxation observed, sometimes, in metallic glasses deformed in supercooled liquid region. We also incorporated an additional work hardening region in the modified deformation map for metallic glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Nanoscale dynamics during self-organized ion beam patterning of Si. II. Kr + bombardment

Understanding the self-organized ion beam nanopatterning of elemental semiconductors, particularly silicon, is of intrinsic scientific and technological interest. This is the second component of a two-part coherent x-ray scattering and x-ray photon correlation spectroscopy (XPCS) investigation of the kinetics and fluctuation dynamics of nanoscale ripple development on silicon during 1 keV Ar + (part I) and Kr + bombardment at 65° polar angle. Here it is found that the ion-enhanced viscous flow relaxation is essentially equal for Kr + -induced patterning as previously found for Ar + patterning despite the difference in ion masses. However, the magnitude of the surface curvature-dependent roughening rate in the early-stage kinetics is larger for Kr + than for Ar + , consistent with expectations that the heavier ion gives an increased mass redistributive contribution to the initial surface instability. As with the Ar + case, fluctuation dynamics in the late stage show a peak in correlation times at the length scale corresponding to the dominant structural feature on the surface—the ripples. Finally, it is shown that speckle motion during the surface evolution can be analyzed to determine spatial inhomogeneities in erosion rate and ripple velocity. This allows the direction and speed of ripple motion to be measured in a real time experiment. In the present case, ripple motion is found to be into the projected direction of the ion source, in contrast to expectations from an existing sputter erosion driven model with parameters derived from binary collision approximation simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A strain density function to analyze particle size effects during high velocity impacts of yttria

Here, the impact behavior of a single yttria (Y 2 O 3 ) nanoparticle onto a Y 2 O 3 substrate is studied as a function of particle velocity (300–1200 m/s) and diameter (12–50 nm) using molecular dynamics simulations. To analyze the results, a strain density function is developed to provide both quantitative and qualitative information about the deformation mechanisms that contribute to the final state of the system. This function provides both clear evidence of shear localization and insight into the conditions for which localization occurs during the impact of Y 2 O 3 particles as they approach experimentally relevant sizes. Further analysis shows that localization and fragmentation only occur for Y 2 O 3 impacts with diameters of at least 50 nm, while particles with diameters of 25 nm and smaller deform primarily by amorphization and viscous flow. Implications for the deposition of films of Y 2 O 3 and other rare-earth oxides via the micro cold spray process are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗
DOE OSTI2020

Shear-induced lift force on spheres in a viscous linear shear flow at finite volume fractions

Several studies have shown a significant increase in drag on a distribution of solid spherical particles within a fluid with increasing particle volume fraction. As a result, many empirical drag laws accounting for the dependence on the Reynolds number and volume fraction can be found in the literature. This study investigates the possibility of a similar effect of the particle volume fraction on the mean hydrodynamic lift force on randomly distributed spherical particles in a linear shear flow. Particle-resolved direct numerical simulations are performed to evaluate the mean lift force, and the results are compared with the case of an isolated particle in a linear shear flow for the same Reynolds number and shear rate. The mean lift force acting on the particles appears to remain nearly the same as that on an isolated particle. However, due to the influence of neighboring particles, there is a substantial force variation in transverse directions on each individual particle, whose magnitude is comparable to the mean drag force. The distribution of drag force in a linear shear flow is shown to be nearly the same as in a uniform flow at the same volume fraction and Reynolds number. A simple stochastic model based on a Gaussian distribution is presented for the lift force variation, and its performance is compared to the prediction of the deterministic pairwise interaction extended point-particle model.

42 ENGINEERING↗
DOE OSTI2022

Bidisperse supension balance model

The suspension balance model (SBM) for viscous Stokes flow has been well studied for the case of identical, monodisperse spherical particles in a channel, however, more work remains to be done to explore the bidisperse and polydisperse SBM. Here, we present a simple extension of the SBM that allows for modeling of suspensions of particles of bidisperse size. The comparison with avaiable experiments and direct simulations is found to be good. Additionally, we present a range of simulations to justify the assumptions made in the polydisperse suspension balance model for non-Brownian, neutrally buoyant particles that vary in size with a bidisperse distribution. The aim is to study the effects of moderate size variation on the rheology of the suspension and the distribution of the particles across the flow, especially in the near-wall regions. It is shown that, at the size ratios considered, bidispersity in the particle size does not significantly affect rheological parameters including the particle volume fraction and relative viscosity in the bulk of a Couette flow. The particle phase stresses are distributed between small and large particles in proportion to the relative volume fractions of each species. In bidisperse suspensions, the large particles in the wall layer develop a spatial structure and form chains in the streamwise direction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Calculation of Hydraulic Conductivity for a Porous Medium with Staggered Array of Solid Grains - 20087

For a porous medium with periodic array of solid grains. the flow of an incompressible viscous fluid in the pore space is considered with an aim of calculating the permeability of the medium. The solid grains are of ellipsoidal shape. For this purpose, a unit-cell boundary-value problem (Stokes problem) is solved numerically. With the assumption of periodicity of the medium structure and the fluid properties (in fact constant), the flow field and pressure distributions are determined numerically with inertial effects ignored. By applying the theory of homogenization to the viscous flow through pore space, a microscale boundary-value problem is defined. After solving the problem, the permeability is determined by taking the micro-cell average of the velocity components, Various porous structures are chosen to allow a wide range of porosity values. The permeability increases mildly for porosity values in the lower range. But it increases sharply for porosity values in the upper range of the porosity. The permeability of the present study agrees with the Kozeny-Carman relation excellently for smaller porosity. Certain discrepancy between the calculated permeability in this study and the Kozeny-Carman relation appears and increases with porosity for larger porosity. It is believed that the discrepancy is due to the replacement of the pores by a bundle of capillaries, which is unrealistic, in the Kozeny-Carman relation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗