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.

SourceResearch connectionFrom yearThrough year

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

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS2002

A Chandra Spectroscopic Survey of Persistent Black Hole Candidates

We present results from observations of persistent black hole candidates with the High-Energy Transmission Gratings on board the Chandra X-Ray Observatory. The sources include LMC X-1, LMC X-3, GRS 1758-258, and Cyg X-1. Along with the published results on 1E 1740.7-2942, we have completed a high-resolution spectroscopic survey of such systems. The observed X-ray spectra of LMC X-1 and LMC X-3 show no prominent discrete features, while absorption edges (Mg K and Si K) are detected in the spectrum of GRS 1758-258. The edges are likely to be of interstellar origin. In most cases, the X-ray continuum can be described well by models that are often adopted in low-resolution studies of black hole candidates: a multitemperature disk spectrum plus a Comptonization component. However, the relative contribution of the two components varies greatly among different sources. For instance, only the disk component is present for LMC X-1 and GRS 1758-258, while the Comptonized component is required for other sources. We discuss general issues related to obtaining disk parameters from modeling X-ray continuum.

Cui, Wei↗

NASA NTRS1987

LTE modeling of inhomogeneous chromospheric structure using high-resolution limb observations

The paper discusses considerations relevant to LTE modeling of rough atmospheres. Particular attention is given to the application of recent high-resolution observations of the solar limb in the far-infrared and radio continuum to the modeling of chromospheric spicules. It is explained how the continuum limb observations can be combined with morphological knowledge of spicule structure to model the physical conditions in chromospheric spicules. This discussion forms the basis for a chromospheric model presented in a parallel publication based on observations ranging from 100 microns to 2.6 mm.

Lindsey, C.↗

NASA NTRS1988

Hydrogen line and continuum emission in young stellar objects. I - Excitation model

Two mechanisms that populate the n = 2 level of hydrogen after the Lyman continuum is depleted are identified. They are ionization of N I from its excited states, followed by charge-exchange between N II and H I, and Ly-beta line wing absorption. Both processes involve absorption of the sub-Lyman continuum between 11 and 13.6 eV. With population in level n = 2 thus maintained, the strong Brackett line fluxes observed are then produced as a result of Balmer photoionization. The magnitude of the sub-Lyman continuum dictates the fraction of Balmer continuum that will be absorbed. Numerical calculations for four young stellar objects with luminosity ranging from 10 to 10,000 solar are performed, and it is concluded that this two-step process of sub-Lyman continuum absorption followed by Balmer photoionization can account for the great majority of observed Brackett line fluxes. The location and mass of the emitting as as determined from the Brackett line fluxes are reported. The Br-alpha luminosity is calculated as a function of the mass loss rate.

Kwan, John↗

NASA NTRS2010

RXTE Observations of the 1A 1118-61 in an Outburst, and the Discovery of a Cyclotron Line

We present the analysis of RXTE monitoring data obtained during the January 2009 outburst of the hard X-ray transient IA 1118-61. Using these observations the broadband (3.5-120 keV) spectrum of the source was measured for the first time ever. We have found that the broadband continuum spectrum of the source is similar to other accreting pulsars and is well described by several conventionally used phenomenological models. We have discovered that regardless of the applied continuum model, a prominent broad absorption feature at approx. 55 keV is observed. We interpret this feature as a Cyclotron Resonance Scattering Feature (CRSF). The observed CRSF energy is one of the highest known and corresponds to a magnetic field of B approx. 4.8 x 10(exp 12) G in the scattering region. Furthermore, our data suggests an iron emission line presence, which was not reported previously for lA 1118-61 as well. Timing properties of the source, including a strong spin-up, were found to be similar to those observed by CGRO/BATSE during the previous outburst, however the broadband capabilities of RXTE reveal a more complicated energy dependency of the pulse-profile.

Doroshenko, V.↗

DOE OSTI2026

Tori, Klein bottles, and modulo 8 parity/time-reversal anomalies of 2+1d staggered fermions

We study the symmetries of lattice staggered fermions in 2+1d. Using the symmetries, we can place the system on any sheared torus or Klein bottle. These different backgrounds provide diagnostics of various ’t Hooft anomalies associated with the crystalline symmetries. We then compare the lattice model to its continuum limit. The symmetries of the lattice system are mapped in a nontrivial way to the symmetries of the continuum theories. Using this map, we match the ’t Hooft anomalies on the lattice and the continuum. Along the way, we develop a general formalism to study Hamiltonian lattice models on nontrivial, compact, flat spaces.

Algorithms and Theoretical Developments↗

DOE OSTI2021

Modeling Lithium Diffusion in Battery Cathodes Considering Chemo-Mechanically Induced Damage

This talk will present a 3D, continuum-level damage model for simulating Lithium diffusion within generated Li_{x}Ni_{0.5}Mn_{0.3}Co_{0.2} (NMC 532) secondary cathode particles. The primary motivation of the particle-level model is to inform cathode-particle design and determine charging profiles that reduce cathode fracture. The model considers NMC 532 secondary particles containing an agglomeration of anisotropic, randomly oriented grains. The model predicts that secondary-particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. The model predicts that small secondary-particles with large grains develop significantly less damage than larger secondary particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first high-rate cycles. This chemo-mechanical “damage saturation” effect indicates that initial secondary-particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

DOE OSTI2020

More exotic field theories in 3+1 dimensions

We continue the exploration of nonstandard continuum field theories related to fractons in 3+1 dimensions. Our theories exhibit exotic global and gauge symmetries, defects with restricted mobility, and interesting dualities. Depending on the model, the defects are the probe limits of either fractonic particles, strings, or strips. One of our models is the continuum limit of the plaquette Ising lattice model, which features an important role in the construction of the X-cube model.

Gorantla, Pranay↗

DOE OSTI2021

fcc lattice, checkerboards, fractons, and quantum field theory

In this work, we consider XY-spin degrees of freedom on an FCC lattice, such that the system respects some subsystem global symmetry. We then gauge this global symmetry and study the corresponding $U(1)$ gauge theory on the FCC lattice. Surprisingly, this $U(1)$ gauge theory is dual to the original spin system. We also analyze a similar $\mathbb{Z}_N$ gauge theory on that lattice. All these systems are fractonic. The $U(1)$ theories are gapless and the $\mathbb{Z}_N$ theories are gapped. We analyze the continuum limits of all these systems and present free continuum Lagrangians for their low-energy physics. Our $\mathbb{Z}_2$ FCC gauge theory is the continuum limit of the well known checkerboard model of fractons. Our continuum analysis leads to a straightforward proof of the known fact that this theory is dual to two copies of the $\mathbb{Z}_2$ X-cube model. We find new models and new relations between known models. The $\mathbb{Z}_N$ FCC gauge theory can be realized by coupling three copies of an anisotropic model of lineons and planons to a certain exotic $\mathbb{Z}_2$ gauge theory. Also, although for $N=2$ this model is dual to two copies of the $\mathbb{Z}_2$ X-cube model, a similar statement is not true for higher $N$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1982

Molecular dynamics simulation of propagating cracks

Steady state crack propagation is investigated numerically using a model consisting of 236 free atoms in two (010) planes of bcc alpha iron. The continuum region is modeled using the finite element method with 175 nodes and 288 elements. The model shows clear (010) plane fracture to the edge of the discrete region at moderate loads. Analysis of the results obtained indicates that models of this type can provide realistic simulation of steady state crack propagation.

Mullins, M.↗

DOE OSTI2026

Active learning of a crystal plasticity flow rule from discrete dislocation dynamics simulations

Continuum-scale material deformation models, such as crystal plasticity (CP), can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e. mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a CP flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization (BO). We demonstrate a semi-autonomous BO loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model CP flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Mesh objective stochastic simulations of quasibrittle fracture

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

42 ENGINEERING↗

DOE OSTI2025

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗

NASA NTRS2018

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

DOE OSTI2024

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Understanding Variability in Mechanical Behavior of Adhesively Bonded Sandwich Joints due to Stochastic Material Damage

A mesoscale progressive damage analysis was applied to the analysis of recently tested adhesively bonded sandwich panels, in which ten nominally identical tension loaded panels failed in three different failure modes, with peakloads ranging from 13.0 –16.1 kips. A parametric finite element model for progressive damage analysis of an adhesively bonded joint (ABJ) considering matrix and fiber damage, and delamination cracks in unidirectional and fabric composite sections of the joint, as well as core crushing, and adhesive fracture was developed. The model was developed using Abaqus/Explicit and utilized continuum damage mechanics material models for intralaminar damage in the facesheet and doubler plies, cohesive elements for delamination and adhesive debonding, and discretized cells considering plasticity to model crushing of the aluminum honeycomb core. Analysis predictions obtained with a pristine model agreed well with the average experimental peak load, strain in the joint, and post-test damage states. Then, various methods to incorporate stochastic material variability into ABJ model were evaluated. These methods included simulation of pre-existing matrix damage, pre-existing adhesive porosity, and uniform modification of adhesive properties. These material and manufacturing defects changed the damage mode predictions in some cases. The analysis of stochastic strength due to the considered defects revealed unintuitive damage mode interaction in the ABJ, explaining some structure-property relationships observed experimentally.

Richard Larson↗

NASA NTRS2006

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗

DOE OSTI2020

Calibration of parallel bond parameters in bonded particle models via physics-informed adaptive moment optimisation

This study proposes an automated calibration procedure for bond parameters in bonded discrete element modelling. By exploring the underlying physical correlations between microscopic parameters of bonds and macroscopic strength parameters of the continuum to be modelled, the microscopic shear strength and tensile strength are identified as independent variables for calibration purpose. Then a physics-informed iterative scheme is proposed to automatically approximate the bond parameters by viewing the micro-macro relation as an implicitly defined mathematical mapping function. As a result of highly non-convex features of this implicit mapping, the adaptive moment estimation (Adam), which is especially suitable for problems with noisy gradients, is adopted as the basic iterative scheme, in conjunction with other numerical techniques to approximately evaluate the partial derivatives involved. The whole procedure offers a simple and effective framework for bond parameter calibration. Finally, a numerical example of SiC ceramic is provided for validation. By compared with some existing calibration methods, the proposed method shows significant advantages in terms of calibration efficiency and accuracy.

42 ENGINEERING↗