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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

A Numerical and Experimental Study of Damage Growth in a Composite Laminate

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

McElroy, Mark↗
NASA NTRS2014

TPSAS-NF1676L-19616-DND

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

Mark McElroy↗
DOE OSTI2024

Precision Mass Measurement of the Proton Dripline Halo Candidate 22 Al

Here, we report the first mass measurement of the proton-halo candidate 22 Al performed with the low energy beam ion trap facility’s 9.4 T Penning trap mass spectrometer at facility for rare isotope beams. This measurement completes the mass information for the lightest remaining proton-dripline nucleus achievable with Penning traps. 22 Al has been the subject of recent interest regarding a possible halo structure from the observation of an exceptionally large isospin asymmetry [J. Lee et al ., Large isospin asymmetry in Si22/O22 Mirror Gamow-Teller transitions reveals the halo structure of 22 Al , Phys. Rev. Lett. 125, 192503 (2020).]. The measured mass excess value of ME=18 092.5⁢(3) keV, corresponding to an exceptionally small proton separation energy of 𝑆𝑝=100.4⁢(8) keV, is compatible with the suggested halo structure. Our result agrees well with predictions from 𝑠⁢𝑑-shell USD Hamiltonians. While USD Hamiltonians predict deformation in the 22 Al ground state with minimal 1⁢𝑠 1/2 occupation in the proton shell, a particle-plus-rotor model in the continuum suggests that a proton halo could form at large quadrupole deformation. These results emphasize the need for a charge radius measurement to conclusively determine the halo nature.

binding energy↗
DOE OSTI2021

Neutron-unbound states in 31 Ne

Background: The Island of Inversion near the N = 20 shell gap is home to nuclei with reordered single-particle energy levels compared to the spherical shell model. Studies of 31 Ne have revealed that its ground state has a halo component, characterized by a valence neutron orbiting a deformed 30 Ne core. This lightly-bound nucleus with a separation energy of only S n = 170 keV is expected to have excited states that are neutron unbound. Purpose: The purpose of this experiment was to investigate the low-lying excited states in 31 Ne that decay by the emission of a single neutron. Methods: An 89 MeV/u 33 Mg beam impinged on a segmented Be reaction target. Neutron unbound states in 31 Ne were populated via a two-proton knockout reaction. The 30 Ne fragment and associated neutron from the decay of 31 Ne were detected by the MoNA-LISA-Sweeper experimental setup at the National Superconducting Cyclotron Laboratory. Invariant mass spectroscopy was used to reconstruct the two-body decay energy ( 30 Ne + n). Results: The two-body decay energy spectrum exhibits two features: a low-lying peak at 0.30 ± 0.17 MeV and a broad enhancement at 1.50 ± 0.33 MeV, each fit with an energy-dependent asymmetric Breit-Wigner line shape representing a resonance in the continuum. Accompanying shell model calculations using the FSU interaction within NuShellX, combined with cross-section cal culations using the eikonal reaction theory, indicate that these peaks in the decay energy spectrum are caused by multiple resonant states in 31 Ne. Conclusions: In this study, excited states in 31 Ne were observed for the first time. Transitions from calculated shell model final states in 31 Ne to bound states in 30 Ne are in good agreement with the measured decay energy spectrum. Cross-section calculations for the two-proton knockout populating 31 Ne states as well as spectroscopic factors pertaining to the decay of 31 Ne into 30 Ne are used to examine the results within the context of the shell model expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Absorption by H2O and H2O-N2 mixtures at 153 GHz

New experimental data on and a theoretical analysis of the absorption coefficient at 153 GHz are presented for pure water vapor and water vapor-nitrogen mixtures. This frequency is 30 GHz lower than the resonant frequency of the nearest strong water line (183 GHz) and complements our previous measurements at 213 GHz. The pressure dependence is observed to be quadratic in the case of pure water vapor, while in the case of mixtures there are both linear and quadratic density components. By fitting our experimental data taken at several temperatures we have obtained the temperature dependence of the absorption. Our experimental data are compared to several theoretical models with and without a continuum contribution, and we find that none of the models is in very good agreement with the data; in the case of pure water vapor, the continuum contribution calculated using the recent theoretical absorption gives the best results. In general, the agreement between the data and the various models is less satisfactory than found previously in the high-frequency wing. The anisotropy in the observed absorption differs from that currently used in atmospheric models.

Bauer, A.↗
NASA NTRS2005

X-Ray Absorbed, Broad-Lined, Red AGN and the Cosmic X-Ray Background

We have obtained XMM spectra for five red, 2MASS AGN, selected from a sample observed by Chandra to be X-ray bright and to cover a range of hardness ratios. Our results confirm the presence of substantial absorbing material in three sources which have optical classifications ranging from Type 1 to Type 2, with an intrinsically flat (hard) power law continuum indicated in the other two. The presence of both X-ray absorption and broad optical emission lines with the usual strength suggests either a small (nuclear) absorber or a favored viewing angle so as to cover the X-ray source but not the broad emission line region (BELR). A soft excess is detected in all three Type 1 sources. We speculate that this soft X-ray emission may arise in an extended region of ionized gas, perhaps linked with the polarized (scattered) light which is a feature of these sources. The spectral complexity revealed by XMM emphasizes the limitations of the low S/N Chandra data. Overall, the new XMM results strengthen our conclusions (Wilkes et al. 2002) that the observed X-ray continua of red AGN are unusually hard at energies greater than 2 keV. Whether due to substantial line-of-sight absorption or to an intrinsically hard or reflection-dominated spectrum, these 'red' AGN have an observed spectral form consistent with contributing significantly to the missing had absorbed population of the Cosmic X-ray Background (CXRB). When absorption and or reflection is taken into account, all these AGN have power law slopes typical of broad-line (Type 1) AGN (Gamma approximately 1.9). This appears to resolve the spectral paradox which for so long has existed between the CXRB and the AGN thought to be the dominant contributors. It also suggests two scenarios whereby Type 1 AGN/QSOs may be responsible for a significant fraction of the CXRB at energies above 2 keV: 1) X-ray absorbed AGN/QSOs with visible broad emission lines; 2) AGN/QSOs with complex spectra whose hardness greater than 2 keV is not detectable in the typically low S/N data of X-ray surveys. Even if absorption is present in only half of the population, the large number of 'red' AGN suggests a development of unification models, where the continuum source is surrounded, over a substantial solid angle, by the wind or atmosphere of an accretion disk/torus. X-ray observations of such AGN not only provide a check on the presence of absorption, but also a unique probe of the absorbing material. Improved information on their space density, in particular as a function of redshift, will soon be provided by Spitzer-Chandra wide area surveys, allowing better estimates of both the importance of red AGN to the full AGN population and their contribution to the CXRB.

Mushotzky, Richard↗
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↗
NASA NTRS1989

On equivalence of discrete-discrete and continuum-discrete design sensitivity analysis

Developments in design sensitivity analysis (DSA) method have been made using two fundamentally different approaches as shown. In the first approach, a discretized structural finite element model is used to carry out DSA. There are three different methods in the discrete DSA approach: finite difference, semi-analytical, and analytical methods. The finite difference method is a popular one due to its simplicity, but a serious shortcoming of the method is the uncertainty in the choice of a perturbation step size of design variables. In the semi-analytical method, the derivatives of stiffness matrix is computed by finite differences, whereas in the analytical method, the derivatives are obtained analytically. For the shape design variable, computation of analytical derivative of stiffness matrix is quite costly. Because of this, the semi-analytical method is a popular choice in discrete shape DSA approach. However, recently, Barthelemy and Haftka presented that the semi-analytical method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. They found that accuracy problems occur even for a simple cantilever beam. In the second approach, a continuum model of the structure is used to carry out DSA.

Choi, Kyung K.↗
NASA NTRS2008

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗
DOE OSTI2025

Geochemical Assessment for Carbon Sequestration in the Conasauga Group, Northwest Georgia, USA

Sedimentary geological formations are known to be great candidates for geological carbon sequestration. Published studies suggest the southeast of the United States contains many formations suitable for carbon storage. The Cassville 1 Stratigraphic Borehole well could act as a potential carbon reservoir for nearby energy resource facilities in Georgia, United States. Although studies have shown that porous formations are adequate for geological carbon sequestration, it is important to understand possible geochemical reactions between CO 2 and the targeted geological formation before injecting any fluids. In this study, a sandstone sample from the Cassville 1 well is being considered for geological carbon sequestration in the Conasauga Group in Northwest Georgia. Here, the collected sandstone sample, consisting of quartz, K-feldspar, micas, kaolinite, and carbonate minerals such as calcite and dolomite, has a 6% porosity. Leveraging the formation composition and porosity, a one-dimensional continuum reactive transport model was built using CrunchFlow to assess possible geochemical reactions between injected CO 2 and the geological formation. Simulation results show that the carbonate minerals, calcite and dolomite, dissolve during the injection period of 10,000 days, increasing formation porosity from 6% to as much as 30%. The rate and extent of carbonate mineral dissolution and resulting porosity increase are highly sensitive to mineral reactive surface area values. No evidence of mineral precipitation was observed, suggesting that dissolution reactions will control porosity evolution during the CO 2 injection period.

42 ENGINEERING↗
DOE OSTI2022

A Square Pulse Thermoreflectance Technique for the Measurement of Thermal Properties

We report on a laser-based square pulse thermoreflectance (SPTR) technique for the measurement of thermal properties for a wide range of materials. SPTR adopts the pump-probe thermoreflectance principle to monitor the evolution of local temperature after square pulse excitation. The technique features a compact setup, high spatial resolution, and fast data collection. By comparing the acquired SPTR signals with a continuum heat transfer model, material thermal properties can be obtained. Taking advantage of various spot sizes and modulation frequencies, SPTR can measure both the thermal diffusivity and thermal conductivity of poorly to moderately conductive materials and the thermal conductivity of conductive materials with satisfactory accuracy, with potential to be applied to more conductive materials. The technique was validated on three materials: fused silica, single crystal CaF2 and single crystal nickel (with conductivities ranging from 1 W·m -1 ·K -1 to 100 W·m -1 ·K -1 ) with typical measurement errors of 5 % to 20 %. The leading sources of error have been identified by Monte Carlo simulations, and the primary limitations of SPTR are discussed. The compact, fiberized platform we describe here will allow instruments based on this methodology to be deployed in complex, multi-analytical environments for the type of high-throughput correlative analyses that are key to materials design and discovery.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Deformation, dislocation evolution and the non-Schmid effect in body-centered-cubic single- and polycrystal tantalum

A physically-informed continuum crystal plasticity model is presented to elucidate deformation mechanisms, dislocation evolution and the non-Schmid effect in body-centered-cubic (bcc) tantalum widely used as a key structural material for mechanical and thermal extremes. We show the unified structural modeling framework informed by mesoscopic dislocation dynamics simulations is capable of capturing salient features of the large inelastic behavior of tantalum at quasi-static (10 -3 s -1 ) to extreme strain rates (5000 s -1 ) and at low (77 K) to high temperatures (873 K) at both single- and polycrystal levels. Here we also present predictive capabilities of the model for microstructural evolution in the material. To this end, we investigate the effects of dislocation interactions on slip activities, instability and the non-Schmid behavior at the single crystal level. Furthermore, ex situ measurements on crystallographic texture evolution and dislocation density growth are carried out for polycrystal tantalum specimens at increasing strains. Numerical simulation results also support that the modeling framework is capable of capturing the main features of the polycrystal behavior over a wide range of strains, strain rates and temperatures. The theoretical, experimental and numerical results at both single- and polycrystal levels provide critical insight into the underlying physical pictures for micro- and macroscopic responses and their relations in this important class of refractory bcc materials undergoing large inelastic deformations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Electrode roughness dependent electrodeposition of sodium at the nanoscale

Na metal is an attractive anode material for rechargeable Na ion batteries, however, the dendritic growth of Na can cause serious safety issues. Along with modifications of solid-electrolyte interphase (SEI), engineering the electrode has been reported to be effective in suppressing Na dendritic growth, likely by reducing localized current density accumulation. However, fundamental understanding of Na growth at the nanoscale is still limited. Here, we report an in-situ study of Na electrodeposition in electrochemical liquid cells with the electrodes in different surface roughness, e.g., flat or sharp curvature. Real time observation using transmission electron microscopy (TEM) reveals the Na electrodeposition with remarkable details. Relatively large Na grains (in the micrometer scale) are achieved on the flat electrode surface. The local SEI thickness variations impact the growth rate, thus the morphology of individual grains. In contrast, small Na grains (in tens of nanometers) grow explosively on the electrode at the point with sharp curvature. The newly formed Na grains preferentially deposit at the base of existing grains close to the electrode. Further studies using continuum-based computational modeling suggest that the growth mode of an alkali metal (e.g. Na) is strongly influenced by the transport properties of SEI. Our direct observation of Na deposition in combination with the theoretical modeling provides insights for comprehensive understanding of electrode roughness and SEI effects on Na electrochemical deposition.

25 ENERGY STORAGE↗
DOE OSTI2020

Computational predictions of porosities, pore size distributions, and conductivities of aerosol deposited particulate films

The structure and resulting physical properties of aerosol deposited particulate films are governed by the size distribution and morphology of the depositing particles as well as the physics governing particle deposition. While particulate film deposition processes are qualitatively understood, the link between particle characteristics and deposition physics to film properties has not been probed systematically. Here we apply a combination of Langevin dynamics deposition simulations, Monte Carlo pore size distribution calculations, and predictions of thermal conductivities to better establish such process parameter-film property relationships. Furthermore, we account for partial coalescence, polydispersity, and aggregate deposition. We establish that the deposition of polydisperse and aggregated particles leads to broadening of the pore size distribution function and an increase in the mode pore size. Via a non-continuum gas conductivity model, we show that aerosol deposited films can achieve porosities and thermal conductivities similar to conventional aerogels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

On the co-nucleation of adjoining twin pairs at grain boundaries in hexagonal close-packed materials

Here this work addresses the formation of adjoining twin pairs (ATPs) at grain boundaries (GBs) in hexagonal close-packed (hcp) metals from the perspective of the co-nucleation (CN) of pairs of deformation twins. A continuum defect mechanics model is proposed to investigate the energetic feasibility of the CN of ATPs resulting from the dissociation of GB dislocations. Application of the model to a set of GBs in Mg reveals that CN is largely preferred over the nucleation of a single twin variant for low misorientation angle GBs, which is consistent with previous experimental analyses. The CN events are subsequently analyzed considering the GB character, and the twin system alignment ($m$') and Schmid factors to obtain better insights into the conditions governing the formation of ATPs via CN. The likelihood of CN even at low $m$' values reveals that CN events could be responsible for the formation of ATP pairs in these scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Pressure Equation of State of 1,3,5-triamino-2,4,6-trinitrobenzene: Insights into the Monoclinic Phase Transition, Hydrogen Bonding, and Anharmonicity

The high-pressure equation of state (EOS) of energetic materials (EMs) is important for continuum and mesoscale models of detonation performance and initiation safety. Additionally, obtaining a high-fidelity EOS of the insensitive EM 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has proven to be difficult because of challenges in experimental characterization at high pressures (HPs). In this work, powder X-ray diffraction patterns were fitted using the recently discovered monoclinic I2/a phase above 4 GPa, which shows that TATB is less compressible than when indexed with the triclinic $P\bar{1}$ phase. First-principles calculations were performed with Perdew–Burke–Ernzerhof (PBE) and PBE0 functionals including thermal effects using the $P\bar{1}$ phase. PBE0 improves the description of hydrogen bonding and thus predicts accurate planar a and b lattice parameters under ambient conditions. However, discrepancies in the predicted lattice parameters above 4–10 GPa compared with experimental measurements indexed with $P\bar{1}$ are further evidence of a structural modification at high pressure. Layer sliding defects are formed during molecular dynamics simulations, which induces an anharmonic effect on the thermal expansion of the c lattice parameter. In short, the results provide several insights into determining high-fidelity EOS parameters for TATB and other molecular crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗