SEARCH · Engineering Papers

Results for “continuum scale modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗

NASA NTRS1986

A comparison of high-resolution radio and optical data for the northeast region of the Cygnus Loop

Radio observations of the northeastern part of the Cygnus Loop obtained with the VLA at 18 cm with a resolution of 4.7 arcsec are presented and discussed. In general, the brightest radio filaments match the brightest optical filaments very well, particularly those which are bright in H-alpha. The agreement in position, morphology, and relative brightness is particularly good where all four optical species are bright. However, a few bright radio features have no bright optical counterparts. The general radio and optical structure of the region can be explained by a model in which shocked interstellar clouds form filaments due to unstable cooling and compression. The features appear to be ropelike, although it is not possible to rule out that they are irregular sheets seen edge-on. A large gradation in cooling on a scale of parsecs is evidenced by the sequence of the optical species in emission, with a corresponding increase in density generally evidenced by the continuum radio emission.

Straka, W. C.↗

DOE OSTI2021

Computational design of microarchitected porous electrodes for redox flow batteries

Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.

25 ENERGY STORAGE↗

DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗

DOE OSTI2022

A parallel discrete dislocation dynamics/kinetic Monte Carlo method to study non-conservative plastic processes

Non-conservative processes play a fundamental role in plasticity and are behind important macroscopic phenomena such as creep, dynamic strain aging, loop raft formation, etc. In the most general case, vacancy-induced dislocation climb is the operating unit mechanism. While dislocation/vacancy interactions have been modeled in the literature using a variety of methods, the approaches developed rely on continuum descriptions of both the vacancy population and its fluxes. However, there are numerous situations in physics where point defect populations display heterogeneous concentrations and/or non-smooth kinetics. Here, a kinetic Monte Carlo (kMC) approach for modeling vacancy transport in response to arbitrary stress fields is used. Vacancies are treated as point particles and are coupled to the dislocation substructure representing a deformed material via an advection term defined by the local stress gradients. The stress fields and the dislocation substructure are evolved using a discrete dislocation dynamics (DDD) module. To extend the coupled model to the treatment of large systems, we have implemented it in the massively-parallel DDD code ParaDiS. To avoid numerical incompatibilities associated with merging deterministic (DDD) and stochastic (kMC) integration algorithms, we cast the entire elasto-plastic-diffusive problem within a single stochastic framework, taking advantage of a parallel kMC algorithm to evolve the system as a single event-driven process. The large-scale implementation enables the study of the evolution of a variety of dislocation-defect scenarios governed by non-conservative transport kinetics. After carrying out an exhaustive numerical and computational analysis of our parallel algorithm, we show results that emphasize situations where inhomogeneous vacancy dynamics are of relevance, and compare discrete kinetics to continuum solutions for several cases.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗

DOE OSTI2023

Habitat quality influences trade-offs in animal movement along the exploration–exploitation continuum

Abstract To successfully establish itself in a novel environment, an animal must make an inherent trade-off between knowledge accumulation and exploitation of knowledge gained (i.e., the exploration–exploitation dilemma). To evaluate how habitat quality affects the spatio-temporal scale of switching between exploration and exploitation during home range establishment, we conducted experimental trials comparing resource selection and space-use of translocated animals to those of reference individuals using reciprocal translocations between habitat types of differing quality. We selected wild pigs ( Sus scrofa ) as a model species to investigate hypotheses related to the movement behavior of translocated individuals because they are globally distributed large mammals that are often translocated within their introduced range to facilitate recreational hunting. Individuals translocated to higher quality habitat (i.e. higher proportions of bottomland hardwood habitats) exhibited smaller exploratory movements and began exploiting resources more quickly than those introduced to lower quality areas, although those in lower-quality areas demonstrated an increased rate of selection for preferred habitat as they gained knowledge of the landscape. Our data demonstrate that habitat quality mediates the spatial and temporal scale at which animals respond behaviorally to novel environments, and how these processes may determine the success of population establishment.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗

NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗

DOE OSTI2021

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2023

Cristae formation is a mechanical buckling event controlled by the inner mitochondrial membrane lipidome

Cristae are high–curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae–shaping proteins have been defined, analogous lipid–based mechanisms have yet to be elucidated. Here, we combine experimental lipidome dissection with multi–scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the inner mitochondrial membrane against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein–mediated curvatures. This model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that cardiolipin is essential in low–oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of cardiolipin is dependent on the surrounding lipid and protein components of the IMM.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2019

The First Swift Intensive AGN Accretion Disk Reverberation Mapping Survey

Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled ∼200–350 times in 0.3–10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these data sets yields three main results: (1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. (2) The UV/optical interband lags are generally consistent with t μ l4 3 as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. (3) The one exception is the U band lags, which are on average a factor of ∼2.2 larger than predicted from the surrounding band data and fits. This excess appears to be due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by JAVELIN versus more standard lag measurement techniques, with JAVELIN reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein.

R. Edelson↗

DOE OSTI2021

Ergodic and nonergodic many-body dynamics in strongly nonlinear lattices

The study of nonlinear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi et al. [Los Alamos Scientific Laboratory Report No. LA-1940, 1955 (unpublished)]. Here, we introduce a family of such systems which consist of chains of N harmonically coupled particles with the nonlinearity introduced by confining the motion of each individual particle to a box or stadium with hard walls. The stadia are arranged on a one-dimensional lattice but they individually do not have to be one dimensional, thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of N . Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh Gordon theory, we do not see any evidence for the kind of nonergodicity famously seen in the work of Fermi et al. Finally, we examine the chain with particles confined to two-dimensional stadia where the individual stadium is already chaotic and find a much more chaotic phase space at small system sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1985

Vertical cloud structure of Jupiter's equatorial plumes

Multiple-scattering radiative transfer calculations were used to deduce the vertical cloud structure (VCS) of Jupiter's equatorial region. The VCS model of the equatorial plumes is obtained through an analysis of Voyager images of the 6190-A methane band and the 6000-A continuum, and ground-based 8900-A methane band images. The VCS of the equatorial plumes is found to be consistent with the hypothesis that the plumes are caused by upwelling at the ammonia condensation level produced by buoyancy due to latent heat release from the condensation of water clouds nearly three scale heights below the plumes.

Stoker, C. R.↗

NASA NTRS2018

3D Mapping of the Neutral X-Ray Absorption in the Local Interstellar Medium: The Gaia and XMM-Newton Synergy

We present a three-dimensional map of the hydrogen density distribution in the Galactic interstellar medium. The hydrogen-equivalent column densities were obtained from the Exploring the X-ray Transient and variable Sky project (EXTRAS) which provides equivalent NH values from X-ray spectral fits of observations within the XMM-Newton Data Release. EXTRAS include multiple fits for each source, allowing an accurate determination of the equivalent column densities, which depends on the continuum modelling of the spectra. A cross-correlation between the EXTRAS catalogue and the first Gaia Data Release was performed in order to obtain accurate parallax and distance measurements. We use a Bayesian method explained in Rezaei Kh. et al. (2017) in order to predict the most probable distribution of the density at any arbitrary point, even for lines of sight along which there are no initial observation. The resulting map shows small-scale density structures which could not have been modelled by using analytic density profiles. In this paper, we present a proof of concept of the kind of science possible with the synergy of these catalogues. However, given the systematic uncertainties connected to the source identification and to the dependence of NH on the spectral model, the present maps should be considered qualitatively at this point.

Gatuzz, Efrain↗

NASA NTRS2021

Leveraging the ALMA Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide

We report the first survey of molecular emission from cometary volatiles using standalone Atacama Compact Array (ACA) observations from the Atacama Large Millimeter/Submillimeter Array (ALMA) toward comet C/2015 ER61 (PanSTARRS) carried out on UT 2017 April 11 and 15, shortly after its April 4 outburst. These measurements of HCN, CS, CH3OH, H2CO, and HNC (along with continuum emission from dust) probed the inner coma of C/2015 ER61, revealing asymmetric outgassing and discerning parent from daughter/distributed source species. This work presents spectrally integrated flux maps, autocorrelation spectra, production rates, and parent scale lengths for each molecule and a stringent upper limit for CO. HCN is consistent with direct nucleus release in C/2015 ER61, whereas CS, H2CO, HNC, and potentially CH3OH are associated with distributed sources in the coma. Adopt- ing a Haser model, parent scale lengths determined for H2CO (Lp ∼ 2200 km) and HNC (Lp ∼ 3300 km) are consistent with previous work in comets, whereas significant extended source production (Lp ∼ 2000 km) is indicated for CS, suggesting production from an unknown parent in the coma. The con- tinuum presents a point-source distribution with a flux density implying an excessively large nucleus, inconsistent with other estimates of the nucleus size. It is best explained by the thermal emission of slowly moving outburst ejectas, with total mass 5–8 × 1010 kg. These results demonstrate the power of the ACA for revealing the abundances, spatial distributions, and locations of molecular production for volatiles in moderately bright comets such as C/2015 ER61.

Molecular spectroscopy (2095) — High resolution sp↗

DOE OSTI2024

Physics-Based Continuum Modeling for an Aqueous Rechargeable Zn/MnO 2 Battery

This study introduces a framework for modeling the aqueous Zn/MnO 2 rechargeable battery. A reaction system and a physics-based continuum model are proposed based on two reaction types, one involving insertion and the second related to dissolution and deposition of a solid reaction product. The model, fitted to empirical data, predicts voltage behavior and capacity limitations during cycling, identifying electrolytic zinc depletion as a limiting mechanism, depending on the original cell construction. The research suggests the need for further material characterization and reaction analysis, which will advance our understanding and facilitate the development of grid-scale energy storage solutions.

Electrochemistry↗

NASA NTRS1985

Metamorphic core complexes: Expression of crustal extension by ductile-brittle shearing of the geologic column

Metamorphic core complexes and detachment fault terranes in the American Southwest are products of stretching of continental crust in the Tertiary. The physical and geometric properties of the structures, fault rocks, and contact relationships that developed as a consequence of the extension are especially well displayed in southeastern Arizona. The structures and fault rocks, as a system, reflect a ductile-through-brittle continuum of deformation, with individual structures and faults rocks showing remarkably coordinated strain and displacement patterns. Careful mapping and analysis of the structural system has led to the realization that strain and displacement were partitioned across a host of structures, through a spectrum of scales, in rocks of progressively changing rheology. By integrating observations made in different parts of the extensional system, especially at different inferred depth levels, it has been possible to construct a descriptive/kinematic model of the progressive deformation that achieved continental crustal extension in general, and the development of metamorphic core complexes in particular.

Davis, G. H.↗