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 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Multi-scale modeling of the electric field assisted sintering process

The electric field assisted sintering (EFAS) process involves tightly coupled physics that influence microstructural evolution in the particles being compacted. It is also an inherently multi-scale phenomenon, with the microstructure of the compact influencing the subsequent engineering-scale response of the sintering system. To improve understanding of how processing parameters influence microstructural evolution, we have developed a multi-scale modeling approach that couples a continuum-level model of the sintering system with a phase-field model for microstructural evolution of particles within the compact. The phase-field model couples the effect of chemical and electrical driving forces on microstructural evolution and includes the effect of charged defect segregation to surfaces and grain boundaries; this segregation leads to enhanced defect transport and heat generation at these interfaces in response to applied electric field. The effect of enhanced heat generation on particle neck growth and the influence of microstructural evolution on the engineering-scale model are demonstrated.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Multi-scale modeling of the electric field assisted sintering process

36 - MATERIALS SCIENCE↗

DOE OSTI2022

Advanced Physics Models for Particle-Particle Interactions (UCD final report)

The LDRD "Advanced Physics Models for Particle-Particle Interactions” collaborative research with Professor Jennifer Sinclair Curtis and PhD student Chanh Nguyen at the University of California, Davis involved meso-scale, discrete element method (DEM) modeling of micron-sized particles impacting a substrate. Using DEM, physical insight into the details of particle response upon impact can be obtained. Such details provide the type of information necessary to enhance existing continuum or macro-scale simulation tools so that modeling for material strength and particle breakage models is based on particle physics versus empirical relationships. While molecular dynamics simulations of atoms can also provide physical insight into the particle response to impact, the scale of particles that are treated is on the order of angstroms, not microns. This study gives evidence to the ability of the DEM simulation technique to provide such fundamental, detailed insight for micron-sized particles as a function of particle properties and impact velocity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

DOE OSTI2023

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗

DOE OSTI2022

Tree hydrodynamic modelling of the soil–plant–atmosphere continuum using FETCH3

Modelling the water transport along the soil–plant–atmosphere continuum is fundamental to estimating and predicting transpiration fluxes. A Finite-difference Ecosystem-scale Tree Crown Hydrodynamics model (FETCH3) for the water fluxes across the soil–plant–atmosphere continuum is presented here. The model combines the water transport pathways into one vertical dimension, and assumes that the water flow through the soil, roots, and above-ground xylem can be approximated as flow in porous media. This results in a system of three partial differential equations, resembling the Richardson–Richards equation, describing the transport of water through the plant system and with additional terms representing sinks and sources for the transfer of water from the soil to the roots and from the leaves to the atmosphere. The numerical scheme, developed in Python 3, was tested against exact analytical solutions for steady state and transient conditions using simplified but realistic model parameterizations. The model was also used to simulate a previously published case study, where observed transpiration rates were available, to evaluate model performance. With the same model setup as the published case study, FETCH3 results were in agreement with observations. Through a rigorous coupling of soil, root xylem, and stem xylem, FETCH3 can account for variable water capacitance, while conserving mass and the continuity of the water potential between these three layers. FETCH3 provides a ready-to-use open access numerical model for the simulation of water fluxes across the soil–plant–atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗

DOE OSTI2020

Moiré quantum chemistry: Charge transfer in transition metal dichalcogenide superlattices

Transition metal dichalcogenide (TMD) bilayers have recently emerged as a robust and tunable moiré system for studying and designing correlated electron physics. In this Rapid Communication, by combining a large-scale first-principles calculation and continuum model approach, we provide an electronic structure theory that maps long-period TMD heterobilayer superlattices onto diatomic crystals with cations and anions. We find that the interplay between the moiré potential and Coulomb interaction leads to filling-dependent charge transfer between different moiré superlattice regions. Further, we show that the insulating state at half filling found in recent experiments on WSe 2 /WS 2 is a charge-transfer insulator rather than a Mott-Hubbard insulator. Our work reveals the richness of simplicity in moiré quantum chemistry.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Electronic structures, charge transfer, and charge order in twisted transition metal dichalcogenide bilayers

Moiré superlattices of transition metal dichalcogenide (TMD) bilayers have been shown to host correlated electronic states, which arise from the interplay of long wavelength moiré potential and long-range Coulomb interaction. Here, we theoretically investigate structural relaxation and single-particle electronic structure of twisted TMD homobilayer. From the large-scale density functional theory calculation and continuum model with layer degrees of freedom, we find that the out-of-plane gating field creates a tunable charge transfer gap at the Dirac point between the first and second moiré valence bands. We further study the charge orders at the fractional band fillings. In the flat band limit, we find from Monte Carlo simulations a series of charge-ordered insulating states at various fillings $n = \frac{1}{4}, \frac{1}{3}, \frac{1}{2}, \frac{2}{3}, 1$. Further, we predict that gating field induces a phase transition between different electron crystals at fixed filling $n = \frac{1}{2}$ or $\frac{2}{3}$. At half-filling $\textit{n} = 1$, the ground state is a Mott insulator with electronically driven ferroelectricity. This work demonstrates that TMD homobilayer provides a powerful platform for the investigation of tunable charge transfer insulator and charge orders.

36 MATERIALS SCIENCE↗

DOE OSTI2022

mesoflow [SWR-22-56]

Mesoflow is a continuum scale simulation tool developed specifically for modeling transport and chemistry at the mesoscale. Our solver utilizes Cartesian block-structured adaptive mesh refinement to resolve complex surface morphologies (of catalysts/biomass particles among others) directly obtained from X-ray tomography data. An immersed boundary based formulation enables rapid representation of complex geometries prevalent in most mesoporous interfaces. The solver is developed on top of open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Our flexible software framework enables integration of complex chemical mechanisms at heterogenous interfaces and time-split algorithms for circumventing highly disparate reaction and flow time-scales. Our current studies indicate a ten-fold performance gain by using graphics-processing-units (GPU) compared to a single processor for representative problem sizes (2 million cell mesh).

Sitaraman, Hariswaran↗

DOE OSTI2023

Data Driven Approach to Dislocation-Based Plasticity Models of Face-Centered Cubic Metals

Dislocation dynamics controls plastic deformation, mechanical strength, and failure of crystalline materials. It also governs fatigue resistance under cyclic loading, creep resistance at elevated-temperature, and radiation resistance for reactor applications. There is a compelling need for understanding fundamental dislocation mechanisms for deformation because virtually all structural metals used in energy systems are fabricated to desired forms and shapes by deformation processes. To date, the most outstanding problem in a physics-based multiscale model of crystal plasticity is the lack of quantitative connections between continuum plasticity (CP) models with the lower scale dislocation models. As a result, existing CP models used in engineering applications are still phenomenological, while evidence continues to mount that they can make inaccurate predictions under realistically complex scenarios. This project takes advantage of the recent advances in high-performance discrete dislocation dynamics (DDD) simulations and data science approaches to establish the first fully connected multiscale plasticity model for pure face-centered cubic (FCC) single crystals.

36 MATERIALS SCIENCE↗

NASA NTRS1994

Pre- and postperihelon abundances of gas and dust in comet Halley

Photometrically calibrated spectra of comet P/Halley (1986 III) were recorded between 1985 September 12-1986 June 10 using the Ohio State University Image Dissector Scanner on the Perkins 72 inch telescope at the Lowell Observatory. Column densities of CN, C3, CH, C2, and NH2 were calculated from measured fluxes in these spectra, and molecular scale lengths were deduced from the radial distribution of CN, C3, C2, and NH2. Production rates were computed using the new scale lengths and a Haser model analysis. Continuum emission at 4260 A was used to derive gas-to-dust ratios. The data indicate than comet Halley was approximately 2-5 times more abundant in gas and dust at postperihelion than preperihelion. On 1986 June 8 we observed the onset of a cometary ourburst which appeared very strong in dust production. The gas-to-dust ratios appeared to be subject to changes as a result of short-term outbursts but otherwise did not exhibit any systematic dependence on heliocentric distance. Reflectivity gradients of the continuum were also measured from the spectra. While most of the continua were red, blue continua were also observed which may be correlated with dust outbursts.

Womack, Maria↗

NASA NTRS2003

Mesoscale Computational Investigation of Shocked Heterogeneous Materials with Application to Large Impact Craters

The propagation of shock waves through target materials is strongly influenced by the presence of small-scale structure, fractures, physical and chemical heterogeneities. Pre-existing fractures often create craters that appear square in outline (e.g. Meteor Crater). Reverberations behind the shock from the presence of physical heterogeneity have been proposed as a mechanism for transient weakening of target materials. Pre-existing fractures can also affect melt generation. In this study, we are attempting to bridge the gap in numerical modeling between the micro-scale and the continuum, the so-called meso-scale. To accomplish this, we are developing a methodology to be used in the shock physics hydrocode (CTH) using Monte-Carlo-type methods to investigate the shock properties of heterogeneous materials. By comparing the results of numerical experiments at the micro-scale with experimental results and by using statistical techniques to evaluate the performance of simple constitutive models, we hope to embed the effect of physical heterogeneity into the field variables (pressure, stress, density, velocity) allowing us to directly imprint the effects of micro-scale heterogeneity at the continuum level without incurring high computational cost.

Crawford, D. A.↗

NASA NTRS1999

Photodissociation Regions in the Interstellar Medium of Galaxies

The interstellar medium of galaxies is the reservoir out of which stars are born and into which stars inject newly created elements as they age. The physical properties of the interstellar medium are governed in part by the radiation emitted by these stars. Far-ultraviolet (6 eV less than h(nu) less than 13.6 eV) photons from massive stars dominate the heating and influence the chemistry of the neutral atomic gas and much of the molecular gas in galaxies. Predominantly neutral regions of the interstellar medium in which the heating and chemistry are regulated by far ultraviolet photons are termed Photo-Dissociation Regions (PDRs). These regions are the origin of most of the non-stellar infrared (IR) and the millimeter and submillimeter CO emission from galaxies. The importance of PDRs has become increasingly apparent with advances in IR and submillimeter astronomy. The IR emission from PDRs includes fine structure lines of C, C+, and O; rovibrational lines of H2, rotational lines of CO; broad middle features of polycyclic aromatic hydrocarbons; and a luminous underlying IR continuum from interstellar dust. The transition of H to H2 and C+ to CO occurs within PDRs. Comparison of observations with theoretical models of PDRs enables one to determine the density and temperature structure, the elemental abundances, the level of ionization, and the radiation field. PDR models have been applied to interstellar clouds near massive stars, planetary nebulae, red giant outflows, photoevaporating planetary disks around newly formed stars, diffuse clouds, the neutral intercloud medium, and molecular clouds in the interstellar radiation field-in summary, much of the interstellar medium in galaxies. Theoretical PDR models explain the observed correlations of the [CII] 158 microns with the COJ = 1-0 emission, the COJ = 1-0 luminosity with the interstellar molecular mass, and the [CII] 158 microns plus [OI] 63 microns luminosity with the IR continuum luminosity. On a more global scale, MR models predict the existence of two stable neutral phases of the interstellar medium, elucidate the formation and destruction of star-forming molecular clouds, and suggest radiation-induced feedback mechanisms that may regulate star formation rates and the column density of gas through giant molecular clouds.

Hollenbach, David J.↗

NASA NTRS2000

Photodissociation Regions

The interstellar medium of galaxies is the reservoir out of which stars are born and into which stars inject newly created elements as they age. The physical properties of the interstellar medium are governed in part by the radiation emitted by these stars. Far-ultraviolet (6 eV< hNu < 13.6 eV) photons from massive stars dominate the heating and influence the chemistry of the neutral atomic gas and much of the molecular gas in galaxies. Predominantly neutral regions of the interstellar medium in which the heating and chemistry are regulated by far ultraviolet photons are termed Photodissociation Regions (PDRs). These regions are the origin of most of the non-stellar infrared (IR) and the millimeter and submillimeter CO emission from galaxies. The importance of PDRs has become increasingly apparent with the advances in IR and submillimeter astronomy. The IR emission from PDRs includes fine structure lines of C, C(+) and O; rovibrational lines of H2; rotational lines of CO; broad mid-IR features of polycyclic aromatic hydrocarbons; and a luminous underlying IR continuum from interstellar dust. The transition of H to H2 and C(+) to CO occurs within PDRs. Comparison of observations with theoretical models of PDRs enables one to determine the density and temperature structure, the elemental abundances, the level of ionization, and the radiation field. PDR models have been applied to interstellar clouds near massive stars, planetary nebulae, red giant outflows, photoevaporating planetary disks around newly formed stars, diffuse clouds, the neutral intercloud medium, and molecular clouds in the interstellar radiation field-in summary, much of the interstellar medium in galaxies. Theoretical PDR models explain the observed correlations of the [CII] 158, micrometers with the CO J=1-0 emission, the CO J=1-0 luminosity with the interstellar molecular mass, and the [CII] 158 micrometers plus [OI] 63 micrometers luminosity with the IR continuum luminosity. On a more global scale, PDR models predict the existence of two stable neutral phases of the interstellar medium, elucidate the formation and destruction of star-forming molecular clouds, and suggest radiation-induced feedback mechanisms that may regulate star formation rates and the column density of gas through giant molecular clouds.

Hollenbach, David J.↗

DOE OSTI2025

Accelerating magnonic simulations with the pseudospectral Landau-Lifshitz equation

The pseudospectral Landau-Lifshitz (PS-LL) model can describe atomic-scale magnetic exchange interactions within a continuum framework. This is achieved by employing a convolution kernel that models the nonlocal interaction in a grid-independent manner. Even though the PS-LL was originally introduced to address atomic exchange, any nonlocal kernel can be modeled. In the field of magnonics, the dipole field is fundamental to describe the dispersion relation of magnons, the quasiparticle representation of angular momentum. Because dipole-dipole interactions are long-range, numerical approaches typically rely on convolutions. Here, we demonstrate that the PS-LL model can be used to perform magnonic simulations with a single convolution kernel derived from analytical solutions. We demonstrate a twofold increase in computational speed compared with the full dipole calculation. This approach is valid insofar as the excitations are linear, which is typically the case for magnons. Our results have the potential to accelerate magnonic research, particularly for the inverse design method, where several simulations must be performed to achieve the desired outcome.

Mathematics and computing↗

DOE OSTI2024

Pseudospectral Landau-Lifshitz description of magnetization dynamics

Magnetic materials host a wealth of nonlinear dynamics, textures, and topological defects. This is possible due to the competition between strong nonlinearity and dispersion, which act at the atomic scale, as well as long-range interactions. However, these features are difficult to study analytically and numerically because of the vastly different temporal and spatial scales involved. Here, in this work, we present a pseudospectral approach for the Landau-Lifshitz equation that invokes energy and momentum conservation embodied in the magnon dispersion relation to accurately describe both atomic and continuum limits. Furthermore, this approach enables analytical study at every scale. We show the applicability of this model in both the continuum and atomic limit by investigating modulational instability and ultrafast evolution of magnetization due to transient grating, respectively, in a one-dimensional ferromagnetic chain with perpendicular magnetic anisotropy. This model provides the possibility of grid-independent multiscale numerical approaches that will enable the description of singularities within a single framework.

1-dimensional spin chains↗