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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Novel Sparse Compositional Technique Reveals Microbial Perturbations

ABSTRACT The central aims of many host or environmental microbiome studies are to elucidate factors associated with microbial community compositions and to relate microbial features to outcomes. However, these aims are often complicated by difficulties stemming from high-dimensionality, non-normality, sparsity, and the compositional nature of microbiome data sets. A key tool in microbiome analysis is beta diversity, defined by the distances between microbial samples. Many different distance metrics have been proposed, all with varying discriminatory power on data with differing characteristics. Here, we propose a compositional beta diversity metric rooted in a centered log-ratio transformation and matrix completion called robust Aitchison PCA. We demonstrate the benefits of compositional transformations upstream of beta diversity calculations through simulations. Additionally, we demonstrate improved effect size, classification accuracy, and robustness to sequencing depth over the current methods on several decreased sample subsets of real microbiome data sets. Finally, we highlight the ability of this new beta diversity metric to retain the feature loadings linked to sample ordinations revealing salient intercommunity niche feature importance. IMPORTANCE By accounting for the sparse compositional nature of microbiome data sets, robust Aitchison PCA can yield high discriminatory power and salient feature ranking between microbial niches. The software to perform this analysis is available under an open-source license and can be obtained at https://github.com/biocore/DEICODE ; additionally, a QIIME 2 plugin is provided to perform this analysis at https://library.qiime2.org/plugins/q2-deicode .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

High radiation pressure on interstellar dust computed by light-scattering simulation on fluffy agglomerates of magnesium-silicate grains with metallic-iron inclusions

Recent space missions have provided information on the physical and chemical properties of interstellar grains such as the ratio β of radiation pressure to gravity acting on the grains in addition to the composition, structure, and size distribution of the grains. Numerical simulation on the trajectories of interstellar grains captured by Stardust and returned to Earth constrained the β ratio for the Stardust samples of interstellar origin. However, recent accurate calculations of radiation pressure cross-sections for model dust grains have given conflicting stories in the β ratio of interstellar grains. The β ratio for model dust grains of so-called “astronomical silicate” in the femto-kilogram range lies below unity, in conflict with β∼1 for the Stardust interstellar grains. Here, I tackle this conundrum by re-evaluating the β ratio of interstellar grains on the assumption that the grains are aggregated particles grown by coagulation and composed of amorphous MgSiO{sub 3} with the inclusion of metallic iron. My model is entirely consistent with the depletion and the correlation of major rock-forming elements in the Local Interstellar Cloud surrounding the Sun and the mineralogical identification of interstellar grains in the Stardust and Cassini missions. I find that my model dust particles fulfill the constraints on the β ratio derived from not only the Stardust mission but also the Ulysses and Cassini missions. My results suggest that iron is not incorporated into silicates but exists as metal, contrary to the majority of interstellar dust models available to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Effectiveness of Mitigation Measures in Reducing Future Primary Particulate Matter Emissions from On-Road Vehicle Exhaust

This work evaluates the effectiveness of on-road primary particulate matter emission reductions that can be achieved by long-term vehicle scrappage and retrofit measures on regional and global levels. Scenario analysis shows that scrappage can provide significant emission reductions as soon as the measures begin, whereas retrofit provides greater emission reductions in later years, when more advanced technologies become available in most regions. Reductions are compared with a baseline that already accounts for implementation of clean vehicle standards. The greatest global emission reductions from a scrappage program occur 5 to 10 years after its introduction and can reach as much as 70%. The greatest reductions with retrofit occur around 2030 and range from 16-31%. Monte Carlo simulations are used to evaluate how uncertainties in the composition of the vehicle fleet affect predicted reductions. Scrappage and retrofit reduce global emissions by 22-60% and 15-31%, respectively, within 95% confidence intervals, under a midrange scenario in the year 2030. The simulations provide guidance about which strategies are most effective for specific regions. Retrofit is preferable for high-income regions. For regions where early emission standards are in place, scrappage is suggested, followed by retrofit after more advanced emission standards are introduced. The early implementation of advanced emission standards is recommended for Western and Eastern Africa

Yan, Fang↗
DOE OSTI2020

Size dependence in chord characteristics from simulated and observed continental shallow cumulus

Abstract. In this study we compare long-term Doppler and Raman lidar observations against a full month of large eddy simulations of continental shallow cumulus clouds. The goal is to evaluate if the simulations can reproduce the mean observed vertical velocity and moisture structure of cumulus clouds and their associated subcloud circulations, as well as to establish if these properties depend on the size of the cloud. We propose methods to compare continuous chords of cloud detected from Doppler and Raman lidars with equivalent chords derived from 1D and 3D model output. While the individual chords are highly variable, composites of thousands of observed and millions of simulated chords contain a clear signal. We find that the simulations underestimate cloud size and fraction but successfully reproduce the observed structure of vertical velocity and moisture perturbations. There is a clear scaling of vertical velocity and moisture anomalies below the chords with chord size, but the moisture anomalies are only 1 %–2 % higher than the horizontal mean values. The differences between the observations and simulations are smaller than the difference in sampling the modeled chords in time or space. The shape of the vertical velocity and moisture anomalies from cloud chords sampled spatially from 3D model snapshots is almost perfectly symmetric. In contrast, the chords sampled temporally from the lidar observations and 1D model output have a marked asymmetry, with stronger updrafts and higher moisture anomalies occurring earlier on.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2021

IACMI Project 4.2: Thermoplastic Composite Development for Wind Turbine Blades

(Section 5.1) Composites made from Arkema’s Elium® thermoplastic resin and Johns Manville fiberglass were researched during this project for applications in wind blade manufacturing. A techno-economic model was developed to model this wind blade manufacturing process using these materials in place of traditional composites made with thermoset resin. This model was based on manufacturing a 61.5-meter wind blade, which showed a 4.7% reduction in wind blade cost as compared traditional thermoset materials. These cost savings were not from the thermoplastic material costing less than traditional thermoset materials, but rather from decreased capital costs, faster cycle times and reduced energy requirements and labor costs. (Section 5.2) An infusion and curing model was developed for thermoplastic composite wind blades using PAM-RTM. The primary goal was to demonstrate the infusion simulation for the Elium® resin system on a 13-meter wind blade. Additionally, the exotherm temperature was predicted and compared to measurements, which showed model results within 10% of actual measurements. (Section 5.3) Composite laminate panels and composite sandwich panels with a balsa core were produced; specimens were cut and characterized. Similar composite specimens were made with Elium® thermoplastic resin and Hexion thermoset epoxy (RIMR135/RIMH1366) to enable comparisons between these resin systems. The static test methods included: tensile, compression, in-plane shear, interlaminar shear, flexural, sandwich core shear flexure, and single cantilever beam tests for sandwich beams. Fatigue testing at room temperature was completed to composite laminate panels at a stress ratio of R=0.1 and R=10. In addition, fatigue testing to laminate panels was completed at -30°C, and at room temperature after conditioning specimens at 70°C and 90% relative humidity. Overall, mechanical test results from Elium® composites are similar to epoxy composites. (Section 5.4) Elium composite panels were produced with intentional defects such as voids and nonwetting of fibers to begin to understand performance sensitivity to defects. A thermal digital image correlation (TDIC) method provides high spatial resolution strain field at elevated temperatures and can be used to identify defective regions within composite panels. Flexural modulus differences of 21% were seen between defect and non-defect panels. Other Elium® composite panels were forced to be defective by boiling the resin after infusion, which created voids throughout the composite laminate. X-ray computed tomography scanning was used to view the internal structure of the defect panels. Defect panels had a significant reduction in fatigue life as compared to baseline panels produced without intentional defects. (Section 5.5) Lap shear specimens were fabricated to compare the lap shear strength of an off-the-shelf adhesive (Plexus MA590) and two new adhesives developed by Arkema (Bostik SAF30 90 and Bostik SAF30 120). ISO standard 4587:2003 was used to standardize the testing method and sample fabrication. Lap shear specimens were made at 1mm, 3mm, and 10mm thicknesses. The Bostik adhesive lap shear test results were similar to Plexus for all thicknesses. (Section 5.6) Fiber-reinforced polymer (FRP) composites are typically used in high-performance applications (e.g., aerospace), and their expansion into high-volume industries (e.g. consumer automotive and wind turbine blade manufacturer or similar) is hindered by their cost and a lack of efficient manufacturing techniques. Monitoring the curing process of these composites during manufacturing can improve the efficiency of the process, and therefore reduce the manufacturing cost. Cure monitoring techniques were developed that use probabilistic estimation methods and surface temperature measurements made using infrared cameras. These techniques enable real-time monitoring of the infusion process to locate manufacturing flaws, and they can, potentially, estimate residual stresses in the part. Their commercialization will help facilitate expansion of FRP composites in high-volume industries. (Section 5.7) A 13-meter composite wind blade was produced with Elium® resin and Johns Manville fiberglass; this blade was made with VARTM processing similar to how megawatt-scale wind blades are currently manufactured, but no post-mold heating was used for this thermoplastic composite blade. The wind blade underwent full-scale validation for static loading (4-different load orientations) and flapwise fatigue loading to simulate 20-years of operational loads. The thermoplastic composite wind blade withstood the loading without any noted issues and performed similar to results from a previous full-scale validation to an equivalent epoxy composite wind blade produced with the same blade molds. (Section 5.8) A study was conducted to determine the feasibility of recycling composite wind turbine blade components fabricated with glass fiber reinforced Elium® thermoplastic resin. Dissolution, which is a process unique to thermoplastic matrices, allows recovery of both the polymer matrix and full-length glass fibers, while maintaining their stiffness and strength throughout the recovery process. The economics of recycling is favorable if 50% of the glass fiber is recovered and resold for a process of $\$$ 0.28/kg, and 90% of the resin is recovered and resold at a price of $\$$ 2.50/kg.(Section 10) Recommendations are outlined for commercializing thermoplastic resin for composite wind blade production, in addition to recommended areas for future research.

17 WIND ENERGY↗
DOE OSTI2014

A thermodynamic model for C-(N-)A-S-H gel: CNASH{sub s}s. Derivation and validation

The main reaction product in Ca-rich alkali-activated cements and hybrid Portland cement (PC)-based materials is a calcium (alkali) aluminosilicate hydrate (C-(N-)A-S-H) gel. Thermodynamic models without explicit definitions of structurally-incorporated Al species have been used in numerous past studies to describe this gel, but offer limited ability to simulate the chemistry of blended PC materials and alkali-activated cements. Here, a thermodynamic model for C-(N-)A-S-H gel is derived and parameterised to describe solubility data for the CaO–(Na{sub 2}O,Al{sub 2}O{sub 3})–SiO{sub 2}–H{sub 2}O systems and alkali-activated slag (AAS) cements, and chemical composition data for C-A-S-H gels. Simulated C-(N-)A-S-H gel densities and molar volumes are consistent with the corresponding values reported for AAS cements, meaning that the model can be used to describe chemical shrinkage in these materials. Therefore, this model can provide insight into the chemistry of AAS cements at advanced ages, which is important for understanding the long-term durability of these materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Simulations of Seismic Wave Propagation on Mars

In this paper, we present global and regional synthetic seismograms computed for 1D and 3D Mars models based on the spectral-element method. For global simulations, we implemented a radially-symmetric Mars model with a 110 km thick crust. For this 1D model, we successfully benchmarked the 3D seismic wave propagation solver SPECFEM3D_GLOBE against the 2D axisymmetric wave propagation solver AxiSEM at periods down to 10 s. We also present higher-resolution body-wave simulations with AxiSEM down to 1 s in a model with a more complex 1D crust, revealing wave propagation effects that would have been difficult to interpret based on ray theory. For 3D global simulations based on SPECFEM3D_GLOBE, we superimposed 3D crustal thickness variations capturing the distinct crustal dichotomy between Mars’ northern and southern hemispheres, as well as topography, ellipticity, gravity, and rotation. The global simulations clearly indicate that the 3D crust speeds up body waves compared to the reference 1D model, whereas it significantly changes surface waveforms and their dispersive character depending on its thickness. We also perform regional simulations with the solver SES3D based on 3D crustal models derived from surface composition, thereby addressing the effects of various distinct crustal features down to 2 s. The regional simulations confirm the strong effects of crustal variations on waveforms. Finally, we conclude that the numerical tools are ready for examining more scenarios, including various other seismic models and sources.

58 GEOSCIENCES↗
DOE OSTI2023

An ab initio molecular dynamics study of varied compositions of the $\text{LiF-NaF-KF}$ molten salt

With increasing interest in molten salt reactors, there becomes a demand for investigations into thermophysical properties of salt systems. The LiF-NaF-KF (FLiNaK) salt system is a primary candidate for use in these reactors. However, the thermophysical properties of compositions outside the eutectic composition are still largely unknown. In this article, properties of ten unique compositions, including four ternary compositions, are investigated using ab initio molecular dynamics simulations across five temperatures between 900 K and 1300 K. The properties of interest are the density, thermal expansion, bulk modulus, compressibility, heat capacity, and enthalpy of mixing. In general, the results were found to be in good agreement with other literature and experimental results. The density and heat capacity had a tendency to be slightly underpredicted. No conclusions could be drawn about the bulk modulus and compressibility in terms of compositional dependence. The thermal expansion had a negative trend with respect to the LiF concentration and no trends were observed for the NaF or KF concentration. The enthalpy of mixing shows minima for the ternary compositions, with the near-equiatomic composition exhibiting the lowest values. Here this work shows the potential for compositional tailoring in the FLiNaK system to optimize thermophysical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multiscale Modeling of Microstructural Evolution in Fused-Coating Additive Manufacturing

A quantitative phase field model with adaptive-computational-domain-expansion (ACDE) method was developed to simulate the microstructural evolution in fused-coating additive manufacturing process. A heat transfer and fluid flow model was established to investigate the solidification behavior of the molten pool in a single deposited layer. The solidification parameters such as temperature gradient and solidification rate were calculated from the macroscopic model and then fed into the phase field model to get the microstructure at different positions within a single deposited layer. Columnar dendrites were formed at the bottom of the deposited layer, and equiaxed dendrites appeared in the upper part of the layer. The ACDE method is used to improve the simulation efficiency. The metallographic experiments and energy-dispersive spectrometer analysis for composition were conducted on the fused-coating specimen, and the results were in good accordance with the simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

Lorentz invariance violation (LIV) is often described by dispersion relations of the form E i 2 = m i 2 +p i 2 +δ i,n E 2+n with delta different based on particle type i , with energy E, momentum p and rest mass m. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients δ i,n tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond 10 19 eV, we obtain δ γ,0 > -10 -21 , δ γ,1 > -10 -40 eV -1 and δ γ,2 > -10 -58 eV -2 in the case of a subdominant proton component up to 10 20 eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as δ had,0 < 10 -19 , δ had,1 < 10 -38 eV -1 and δ had,2 < 10 -57 eV -2 at 5σ CL.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-T-455: Characterization of 3D Printed Materials for Proton Beam Therapy

Purpose: The widespread availability of low cost 3D printing technologies provides an alternative fabrication method for customized proton range modifying accessories such as compensators and boluses. However the material properties of the printed object are dependent on the printing technology used. In order to facilitate the application of 3D printing in proton therapy, this study investigated the stopping power of several printed materials using both proton pencil beam measurements and Monte Carlo simulations. Methods: Five 3–4 cm cubes fabricated using three 3D printing technologies (selective laser sintering, fused-deposition modeling and stereolithography) from five printers were investigated. The cubes were scanned on a CT scanner and the depth dose curves for a mono-energetic pencil beam passing through the material were measured using a large parallel plate ion chamber in a water tank. Each cube was measured from two directions (perpendicular and parallel to printing plane) to evaluate the effects of the anisotropic material layout. The results were compared with GEANT4 Monte Carlo simulation using the manufacturer specified material density and chemical composition data. Results: Compared with water, the differences from the range pull back by the printed blocks varied and corresponded well with the material CT Hounsfield unit. The measurement results were in agreement with Monte Carlo simulation. However, depending on the technology, inhomogeneity existed in the printed cubes evidenced from CT images. The effect of such inhomogeneity on the proton beam is to be investigated. Conclusion: Printed blocks by three different 3D printing technologies were characterized for proton beam with measurements and Monte Carlo simulation. The effects of the printing technologies in proton range and stopping power were studied. The derived results can be applied when specific devices are used in proton radiotherapy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Temperature and composition dependent screw dislocation mobility in austenitic stainless steels from large-scale molecular dynamics

Abstract Extensive molecular dynamics simulations are performed to determine screw dislocation mobility in austenitic Fe 0.7 Ni x Cr 0.3-x stainless steels as a function of temperature ranging from 100 to 1300 K, resolved shear stress from 30 to 140 MPa, and Ni composition from 0.0 to 30.0 at%. These mobility data are fitted to a linear mobility law with a nonzero stress offset, referred to as the threshold stress. We find that both the linear drag coefficient and the threshold stress increase with Ni composition. The drag coefficient increases with temperature, whereas the threshold stress decreases with temperature. Based on these calculations, we determine fitting functions for the linear solute drag coefficient as a function of temperature and composition. The mobility laws determined in this study may serve to inform dislocation dynamics simulations pertinent to dislocation network evolution at elevated temperatures for a wide composition range of austenitic stainless steels.

Chu, Kevin (ORCID:000000031929710X)↗
DOE OSTI2019

Simulations of particle interactions in a high-current RFQ

High-intensity beams of composite ions are susceptible to particle interactions that do not affect single ions, such as stripping of electrons and charge exchange. Beam dynamics simulation techniques used in accelerator design do not often include these particle interactions, and so the transmission of the real beam can be significantly different from the design simulations. We have modified the IMPACT-T code to include particle interactions as part of the beam dynamics simulations, including space-charge effects and beam losses. We have carried out validation simulations of interactions in a molecular hydrogen ion beam through an existing high-intensity deuteron RFQ design, as beam commissioning is often carried out using molecular hydrogen beams, to avoid deuteron-deuteron interactions in the accelerator leading to neutron activation. The new code framework successfully models the interactions within the beam, tracks the different source and product ions through the accelerator, and calculates the collective effects such as space charge from all particles together. Finally, we also discuss other simulation codes that could be modified to include similar interaction effects, and possible applications to other types of accelerators.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Changes in Tropical Cyclones Undergoing Extratropical Transition in a Warming Climate: Quasi-Idealized Numerical Experiments of North Atlantic Landfalling Events

The current study extends earlier work that demonstrated future extratropical transition (ET) events will feature greater intensity and heavier precipitation to specifically consider potential changes in the impacts of landfalling ET events in a warming climate. A quasi-idealized modeling framework allows comparison of highly similar present-day and future event simulations; the model initial conditions are based on observational composites, increasing representativeness of the results. The future composite ET event features substantially more impactful weather conditions in coastal areas, with heavier precipitation and greater storm intensity. Specifically, a Category 2 present-day storm attained Category 4 Saffir-Simpson intensity in the future simulation and maintained greater intensity throughout the entire life cycle, although the storm undergoes less reintensification during the post-ET process, a result of reduced baroclinic conversion. These findings suggest increased potential for coastal hazards due to stronger tropical cyclone winds and heavier rainfall, leading to more severe coastal flooding and storm surge.

54 ENVIRONMENTAL SCIENCES↗