Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum 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.

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

At least 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗
DOE OSTI2026

Ab initio symmetry-adapted approaches to nuclear reactions

In this review, we discuss recent applications of the ab initio symmetry-adapted no-core shell-model (SA-NCSM) theory for study and prediction of structure and reactions of stable and unstable nuclei from light to medium mass range. We explore structure properties of neutron-rich He, Li, and Mg isotopes, with a focus on nuclear collectivity, clustering, and spectroscopic factors, as well as multi-particle excitations of utmost significance in the proximity of the drip lines. In addition, we present extensions of the SA-NCSM with continuum for determining the microscopic structure of reaction fragments, which enables calculations of reaction cross sections for targets from the lightest He to 40Ca, rooted in first principles. We illustrate this for neutron and proton elastic scattering, deuteron and alpha capture reactions, and alpha knock-out reactions. Furthermore, we discuss microscopic optical potentials with uncertainty quantification, a critical ingredient in many reaction models. As a result, we also discuss the impact of alpha clustering on reactions of significance to nuclear astrophysics, as well as on beta decays and beyond-the-standard-model physics.

Ab initio approaches to nuclear reactions↗
DOE OSTI2025

Wino and Real Minimal Dark Matter Excluded by Fermi Gamma-Ray Observations

We show that minimal, fermionic dark matter (DM) models in the $n$-dimensional representation of the weak force with zero hypercharge that make up 100% of the DM under the standard cosmological history are strongly excluded for $n < 9$. This includes the thermal wino, which we show is ruled out even allowing for DM core sizes up to $\sim$6.7 kpc with the preferred local DM density or $\sim$3.7 kpc in addition to the local DM density being half the preferred value, at less than $0.2$ GeV/cm$^3$. We reach these conclusions through dedicated searches with 14 years of Fermi gamma-ray data in the inner Galaxy between 30 GeV and 2 TeV for the continuum gamma-rays produced in the decays of unstable particles produced in DM annihilation and bound-state formation processes. We consider a variety of Milky Way DM profiles in our analyses, including those motivated by modern hydrodynamic cosmological simulations, and show that all the $n < 9$ minimal DM models are disfavored even under the most conservative assumptions for these density profiles. While wino, quintuplet ($n=5$), and $n = 7$ DM models are strongly disfavored by our analyses under the standard cosmology, we discuss how non-standard cosmological histories or DM sub-fractions could still allow for these particles to be realized in nature, with discovery opportunities at next-generation particle colliders and gamma-ray telescopes.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2025

A novel peridynamics-based approach to predict pharmaceutical tablet robustness

The pharmaceutical drug product development process can be greatly accelerated through the use of modeling and simulation techniques to predict the manufacturability and performance of a given formulation. The anticipation and possible mitigation of tablet damage due to manufacturing stresses represents a specific area of interest in the pharmaceutical industry for predicting formulation and tableting performance. While the finite element method (FEM) has been extensively used for predicting the mechanical behavior of powder material in the compaction processes, a shortcoming of the approach is the inherent difficulty to predict discontinuities (e.g., damage or cracking) within a tablet as FEM is a continuum-based approach. In this work, we propose a novel method utilizing peridynamics (PD), a numerical method that can capture discontinuities such as tablet fracture, to predict the evolution of damage and breakage in pharmaceutical tablets. The approach links (1) the finite element method – to elucidate the behavior of powders during die compaction – with (2) the peridynamics modeling technique – to model the discontinuous nature of damage and predict tablet breakage during the critical stages of unloading and ejection from the compression die. This short communication presents a proof of concept including a workflow to calibrate the linked FEM-PD simulation models. Further, it demonstrates promising results from a preliminary experimental validation of the approach. Following further development, this approach could be used to guide the optimization of compression processes through targeted changes to formulation material properties, compression process conditions, and/or tooling geometries to deliver improved process efficiency and tablet robustness.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-invertible symmetries and LSM-type constraints on a tensor product Hilbert space

We discuss the exact non-invertible Kramers-Wannier symmetry of 1+1d lattice models on a tensor product Hilbert space of qubits. This symmetry is associated with a topological defect and a conserved operator, and the latter can be presented as a matrix product operator. Importantly, unlike its continuum counterpart, the symmetry algebra involves lattice translations. Consequently, it is not described by a fusion category. In the presence of this defect, the symmetry algebra involving parity/time-reversal is realized projectively, which is reminiscent of an anomaly. Different Hamiltonians with the same lattice non-invertible symmetry can flow in their continuum limits to infinitely many different fusion categories (with different Frobenius-Schur indicators), including, as a special case, the Ising CFT. The non-invertible symmetry leads to a constraint similar to that of Lieb-Schultz-Mattis, implying that the system cannot have a unique gapped ground state. It is either in a gapless phase or in a gapped phase with three (or a multiple of three) ground states, associated with the spontaneous breaking of the lattice non-invertible symmetry.

Seiberg, Nathan (ORCID:000000033897046X)↗
DOE OSTI2026

Effects of screening and pressure ionization on the electron broadening of spectral lines in dense plasmas

Collisions between electrons and radiating atoms broaden spectral absorption and emission lines in dense plasmas. High densities also introduce screening and pressure ionization effects that distort the wave functions of both bound and free electrons. In order to study how dense plasma effects influence the electron broadening of spectral lines, this paper incorporates electron wave functions from an average-atom (AA) model to calculate the linewidth of the B III 2⁢𝑝−2⁢𝑠 transition at 𝑇 = 10 eV for mass densities ranging from 𝜌 = 10 −4 to 0.4⁢ g⁡/cc. The calculation method uses the impact approximation, allowing the linewidth to be written in terms of electron-collision cross sections and an interference term. Compared to an otherwise identical calculation that uses Coulomb free wave functions, the AA method is found to modify both the cross sections and the resulting linewidth at sufficiently high density by introducing screening and pressure-ionized bound states. Screening lowers the cross sections at low energies and near electron excitation thresholds, while pressure-ionized bound states introduce resonances into the continuum. Thus, as the density increases, the relative linewidth between the AA and Coulomb calculations follows a general decrease because of screening, with sharp increases at various intervals due to pressure ionization. Finally, the AA results are also compared with a common approach to introduce screening through the interaction potential and reduced models that use the Bethe formula for the inelastic electron-collision cross sections.

electronic excitation & ionization↗
DOE OSTI2020

Clues on the history of early-type galaxies from SDSS spectra and GALEX photometry

ABSTRACT Stellar population studies of early-type galaxies (ETGs) based on their optical stellar continuum suggest that these are quiescent systems. However, emission lines and ultraviolet photometry reveal a diverse population. We use a new version of the starlight spectral synthesis code and state-of-the-art stellar population models to simultaneously fit Sloan Digital Sky Survey spectra and Galaxy Evolution Explorer photometry for a sample of 3453 galaxies at z < 0.1 with near ultraviolet (NUV) − r > 5 that are classified as elliptical by Galaxy Zoo. We reproduce far ultraviolet (FUV) magnitudes for 80 per cent of UV upturn galaxies selected using criteria from the literature, suggesting that additional stellar population ingredients such as binaries and extreme horizontal branch stars may have a limited contribution to the UV upturn. The addition of ultraviolet data leads to a broadening of the distributions of mean stellar ages, metallicities, and attenuation. Stellar populations younger than $1\,$ Gyr are required to reproduce the ultraviolet emission in 17 per cent of our sample. These systems represent 43 per cent of the sample at 5 < NUV − r < 5.5 and span the same stellar mass range as other ETGs in our sample. ETGs with young stellar components have larger H α equivalent widths (WH α) and larger dust attenuation. Emission line ratios and WH α indicate that the ionizing source in these systems is a mixture of young and old stellar populations. Their young stellar populations are metal-poor, especially for high-mass galaxies, indicating recent star formation associated with rejuvenation events triggered by external processes, such as minor mergers.

Werle, A.↗
DOE OSTI2022

X 17 boson and the H 3 ( p , e + e – ) He 4 and He 3 ( n , e + e – ) He 4 processes: A theoretical analysis

The present work deals with e + – e – pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ( χ EFT ) . Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross section for such a process. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The ab initio calculations use exact hyperspherical-harmonics methods to describe the bound state and low-energy spectrum of the A = 4 continuum, and they fully account for initial state interaction effects in the 3 + 1 clusters. While electromagnetic interactions are treated to high orders in the chiral expansion, the interactions of the hypothetical boson with nucleons are modeled in leading-order χ EFT (albeit, in some instances, selected subleading contributions are also accounted for). We also provide an overview of possible future experiments probing pair production in the A = 4 system at a number of candidate facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Final Report for Tech-X Corporations’s contribution to Enhancing Understanding of High Energy Density Plasmas from Wire Array and Solid Liner Implosions Using Fluid Modeling with Kinetic Closures

Recent results from experiments and simulation of magnetically driven pulsed power liners have explored the role of the early-time electrothermal instability in the evolution of the magneto-Rayleigh-Taylor instability. Our focus will be on understanding the development of such instabilities and the potential stabilization mechanisms via electron thermal conduction, viscosity, and large magnetic fields which we expect could play a significant role in supporting the success of the MagLIF program. Experiments have shown that studies of high-energy density plasmas from wire-array implosions require physics modeling that goes well beyond simple models such as ideal magnetohydrodynamics. The goal of this work is to provide increased understanding of these experiments by employing simulations with a multi-fluid extended-MHD model which uses kinetic closures for thermal conduction, resistivity and viscosity. We will use codes easily available to the wider research community, including university students, with a secondary goal of providing the community with well-benchmarked tools capable of advanced modeling of high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Computationally efficient method for determining limiting velocities of edge dislocations in anisotropic crystals

The continuum-limit theory of dislocations in crystals predicts divergences in the elastic energy at crystal-geometry dependent limiting velocities vL, which separate subsonic, transsonic, and supersonic dislocation glide regimes and are therefore import for material strength models at high strain rates. Although it is known how to calculate those limiting velocities, there is one special case - edge dislocations with reflection symmetry, but non-vanishing elastic constants c16 or c26 - where previous methods have been notoriously slow. In this letter, we address this deficiency by deriving a computationally efficient method for determining the limiting velocities of edge dislocations with reflection symmetry which is two orders of magnitude faster than the previous method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Changing-look Blazar B2 1420+32

Abstract Blazars are active galactic nuclei with their relativistic jets pointing toward the observer, comprising two major subclasses, flat-spectrum radio quasars (FSRQs) and BL Lac objects. We present multiwavelength photometric and spectroscopic monitoring observations of the blazar B2 1420+32, focusing on its outbursts in 2018–2020. Multiepoch spectra show that the blazar exhibited large-scale spectral variability in both its continuum and line emission, accompanied by dramatic gamma-ray and optical variability by factors of up to 40 and 15, respectively, on week to month timescales. Over the last decade, the gamma-ray and optical fluxes increased by factors of 1500 and 100, respectively. B2 1420+32 was an FSRQ with broad emission lines in 1995. Following a series of flares starting in 2018, it transitioned between BL Lac and FSRQ states multiple times, with the emergence of a strong Fe pseudocontinuum. Two spectra also contain components that can be modeled as single-temperature blackbodies of 12,000 and 5200 K. Such a collection of “changing-look” features has never been observed previously in a blazar. We measure gamma-ray–optical and interband optical lags implying emission-region separations of less than 800 and 130 gravitational radii, respectively. Since most emission-line flux variations, except the Fe continuum, are within a factor of 2–3, the transitions between FSRQ and BL Lac classifications are mainly caused by the continuum variability. The large Fe continuum flux increase suggests the occurrence of dust sublimation releasing more Fe ions in the central engine and an energy transfer from the relativistic jet to subrelativistic emission components.

Astronomy & Astrophysics↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

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

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗
DOE OSTI2025

Electroexcitation of Nucleon Resonances and Emergence of Hadron Mass

Developing an understanding of phenomena driven by the emergence of hadron mass (EHM) is one of the most challenging problems in the Standard Model. This discussion focuses on the impact of results on nucleon resonance (𝑁 ∗ ) electroexcitation amplitudes (or 𝛾 𝑣 𝑝𝑁 ∗ electrocouplings) obtained from experiments during the 6 GeV era in Hall B at Jefferson Lab on understanding EHM. Analyzed using continuum Schwinger function methods (CSMs), these results have revealed new pathways for the elucidation of EHM. A good description of the Δ(1232)3/2 + , 𝑁(1440)1/2 + , and Δ(1600)3/2 + electrocouplings, achieved by CSM analyses that express a realistic dressed quark mass function, sheds light on the strong interaction dynamics underlying EHM. Extensions to 𝑁 ∗ studies for higher-mass states are outlined, as well as experimental results anticipated in the 12 GeV era at Jefferson Lab and those that would be enabled by a further increase in the beam energy to 22 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A VLA Survey of Late-time Radio Emission from Superluminous Supernovae and the Host Galaxies

We present the results of 3 GHz radio continuum observations of 23 superluminous supernovae (SLSNe) and their host galaxies by using the Karl G. Jansky Very Large Array conducted 5–21 yr after the explosions. The sample consists of 15 Type I and 8 Type II SLSNe at z < 0.3, providing one of the largest samples of SLSNe with late-time radio data. We detected radio emission from one SLSN (PTF10hgi) and five hosts with a significance of >5σ. No time variability is found in late-time radio light curves of the radio-detected sources in a timescale of years except for PTF10hgi, whose variability is reported in a separate study. Comparison of star formation rates (SFRs) derived from the 3 GHz flux densities with those derived from SED modeling based on UV–NIR data shows that four hosts have an excess of radio SFRs, suggesting obscured star formation. Upper limits for undetected hosts and stacked results show that the majority of the SLSN hosts do not have a significant obscured star formation. By using the 3 GHz upper limits, we constrain the parameters for afterglows arising from interaction between initially off-axis jets and circumstellar medium (CSM). Here, we found that the models with higher energies ( E iso ≳ several × 10 53 erg) and CSM densities ( n ≳ 0.01 cm –3 ) are excluded, but lower energies or CSM densities are not excluded with the current data. We also constrained the models of pulsar wind nebulae powered by a newly born magnetar for a subsample of SLSNe with model predictions in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗