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 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach

With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multi-modal Approach to Modeling Creep Deformation in Nickel-base Superalloy

This study developed the multi-modal modeling approaches to elucidate the physical and overall mechanical properties of two complex alloy systems of Ni-based Superalloys, namely Haynes 282 and Inconel 740. The overall strategies employed systematically in different scales of modeling are based on the common structural framework of the multi phases that made up the Ni-based Superalloys, namely the matrix phase of the Gamma ($\gamma$) phase which is embedded with ordered Gamma prime ($\gamma$’) precipitates and decorated with the dispersions of mostly carbide precipitates. With this in mind, we had set up the modeling platforms to incorporate a wide range of scales, starting from electronic structures to the continuum level via the crystal plasticity model.

20 FOSSIL-FUELED POWER PLANTS↗

NASA NTRS1975

Energy transfer by radiation in non-grey atomic gases in isothermal and non-isothermal slabs

A multiband model for the absorption coefficient of atomic hydrogen-helium plasmas is constructed which includes continuum and line contributions. Emission from 28 stronger lines of 106 that have been screened is considered, of which 21 are from hydrogen and 7 belong to helium, with reabsorption due to line-line, line-continuum overlap accurately accounted for. The model is utilized in the computation of intensities and fluxes from shock-heated slabs of 85% H2-15% He mixtures for slab thicknesses from 1 to 30 cm, temperature from 10,000 to 20,000 K, and for different densities. In conjunction with the multiband model, simple numerical schemes have been devised which provide a quick and comprehensive way of computing radiative energy transfer in nonisothermal and nongrey gases.

Poon, P. T. Y.↗

DOE OSTI2020

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

DOE OSTI2025

exoALMA. XI. ALMA Observations and Hydrodynamic Models of LkCa 15: Implications for Planetary Mass Companions in the Dust Continuum Cavity

In the past decade, the Atacama Large Millimeter/submillimeter Array (ALMA) has revealed a plethora of substructures in the disks surrounding young stars. These substructures have several proposed formation mechanisms, with one leading theory being the interaction between the disk and newly formed planets. In this Letter, we present high angular resolution ALMA observations of LkCa 15’s disk that reveal a striking difference in dust and CO emission morphology. The dust continuum emission shows a ringlike structure characterized by a dust-depleted inner region of ~40 au in radius. Conversely, the CO emission is radially smoother and shows no sign of gas depletion within the dust cavity. We compare the observations with models for the disk–planet interaction, including radiative transfer calculation in the dust and CO emission. This source is particularly interesting, as the presence of massive planets within the dust cavity has been suggested based on previous near-IR observations. We find that the level of CO emission observed within the dust cavity is inconsistent with the presence of planets more massive than Jupiter orbiting between 10 and 40 au. Instead, we argue that the LkCa 15 innermost dust cavity might be created either by a chain of low-mass planets or by other processes that do not require the presence of planets.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

An Automated Meshing Framework for Progressive Damage Analysis of Fabrics Using CompDam

A meshing and model construction framework for progressive damage analysis of fabric composites is proposed. A mesoscale idealization is used such that the fiber tows are considered homogeneous and transversely isotropic continua, and the matrix is a homogeneous isotropic continuum. By using this idealization, relatively mature damage modeling techniques developed for laminates of unidirectional plies are transferrable to fabrics. Cohesive interface elements and a continuum damage mechanics (CDM) model implemented in the software CompDam are used to account for the key damage modes including cracking at the tow/matrix interface, transverse matrix cracking in the tows, and fiber breakage in the tows. Matrix degradation is accounted for by plasticity. The modeling approach is demonstrated for a 3D orthogonal woven fabric. The main novelty of the proposed model construction framework is in its flexibility resulting from a formulation that is agnostic to the fabric architecture and implementation using mesh based operations. The algorithm does not make any assumptions about the tow sizes, shape, path, or proximity to other tows, thus it is applicable to a broad range of fabric architectures. The matrix region mesh surrounding the tows is constructed from domain boundary and tow surface meshes (i.e. through a series of mesh-based operations) such that the model domain can be arbitrarily shaped.

Andrew Cole Bergan↗

DOE OSTI2024

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Precision test of gauge/gravity duality in D0-brane matrix model at low temperature

We test the gauge/gravity duality between the matrix model and type IIA string theory at low temperatures with unprecedented accuracy. To this end, we perform lattice Monte Carlo simulations of the Berenstein-Maldacena-Nastase (BMN) matrix model, which is the one-parameter deformation of the Banks-Fischler-Shenker-Susskind (BFSS) matrix model, taking both the large N and continuum limits. We leverage the fact that sufficiently small flux parameters in the BMN matrix model have a negligible impact on the energy of the system while stabilizing the flat directions so that simulations at smaller N than in the BFSS matrix model are possible. Hence, we can perform a precision measurement of the large N continuum energy at the lowest temperatures to date. The energy is in perfect agreement with supergravity predictions including estimations of α'-corrections from previous simulations. At the lowest temperature where we can simulate efficiently (T = 0.25λ 1/3 , where λ is the ’t Hooft coupling), the difference in energy to the pure supergravity prediction is less than 10%. Furthermore, we can extract the coefficient of the 1/N 4 corrections at a fixed temperature with good accuracy, which was previously unknown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1993

A numerical analysis of high-temperature heat pipe startup from the frozen state

Continuum and rarefied vapor flows co-exist along the heat pipe length for most of the startup period. A two-region model is proposed in which the vapor flow in the continuum region is modeled by the compressible Navier-Stokes equations, and the vapor flow in the rarefied region is simulated by a self-diffusion model. The two vapor regions are linked with appropriate boundary conditions, and heat pipe wail, wick, and vapor flow are solved as a conjugate problem. The numerical solutions for the entire heat pipe startup process from the frozen state are compared with the corresponding experimental data with good agreement.

Cao, Y.↗

DOE OSTI2023

Modeling support for the development of material surveillance specimens and procedures

This report describes modeling and simulation activities performed at Argonne National Laboratory supporting the development of passively actuated mechanical test articles for material surveillance in Molten Salt Reactors (MSRs). These test articles are a critical technology in formulating material surveillance programs for future MSRs to monitor the degradation in the structural properties of the materials in critical plant components. The main activity described in this report is the development of a method for inferring the amount of mechanical damage a test article has experienced during some duration of exposure to plant thermal and environmental conditions, using only mechanical test data collected from the test articles before and after exposure. The basic approach is to develop a model of the test article, including a description of mechanical degradation through continuum damage mechanics, and then use this model to cast the problem of inferring mechanical degradation in the test article materials into a shooting problem for a set of ordinary differential equations. We can then solve the shooting problem to determine the amount of damage accumulated in the sample. The report also details a few miscellaneous simulation studies completed at Argonne to support the development of the test articles themselves at Idaho National Laboratory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Active Galactic Nuclei Continuum Reverberation Mapping Based on Zwicky Transient Facility Light Curves

We perform a systematic survey of active galactic nuclei (AGNs) continuum lags using ~3 days cadence gri-band light curves from the Zwicky Transient Facility. We select a sample of 94 type 1 AGNs at z < 0.8 with significant and consistent inter-band lags based on the interpolated cross-correlation function method and the Bayesian method JAVELIN. Within the framework of the "lamp-post" reprocessing model, our findings are: (1) The continuum emission (CE) sizes inferred from the data are larger than the disk sizes predicted by the standard thin-disk model. (2) For a subset of the sample, the CE size exceeds the theoretical limit of the self-gravity radius (12 lt-days) for geometrically thin disks. (3) The CE size scales with continuum luminosity as R CE ∝ L 0.48±0.04 with a scatter of 0.2 dex, analogous to the well-known radius–luminosity relation of broad Hβ. These findings suggest a significant contribution of diffuse continuum emission from the broad-line region (BLR) to AGN continuum lags. We find that the R CE –L relation can be explained by a photoionization model that assumes ~23% of the total flux comes from the diffuse BLR emission. In addition, the ratio of the CE size and model-predicted disk size anticorrelates with the continuum luminosity, which is indicative of a potential nondisk BLR lag contribution evolving with the luminosity. Finally, a robust positive correlation between the CE size and black hole mass is detected.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1977

Spectrophotometry of the transient X-ray source A0620-00

Optical observations of the transient X-ray source A0620-00 over the wavelength range from about 3200 to 6600 A are analyzed which were made with the 5-m Hale telescope using a multichannel spectrometer, a digital spectrograph with a cooled integrating SIT vidicon, and an image photon-counting system attached to the camera of the coude spectrograph. Infrared observations at 1.65, 2.2, and 3.5 microns are also examined. Individual spectral lines of special interest are discussed, the E(B-V) color excess is estimated to be about 0.40, and a distance of 1200 pc is adopted. Analysis of the continuous spectrum yields several simple binary models. A specific model is proposed in which one of the continuum sources is an optically thick object characterized by a model atmosphere at a temperature of 25,000 to 30,000 K and the second source is either a cool black body at 4000 K or a partially optically thick source of bound-free plus free-free radiation from a hydrogen gas with an electron temperature of at least 100,000 K. Some possible parameters of these sources and of the accretion disk are considered.

Oke, J. B.↗

NASA NTRS1985

Photodissociation regions. I - Basic model. II - A model for the Orion photodissociation region

A theoretical parameter study of the temperature and chemical structure of dense photodissociation regions and their resultant spectrum is presented. Models are discussed which are relevant not only to the dust and gas between molecular clouds and H II regions, but also apply to any neutral cloud illuminated by intense FUV fluxes. The models relate observed line and continuum emission from these regions to physical parameters such as the gas density and temperature, the elemental and chemical abundances, the local radiation field, and the grain properties. The results are applied to observational data from the OMC-1 region. The model shows that the observed high brightness temperature of the C I 609 microns line can be explained by emission from the C(+)/C/CO transition region. This difference with previous chemical models is due to a higher gas phase elemental abundance of carbon, to the charge exchange reactions of C(+) with S and SiO, and to carbon self-shielding.

Tielens, A. G. G. M.↗

DOE OSTI2021

Low-energy limit of some exotic lattice theories and UV/IR mixing

We continue our exploration of exotic, gapless lattice and continuum field theories with subsystem global symmetries. In an earlier paper, we presented free lattice models enjoying all the global symmetries (except continuous translations), dualities, and anomalies of the continuum theories. Here, we study in detail the relation between the lattice models and the corresponding continuum theories. We do that by analyzing the spectrum of the theories and several correlation functions. These lead us to uncover interesting subtleties in the way the continuum limit can be taken. In particular, in some cases, the infinite volume limit and the continuum limit do not commute. Furthermore, this signals a surprising UV/IR mixing, i.e., long distance sensitivity to short distance details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1995

Ultraviolet imaging telescope and optical emission-line observations of H II regions in M81

Images of the type Sab spiral galaxy M81 were obtained in far-UV and near-UV bands by the Ultraviolet Imaging Telescope (UIT) during the Astro-1 Spacelab mission of 1990 December. Magnitudes in the two UV bands are determined for 52 H II regions from the catalog of Petit, Sivan, & Karachentsev (1988). Fluxes of the H-alpha and H-beta emission lines are determined from CCD images. Extinctions for the brightest H II regions are determined from observed Balmer decrements. Fainter H II regions are assigned the average of published radio-H-alpha extinctions for several bright H II regions. The radiative transfer models of Witt, Thronson, & Capuano (1992) are shown to predict a relationship between Balmer Decrement and H-alpha extinction consistent with observed line and radio fluxes for the brightest 7 H II regions and are used to estimate the UV extinction. Ratios of Lyman continuum with ratios predicted by model spectra computed for initial mass function (IMF) slope equal to -1.0 and stellar masses ranging from 5 to 120 solar mass. Ages and masses are estimated by comparing the H-alpha and far-UV fluxes and their ratio with the models. The total of the estimated stellar masses for the 52 H II regions is 1.4 x 10(exp 5) solar mass. The star-formation rate inferred for M81 from the observed UV and H-alpha fluxes is low for a spiral galaxy at approximately 0.13 solar mass/yr, but consistent with the low star-formation rates obtained by Kennicutt (1983) and Caldwell et al. (1991) for early-type spirals.

Hill, Jesse K.↗

NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗

DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗