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.

SourceResearch connectionFrom yearThrough year

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

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Computational prediction of dielectric breakdown strength of a transformer paper in oil with uncertainty quantification

The determination of the dielectric breakdown strengths of microstructurally heterogeneous materials has been a primarily experimental endeavor. We report the development of a microstructure-level model for computationally predicting the breakdown strength and analyzing the interactions between electromagnetic pulses (EMP) and the constituents in a composite of cellulose-based paper and mineral oil found in electrical transformers. The model allows explicit simulation of the material breakdown process by tracking the transition of dielectric constituents from non-conductive to conductive states. The focus is on the electric fields induced in the materials and the overall conditions for dielectric breakdown (defined as the onset of avalanche) caused by the electric field induced in the composite. Responses to three distinct pulse shapes, i.e., Steep Front (SF), Lightning (L), and AC with spectra spanning 60–9 × 105 Hz are considered. It is found that the breakdown strength of the material is significantly affected by microstructure heterogeneities, the spatial variations of the constituent properties, and the pulse shapes. A probabilistic characterization of the breakdown strength is computationally obtained and compared with experimental measurements. Although one particular material is analyzed, the model and approach are applicable to other heterogeneous materials as well.

breakdowns↗

DOE OSTI2020

Temperature drives local contributions to beta diversity in mountain streams: stochastic and deterministic processes

Community variation (i.e. beta diversity) in geographical coenoclines is of substantial interest in ecology, yet the elevational patterns in beta diversity are still understudied, and no general conclusions have been reached. Here, we explored local contributions to beta diversity (LCBD) of stream bacteria, diatoms and macroinvertebrates on mountainsides. Across all taxa, we revealed a general U-shaped LCBD-elevation relationships (LCBDer), suggesting higher uniqueness of community composition at the two ends of the elevational gradients. This pattern was confirmed by coenocline simulations based on observed communities, and supported by the dominance of species replacement in total beta diversity. Across mountains and taxonomic groups, temperature was the main environmental factor underlying elevational patterns in LCBD and its metric relatives, which were also mediated by local environmental variables. The observed magnitudes of U-shaped LCBDer across mountains and taxonomic groups were specifically related to the heterogeneity of stream morphology, and the mean annual temperature on each mountain. U-shaped LCBDer were strengthened by the occurrence of generalist species with large elevational ranges, while rare species had a lesser effect on the patterns. Raup-Crick beta diversity indicated that both stochastic and deterministic processes drove the community assembly, and the latter one was positively related to the magnitudes of LCBDer. Our synthesis across mountains and taxonomic groups clearly shows that there are consistent elevational patterns in LCBD among taxonomic groups, and that these patterns are explained by similar ecological mechanisms, producing a more complete picture in understanding and bridging the spatial variation in biodiversity under changing climate.

Wang, Jianjun↗

DOE OSTI2013

Biologically inspired growth of hydroxyapatite crystals on bio-organics-defined scaffolds

Graphical abstract: Petal-like crystals were observed to form on the surface of the BG/COL/ChS scaffolds. Highlights: ► Porous scaffolds were prepared using bioglass, collagen and chondroitin sulfate. ► Highly oriented HA crystals were grown on scaffolds using simulated body fluids ► The microstructure and orientation of HA were explained by molecular configuration. - Abstract: Several bio-organics-defined composite scaffolds were prepared using 58s-bioglass (BG), collagen (Col) and chondroitin sulfate (ChS). These scaffolds possess highly porous structure. X-ray diffraction of these scaffolds strongly indicated that hydroxyapatite (HA) crystals formed on their surfaces in simulated body fluids within 3 d, and similar formation process of crystals could be obtained on BG/Col and BG/Col/ChS scaffolds. The morphology and structure of the crystals were further examined by scanning electron microscopy. The results obtained indicate that an apatite with petal-like structure similar to that found on BG/Col scaffolds can be produced on BG/Col/ChS scaffolds through biomimetic synthesis, while that on BG/ChS scaffolds took place differently. The differences could be explained by self-assembly processes and the different macromolecular configurations of the Col and ChS fibrils which self-assemble spontaneously into their fibers. On the other hand, the bio-organics-defined composites have good cell biocompability. The results may be applicable to develop tailored biomaterials for peculiar bone substitute.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Modeling of high composition AlGaN channel high electron mobility transistors with large threshold voltage

We report on the potential of high electron mobility transistors (HEMTs) consisting of high composition AlGaN channel and barrier layers for power switching applications. Detailed two-dimensional (2D) simulations show that threshold voltages in excess of 3 V can be achieved through the use of AlGaN channel layers. We also calculate the 2D electron gas mobility in AlGaN channel HEMTs and evaluate their power figures of merit as a function of device operating temperature and Al mole fraction in the channel. Our models show that power switching transistors with AlGaN channels would have comparable on-resistance to GaN-channel based transistors for the same operation voltage. The modeling in this paper shows the potential of high composition AlGaN as a channel material for future high threshold enhancement mode transistors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Nonadiabatic Flamelet Formulation for Predicting Wall Heat Transfer in Rocket Engines

A flamelet-based combustion model is introduced for the prediction of wall heat transfer in rocket engines and confined combustion systems. To account for convective heat loss due to the interaction of the flame with the wall, a permeable thermal boundary condition is introduced in the counterflow diffusion flame configuration. The solution of the resulting nonadiabatic flame structure forms a three-dimensional manifold, which is parameterized in terms of mixture fraction, progress variable, and temperature. The performance of the model is first assessed through a direct numerical simulation analysis of an H 2 /O 2 diffusion flame that is stabilized at an inert isothermal wall. The developed nonadiabatic flamelet model is shown to accurately predict the temperature, chemical composition, and wall heat transfer. Combined with a presumed probability density function closure, the model is then applied to large-eddy simulation of a single-injector rocket combustor to examine effects of heat transfer on the turbulent flame structure in rocket engines.

42 ENGINEERING↗

DOE OSTI2022

Composition Measurements of the LAW HPVR Glasses

This report provides the results from the chemical analyses of the glass compositions of the Low-Activity Waste High PCT and VHT Response study glasses, a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory. These data will be used in the development, validation, and implementation of enhanced property/composition models for waste glass vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the quenched glasses to allow for comparisons with targeted compositions. The relative differences between the targeted and measured concentrations of Cl - , K 2 O, Na 2 O, and ZrO 2 for several of the glasses were greater than 10%. These results can be used in further characterization of this series of glasses, including the normalization of Product Consistency Test results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Simulation-driven design optimization of reaction injection molding (RIM) process for polydicyclopentadiene (pDCPD): Minimizing cycle time, defects, and warpage

Replacing metal components in trucks, trailers, and buses with lightweight polymer composites is challenging due to high temperatures and complex manufacturing. The Reaction Injection Molding (RIM) process using Dicyclopentadiene (DCPD) resin offers a solution by producing robust parts with excellent stiffness, impact strength, and resistance properties. Simulations are essential for optimizing this process, predicting defects, and improving quality. However, most commercial software is tailored for thermoplastics, requiring thermoset users to generate their own datasets. In this study, a material data card for DCPD was developed to perform RIM simulations. Design of Experiments (DOE) was used to identify key factors affecting filling, curing, and warpage, aiming to minimize cycle time and defects. The simulations explored varying injection gate parameters (size, location, number) and process conditions (mold/resin temperature, injection/curing pressure). Results showed that gate design significantly impacts filling behavior and defects. A single central gate provided balanced flow with fewer defects, while two corner gates led to more defects. Additionally, lower injection pressure increased filling time, while higher mold temperature accelerated curing but led to more warpage. In conclusion, this optimization framework aims to enhance DCPD part performance and promote sustainable manufacturing by reducing waste and energy consumption.

42 ENGINEERING↗

DOE OSTI2015

Virus Evolution Serial

This is the first code, designed to run on a desktop, which models the intracellular replication and the cell-to-cell infection and demonstrates virus evolution at the molecular level. This code simulates the infection of a population of "idealized biological cells" (represented as objects that do not divide or have metabolism) with "virus" (represented by its genetic sequence), the replication and simultaneous mutation of the virus which leads to evolution of the population of genetically diverse viruses. The code is built to simulate single-stranded RNA viruses. The input for the code is 1. the number of biological cells in the culture, 2. the initial composition of the virus population, 3. the reference genome of the RNA virus, 4. the coordinates of the genome regions and their significance and, 5. parameters determining the dynamics of virus replication, such as the mutation rate. The simulation ends when all cells have been infected or when no more infections occurs after a given number of attempts. The code has the ability to simulate the evolution of the virus in serial passage of cell "cultures", i.e. after the end of a simulation, a new one is immediately scheduled with a new culture of infected cells. The code outputs characteristics of the resulting virus population dynamics and genetic composition of the virus population, such as the top dominant genomes, percentage of a genome with specific characteristics.

Vassilevska, Tanya↗

DOE OSTI2019

Active Mixing of Reactive Materials for 3D Printing

The mixing of materials during additive manufacturing is a major benefit which allows one to compositionally and spatially tailor material properties, for example to locally control the reactivity in fuel–oxidizer systems known as thermites. This work characterizes an active mixing printhead used in conjunction with a 3D printing process known as Direct Ink Writing. Besides compositional control, a major benefit of this approach is that it offers a safe method for working with these materials, which can otherwise be hazardous once mixed. Custom fuel and oxidizer inks are fed at fixed volumetric rates into an active mixing printhead, and both the rotational speed of the mixing impeller as well as the fuel–oxidizer ratio are varied. Upon ignition, the propagation speed increases with the rotational speed of the mixer and plateaus above a critical value of ≈750 rpm. The critical mixing speed is corroborated by computational fluid simulations and an analytical expression that considers the inks’ complex fluid behavior. Additionally, varying the composition results in a wide range of propagation speeds with peak reactivity corresponding to a fuel–rich formulation ( Φ = 1.5). As a result, a test article incorporating a fast– and slow–burning region demonstrates how spatial composition can manipulate the reactivity.

3D printing↗

DOE OSTI2015

SU-E-T-558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry

Purpose: Monte Carlo simulation on GPU has experienced rapid advancements over the past a few years and tremendous accelerations have been achieved. Yet existing packages were developed only in voxelized geometry. In some applications, e.g. radioactive seed modeling, simulations in more complicated geometry are needed. This abstract reports our initial efforts towards developing a quadric geometry module aiming at expanding the application scope of GPU-based MC simulations. Methods: We defined the simulation geometry consisting of a number of homogeneous bodies, each specified by its material composition and limiting surfaces characterized by quadric functions. A tree data structure was utilized to define geometric relationship between different bodies. We modified our GPU-based photon MC transport package to incorporate this geometry. Specifically, geometry parameters were loaded into GPU’s shared memory for fast access. Geometry functions were rewritten to enable the identification of the body that contains the current particle location via a fast searching algorithm based on the tree data structure. Results: We tested our package in an example problem of HDR-brachytherapy dose calculation for shielded cylinder. The dose under the quadric geometry and that under the voxelized geometry agreed in 94.2% of total voxels within 20% isodose line based on a statistical t-test (95% confidence level), where the reference dose was defined to be the one at 0.5cm away from the cylinder surface. It took 243sec to transport 100million source photons under this quadric geometry on an NVidia Titan GPU card. Compared with simulation time of 99.6sec in the voxelized geometry, including quadric geometry reduced efficiency due to the complicated geometry-related computations. Conclusion: Our GPU-based MC package has been extended to support photon transport simulation in quadric geometry. Satisfactory accuracy was observed with a reduced efficiency. Developments for charged particle transport in this geometry are currently in progress.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2016

Final report for the DOE Early Career Award #DE-SC0003912

This DoE supported early career project was aimed at developing computational models, theory and simulation methods that would be then be used to predict assembly and morphology in polymer nanocomposites. In particular, the focus was on composites in active layers of devices, containing conducting polymers that act as electron donors and nanoscale additives that act as electron acceptors. During the course this work, we developed the first of its kind molecular models to represent conducting polymers enabling simulations at the experimentally relevant length and time scales. By comparison with experimentally observed morphologies we validated these models. Furthermore, using these models and molecular dynamics simulations on graphical processing units (GPUs) we predicted the molecular level design features in polymers and additive that lead to morphologies with optimal features for charge carrier behavior in solar cells. Additionally, we also predicted computationally new design rules for better dispersion of additives in polymers that have been confirmed through experiments. Achieving dispersion in polymer nanocomposites is valuable to achieve controlled macroscopic properties of the composite. The results obtained during the course of this DOE funded project enables optimal design of higher efficiency organic electronic and photovoltaic devices and improve every day life with engineering of these higher efficiency devices.

14 SOLAR ENERGY↗

DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Performance limits of ion extraction systems with non-circular apertures

A three-dimensional computer simulation is used to determine the perveance limitations of ion extraction systems with non-circular apertures. The objective of the study is to analyze the possibilities to improve mechanical strength of the ion optics made of carbon-carbon composite materials. Non-circular grid apertures are better suited to the physical structure of carbon-carbon composite materials, than conventionally used circular holes in a hexagonal pattern, because they allow a fewer number of cut fibers. However, the slit-type accelerating systems, usually regarded as the main alternative to the conventional ion optics, have an intolerably narrow range of operating perveance values at which there is no direct ion impingement on the acceleration grid. This paper presents results of comparative analysis of a number of different ion optical systems with non-circular apertures and conventional ion optical systems with circular apertures. It has been revealed that a relatively wide perveance range without direct ion impingement may be obtained with apertures shaped as a square with rounded corners. Numerical simulations show that this geometry may have equivalent perveance range as the traditional geometry with circular apertures while being more mechanically robust. In addition, such important characteristics, as the effective transparency for both the ions and the neutral atoms, the height of the potential barrier reflecting the downstream plasma electrons and the angular divergence of the beamlet also can be very close to these parameters for the optics with circular apertures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2019

Short- and Intermediate-Range Structure and Dynamics of Fe-Ni-C Liquid Under Compression

Properties of liquid Fe alloys under high-pressure conditions are crucial for understanding the composition, thermal state, and dynamics of Earth’s core. Experiments on such liquids, however, are often performed under pressures far below those of the outer core, necessitating long extrapolations of experimental results to core conditions. Such estimates can be complicated by light elements possibly forming pressure-dependent molecular clusters that can significantly affect the physical properties of liquids as core conditions are approached. First-principles molecular dynamics simulations were employed to compute the properties of an Fe-Ni-C liquid with a composition of Fe 3.7 Ni 0.37 C at 1673 K and pressures from 0 to 67 GPa to benchmark computational methods on pressure effects on the structure and properties of the liquid relative to low pressure experimental results. The short-range structure is manifested by the coordination number (CN) of Fe/Ni-Fe/Ni being around 12, indicative of a nearly close-packed structure in the pressure range, and the CN of C-Fe/Ni gradual increasing from 6.5 to 8.5, indicative of an approximately octahedral to cubic transition as pressure increases. The Fe/Ni-Fe/Ni bond distance, however, is found to be 10 times more compressible than the C-Fe/Ni distance. The intermediate-range structure of Fe/Ni-Fe/Ni and C-Fe/Ni subsystems, described by a partial configurationally decomposed distribution function, undergoes substantial changes, characterized by a significant increase of the number of polyhedra that share 3 atoms with each other. Such dense configurations are related to an increased bulk modulus, decreased diffusion coefficient, decreased activation volume for diffusion, and increased shear viscosity. Reproducing the experimental observations at low pressures provides important support for modeling the liquid under the conditions relevant to the outer core.

58 GEOSCIENCES↗

DOE OSTI2020

Molecular-Level Behavior of Imidazolium-Based Ionic Liquid Mixtures

Understanding the molecular-level organization of anions and cations in a mixture of different ionic liquids (ILs) is challenging, but it can be very important for designing solvents and IL-based membrane materials. In this work, we perform all-atom molecular dynamics (MD) simulations to understand the intermolecular interactions and structuring of both neat and mixed IL compositions containing 1-benzyl-3-methylimidazolium [ZMIM + ], 1-butyl-3-methylimidazolium [BMIM + ], bistriflimide [Tf 2 N – ], and tetrafluoroborate [BF 4 – ]. We observe very strong underlying molecular-level migrations in these mixtures, and this tends to mainly involve the anion species. As compared to the systems with a single anion, the mixed anions lead to large increases/decreases in the radial distribution function peak heights, whereas the cation mixing is relatively benign. The [BF 4 – ] anions have a strong tendency to displace the [Tf 2 N – ] anions in the mixtures, and this is seen with respect to both [BMIM + ] and [ZMIM + ].

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Theoretical Study of Doped and Co-doped Gallium Nitride with (Cr, Ni): an LDA-SIC and Monte Carlo Study

Using a composition of first-principle density functional theory (DFT) calculations and Monte Carlo simulation with the calculation within the self-interaction-correction local-density approximation, the electronic, magnetic, and optical properties of doped and co-doped GaN (Ga{sub 1−x}Cr{sub x}N, Ga{sub 1−x}Ni{sub x}N) and GaN (Ga{sub 1−2x}Cr{sub x}Ni{sub x}N) (x = 0.03 and 0.06), respectively, have been investigated. Our results have proven that the half-metallic ferromagnetic state still persists of co-doped GaN. In comparing the total energies, ab initio calculations certify the magnetic-state stability of the ferromagnetic phase compared with the spin-glass state. The exchange interactions obtained from the ab initio calculations were used as input parameters in a classical Ising model by Monte Carlo simulation to confirm the half-metallic ferromagnetic states with high Neel temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Morphologies of dealloying corrosion attack at grain boundaries

Dealloying corrosion at grain boundaries severely compromises the performance of polycrystalline materials across a wide variety of technological applications. The impact of this phenomenon depends upon the morphology and rate of intergranular dealloying, which can range from planar to wormhole-like patterns that rapidly advance into the alloy. Using 2D and 3D multi-phase field simulations, we reveal how diverse microstructures result from a fundamental interplay between alloy composition and a grain boundary migration mechanism that alters diffusional pathways of dealloying. Inside intergranular dealloying channels, corrosion product buildup can spawn new channels that branch into grain interiors, such that alloys can be degraded from the inside out. These processes further lead to an atypical coarsening mode assisted by diffusion in the dealloying agent. We summarize a unifying explanation for dealloying morphology selection in polycrystalline alloys, which provides an important step towards their optimization for dealloying corrosion environments.

Corrosion↗

DOE OSTI2017

Measurement of key resonance states for the P 30 ( p , γ ) S 31 reaction rate, and the production of intermediate-mass elements in nova explosions

We report the first experimental constraints on spectroscopic factors and strengths of key resonances in the P-30(p, gamma)S-31 reaction critical for determining the production of intermediate-mass elements up to Ca in nova ejecta. The P-30(d,n)S-31 reaction was studied in inverse kinematics using the GRETINA gamma-ray array to measure the angle-integrated cross-sections of states above the proton threshold. In general, negative parity states are found to be most strongly produced but the absolute values of spectroscopic factors are typically an order of magnitude lower than predicted by the shell-model calculations employing WBP Hamiltonian for the negative-parity states. The results clearly indicate the dominance of a single 3/2(-) resonance state at 196 keV in the region of nova burning T approximate to 0.10-0.17 GM, well within the region of interest for nova nucleosynthesis. Hydrodynamic simulations of nova explosions have been performed to demonstrate the effect on the composition of nova ejecta.

Kankainen, A. (ORCID:0000000310827602)↗