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 199 records · Page 11

Surface Segregation Studies in Ternary Noble Metal Alloys: Comparing DFT and Machine Learning with Experimental Data

Abstract Surface segregation, whereby the surface composition of an alloy differs systematically from the bulk, has historically been hard to study, because it requires experimental and modeling methods that span alloy composition space. In this work, we study surface segregation in catalytically relevant noble and platinum‐group metal alloys with a focus on three ternary systems: AgAuCu, AuCuPd, and CuPdPt. We develop a data set of 2478 fcc slabs with those compositions including all three low‐index crystallographic orientations relaxed with Density Functional Theory using the PBEsol functional with D3 dispersion corrections. We fine‐tune a machine learning model on this data and use the model in a series of 1800 Monte Carlo simulations spanning ternary composition space for each surface orientation and ternary chemical system. The results of these simulations are validated against prior experimental surface segregation data collected using composition spread alloy films for AgAuCu and AuCuPd. Our findings reveal that simulations conducted using the (110) orientation most closely match experimentally observed surface segregation trends, and while predicted trends qualitatively match observation, biases in the PBEsol functional limit numeric accuracy. This study advances understanding of surface segregation and the utility of computational studies and highlights the need for further improvements in simulation accuracy.

Chemistry↗

Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites

This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021

Materials Science↗

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗

Optimization of Type 4 composite pressure vessels using genetic algorithms and simulated annealing

We propose and demonstrate two methodologies for weight minimization of Type 4 (plastic lined, fiber wrapped) compressed hydrogen pressure vessels: genetic algorithms and simulated annealing. We consider 70 MPa vessels with a safety factor of 2.25, and analyze the vessels with classical laminate theory. We propose an objective function based on Tsai-Wu criterion, composite thickness, a safety factor, and a penalization factor. Here, the optimum results are analyzed and compared by a high resolution finite element model. Computer simulations show that the proposed methodology produces more efficient designs by reducing the weight by up to 9.8% and 11.2% when compared to previously published vessel optimization research.

08 HYDROGEN↗

Hanford Composite Analysis Special Studies FY2025 - Groundwater Flow Simulation for use in a Comparison of Simulated Concentrations in the Saturated Zone Estimated by the Plateau-to-River Model Versions 8.3 and 9.1

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This Environmental Calculation File (ECF) details the application of the Plateau-to-River (P2R) Model (CP-57037, Rev. 3, Model Package Report for the P2R Model: Version 9.1) to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support a special analysis associated with the Hanford Site CA. The simulated flow field will support the simulation of fate and transport of contaminants for comparison to the results obtained as part of the CA to evaluate impact of the updated version of the P2R Model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effects of 2D filler on rheology of additive manufacturing polymers: Simulation and experiment on polyetherketoneketone-mica composites

Polymeric composites have attracted increasing interest in fundamental studies and industrial applications due to their often-superior characteristics over conventional pure polymers. In this work, we studied the relationship between the rheological properties, critical in processing methods such as additive manufacturing, and the structure of a polymer-2D filler system by experiment and simulation methods. Considering the high performance polymer polyetherketoneketone (PEKK) with 2D filler (mica) composite as an example, the effects of filler size, filler content, surface chemistry and shear rate were systematically investigated. Here, it is demonstrated that the interaction between the polymer and 2D filler brings two different consequences on tuning the viscosity of polymer-filler composite. The stronger interaction increases the shear stress of the polymer-2D filler composite system, leading to higher viscosity. On the other hand, the stronger interaction can promote more absorption of polymer chains on the surface, which makes the local polymer network disentangled and decreases the continuous viscosity from the polymer network. This work provides a fundamental understanding on the relationship between macroscopic rheological properties and the structural properties and energetics at molecular level of polymer-2D filler composite.

36 MATERIALS SCIENCE↗

Beam Current Effect on Microstructure and Properties of Electron-Beam-Melted Ti-6Al-4V Alloy

In this study, a noncommercial 3D printing machine was used to fabricate Ti-6Al-4V alloy by electron-beam melting (EBM). The influence of beam current on the microstructure, phase composition and mechanical properties of electron-beam-melted Ti-6Al-4V alloy was investigated. Numerical simulation is implemented to evaluate thermal fields during electron-beam melting of Ti-6Al-4V powder. The decrease in beam current from 3.5 to 2.5 mA leads to refinement of microstructure: The average width of α plates decreases from 10 down to 6 μm. The formation of finer microstructure is attributed to higher cooling rate at lower beam current confirmed by simulation. The phase composition of EBM Ti-6Al-4V indirectly depends on the beam current. High content of β phase (7%) was achieved at the beam current of 3 mA. The produced Ti-6Al-4V samples are characterized by high microhardness (470-520 HV).

36 MATERIALS SCIENCE↗

Predictive Contaminant Transport Simulation with the P2R Model for the Composite Analysis Base Case

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1 Chg. 1, Radioactive Waste Management, a composite analysis (CA) must be completed for the Central Plateau at the Hanford Site. The updated Hanford Site CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within and downgradient of the Central Plateau. Numerical simulations using the MT3D-MST computer code were conducted to evaluate the fate and transport of existing distributions of 16 radionuclides as identified in CP-62184, Hanford Site Composite Analysis: Radionuclide Selection for Groundwater Pathway Evaluation: tritium (H-3), carbon-14 (C-14), chloride- 36 (Cl-36), strontium-90 (Sr-90), technetium-99 (Tc-99), rhenium-187 (Re-187), radium-226 (Ra-226), thorium-230 (Th-230), uranium-232 (U-232), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), neptunium-237 (Np-237), and uranium-238 (U- 238). This document details the application of the P2R Model version 8.3 to predict the fate and transport of contaminants in groundwater on the Central Plateau for the 10,052 year simulation to support the CA. The simulation of fate and transport of contaminants reported in this environmental calculation file (ECF) will support dose predictions as part of the updated Hanford Site CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Metal-ceramic materials. Study and prediction of effective mechanical properties

Mechanical behavior of stochastic metal-ceramic composite materials was numerically simulated on mesoscopic scale level. Deformation of mesoscopic volumes of composites, whose structure consists of a metal matrix and randomly distributed ceramic inclusions, was numerically simulated. The results of the numerical simulation were used for evaluation of the effective elastic and strength properties of metal-ceramic materials with different parameters of the structure. The values of the effective mechanical properties of investigated materials were obtained, and the character of the dependence of the effective elastic and strength properties on the structure parameters of composites was determined.

36 MATERIALS SCIENCE↗

Black Hole to Photosphere: 3D GRMHD Simulations of Collapsars Reveal Wobbling and Hybrid Composition Jets

Long-duration γ-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth via black hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over six orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important, and launches bipolar jets. The jets reach the photosphere at ~10 12 cm with an opening angle θ j ~ 6° and a Lorentz factor Γ j ≲ 30, unbinding ≳90% of the star. We find that (i) the disk–jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle θ t ~ 12°, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection θ j + θ t suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs are rarer than currently thought and emerge in only ~0.1% of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where ~10% of the gas maintains magnetization σ ≳ 0.1. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.

79 ASTRONOMY AND ASTROPHYSICS↗

A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity (κ C ) and area compressibility (K A ). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate K A values, current approaches suffer significant technical limitations. Here, in this study, we present a novel, to our knowledge, computational framework that allows for a direct estimation of K A values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives K A from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model’s predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.

59 BASIC BIOLOGICAL SCIENCES↗

Assessing modifications to the Abdul-Razzak and Ghan aerosol activation parameterization (version ARG2000) to improve simulated aerosol–cloud radiative effects in the UK Met Office Unified Model (UM version 13.0)

The representation of aerosol activation is a key source of uncertainty in global composition-climate model simulations of aerosol–cloud interactions. The Abdul-Razzak and Ghan (ARG) activation parameterization is used in several global and regional models that employ modal aerosol microphysics schemes. In this study, we investigate the ability of the ARG parameterization to reproduce simulations with a cloud parcel model and find its performance is sensitive to the geometric standard deviations (widths) of the lognormal aerosol modes. We recommend adjustments to three constant parameters in the ARG equations, which improve the performance of the parameterization for small mode widths and its ability to simulate activation in polluted conditions. For the accumulation mode width of 1.4 used in the Met Office Unified Model (UM), the modifications decrease the mean bias in the activated fraction of aerosols compared to a cloud parcel model from −6.6 % to +1.2 %. We implemented the improvements in the UM and compared simulated global cloud droplet concentrations with satellite observations. The simulated cloud radiative effect changes by −1.43 W m −2 (6 %) and aerosol indirect radiative forcing over the industrial period changes by −0.10 W m −2 (10 %).

Ghosh, Pratapaditya [Carnegie Mellon University, P↗

American Society of Composites, 32nd Technical Conference

Resin transfer molding (RTM) has become increasingly popular for the manufacturing of composite parts. To enable high volume manufacturing and obtain good quality parts at an acceptable cost to automotive industry, accurate process simulation tools are necessary to optimize the process conditions. Towards that goal, General Motors and the ESI-group are involved in developing a state of the art process simulation tool for composite manufacturing in a project supported by the Department of Energy. This paper describes the modeling of various stages in resin transfer molding such as resin injection, resin curing, and part distortion. An instrumented RTM system located at the General Motors Research and Development center was used to perform flat plaque molding experiments. The experimental measurements of fill time, in-mold pressure versus time, cure variation with time, and part deformation were compared with the model predictions and very good correlations were observed.

42 ENGINEERING↗

Predicting the mineral composition of dust aerosols – Part 2: Model evaluation and identification of key processes with observations

A global compilation of nearly sixty measurement studies is used to evaluate two methods of simulating the mineral composition of dust aerosols in an Earth system model. Both methods are based upon a Mean Mineralogical Table (MMT) that relates the soil mineral fractions to a global atlas of arid soil type. The Soil Mineral Fraction (SMF) method assumes that the aerosol mineral fractions match the fractions of the soil. The MMT is based upon soil measurements after wet sieving, a process that destroys aggregates of soil particles that would have been emitted from the original, undisturbed soil. The second method approximately reconstructs the emitted aggregates. This model is referred to as the Aerosol Mineral Fraction (AMF) method because the mineral fractions of the aerosols differ from those of the wet-sieved parent soil, partly due to reaggregation. The AMF method remedies some of the deficiencies of the SMF method in comparison to observations. Only the AMF method exhibits phyllosilicate mass at silt sizes, where they are abundant according to observations. In addition, the AMF quartz fraction of silt particles is in better agreement with measured values, in contrast to the overestimated SMF fraction. Measurements at distinct clay and silt particle sizes are shown to be more useful for evaluation of the models, in contrast to the sum over all particles sizes that is susceptible to compensating errors, as illustrated by the SMF experiment. Model errors suggest that allocation of the emitted silt fraction of each mineral into the corresponding transported size categories is an important remaining source of uncertainty. Evaluation of both models and the MMT is hindered by the limited number of size-resolved measurements of mineral content that sparsely sample aerosols from the major dust sources. In conclusion, the importance of climate processes dependent upon aerosol mineral composition shows the need for global and routine mineral measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Novel CAE Method for Compression Molding Simulation of Carbon Fiber-Reinforced Thermoplastic Composite Sheet Materials

Its high-specific strength and stiffness with lower cost make discontinuous fiber-reinforced thermoplastic (FRT) materials an ideal choice for lightweight applications in the automotive industry. Compression molding is one of the preferred manufacturing processes for such materials as it offers the opportunity to maintain a longer fiber length and higher volume production. In the past, we have demonstrated that compression molding of FRT in bulk form can be simulated by treating melt flow as a continuum using the conservation of mass and momentum equations. However, the compression molding of such materials in sheet form using a similar approach does not work well. The assumption of melt flow as a continuum does not hold for such deformation processes. To address this challenge, we have developed a novel simulation approach. First, the draping of the sheet was simulated as a structural deformation using the explicit finite element approach. Next, the draped shape was compressed using fluid mechanics equations. The proposed method was verified by building a physical part and comparing the predicted fiber orientation and warpage measurements performed on the physical parts. The developed method and tools are expected to help in expediting the development of FRT parts, which will help achieve lightweight targets in the automotive industry.

36 MATERIALS SCIENCE↗

Roadmap for thermal property measurements of Molten Salt Reactor systems

Thermal properties of molten salts, both with and without fissile and fertile material, are fundamental to modeling and simulation that supports design, construction, operation, and accident progression evaluation of molten salt reactors. This roadmap provides guidance to the research community on the types and amounts of required data for base, or carrier, salt systems, recommends simulated burnup fuel salt compositions to evaluate, and provides a technical foundation for a focused salt property evaluation matrix. The property values generated from both computation and experimental methods are used to develop models that constitute the Molten Salt Thermal Properties Database (MSTDB). While the building blocks of MSTDB are fundamental in nature, its data facilitates constructing engineering models that allow stakeholders to investigate the behavior of specific compositions of interest and supports broader modeling and simulation of MSRs through coupling to multi-physics, multi-scale mass accountancy tools.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1 Chg. 1, Radioactive Waste Management, a Composite Analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This document details the application of the P2R Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the CA. The simulated flow field will support the simulation of fate and transport of contaminants for use in predicting dose as part of the CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case

In order to meet the requirements of U.S. Department of Energy (DOE) in DOE O 435.1 Chg.1, Radioactive Waste Management, a Composite Analysis (CA) must be completed for the Central Plateau at the Hanford Site. The CA requires estimates of fate and transport of radionuclides in the groundwater from multiple sources within the modeling domain. This document details the application of the P2R Model version 8.3 to predict the flow of groundwater on the Central Plateau for the 10,000-year simulation to support the CA. The simulated flow field will support the simulation of fate and transport of contaminants for use in predicting dose as part of the CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗