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 91 records · Page 5

Formation and Evolution Simulations of Saturn, Including Composition Gradients and Helium Immiscibility

In this paper, the formation of Saturn is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. Previous models are enhanced to include the dissolution of accreting planetesimals, composed of water ice, rock, and iron, in the gaseous envelope of the planet, leading to a nonuniform composition with depth. The immiscibility of helium in metallic hydrogen layers is also considered. The calculations start at a mass of 0.5 Earth masses and are extended to the present day. At 4.57 Gyr, the model, proceeding outward, has the following structure: (i) a central core composed of 100% heavy elements and molecules, (ii) a region with a decreasing heavy-element mass fraction, down to a value of 0.1, (iii) a layer of uniform composition with the helium mass fraction Y enhanced over the primordial value, (iv) a helium rain region with a gradient in Y, (v) an outer convective, adiabatic region with uniform composition in which Y is reduced from the primordial value, and (vi) the very outer layers where cloud condensation of the heavy elements occurs. Models of the distribution of heavy elements as a function of radius are compared with those derived to fit the observations of the Cassini mission, with rough qualitative agreement. The helium mass fraction in Saturn’s outer layers is estimated to be around 20%. Models are found that provide good agreement with Saturn’s intrinsic luminosity and radius.

79 ASTRONOMY AND ASTROPHYSICS↗

Molecular dynamics study of mechanical properties of carbon nanotube reinforced aluminum composites

Atomistic simulations were conducted to estimate the effect of the carbon nanotube (CNT) reinforcement on the mechanical behavior of CNT-reinforced aluminum (Al) nanocomposite. The periodic system of CNT-Al nanocomposite was built and simulated using molecular dynamics (MD) software LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). The mechanical properties of the nanocomposite were investigated by the application of uniaxial load on one end of the representative volume element (RVE) and fixing the other end. The interactions between the atoms of Al were modeled using embedded atom method (EAM) potentials, whereas Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used for the interactions among carbon atoms and these pair potentials are coupled with the Lennard-Jones (LJ) potential. The results show that the incorporation of CNT into the Al matrix can increase the Young’s modulus of the nanocomposite substantially. In the present case, i.e. for approximately 9 with % reinforcement of CNT can increase the axial Young’s modulus of the Al matrix up to 77 % as compared to pure Al.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Effect of Reaction Kinetics on Graphene-Edge Morphology and Composition

New simulations of graphene growth in flame environments are presented. The simulations employed a kinetic Monte Carlo (KMC) algorithm coupled to molecular mechanics (MM) geometry optimization to track individual graphenic species as they evolve. Focus was given to incorporation of five-member rings and resulting curvature and edge defects. The model code was re-written to be more computationally efficient enabling a larger set of simulations to be run, decreasing stochastic fluctuations in the averaged results. The model also included updated rate coefficients for graphene edge reactions recently published in the literature. Sensitivity analysis, enabled by the more efficient coding, confirmed the importance of the updated rates as well as providing further insight on key reactions controlling layer growth and morphology. The new KMC code was applied to simulate graphene-layer growth in the environments of laminar premixed flames. Finally, in these flame simulations, C-H ratios of evolving structures were computed and compared to those from experiment and from an alternate model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and numerical investigation of hydrate-guest molecule exchange kinetics

The 2012 Ignik Sikumi #1 Field Trial in Alaska was designed to test at field scale a methane hydrate production methodology that involves injecting carbon dioxide (CO 2 ) in situ to exchange with methane (CH 4 ) within a hydrate structure, thus releasing the methane for production. Since the completion of the field trial, experimental and numerical investigations have sought to better understand observations in terms of CH 4 production and changes to the hydrate bearing formation. Collectively, these insights may lead to more effective strategies for producing CH 4 from naturally occurring hydrates. This study presents results from a laboratory experiment performed at conditions experienced during the field trial (8.27 MPa and 275.15 K) to collect data relating to the CH 4 molecular exchange process. The experiment intentionally isolated the kinetic guest molecule exchange process from additional hydrate formation and mechanical changes to the hydrate bearing sand. Data were used to inform development of numerical models that were then used to fit key hydrate parameters (bound water saturation, kinetic formation constant, kinetic exchange constant, and preferential exchange weighting factors) and analyze experimental results. Results were in agreement with observations of preferential CO 2 /CH 4 guest molecule exchange at Ignik Sikumi #1 and from previous experimental studies while providing quantitative estimates of changing hydrate compositions. Additionally, simulations confirmed hydrate behavior and composition profiles otherwise indirectly evidenced by experimental results.

03 NATURAL GAS↗

Mapping Composition Evolution through Synthesis, Purification, and Depolymerization of Random Heteropolymers

Random heteropolymers (RHPs) consisting of three or more comonomers have been routinely used to synthesize functional materials. While increasing the monomer variety diversifies the side-chain chemistry, this substantially expands the sequence space and leads to ensemble-level sequence heterogeneity. Most studies have relied on monomer composition and simulated sequences to design RHPs, but the questions remain unanswered regarding heterogeneities within each RHP ensemble and how closely these simulated sequences reflect the experimental outcomes. Here, we quantitatively mapped out the evolution of monomer compositions in four-monomer-based RHPs throughout a design-synthesis-purification-depolymerization process. By adopting a Jaacks method, we first determined 12 reactivity ratios directly from quaternary methacrylate RAFT copolymerization experiments to account for the influences of competitive monomer addition and the reversible activation/deactivation equilibria. The reliability of in silico analysis was affirmed by a quantitative agreement (<4% difference) between the simulated RHP compositions and the experimental results. Furthermore, we mapped out the conformation distribution within each ensemble in different solvents as a function of monomer chemistry, composition, and segmental characteristics via high-throughput computation based on self-consistent field theory (SCFT). These comprehensive studies confirmed monomer composition as a viable design parameter to engineer RHP-based functional materials as long as the reactivity ratios are accurately determined and the livingness of RHP synthesis is ensured.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Element Composition Variation in Lead Zirconate Titanate upon the Ion-Plasma Deposition: Experiment and Simulation

The element composition modification is experimentally studied in PZT thin films obtained via the RF magnetron sputtering with variable working gas pressure. The experimental data were interpreted through the statistic simulation of thermalization and diffusion of sprayed Pb, Zr, and Ti atoms during the ion-plasma deposition of PZT thin films. The simulation data are shown to adequately describe (within the limits of 5%) the change in element composition at various pressures of the gas environment. The study conducted enables one to produce PZT thin films with a set element composition, which is essential for optimizing the electric and physical properties of solid solutions at the morphotropic phase boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

INTEGRATION OF COST MODELS AND PROCESS SIMULATION TOOLS FOR OPTIMUM COMPOSITE MANUFACTURING PROCESS

Manufacturing cost of resin transfer molded composite parts is significantly influenced by the cycle time, which is strongly related to the time for both filling and curing of the resin in the mold. The time for filling can be optimized by various injection strategies, and by suitably reducing the length of the resin flow distance during the injection. The curing time can be reduced by the usage of faster curing resins, but it requires a high pressure injection equipment, which is capital intensive. Predictive manufacturing simulation tools that are being developed recently for composite materials are able to provide various scenarios of processing conditions virtually well in advance of manufacturing the parts. In the present study, we integrate the cost models with process simulation tools to study the influence of various parameters such as injection strategies, injection pressure, compression control to minimize high pressure injection, resin curing rate, and demold time on the manufacturing cost as affected by the annual part volume. A representative automotive component was selected for the study and the results are presented in this paper

42 ENGINEERING↗

Uncertainty quantification in multiscale simulation of woven fiber composites

We report woven fiber composites have been increasingly employed as light-weight materials in aerospace, construction, and transportation industries due to their superior properties. These materials possess a hierarchical structure that necessitates the use of multiscale simulations in their modeling. To account for the inherent uncertainty in materials, such simulations must be integrated with statistical uncertainty quantification (UQ) and propagation (UP) methods. However, limited advancement has been made in this regard due to the significant computational costs and complexities in modeling spatially correlated structural variations coupled at different scales. In this work, a non-intrusive approach is proposed for multiscale UQ and UP to address these limitations. We introduce the top-down sampling method that allows to model non-stationary and continuous (but not differentiable) spatial variations of uncertainty sources by creating nested random fields (RFs) where the hyperparameters of an ensemble of RFs is characterized by yet another RF. We employ multi-response Gaussian RFs in top-down sampling and leverage statistical techniques (such as metamodeling and dimensionality reduction) to address the considerable computational costs of multiscale simulations. We apply our approach to quantify the uncertainty in a cured woven composite due to spatial variations of yarn angle, fiber volume fraction, and fiber misalignment angle. Finally, our results indicate that, even in linear analysis, the effect of uncertainty sources on the material’s response could be significant.

42 ENGINEERING↗

Constraining Aging Processes of Black Carbon in the Community Atmosphere Model Using Environmental Chamber Measurements

Abstract The direct radiative forcing of black carbon aerosol (BC) on the Earth system remains unsettled, largely due to the uncertainty with physical properties of BC throughout their lifecycle. Here we show that ambient chamber measurements of BC properties provide a novel constraint on the crude BC aging representation in climate models. Observational evidence for significant absorption enhancement of BC can be reproduced when the aging processes in the four‐mode version of the Modal Aerosol Module (MAM4) aerosol scheme in the Community Atmosphere Model version 5 are calibrated by the recent in situ chamber measurements. An observation‐based scaling method is developed in the aging timescale calculation to alleviate the influence of biases in the simulated model chemical composition. Model sensitivity simulations suggest that the different monolayer settings in the BC aging parameterization of MAM4 can cause as large as 26% and 24% differences in BC burden and radiative forcing, respectively. We also find that an increase in coating materials (e.g., sulfate and secondary organic aerosols) reduces BC lifetime by increasing the hygroscopicity of the mixture but enhances its absorption, resulting in a net increase in BC direct radiative forcing. Our results suggest that accurate simulations of BC aging processes as well as other aerosol species are equally important in reducing the uncertainty of BC forcing estimation.

54 ENVIRONMENTAL SCIENCES↗

Oxidation of Glycolate in the Defense Waste Processing Facility (DWPF) Recycle Collection Tank - 20305

The Savannah River Site's Defense Waste Processing Facility (DWPF) operations are being upgraded with the introduction of the Nitric-Glycolic Flowsheet. Glycolic acid has been shown superior to formic acid as the reducing acid used during chemical processing. The new flowsheet improves or maintains necessary parameters such as 1) reduction of mercury, 2) adjustment of feed rheology and 3) adjustment of melter oxidation/reduction potential. Further, the potential for catalytic hydrogen generation in DWPF processing is virtually eliminated. DWPF process condensates are collected and returned to the SRS Concentration, Storage and Transfer Facilities (CSTF). The Recycle Collection Tank (RCT) collects off-gas condensate during chemical processing, vitrification, and other unit operations performed in DWPF and is the singular return vessel delivering recycle effluent back to CSTF. Each batch of recycle may contain a small amount of glycolate from chemical processing and melter off-gas condensates. To avoid potential flammability issues due to thermolysis of glycolate in the CSTF, chemical oxidation within the RCT has been investigated as an option for mitigating the transfer of glycolate. Sodium permanganate has been down-selected as the best option for oxidation of glycolate. Testing was performed using both 2-L and 22-L reactors (16,800:1 and 1,530:1 scale by volume) with non-radioactive waste simulants to approximate the expected RCT compositions. RCT simulants were evaluated at various process pH and temperature conditions. Also, RCT operations, namely the sequence of addition of corrosion inhibitors (NaOH and NaNO{sub 2}) versus a permanganate strike, were evaluated. Glycolate was introduced via a sludge simulant to mimic both expected entrainment and abnormal process foam-over conditions - the range being between 68 and 5100 mg/kg glycolate. Glycolate destruction was monitored by ion chromatography (IC). The corresponding manganese behavior was monitored in real-time using in situ ultraviolet-visible (UV-Vis) spectroscopy. RCT glycolate content can be reduced to below the IC detection limit within 90 minutes for all concentrations investigated. Ion Chromatography analysis revealed that under alkaline conditions, glycolate is primarily oxidized to oxalate with no significant formation of CO{sub 2} or carbonate, and nitrite is not oxidized to nitrate. Initially, complete oxidation of organics species and nitrite was assumed. Determination of the mechanistic chemical reaction has allowed the required amount of permanganate to be more accurately predicted and the total addition to be significantly reduced. UV-Vis measurements reveal that permanganate (Mn{sup 7+}) is reduced to manganate (Mn{sup 6+}) in the RCT. The oxidant stoichiometry is defined by using the initial permanganate to glycolate (P/G) molar ratio. At low initial glycolate concentration (68 and 140 mg/kg), the minimum required initial permanganate to glycolate (P/G) molar ratio was found to be 5-6. With high initial glycolate concentrations (5100 mg/kg) a lower (P/G) molar ratio of ∼2.5 was needed. The final portion of this effort supporting the nitric/glycolic flowsheet will be to test actual (fully radioactive) RCT samples as per the above simulant tests. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Broad Grid of 2D Kilonova Emission Models

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M ⊙ and 0.05 to 0.3c, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (Y e ) ejecta with a robust r-process composition and either a spherical or lobed high-Y e ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗

Measuring the composition-curvature coupling in binary lipid membranes by computer simulations

The coupling between local composition fluctuations in binary lipid membranes and curvature affects the lateral membrane structure. We propose an efficient method to compute the composition-curvature coupling in molecular simulations and apply it to two coarse-grained membrane models—a minimal, implicit-solvent model and the MARTINI model. Both the weak-curvature behavior that is typical for thermal fluctuations of planar bilayer membranes as well as the strong-curvature regime corresponding to narrow cylindrical membrane tubes are studied by molecular dynamics simulation. The simulation results are analyzed by using a phenomenological model of the thermodynamics of curved, mixed bilayer membranes that accounts for the change of the monolayer area upon bending. Additionally the role of thermodynamic characteristics such as the incompatibility between the two lipid species and asymmetry of composition are investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

Supplementary Data for "Evaluation of the Economic Implications of Varied Pressure Drawdown Strategies Generated Using a Real-time, Rapid Predictive, Multi-fidelity Model for Unconventional Oil and Gas Wells" by Bello, K., Vikara, D., Sheriff, A., Viswanathan, H., Carr, T., Sweeney, M., O'Malley, D., Marquis, M., Vactor, R.T., and Cunha, L.

The Bello et al. study evaluates the impact of contrasting pressure drawdown on gas productivity and the resulting economics of a well in the Marcellus Shale of the Appalachian Basin. This research applies a techno-economic analysis approach to help identify potential ways pressure management strategies can be used to improve cumulative recovery of hydraulically fractured horizontal wells while maintaining project profitability. Gas production forecast outlook scenarios of the Marcellus Shale Energy and Environment Laboratory Laboratory's MIP-3H well were generated under varying pressure drawdowns using two approaches: 1) a novel physics-informed machine learning (PIML) workflow and 2) via traditional reservoir simulation in Computer Modeling Group’s (CMG) GEM Compositional & Unconventional Simulator. Cash flow and other economic metrics of interest were compiled on the production outlook using the U.S. Department of Energy's (DOE) National Energy Technology Laboratory (NETL) Unconventional Shale Well Economic Model (UShWEM).The sheets within this Microsoft ExcelTM workbook provide the economic metric outputs for the baseline condition and the one-at-a-time (OAT) sensitivity analysis of UShWEM's input parameters for each of the production scenarios evaluated.

Fracture Network Model↗