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

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimating outflow masses and velocities in merger simulations: Impact of r -process heating and neutrino cooling

Here, the determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fall-back timescale of any bound ejecta. We derive a model for the asymptotic velocity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full three-dimensional fluid simulations of the outflows, as nonlocal effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long three-dimensional simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

Dim, Odera↗

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of five subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the five pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Temperature- and composition-dependent hydrogen diffusivity in palladium from statistically-averaged molecular dynamics

Solid-state hydrogen storage materials undergo complex phase transformations whose kinetics is often limited by hydrogen diffusion. Among metal hydrides, palladium hydride undergoes a diffusional phase transformation upon hydrogen uptake, during which the hydrogen diffusivity varies with hydrogen composition and temperature. Here we perform robust statistically-averaged molecular dynamics simulations to obtain a well-converged analytical expression for hydrogen diffusivity in bulk palladium that is valid throughout all stages of the reaction. Our studies confirm significant dependence of the diffusivity on composition and temperature that elucidate key trends in the available experimental measurements. Whereas at low hydrogen compositions, a single process dominates, at high hydrogen compositions, diffusion is found to exhibit behavior consistent with multiple hopping barriers. Further analysis, supported by nudged elastic band computations, suggests that the multi-barrier diffusion can be interpreted as two distinct mechanisms corresponding to hydrogen-rich and hydrogen-poor local environments.

36 MATERIALS SCIENCE↗

A hybrid origin of the Martian crustal dichotomy: Degree-1 convection antipodal to a giant impact

The Martian crustal dichotomy is the stark ~5 km difference in surface elevation and ~26 km difference in crustal thickness between the northern lowlands and southern highlands that originated within 100s of Myr of Mars’ formation. The origin of the dichotomy has broad implications for the geodynamic history of Mars, but purely exogenic or endogenic theories so far cannot explain all of the large scale geophysical observations associated with dichotomy formation. A giant impact can produce the shape and slope of the dichotomy boundary, but struggles to explain Mars’ remanent crustal magnetic signatures and the ultimate formation of Tharsis. Degree-1 mantle convection can relate the crustal dichotomy to the formation of Tharsis, but does not explain the elliptical dichotomy shape and must be initiated by a large pre-existing viscosity jump in the mantle. We propose a hybrid model of dichotomy formation in which a giant impact induces degree-1 convection with an upwelling antipodal to the impact site. In this scenario, a giant impact in the northern hemisphere excavates crust, creating an initial difference in crustal thickness and possibly composition between the two hemispheres. Over 10s to 100s of Myr, the dominant upwelling(s) would migrate to be under the thicker, insulating crust in the southern hemisphere, generating melt that further thickens the southern crust. We examine this process using 3-D mantle convection simulations, and find that a hemispherical difference in crustal thickness and composition caused by a giant impact can induce degree-1 convection with the upwelling(s) antipodal to the impact site in <100 Myr.

58 GEOSCIENCES↗

Exploration of oxidative chemistry and secondary organic aerosol formation in the Amazon during the wet season: explicit modeling of the Manaus urban plume with GECKO-A

The GoAmazon 2014/5 field campaign took place in Manaus, Brazil, and allowed the investigation of the interaction between background-level biogenic air masses and anthropogenic plumes. We present in this work a box model built to simulate the impact of urban chemistry on biogenic secondary organic aerosol (SOA) formation and composition. An organic chemistry mechanism is generated with the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A) to simulate the explicit oxidation of biogenic and anthropogenic compounds. A parameterization is also included to account for the reactive uptake of isoprene oxidation products on aqueous particles. The biogenic emissions estimated from existing emission inventories had to be reduced to match measurements. The model is able to reproduce ozone and NO x for clean and polluted situations. The explicit model is able to reproduce background case SOA mass concentrations but does not capture the enhancement observed in the urban plume. The oxidation of biogenic compounds is the major contributor to SOA mass. A volatility basis set (VBS) parameterization applied to the same cases obtains better results than GECKO-A for predicting SOA mass in the box model. The explicit mechanism may be missing SOA-formation processes related to the oxidation of monoterpenes that could be implicitly accounted for in the VBS parameterization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Techno-Economic Wind Blade Manufacturing Model to Identify Opportunities for Cost Improvements Phase II IACMI Project 4.6/4.8

In IACMI Project 4.6 and IACMI Project 4.8, an Excel-based Techno-Economic Model (TEM) of the manufacturing process for composite wind turbine blades and a DELMIA Factory Flow Simulation of a generic wind blade manufacturing facility was developed. Together, these two tools provide a combined economic modeling capability that accounts for the material, labor, overhead and full-lifecycle operating costs associated with wind blade manufacturing as well as the impact of process flow and factory layout on overall manufacturing efficiency. The tools provide a novel means of detailed comparative analysis of the economic feasibility of proposed technologies and process changes for blade manufacturing. The modeling tools were developed with close support from members of industry and visits to multiple blade manufacturing facilities. With industry oversight, a detailed generalized manufacturing process plan and facility layout were developed with manufacturing parameters, material costs and economic factors based on historical data. Dassault Systèmes and the University of Texas at Dallas (UTD) contributed to the development of the Techno-Economic Model by providing macros to enable the generation of Bill of Material (BOM) data from a 3D blade design in either CATIA or NuMAD format, respectively. The TEM was built with the capability to directly import a Bill of Materials for economic analysis, and with the addition of the macros provided by Dassault and UTD, the TEM can directly import blade designs from both CATIA and NuMAD file formats. The modeling tools developed in Project 4.6 were used to investigate four wind blade manufacturing concepts in detail and select one to explore with laboratory-scale experimentation in Project 4.8. The four manufacturing concepts that were investigated were down-selected by the full project team from a larger list of concepts. The selections were made based on a number of criteria ranking viability and level of interest for each concept. The ‘One-Step Close’ manufacturing concept was ultimately selected for investigation in Project 4.8 and the demonstration was performed at the NREL CoMET facility. The TPI advanced manufacturing facility in Warren, RI contributed the production of several prototype components, the designs for which were developed by Janicki Industries. The demonstration project provided clear indication of the viability of the One-Step Close manufacturing concept for blade manufacturing and good validation of the Techno-Economic Model’s prediction of its economic impact.

17 WIND ENERGY↗

Cs retention and diffusion in C-S-H at different Ca/Si ratio

Cement and concrete are commonly used in the construction of repository sites for radioactive wastes. The correct isolation of those contaminants requires good adsorption and low diffusion rates. Both parameters are highly affected by many factors, such as pH, temperature or composition. The large variability of experimental conditions and formulations makes it extraordinarily difficult to tackle the influence of each of them independently in experimental samples. To this effect, molecular dynamics simulations have been employed in this study to investigate the role of the composition in the capacity to retain Cs and diffusivity of these ions in calcium silicate hydrate (C-S-H) gel pores. The results indicate that the adsorption of Cs ions is worsened at high Ca/Si ratios due to a lower interaction of the cations with the C-S-H surface, while hydrophilicity of the C-S-H nanopore rises, resulting in higher long-range ordering and lower diffusion coefficients of Cs ions.

36 MATERIALS SCIENCE↗

Chemical Imaging of Atmospheric Particles

Airborne particles are very dynamic and highly reactive components of the Earth’s atmosphere. Their high surface area and water content provide a unique reaction environment for multiphase chemistry that continually modifies particle composition and properties that consequently impact air quality as well as concentrations of gas-phase species. By absorbing and scattering solar and terrestrial radiation, particles directly influence the planet’s radiative balance. Their indirect effects include modifying the nucleation, lifetime, and physical properties of clouds. Due to the sensitivity of the atmospheric environment to all these variables, fundamental studies of chemical transformations of atmospheric particles, their sources, continuously evolving composition, and physical properties are of highest research priority. Accurate descriptions of particles and their effects in the atmosphere require comprehensive information not only on the particle-type populations and their size distributions and concentrations, but also on the diversity and the spatial heterogeneity of chemical components within individual particles. Developments and applications of modern chemical imaging approaches for off-line characterization of atmospheric particles have been at the forefront of modern experimental studies and have resulted in a transformative impact in atmospheric chemistry and physics. This Account presents a synopsis of recent advances in chemical imaging of atmospheric particles collected on substrates during field and laboratory experiments. The unique advantage of chemical imaging methods is that they simultaneously provide two analytical measurements: imaging of particles to assess variability in their individual sizes and morphology, as well as particle specific speciation of their composition and spatial heterogeneity of different chemical components within individual particles. We also highlight analytical chemistry approaches that enable chemical imaging of particles with different levels of elemental and molecular specificity, including applications of multimodal methodologies where the same or similar groups of particles are probed by two or more complementary techniques. These approaches provide unique experimental insights on the nature and sources of particles, understanding their physical properties, atmospheric reactivity, and transformations. Chemical imaging data provide unique experimental input for atmospheric models that simulate aging and changes in particle-type populations, internal composition, and their associated optical and cloud forming properties. Here, we highlight applications of chemical imaging in selected recent studies, discuss their existing limitations, and forecast future research directions for this area.

54 ENVIRONMENTAL SCIENCES↗

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau-to-River (P2R) Model is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by Central Plateau Cleanup Company (CPCCo) at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the composite analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Meso-modeling of Carbon Fiber Composite for Crash Safety Analysis

In the conventional approach, the material properties for crash safety simulations are typically obtained from standard coupon tests, where the test results only provide single layer material properties used in crash simulations. However, the lay-up effects for the failure behaviors of the real structure were not considered in numerical simulations. Hence, there was discrepancy between the crash simulations and experimental tests. Consequently, an intermediate stage is required for accurate predictions. Some component tests are required to correlate the material models in the intermediate stage. In this paper, a Mazda Tube under high-impact velocity is chosen as an example for the crash safety analysis. The tube consists of 24 layers of uni-directional (UD) carbon fiber composite materials, in which 4 layers are perpendicular to, while the other layers are parallel to the impact direction. An LS-DYNA meso-model was constructed with orthotropic material models counting for the single-layer material behaviors. Between layers, a node-based tie-break contact was used for modeling the delamination of the composite material. Since fiber directions are not single-oriented, the lay-up effects could be an important effect. From the first numerical trial, premature material failure occurred due to the use of material parameters obtained directly from the coupon tests. Some parametric studies were conducted to identify the cause of the numerical instability. The finding is that the material failure strength used in the numerical model needs to be enlarged to stabilize the numerical model. Some hypothesis was made to provide the foundation for enlarging the failure strength and the corresponding experiments will be conducted to validate the hypothesis.

Lin, Shih-Po↗

Mechanical properties, failure mechanisms, and scaling laws of bicontinuous nanoporous metallic glasses

Molecular dynamics simulations are employed to study the mechanical properties of nanoporous Cu x Zr 1-x metallic glasses (MGs) with five different compositions, x = 0.28, 0.36, 0.50, 0.64, and 0.72, and porosity in the range 0.1 < Φ < 0.7. Results from tensile loading simulations indicate a strong dependence of Young's modulus, E, and Ultimate Tensile Strength (UTS) on porosity and composition. By increasing the porosity from Φ = 0.1 to Φ = 0.7, the topology of the nanoporous MG shifts from closed cell to open-cell bicontinuous. The change in nanoporous topology enables a brittle-to-ductile transition in deformation and failure mechanisms from a single critical shear band to necking and rupture of ligaments. Genetic Programming (GP) is employed to find scaling laws for E and UTS as a function of porosity and composition. Here, a comparison of the GP-derived scaling laws against existing relationships shows that the GP method is able to uncover expressions that can predict accurately both the values of E and UTS in the whole range of porosity and compositions considered.

36 MATERIALS SCIENCE↗

Chemical composition analysis and product consistency tests supporting refinement of the Nepheline model for the high aluminum Hanford Glass composition region

In this report, SRNL provides chemical analyses and Product Consistency Test (PCT) results for a series of simulated HLW glasses fabricated by Pacific Northwest National Laboratory (PNNL) as part of an ongoing nepheline crystallization study. The results of these analyses will be used to improve the ability to predict crystallization of nepheline as a function of composition and heat treatment for glasses formulated at high alumina concentrations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Chemical composition measurements of the low activity waste (LAW) EPA-Series glasses

In this report, the Savannah River National Laboratory provides chemical analysis results for a series of simulated low activity waste glasses provided by Pacific Northwest National Laboratory as part of an ongoing development task. The measured chemical composition data are reported and compared with the targeted values for each component for each glass. A detailed review showed no indications of errors in the preparation or measurement of the study glasses. All of the measured sums of oxides for the study glasses fell within the interval of 100.2 to 100.8 wt %, indicating recovery of all components. Comparisons of the targeted and measured chemical compositions showed that the measured values for the glasses met the targeted concentrations within 10% for those components present at more than 5 wt %.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Nonorthogonal Cascade Catalysis in Multicompartment Micelles

Multicompartment micelles (MCMs) containing acid and base sites in discrete domains are prepared from poly(norbornene)-based amphiphilic bottlebrush copolymers in aqueous media. Here, the acid and base sites are localized in different compartments of the micelle, enabling the nonorthogonal reaction sequence: deacetalization – Knoevenagel condensation – Michael addition of acetals to 2-amino chromene derivatives. Computational simulations using dissipative particle dynamics (DPD) elucidated the bottlebrush composition required to effectively site-isolate the nonorthogonal catalysts. This contribution presents MCMs as a new class of nanostructures for one-pot multistep nonorthogonal cascade catalysis, laying the groundwork for the isolation of three or more incompatible catalysts to synthesize value-added compounds in a single reaction vessel, in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Printed assemblies of microscale triple–junction inverted metamorphic GaInP/GaAs/InGaAs solar cells

Inverted metamorphic (IMM) multijunction solar cells represent a promising material platform for ultrahigh efficiency photovoltaic systems (UHPVs) with a clear pathway to beyond 50% efficiency. The conventional device processing of IMM solar cells, however, typically involves wafer bonding of a centimeter-scale die and destructive substrate removal, thereby imposing severe restrictions in achievable cell size, type of module substrate, spatial layout, as well as cost effectiveness. Here, we report material design and fabrication strategies for microscale triple-junction IMM (3J IMM) Ga 0.51 In 0.49 P/GaAs/In 0.26 Ga 0.74 As solar cells that can overcome these difficulties. Specialized schemes of delineation and undercut etching enable the defect-free release of microscale IMM solar cells and printed assemblies on a glass substrate in a manner that preserves the growth substrate, where efficiencies of 27.3% and 33.9% are demonstrated at simulated AM1.5D one- and 351 sun illumination, respectively. A composite carrier substrate where released IMM microcells are formed in fully functional, print-ready configurations allows high-throughput transfer printing of individual IMM microcells in a programmable spatial layout on versatile choices of module substrate, all desired for CPV applications.

14 SOLAR ENERGY↗