Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONTINUUM”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interplay between hydrogen, temperature, and character angle on dissociated dislocation energies in Fe–Ni–Cr austenitic stainless steels

Dislocation energy has an important role in the mechanical performance of structural metals. While dislocation energies cannot be fully obtained from continuum theories due to the contribution of the dislocation core, they have been calculated via atomistic simulations in elemental metals. However, constraints on the local atomic environments have prevented the use of such approaches in systems that incorporate alloying or interstitial solutes. In this work, we develop robust molecular dynamics methods to resolve these issues through a geometric construction of dislocation dipoles and the calculation of time-averaged energies. Furthermore, we apply these methods to calculate dislocation energies (including core energies) in an Fe 70 Ni 11 Cr 19 austenitic steel at a variety of character angles, hydrogen concentrations, temperatures, and dipole spacings. The resulting highly converged energies show an excellent agreement with continuum expressions. Overall, hydrogen concentrations up to 1.0 % do not have a significant effect on the elastic parameters and dislocation energy. The methods and insights derived in this work have the potential to facilitate the calculation of dislocation energies in a wide range of systems, and to guide our understanding of hydrogen embrittlement.

Alloyed systems↗
DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mesoscale informed parameter estimation through machine learning: A case-study in fracture modeling

Scale bridging is a critical need in computational sciences, where the modeling community has developed accurate physics models from first principles, of processes at lower length and time scales that influence the behavior at the higher scales of interest. However, it is not computationally feasible to incorporate all of the lower length scale physics directly into upscaled models. This is an area where machine learning has shown promise in building emulators of the lower length scale models, which incur a mere fraction of the computational cost of the original higher fidelity models. We demonstrate the use of machine learning using an example in materials science estimating continuum scale parameters by emulating, with uncertainties, complicated mesoscale physics. Additionally, we describe a new framework to emulate the fine scale physics, especially in the presence of microstructures, using machine learning, and showcase its usefulness by providing an example from modeling fracture propagation. Our approach can be thought of as a data-driven dimension reduction technique that yields probabilistic emulators. Our results show well-calibrated predictions for the quantities of interests in a low-strain simulation of fracture propagation at the mesoscale level. Furthermore, on average, we achieve ~10% relative errors on time-varying quantities like total damage and maximum stresses. Successfully replicating mesoscale scale physics within the continuum models is a crucial step towards predictive capability in multi-scale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗
DOE OSTI2022

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A improved equation of state for Xe gas bubbles in gamma U-Mo fuels

A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.

Beeler, Benjamin↗
DOE OSTI2021

Temperature-dependent model of helium bubble expansion and bursting in tungsten

Continuum-scale models that can reliably predict the behavior of helium in tungsten are of interest to the fusion community due to the projected impact of these materials on fusion reactor operation. In this work, we perform molecular dynamics simulations of spherical helium bubbles of various sizes in tungsten at different temperatures and depths with the goal of determining a mathematical model of the pressure and density at which the bubbles initially expand or burst as a function of depth, size, temperature, and surface orientation. The bubbles are small enough that their loop-punching pressures cannot be accurately predicted with continuum mechanics, and their expansion behavior is important, as it appears to cause many of the features observed on helium-irradiated tungsten surfaces. We vary the temperature, bubble size, bubble depth, and surface orientation in each case, recording the bubble pressure and density that result in bubble expansion. An exponential function with three adjustable parameters is found to fit the results well; the parameters that best fit our results are provided.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Resonant ultrasound spectroscopy measurement and modeling of additively manufactured octet truss lattice cubes

In this work, the elastic characteristics of additively manufactured Ti-5553 octet truss lattice cubes were examined using resonant ultrasound spectroscopy (RUS) techniques and correlated to X-ray computed tomography measurements of strut continuity. RUS models and measurements identified an outlier lattice possessing lower effective modulus despite larger effective density. Here, experimental measurements of density and effective Young's modulus were treated as input parameters in a homogenous anisotropic continuum model. The continuum model exhibited good agreement with experimental RUS measurements, thus confirming the potential for a simplified approximation of the octet truss lattice.

42 ENGINEERING↗
DOE OSTI2022

Bound states and electromagnetic radiation of relativistically rotating cylindrical wells

Here, we compute the effect of rigid rotation on the non-relativistic bound states. The energy levels of the bound states increase with the angular velocity of rotation until at certain value of the angular velocity they are completely pushed out into the continuum which corresponds to dissociation of the bound states. When the angular velocity exceeds the critical value at which the ground state disappears into the continuum, no bound state is possible. This effect should have important consequences for the phenomenology of the quark-gluon plasma. One of the ways to study it experimentally is to observe the electromagnetic radiation emitted by a rotating bound state. We compute the corresponding intensity of electromagnetic radiation and show that it strongly depends on the angular velocity of rotation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Ab initio calculations of scattering and fusion reactions

We review recent advancements achieved within the ab initio no-core shell model with continuum (NCSMC), a unified first-principle (or, ab initio) approach to nuclear bound and continuum states. We highlight its application to reactions of astrophysical interests, including the 7 Be(p, γ) 8 B solar fusion and α-α scattering, and to the structure of the exotic 12 Be nucleus. Furthermore, our findings provide a solid foundation for integrating ab initio calculations with experimental measurements, allowing for more precise evaluations of key thermonuclear reaction rates at astrophysical energies, and offer new insights into the factors that contribute to the emergence of clustering.

Quaglioni, Sofia [Lawrence Livermore National Labo↗
DOE OSTI2021

Non-perturbative BRST symmetry and the spectral structure of the ghost propagator

In BRST-quantised Yang-Mills theory the existence of BRST symmetry imposes significant constraints on the analytic structure of the continuum theory. In particular, the presence of this symmetry in the non-perturbative regime implies that any on-shell state with vanishing norm must have an associated partner state with identical mass, but negative inner product. In this work we demonstrate that for quantum chromodynamics (QCD) this constraint gives rise to an interconnection between the ghost and gluon spectra, and in doing so provides a non-trivial test of whether BRST symmetry is realised non-perturbatively. By analysing infrared lattice data for the minimal Landau gauge ghost propagator in pure SU(3) Yang-Mills theory, and comparing this with previous results for the gluon propagator, we show that this interconnection is violated, and hence conclude that continuum and current lattice formulations of Yang-Mills theory in Landau gauge represent two distinct realisations of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The effect of hydrogen-bonding on flavin’s infrared absorption spectrum

Cluster and continuum solvation computational models are employed to model the effect of hydrogen bonding interactions on the vibrational modes of lumiflavin. Calculated spectra were compared to experimental Fourier-transform infrared (FTIR) spectra in the diagnostic 1450–1800 cm –1 range, where intense V C=C , V C=N , V C 2 =O , and V C 4 =O stretching modes of flavin’s isoalloxazine ring are found. Local mode analysis is used to describe the strength of hydrogen-bonding in cluster models. The computations indicate that V C=C and V C=N mode frequencies are relatively insensitive to intermolecular interactions while the V C 2 =O and V C 4 =O modes are sensitive to direct (and also indirect for V C 2 =O ) hydrogen-bonding interactions. Although flavin is neutral, basis sets without the diffuse functions provide incorrect relative frequencies and intensities. The 6-31+G* basis set is found to be adequate for this system, and there is limited benefit to considering larger basis sets. Calculated vibrational mode frequencies agree with experimentally determined frequencies in solution when cluster models with multiple water molecules are used. Accurate simulation of relative FTIR band intensities, on the other hand, requires a continuum (or possibly quantum mechanical/molecular mechanical) model that accounts for long-range electrostatic effects. Lastly, an experimental peak at ca. 1624 cm –1 that is typically assigned to the V C 2 =O vibrational stretching mode has a complicated shape that suggests multiple underlying contributions. Our calculations show that this band has contributions from both the C 6 –C 7 and C 2 = O stretching vibrations.

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

Heat transfer coefficients of moving particle beds from flow-dependent thermal conductivity and near-wall resistance

Accurate determination of heat transfer coefficients for flowing packed particle beds is essential to the design of particle heat exchangers and other thermal and thermochemical equipment. While such dense granular flows mostly fall into the well-known plug-flow regime, the discrete nature of granular materials alters the thermal transport processes in both the near-wall and bulk regions of flowing particle beds from their stationary counterparts. As a result, heat transfer correlations based on the stationary particle bed thermal conductivity could be inadequate for flowing particles in a heat exchanger. Most earlier works have achieved a reasonable agreement with experiments by treating granular heat transfer media as a plug-flow continuum with a near-wall thermal resistance in series. However, the thermal conductivity values of the continuum were often obtained from measurements on stationary beds owing to the difficulty of flowing bed measurements. In this work, it was found that the properties of a stationary bed are highly sensitive to the method of particle packing and there is a decrease in the particle bed thermal conductivity and increase in the near-wall thermal resistance, measured as an effective air gap thickness, on the onset of particle flow. These variations in thermal conductivity of stationary and flowing particle beds can lead to errors in heat transfer coefficient calculations. Therefore, the heat transfer coefficients for granular flows were calculated using experimentally determined flowing particle bed thermal conductivity and near-wall air gap for ceramic particles – CARBO CP 40/100 (mean diameter = 275 µm), HSP 40/70 (404 µm) and HSP 16/30 (956 µm); at velocities of 5–15 mm·s –1 ; and temperatures of 300–650 °C. The thermal conductivity and air gap values for CP 40/100 and HSP 40/70 were further used to calculate heat transfer coefficients across different particle bed temperatures and velocities for different parallel-plate heat exchanger dimensions. Furthermore, these calculations, which show good agreement with measured HTC values reported in literature, can be used as a guide for heat exchanger designs. Graphical abstract

14 SOLAR ENERGY↗
DOE OSTI2021

Establishing a unified framework for ion solvation and transport in liquid and solid electrolytes

Electrolytes used in rechargeable batteries must enable rapid translation of the working ion between macroscopically separated electrodes. These electrolytes are, however, usually designed and synthesized using atomic-level insights. Whether the ideal electrolyte for a particular battery is a solid or a liquid remains an important unresolved question, especially as solids with conductivities comparable with liquids are discovered. Here, to help resolve such questions, we present the first steps toward a unified framework for relating atomic and continuum scale phenomena. Solvation shells in liquids are entities that translate with the working ion for a short while before they break up due to Brownian motion. By contrast, solvation cages in classical solids and polymers cannot not translate with the working ion. Mobility of the entities that make up the cages and shells, which is quantified by an order parameter, is shown to influence translation of the working ion on continuum length scales.

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

Differential localization patterns of pyruvate kinase isoforms in murine naïve, formative, and primed pluripotent states

Highlights: • PKM1/2 protein abundance is greater in formative mEpiLCs compared to naïve mESCs or primed mEpiSCs. The ratio of PKM1/2 is maintained. • PKM1/2 are localized in both nuclear and cytoplasmic regions relative to GAPDH and OCT4 across the pluripotent continuum. • PKM1 localization is strongly correlated to OCT4 localization, and moderately to GAPDH in formative mEpiLCs. • PKM1/2 localization is moderately correlated to OCT4 and GAPDH localization in naïve mESCs. • PKM1/2 localization are moderately correlated to GAPDH localization in primed mEpiSCs. Mouse embryonic stem cells (mESCs) and mouse epiblast stem cells (mEpiSCs) represent opposite ends of the pluripotency continuum, referred to as naïve and primed pluripotent states, respectively. These divergent pluripotent states differ in several ways, including growth factor requirements, transcription factor expression, DNA methylation patterns, and metabolic profiles. Naïve cells employ both glycolysis and oxidative phosphorylation (OXPHOS), whereas primed cells preferentially utilize aerobic glycolysis, a trait shared with cancer cells referred to as the Warburg Effect. Until recently, metabolism has been regarded as a by-product of cell fate, however, evidence now supports metabolism as being a driver of stem cell state and fate decisions. Pyruvate kinase muscle isoforms (PKM1 and PKM2) are important for generating and maintaining pluripotent stem cells (PSCs) and mediating the Warburg Effect. Both isoforms catalyze the final, rate limiting step of glycolysis, generating adenosine triphosphate and pyruvate, however, the precise role(s) of PKM1/2 in naïve and primed pluripotency is not well understood. The primary objective of this study was to characterize the cellular expression and localization patterns of PKM1 and PKM2 in mESCs, chemically transitioned epiblast-like cells (mEpiLCs) representing formative pluripotency, and mEpiSCs using immunoblotting and confocal microscopy. The results indicate that PKM1 and PKM2 are not only localized to the cytoplasm, but also accumulate in differential subnuclear regions of mESC, mEpiLCs, and mEpiSCs as determined by a quantitative confocal microscopy employing orthogonal projections and airyscan processing. Importantly, we discovered that the subnuclear localization of PKM1/2 changes during the transition from mESCs, mEpiLCs, and mEpiSCs. Finally, we have comprehensively validated the appropriateness and power of the Pearson's correlation coefficient and Manders's overlap coefficient for assessing nuclear and cytoplasmic protein colocalization in PSCs by immunofluorescence confocal microscopy. We propose that nuclear PKM1/2 may assist with distinct pluripotency state maintenance and lineage priming by non-canonical mechanisms. These results advance our understanding of the overall mechanisms controlling naïve, formative, and primed pluripotency.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A seamless multiscale operator neural network for inferring bubble dynamics

Modelling multiscale systems from nanoscale to macroscale requires the use of atomistic and continuum methods and, correspondingly, different computer codes. Here, we develop a seamless method based on DeepONet, which is a composite deep neural network (a branch and a trunk network) for regressing operators. In particular, we consider bubble growth dynamics, and we model tiny bubbles of initial size from 100 nm to 10 $\mathrm {\mu }\textrm {m}$ , modelled by the Rayleigh–Plesset equation in the continuum regime above 1 $\mathrm {\mu }\textrm {m}$ and the dissipative particle dynamics method for bubbles below 1 $\mathrm {\mu }\textrm {m}$ in the atomistic regime. After an offline training based on data from both regimes, DeepONet can make accurate predictions of bubble growth on-the-fly (within a fraction of a second) across four orders of magnitude difference in spatial scales and two orders of magnitude in temporal scales. The framework of DeepONet is general and can be used for unifying physical models of different scales in diverse multiscale applications.

Mechanics↗
DOE OSTI2022

Simulating ice-shelf extent using damage mechanics

Abstract Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.

54 ENVIRONMENTAL SCIENCES↗