Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 865 records · Page 48

Better calibration of cloud parameterizations and subgrid effects increases the fidelity of the E3SM Atmosphere Model version 1

Abstract. Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

54 ENVIRONMENTAL SCIENCES↗

Composition Measurements of the HLW APPS Glasses

This report provides the results from the chemical analyses of the glass compositions of the High Activity Waste Aspen Process Performance Simulation (HLW APPS) study glasses, a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory. These data will be used in the development, validation, and implementation of enhanced property/composition models for waste glass vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the quenched glasses to allow for comparisons with targeted compositions. The relative differences between the targeted and measured concentrations of F- for three glasses and B 2 O 3 for five glasses were greater than 10%. These results can be used in further characterization of this series of glasses, including the normalization of Product Consistency Test results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Inclusion Characterization and Formation Mechanisms in Spring Steel Deoxidized by Silicon

To elucidate the mechanism of formation of inclusions in microalloyed spring steel deoxidized by Si, their number, morphology, and chemical compositions were analyzed by electron probe X-ray microanalysis. Wavelength-dispersive spectrometry and Monte Carlo simulation were used to confirm the existence of inherent FeO in the inclusions. Additions of Nb and V did not have an observable influence on the formation of inclusions. Based on the major chemical compositions, the inclusions in spring steel were classified into four types: CaO-SiO{sub 2}, Al{sub 2}O{sub 3}-SiO{sub 2}, Al{sub 2}O{sub 3}-SiO{sub 2}-CaO, and MnS. Most CaO-SiO{sub 2} inclusions were around 10 μm in diameter, with their composition close to that of the refining slag and liquidus temperatures all below 1400 °C. The Al{sub 2}O{sub 3}-SiO{sub 2} inclusions were mainly attributed to deoxidization products. They were usually multiphased with liquidus temperatures higher than 1500 °C and smaller than 10 μm in diameter. Most Al{sub 2}O{sub 3}-SiO{sub 2}-CaO inclusions with high Al{sub 2}O{sub 3} and SiO{sub 2} contents originated from the coalescence of CaO-SiO{sub 2} and Al{sub 2}O{sub 3}-SiO{sub 2} inclusions; those with high CaO and SiO{sub 2} contents were considered to form via reduction of CaO-SiO{sub 2} inclusions by Al dissolved in the steel.

36 MATERIALS SCIENCE↗

Investigation of precipitate refinement in Mg alloys by an analytical composite failure model

An analytical model is developed to simulate precipitate refinement in second phase strengthened magnesium alloys. The model is developed based on determination of the stress fields inside elliptical precipitates embedded in a rate dependent inelastic matrix. The stress fields are utilized to determine the failure mode that governs the refinement behavior. Using an AZ31 Mg alloy as an example, the effects the applied load, aspect ratio and orientation of the particle is studied on the macroscopic failure of a single α-Mg17Al12 precipitate. Additionally, a temperature dependent version of the corresponding constitutive law is used to incorporate the effects of temperature. In plane strain compression, an extensional failure mode always fragments the precipitates. The critical strain rate at which the precipitates start to fail strongly depends on the orientation of the precipitate with respect to loading direction. The results show that the higher the aspect ratio is, the easier the precipitate fractures. Precipitate shape is another factor influencing the failure response. In contrast to elliptical precipitates with high aspect ratio, spherical precipitates are strongly resistant to sectioning. In pure shear loading, in addition to the extensional mode of precipitate failure, a shearing mode may get activated depending on orientation and aspect ratio of the precipitate. The effect of temperature in relation to strain rate was also verified for plane strain compression and pure shear loading cases.

Tabei, Ali↗

Molecular Dynamics Simulation on the Effect of MgO/Al{sub 2}O{sub 3} Ratio on Structure and Properties of Blast Furnace Slag Under Different Basicity Conditions

The SiO{sub 2}-Al{sub 2}O{sub 3}-CaO-MgO is the basic structural system of blast furnace slag and the composition directly affects the performance of the slag. Molecular dynamics simulations were carried out to analyze the local structure, structural unit, bond angle, transport properties, and enthalpy of the slag with the increase of MgO/Al{sub 2}O{sub 3} mass ratio under different basicity conditions. It was found that the change of MgO/Al{sub 2}O{sub 3} ratio does not affect the short-range order structure, but it will reduce the overall stability of the network structure. Through the analysis of the structural unit of the slag, it was found that the polymerization degree of the system decreases with the increase of MgO/Al{sub 2}O{sub 3} ratio, indicating that the complex network structure of the system is partially depolymerized. Besides, the self-diffusion coefficients of each ion were concluded and the magnitudes were observed to be in the following order: Mg{sup 2+}>Ca{sup 2+}>Al{sup 3+}>O{sup 2-}>Si{sup 4+}. The diffusivity of the slag increases with the increase of MgO/Al{sub 2}O{sub 3} ratio and the corresponding viscosity decreases. And the enthalpy of the slag is increased with the increase of MgO/Al{sub 2}O{sub 3} ratio or basicity, therefore the fuel consumption should also be properly adjusted during the operation of blast furnace to ensure the target temperature.

36 MATERIALS SCIENCE↗

Constitutive Modeling for Predicting High-Temperature Flow Behavior in Aluminum 5083+10 Wt Pct SiC{sub p} Composite

A constitutive model capable of predicting material flow behavior with high precision is essential for optimizing secondary processing parameters through simulation techniques. In this study, the hot deformation behavior of aluminum 5083+10 wt pct SiC particulate composite was predicted using constitutive equations based on the modified Johnson–Cook (JC), modified Zerilli–Armstrong (ZA), and strain-compensated Arrhenius models. The models were established on the basis of the true stress–strain values obtained from an isothermal hot compression test conducted on the INSTRON 8801 universal tensile testing machine under a temperature range of 473 K to 773 K and strain rate of 0.01 to 10 s{sup −1}. The prediction ability of the modified models was compared by calculating the correlation coefficient (R), average absolute relative error, and relative error. All the models could precisely predict the hot flow behavior of the composite. The modified ZA model had the highest accuracy. The JC model required the least number of material constants and had the lowest calculation time, followed by the modified ZA and strain-compensated Arrhenius models.

36 MATERIALS SCIENCE↗

High-throughput ion irradiation of additively manufactured compositionally complex alloys

Several advanced nuclear reactor designs promise efficiency and safety improvements over the current reactor fleet but are limited by the current set of ASME code-qualified materials. Novel alloys including high-entropy alloys (HEAs), and more broadly compositionally complex alloys (CCAs), have shown promising irradiation-tolerance. However, the vast range of alloy compositions adds to an already time-consuming alloy development process. In this study, to accelerate the development of novel alloys for nuclear applications, a high-throughput (HTP) methodology has been employed. Additive manufacturing has been used to produce a compositional array of unary, binary, ternary, and quaternary alloys, including several CCAs, which span the Cr-Fe-Mn-Ni composition space. Additionally, the compositional array was homogenized at 1000°C for 24 hours and each sample was irradiated using 4-MeV Ni2+ ions at room temperature to a peak damage of 50 dpa, as estimated using SRIM, at the University of Wisconsin Ion Beam Laboratory. A custom XY stage was built to accommodate the large compositional array and half of each sample was masked during irradiation enabling both the irradiated and unirradiated properties of each alloy to be characterized side-by-side. Each alloy was characterized using X-ray fluorescence (XRF), X-ray diffraction (XRD), scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS), and nanoindentation. CALPHAD simulations spanning the entire Cr-Fe-Mn-Ni composition space at 1000°C were performed to compare predicted equilibrium phases with phases identified experimentally from the unirradiated regions of each alloy. Nanoindentation measurements indicate radiation-induced hardening ranging from ~1-1.5 GPa in each of Cr-Fe-Mn-Ni CCAs, which is relatively insensitive to modest changes in alloy composition and comparable to hardening observed in neutron irradiated Cr-Fe-Mn-Ni CCAs in the literature. Overall, a substantial time savings was realized by employing HTP synthesis, irradiation, and characterization in this study compared to conventional techniques, the implications of which are discussed.

36 MATERIALS SCIENCE↗

An ICME Modeling Framework for Titanium/Tungsten-Carbide Metal Matrix Composites

This report describes a collaborative project to develop a validated, predictive model for the high temperature mechanical properties of a titanium-matrix, tungsten-carbide/cobalt-reinforced metal matrix composite. The modeling approach was to first develop a detailed, microstructural model linking the material structure and the interfacial debonding properties to the effective properties of the material. The project then completed a throughput simulation campaign to generate a large number of simulations for discrete microstructures and different debonding parameters. Finally, the project trained a fast, Gaussian process surrogate model against this simulation database to provide a quick model linking the material compositions, structure, and processing parameters to the resulting material properties. This model was validated against high temperature tensile test data on a few particular composite compositions. The tests validate the model predictions for ultimate tensile strength and uniform elongation/ductility, meaning the final surrogate model can now be used to tune the material composition and processing parameters to identify optimal composite compositions for particular applications.

36 MATERIALS SCIENCE↗

Selectivity trend of gas separation through nanoporous graphene

By means of molecular dynamics (MD) simulations, we demonstrate that porous graphene can efficiently separate gases according to their molecular sizes. The flux sequence from the classical MD simulation is H{sub 2}>CO{sub 2}≫N{sub 2}>Ar>CH{sub 4}, which generally follows the trend in the kinetic diameters. This trend is also confirmed from the fluxes based on the computed free energy barriers for gas permeation using the umbrella sampling method and kinetic theory of gases. Both brute-force MD simulations and free-energy calcualtions lead to the flux trend consistent with experiments. Case studies of two compositions of CO{sub 2}/N{sub 2} mixtures further demonstrate the separation capability of nanoporous graphene. - Graphical abstract: Classical molecular dynamics simulations show the flux trend of H{sub 2}>CO{sub 2}≫N{sub 2}>Ar>CH{sub 4} for their permeation through a porous graphene, in excellent agreement with a recent experiment. - Highlights: • Classical MD simulations show the flux trend of H{sub 2}>CO{sub 2}≫N{sub 2}>Ar>CH{sub 4} for their permeation through a porous graphene. • Free energy calculations yield permeation barriers for those gases. • Selectivities for several gas pairs are estimated from the free-energy barriers and the kinetic theory of gases. • The selectivity trend is in excellent agreement with a recent experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adding stable carbon isotopes improves model representation of the role of microbial communities in peatland methane cycling

Abstract Climate change is expected to have significant and uncertain impacts on methane (CH 4 ) emissions from northern peatlands. Biogeochemical models can extrapolate site‐specificCH 4 measurements to larger scales and predict responses of CH 4 emissions to environmental changes. However, these models include considerable uncertainties and limitations in representing CH 4 production, consumption, and transport processes. To improve predictions of CH 4 transformations, we incorporated acetate and stable carbon (C) isotopic dynamics associated with CH 4 cycling into a biogeochemistry model, DNDC. By including these new features, DNDC explicitly simulates acetate dynamics and the relative contribution of acetotrophic and hydrogenotrophic methanogenesis (AM and HM) to CH 4 production, and predicts the C isotopic signature (δ 13 C) in soil C pools and emitted gases. When tested against biogeochemical and microbial community observations at two sites in a zone of thawing permafrost in a subarctic peatland in Sweden, the new formulation substantially improved agreement with CH 4 production pathways and δ 13 C in emitted CH 4 (δ 13 C‐CH 4 ), a measure of the integrated effects of microbial production and consumption, and of physical transport. We also investigated the sensitivity of simulated δ 13 C‐CH 4 to C isotopic composition of substrates and, to fractionation factors for CH 4 production (α AM and α HM ), CH 4 oxidation (α MO ), and plant‐mediated CH 4 transport (α TP ). The sensitivity analysis indicated that the δ 13 C‐CH 4 is highly sensitive to the factors associated with microbial metabolism (α AM , α HM , and α MO ). The model framework simulating stable C isotopic dynamics provides a robust basis for better constraining and testing microbial mechanisms in predicting CH 4 cycling in peatlands.

54 ENVIRONMENTAL SCIENCES↗

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗

DETECTING THE RAPIDLY EXPANDING OUTER SHELL OF THE CRAB NEBULA: WHERE TO LOOK

We present a range of steady-state photoionization simulations, corresponding to different assumed shell geometries and compositions, of the unseen postulated rapidly expanding outer shell to the Crab Nebula. The properties of the shell are constrained by the mass that must lie within it, and by limits to the intensities of hydrogen recombination lines. In all cases the photoionization models predict very strong emissions from high ionization lines that will not be emitted by the Crab's filaments, alleviating problems with detecting these lines in the presence of light scattered from brighter parts of the Crab. The near-NIR [Ne VI] {lambda}7.652 {mu}m line is a particularly good case; it should be dramatically brighter than the optical lines commonly used in searches. The C IV {lambda}1549 doublet is predicted to be the strongest absorption line from the shell, which is in agreement with Hubble Space Telescope observations. We show that the cooling timescale for the outer shell is much longer than the age of the Crab, due to the low density. This means that the temperature of the shell will actually ''remember'' its initial conditions. However, the recombination time is much shorter than the age of the Crab, so the predicted level of ionization should approximate the real ionization. In any case, it is clear that IR observations present the best opportunity to detect the outer shell and so guide future models that will constrain early events in the original explosion.

79 ASTRONOMY AND ASTROPHYSICS↗

Architectural Improvements and New Processing Tools for the Open XAL Online Model

The online model is the component of Open XAL providing accelerator modeling, simulation, and dynamic synchronization to live hardware. Significant architectural changes and feature additions have been recently made in two separate areas: 1) the managing and processing of simulation data, and 2) the modeling of RF cavities. Simulation data and data processing have been completely decoupled. A single class manages all simulation data while standard tools were developed for processing the simulation results. RF accelerating cavities are now modeled as composite structures where parameter and dynamics computations are distributed. The beam and hardware models both maintain their relative phase information, which allows for dynamic phase slip and elapsed time computation.

Accelerator Control↗

Modeling and Depletion Simulations for a High Flux Isotope Reactor Cycle with a Representative Experiment Loading

The purpose of this report is to document a high-fidelity VESTA/MCNP High Flux Isotope Reactor (HFIR) core model that features a new, representative experiment loading. This model, which represents the current, high-enriched uranium fuel core, will serve as a reference for low-enriched uranium conversion studies, safety-basis calculations, and other research activities. A new experiment loading model was developed to better represent current, typical experiment loadings, in comparison to the experiment loading included in the model for Cycle 400 (operated in 2004). The new experiment loading model for the flux trap target region includes full length 252 Cf production targets, 75 Se production capsules, 63 Ni production capsules, a 188 W production capsule, and various materials irradiation targets. Fully loaded 238 Pu production targets are modeled in eleven vertical experiment facilities located in the beryllium reflector. Other changes compared to the Cycle 400 model are the high-fidelity modeling of the fuel element side plates and the material composition of the control elements. Results obtained from the depletion simulations with the new model are presented, with a focus on time-dependent isotopic composition of irradiated fuel and single cycle isotope production metrics.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas

Carbon-hydrogen plasmas and hydrocarbon materials are of broad interest to laser shock experimentalists, high energy density physicists, and astrophysicists. Accurate equations of state (EOSs) of hydrocarbons are valuable for various studies from inertial confinement fusion to planetary science. By combining path integral Monte Carlo (PIMC) results at high temperatures and density functional theory molecular dynamics results at lower temperatures, we compute the EOSs for hydrocarbons from simulations performed at 1473 separate (ρ, T)-points distributed over a range of compositions. These methods accurately treat electronic excitation effects with neither adjustable parameter nor experimental input. PIMC is also an accurate simulation method that is capable of treating many-body interaction and nuclear quantum effects at finite temperatures. These methods therefore provide a benchmark-quality EOS that surpasses that of semi-empirical and Thomas-Fermi-based methods in the warm dense matter regime. By comparing our first-principles EOS to the LEOS 5112 model for CH, we validate the specific heat assumptions in this model but suggest that the Grüneisen parameter is too large at low temperatures. Based on our first-principles EOSs, we predict the principal Hugoniot curve of polystyrene to be 2%-5% softer at maximum shock compression than that predicted by orbital-free density functional theory and SESAME 7593. By investigating the atomic structure and chemical bonding of hydrocarbons, we show a drastic decrease in the lifetime of chemical bonds in the pressure interval from 0.4 to 4 megabar. We find the assumption of linear mixing to be valid for describing the EOS and the shock Hugoniot curve of hydrocarbons in the regime of partially ionized atomic liquids. We make predictions of the shock compression of glow-discharge polymers and investigate the effects of oxygen content and C:H ratio on its Hugoniot curve. In conclusion, our full suite of first-principles simulation results may be used to benchmark future theoretical investigations pertaining to hydrocarbon EOSs and should be helpful in guiding the design of future experiments on hydrocarbons in the gigabar regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preliminary Formation Analysis for Compressed Air Energy Storage in Depleted Natural Gas Reservoirs: A Study for the DOE Energy Storage Systems Program

The purpose of this study is to develop an engineering and operational understanding of CAES performance for a depleted natural gas reservoir by evaluation of relative permeability effects of air, water and natural gas in depleted natural gas reservoirs as a reservoir is initially depleted, an air bubble is created, and as air is initially cycled. The composition of produced gases will be evaluated as the three phase flow of methane, nitrogen and brine are modeled. The effects of a methane gas phase on the relative permeability of air in a formation are investigated and the composition of the produced fluid, which consists primarily of the amount of natural gas in the produced air are determined. Simulations of compressed air energy storage (CAES) in depleted natural gas reservoirs were carried out to assess the effect of formation permeability on the design of a simple CAES system. The injection of N 2 (as a proxy to air), and the extraction of the resulting gas mixture in a depleted natural gas reservoir were modeled using the TOUGH2 reservoir simulator with the EOS7c equation of state. The optimal borehole spacing was determined as a function of the formation scale intrinsic permeability. Natural gas reservoir results are similar to those for an aquifer. Borehole spacing is dependent upon the intrinsic permeability of the formation. Higher permeability allows increased injection and extraction rates which is equivalent to more power per borehole for a given screen length. The number of boreholes per 100 MW for a given intrinsic permeability in a depleted natural gas reservoir is essentially identical to that determined for a simple aquifer of identical properties. During bubble formation methane is displaced and a sharp N 2 -methane boundary is formed with an almost pure N 2 gas phase in the bubble near the borehole. During cycling mixing of methane and air occurs along the boundary as the air bubble boundary moves. The extracted gas mixture changes as a function of time and proximity of the bubble boundary to the well. For all simulations reported here, with a formation radius above 50 m the maximum methane composition in the produced gas phase was less than 0.5%. This report provides an initial investigation of CAES in a depleted natural gas reservoir, and the results will provide useful guidance in CAES system investigation and design in the future.

Gardner, William Payton↗

Composite entanglement topology and extensional rheology of symmetric ring-linear polymer blends

Extensive molecular simulations are applied to characterize the equilibrium dynamics, entanglement topology, and nonlinear extensional rheology of symmetric ring-linear polymer blends with systematically varied ring fraction ϕR. Chains with degree of entanglement Z≈14 are mixed to produce 10 well-entangled systems with ϕR varying from neat linear to neat ring melts. Z is large enough that except for very large ϕR, the rings are threaded by multiple linear chains in equilibrium. Primitive path analysis is used to visualize and quantify the structure of the composite ring-linear entanglement network. We measure the quantity of ring-linear threading and linear-linear entanglement as a function of ϕR and identify with simple arguments a ring fraction ϕR≈0.4 where the topological constraints of the composite entanglement network are maximized. These topological analyses are used to rationalize the ϕR-dependence of ring and linear chain dynamics, conformations, and rheology. Simulations of startup uniaxial elongation flows demonstrate the extensional stress overshoot observed in recent filament stretching experiments and characterize how it depends on the blend composition and entanglement topology. The overshoot is driven by an overstretching and recoil of ring polymers due to the convective unthreading of rings from linear chains.

O’Connor, Thomas C.↗

Transport properties of nanocomposite and its simulation with L-R-C circuit

The nano particles are represented in this communication by L-R-C equivalent circuit. The dc current voltage characteristics (CVC) of the proposed circuit have simulated using Circuit-Maker ® 2000. Experimental investigation on ZnO nano-composite with capping material gum acacia shows similar CVC. NPs are represented by C-R combinations to manifest the Coulomb blockade effect of a quantum dot. The capping material is represented by an inductor along with a resistance in series. Nine NPs with capping matrix are simulated. The dc current voltage characteristics (CVC) and gross feature of polarization nature obtained by experiment and simulation study are consistent.

77 NANOSCIENCE AND NANOTECHNOLOGY↗