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 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Modeling Chemo-Mechanics with Electrolyte Infiltration to Quantify Degradation of Cathode Particles

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. To that end, this talk investigates the formation of stress-induced fracture within polycrystalline cathode particles and the impact on capacity loss. Physically based cathode aging dynamics is simulated in a 3D, continuum-level chemo-mechanical model. The model captures anisotropic Li diffusion within a single polycrystalline particle comprised of hundreds to thousands of randomly oriented grains. A recent addition to this model includes electrolyte infiltration, which occurs when the electrolyte seeps into surface cracks within the particle. The model predicts that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains are misaligned, they expand a different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicts that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicates that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. The principle contribution of this research is the use of an anisotropic chemo-mechanical model to test how particle geometry affect capacity fade, which predicts that particle size has a stronger effect on capacity fade than grain size and ultimately that small particles with large grains have the least capacity fade.

cathode↗

DOE OSTI2020

The charge reduction rate for multiply charged polymer ions via ion–ion recombination at atmospheric pressure

The charge reduction of multiply charged macromolecular ions via recombination with small ions in the gas phase is commonly employed to modulate the charge on macromolecules prior to mass spectrometric and mobility analyses. We employ a recently developed continuum-Molecular Dynamics (MD) calculation approach to determine the recombination rate coefficient of multiply charged (1 to 7 excess positive charged) polyethylene glycol ions (mass of 4600 Da) with smaller singly charged anions, modeled as NO 2 – ions. The continuum-MD approach accounts explicitly for the influence of the background gas on the recombination process, accounts explicitly for ion translational, vibrational, and rotational motion, and enables recombination rate coefficient calculation in nitrogen near atmospheric pressure, wherein neither low pressure nor high pressure recombination theories are strictly applicable. Continuum-MD simulations yield recombination rate coefficients near 3.9 × 10 –14 m 3 s –1 for singly charged ions, increasing to 3.0 × 10 –11 m 3 s –1 for the +7 ion. Pre-existing collision rate coefficient expressions for rigid ions are found to be within a factor of 2–5 of calculations for all charge states, but their use requires knowledge of an appropriate collision distance, which is not well-defined for flexible polymer ions. Continuum-MD-inferred rate coefficients are incorporated into a model of charge reduction, and the charge state distribution versus anion concentration determined with it is compared to charge reduction measurements made via atmospheric pressure differential mobility analysis. Furthermore, good agreement is observed between simulations and experiments; although results are extremely sensitive to the recombination rate coefficients, experimental results are bound by models utilizing rates within a factor of 2 (0.5–2.0×) of the continuum-MD rates.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS1973

Direct integration transmittance model

A transmittance model was developed for the 200-2000/cm region for interpretation of high spectral resolution measurements of laboratory absorption and of planetary thermal emission. The high spectral resolution requires transmittances to be computed monochromatically by summing the contribution of individual molecular absorption lines. A magnetic tape atlas of H2O,O3, and CO2 molecular line parameters serves as input to the transmittance model with simple empirical representations used for continuum regions wherever suitable laboratory data exist. The theoretical formulation of the transmittance model and the computational procedures used for the evaluation of the transmittances are discussed. Application is demonstrated of the model to several homogenous path laboratory absorption examples.

Kunde, V. G.↗

NASA NTRS1974

Direct integration transmittance model

A transmittance model has been developed for interpretation of high spectral resolution measurements of laboratory absorption and of planetary thermal emission. The high spectral resolution requires transmittances to be computed monochromatically by summing the contribution of individual molecular absorption lines. A magnetic tape atlas of H2O, O3, and CO2 molecular line parameters serves as input to the transmittance model with simple empirical representations used for continuum regions wherever suitable laboratory data exist. The theoretical formulation of the transmittance model and the computational procedures used for the evaluation of the transmittances are discussed, and application of the model to several homogeneous-path laboratory absorption examples is demonstrated.

Kunde, V. G.↗

NASA NTRS1992

Reaction zone structure for strong, weak overdriven, and weak underdriven oblique detonations

A simple dynamic systems analysis is used to give examples of strong, weak overdriven, and weak underdriven oblique detonations. Steady oblique detonations consisting of a straight lead shock attached to a solid wedge followed by a resolved reaction zone structure are admitted as solutions to the reactive Euler equations. This is demonstrated for a fluid that is taken to be an inviscid, calorically perfect ideal gas that undergoes a two-step irreversible reaction with the first step exothermic and the second step endothermic. This model admits solutions for a continuum of shock wave angles for two classes of solutions identified by a Rankine-Hugoniot analysis: strong and weak overdriven waves. The other class, weak underdriven, is admitted for eigenvalue shock-wave angles. Chapman-Jouguet waves, however, are not admitted. These results contrast those for a corresponding onestep model that, for detonations with a straight lead shock, only admits strong, weak overdriven, and Chapman-Jouguet solutions.

Powers, Joseph M.↗

NASA NTRS2001

CFD Simulation of Liquid Rocket Engine Injectors

Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.

Farmer, Richard↗

NASA NTRS2008

Modeling the Non-Linear Response of Fiber-Reinforced Laminates Using a Combined Damage/Plasticity Model

The present work is concerned with modeling the non-linear response of fiber reinforced polymer laminates. Recent experimental data suggests that the non-linearity is not only caused by matrix cracking but also by matrix plasticity due to shear stresses. To capture the effects of those two mechanisms, a model combining a plasticity formulation with continuum damage has been developed to simulate the non-linear response of laminates under plane stress states. The model is used to compare the predicted behavior of various laminate lay-ups to experimental data from the literature by looking at the degradation of axial modulus and Poisson s ratio of the laminates. The influence of residual curing stresses and in-situ effect on the predicted response is also investigated. It is shown that predictions of the combined damage/plasticity model, in general, correlate well with the experimental data. The test data shows that there are two different mechanisms that can have opposite effects on the degradation of the laminate Poisson s ratio which is captured correctly by the damage/plasticity model. Residual curing stresses are found to have a minor influence on the predicted response for the cases considered here. Some open questions remain regarding the prediction of damage onset.

Schuecker, Clara↗

NASA NTRS1989

Flare model sensitivity of the Balmer spectrum

Careful studies of various chromospheric spectral signatures are very important in order to explore their possible sensitivity to the modifications of the thermodynamic quantities produced by the flare occurrence. Pioneer work of Canfield and co-workers have shown how the H alpha behavior is able to indicate different changes in the atmospheric parameters structure associated to the flare event. It was decided to study the behavior of the highest Balmer lines and of the Balmer continuum in different solar flare model atmospheres. These spectral features, originating in the deep photosphere in a quiet area, may have a sensitivity different from H alpha to the modification of a flare atmosphere. The details of the method used to compute the Stark profile of the higher Balmer line (n is greater than or equal to 6) and their merging were extensively given elsewhere (Donati-Falchi et al., 1985; Falchi et al., 1989). The models used were developed by Ricchiazzi in his thesis (1982) evaluating the chromospheric response to both the nonthermal electron flux, for energy greater than 20 kev, (F sub 20) and to the thermal conduction, (F sub c). The effect of the coronal pressure values (P sub O) at the apex of the flare loop is also included.

Falchi, A.↗

NASA NTRS1976

Helium emission from model flare layers

The emission of visible and UV He I and He II line radiation from a plane-parallel model flare layer characterized by electron temperatures of 10,000 to 50,000 K and electron densities of 10 to the 10th power to 10 to the 15th power per cu cm is analyzed by solving the statistical-equilibrium equations for a 30-level He I-II-III system, using parametric representations of the line and continuum radiation fields. The atomic model was chosen to provide accurate solutions for the first two resonance lines of He I and He II as well as for the D3 and 10,830-A lines of He I. Reaction rates are discussed, and sample solutions to the steady-state population equations are given for a generally optically thin gas assumed to be irradiated over 2pi sr by a blackbody spectrum at 6000 K. Specific results are examined for ionization equilibrium, level populations, approximate optical depths of a 1000-km-thick flare layer, line intensities, and upper-level population rates.

Kulander, J. I.↗

NASA NTRS1989

The Lyman edge test of the quasar emission mechanism

The Lyman edge region in the spectra of 11 high-redshift quasars is studied using data designed to detect broadened Lyman edge absorption predicted by thermal models of the 'big blue bump' continuum component, especially thin accretion disks. Three objects have broadened, partial edges near the systemic redshifts, as expected for thermal models. In two cases, however, narrow Lyamn lines with the edge absorption, ruling out the disk interpretation. Thus, the data are negative for thermal models, at least for opaque thin accretion disks.

Antonucci, R. R. J.↗

DOE OSTI2023

A superconfiguration calculation of opacity with consistent bound and continuum electron treatments using green’s functions

One of the challenges in calculating the opacity of dense plasmas is the difficulty in consistently modeling electrons bound to nuclei and those that exist within the continuum of free states in electronic structure models. We address this issue by adapting the green's function approach, originally developed for use in average atom calculations, to the determination of superconfiguration (SC) electronic structure. The spectra created using these SCs indicate that a consistent treatment of continuum electronic structure is important for phenomena involving electrons near ionization thresholds, such as the pressure ionization of bound states and the opacity due to transitions near bound-free edges. Though important for dense plasmas, the detailed incorporation of continuum electrons into structure calculations does not have significant impact on the recent discrepancies between the predicted and measured opacity of hot, dense iron (Bailey et al 2015 Nature517 56). Here we find that the inclusion of plasma effects through an ion-sphere model along with our treatment of continuum electronic states gives a description of pressure ionization in hot, dense aluminum that is in better agreement with experiment than methods that rely on perturbative descriptions of the plasma environment (Hoarty et al 2013 Phys. Rev. Lett.110 265003).

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Tardigrade-examples V0.1.0

Tardigrade-examples is a repository of computational workflows that exercise the Tardigrade software package. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling. These capabilities include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This repository investigates continuum upscaling of various direct numerical simulations (DNS) conducted in Abaqus finite element (FE), Ratel FE, and GEOS material point method (MPM) software. Verification of the upscaling workflow is first investigated by considering DNS of trivial stress states for homogeneous materials, results of which indicate that classical continuum behavior is recovered as expected. DNS of heterogeneous materials are then considered.

Allard, Thomas↗

DOE OSTI2024

Lattice regularization of reduced Kähler-Dirac fermions and connections to chiral fermions

We show how a path integral for reduced Kähler-Dirac fermions suffers from a phase ambiguity associated with the fermion measure that is an analog of the measure problem seen for chiral fermions. However, unlike the case of chiral symmetry, a doubler free lattice action exists which is invariant under the corresponding onsite symmetry. This allows for a clear diagnosis and solution to the problem using mirror fermions resulting in a unique gauge invariant measure. By introducing an appropriate set of Yukawa interactions which are consistent with ’t Hooft anomaly cancellation we conjecture the mirrors can be decoupled from low energy physics. Moreover, the minimal such Kähler-Dirac mirror model yields a light sector which corresponds, in the flat space continuum limit, to the Pati-Salam GUT model.

Physics↗

DOE OSTI2021

Climate-Relevant Gas Absorption Properties from AWARE and Other ARM Spectral Measurements (Final Report)

The objective of this project has been to use measurements from the Atmospheric Emitted Radiance Interferometer (AERI) deployed at ARM sites to improve our knowledge of uncertain infrared spectroscopic parameters of importance to climate, remote sensing, and data assimilation. A primary focus has been the water vapor continuum in the infrared atmospheric window. Radiation codes used to predict climate and weather base their representation of the water vapor continuum in this window on the MT_CKD model, but in this spectral region the MT_CKD water vapor continuum absorption coefficients were derived almost two decades ago by an analysis of AERI measurements for conditions with a limited range of precipitable water vapor (PWV) and temperature values, leading to uncertainty in the derived coefficients. This project has included a new, comprehensive analysis of all aspects of the water vapor continuum in the atmospheric window -- the self continuum, the foreign continuum, and the self continuum temperature dependence, all resolved spectrally. Having a wide range of PWV amounts and temperatures from deployments of the ARM Mobile Facility for GoAmazon as well as the more than a decade of AERI measurements at SGP has provided an appropriate foundation for such a comprehensive analysis.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2018

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗

DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗

NASA NTRS1995

EUVE spectroscopy of epsilon Canis Majoris (B2 II) from 70 to 730 A

We present spectra of the brightest stellar source of extreme ultraviolet (EUV) radiation longward of 400 A, the B2 II star, epsilon CMa. These data were taken with the three spectrometers aboard the NASA Extreme Ultraviolet Explorer satellite (EUVE) during the first cycle of pointed observations. We report on our initial studies of the continuum and line spectrum of the stellar photosphere in the 320 to 730 A region, and on the wind emission lines observed in the 170-375 A region. This is the first EUV spectrum of an early-type star, and thus makes epsilon CMa the most comprehensively observed B star from the X-ray to infrared regimes. The radiation in both the H Lyman continuum and He I continuum (shortward of 504 A) are found to be significantly greater than predicted by both Local Thermodynamic Equilibrium (LTE) and non-LTE model atmospheres. Since epsilon CMa also exhibits a mid-infrared excess, this points to the outer layers being warmer than the models indicate. The anomalously large Lyman continuum flux, combined with the very low column density measured in the direction toward this star implies that it is the dominant source of hydrogen ionization of the local interstellar medium in the immediate vicinity of the sun. All of the lines predicted to be strong from model atmospheres are present and several wind absorption features are also identified. We have detected emission lines from highly ionized iron that are consistent with the ROSAT Position Sensitive Proportional Counter (PSPC) observations if a multi-temperature emission model is used, and the assumption is made that there is significant absorption beyond that of the neutral phase of the ISM. The spectrum shows strong O III 374 A line emission produced by the Bowen flourescence mechanism, which has not previously been observed in the spectra of hot stars.

Cassinelli, J. P.↗

NASA NTRS2011

The Extreme Spin of the Black Hole in Cygnus X-1

The compact primary in the X-ray binary Cygnus X-1 was the first black hole to be established via dynamical observations. We have recently determined accurate values for its mass and distance, and for the orbital inclination angle of the binary. Building on these results, which are based on our favored (asynchronous) dynamical model, we have measured the radius of the inner edge of the black hole s accretion disk by fitting its thermal continuum spectrum to a fully relativistic model of a thin accretion disk. Assuming that the spin axis of the black hole is aligned with the orbital angular momentum vector, we have determined that Cygnus X-1 contains a near-extreme Kerr black hole with a spin parameter a* > 0.95 (3(sigma)). For a less probable (synchronous) dynamical model, we find a. > 0.92 (3 ). In our analysis, we include the uncertainties in black hole mass, orbital inclination angle, and distance, and we also include the uncertainty in the calibration of the absolute flux via the Crab. These four sources of uncertainty totally dominate the error budget. The uncertainties introduced by the thin-disk model we employ are particularly small in this case given the extreme spin of the black hole and the disk s low luminosity.

Gou, Lijun↗