Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer 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 1,027 records · Page 57

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗

Analysis of chemical stress and the propensity for cracking during the vertical Bridgman growth of BaBrCl:Eu

Computational models are employed to analyze residual chemical stresses arising from compositional variations in europium-doped BaBrCl crystals grown by the vertical Bridgman method. We find that significant chemical stress is produced by radial segregation of Eu in this system. In particular, the distribution of normal stresses is set by the radial concentration gradient, whose changing sign produces surface states in tension or compression. Additionally, crack opening from surface flaws will be promoted or suppressed by tensile or compressive surface stresses, respectively. Thus, crystal growth processing strategies that change the radial dopant concentration gradients are posited to affect the propensity for cracking. For this system, surface stresses are changed from states of tension to compression when the growth rate is increased, thus improving the chances to avoid cracking—a strategy that defies classical wisdom that dictates slower growth to improve outcomes. Similar strategies affecting segregation may prove beneficial to tailor chemical stress fields to reduce cracking in other crystal growth systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sample environment effects on synchrotron-measured temperature profiles in an approximant of optical floating zone crystal growth

Even though the growth of crystals using optical floating zone furnaces has had an immense scientific impact, the implementation of this method remains more of an art than a science due to the difficulty of obtaining quantitative information about the sample thermal profile during crystal growth. Building on recent work demonstrating that in situ synchrotron studies can be used to map sample rod temperatures during heating, investigations were carried out to better understand how the sample environment affects the sample temperature profile. Through a combination of experimental studies and modeling efforts, it is shown that the environment in the furnace can strongly influence the sample temperature at the lamp focus, the steepness of the vertical temperature gradient, and the timescale required for sample heating and cooling - effects which can combine to produce a strong history-dependence to sample temperature profiles. It is demonstrated that the furnace effects can be effectively captured in thermal models, allowing both the steady-state and time-dependent behavior of the sample to be accurately reproduced with predictive models, and providing a launching point for improved furnace designs that can more readily deliver desired thermal profiles.

36 MATERIALS SCIENCE↗

Optimizing $\mathrm{ACRT}$ to reduce inclusion formation during the $\mathrm{VGF}$ growth of cadmium zinc telluride: II. Application to experiments

Here two different accelerated crucible rotation technique (ACRT) rotation schedules are assessed, via theoretical computations and growth experiments, according to their ability to reduce tellurium-rich inclusions during the vertical gradient freeze (VGF) of cadmium zinc telluride (CZT). A vigorous rotation schedule based on classical ACRT guidelines produces a well-mixed melt that is expected to reduce tellurium inclusions compared to growth without crucible rotation. In contrast, a new ACRT schedule, derived from model-based optimization for interface stability, employs slower rotation with longer periods and is predicted to further reduce inclusions. Growth experiments corroborate these expectations. Namely, both ACRT schedules result in crystals with inclusion size and volume significantly decreased from levels found in material grown with no rotation, and material grown using the computationally optimized ACRT parameters exhibits a median inclusion size that is smaller and with a sharper distribution than in material grown via classical ACRT.

36 MATERIALS SCIENCE↗

Optimizing $\mathrm{ACRT}$ to reduce inclusion formation during the $\mathrm{VGF}$ growth of cadmium zinc telluride: I. Computational approach

Theoretical computations reveal insight into parameters that govern the accelerated crucible rotation technique (ACRT) applied to the gradient freeze (VGF) growth of cadmium zinc telluride (CZT). A metric based on the classic Mullins and Sekerka instability criterion is put forth that provides a quantitative means of assessing the impact of different ACRT rotation schedules on morphological stability and inclusion formation during CZT growth. Here this metric is employed in a 2-k factorial design of experiments that identifies the maximum rotation rate as the most important ACRT parameter, followed by the acceleration period and the rest period. Subsequent optimization of rotation schedule is performed for a specific VGF growth system, and mechanistic explanations are provided for the preferred ACRT schedule, which improves interface stability via slower rotation rates and longer acceleration periods than would be applied using classical ACRT schedules.

36 MATERIALS SCIENCE↗

The effects of ACRT on melt undercooling during the growth of CZT via the traveling heater method: Ekman versus Taylor-Görtler flows

For the first time, a comprehensive model for the traveling heater method, including detailed species transport and phase-change thermodynamics, is applied to understand the effects of accelerated crucible rotation (ACRT) on the growth of cadmium zinc telluride. Here, a rotation schedule that is designed from classical considerations does not completely mix the melt, rather undercooled regions of localized melt are swept periodically across the growth interface. Even without melt homogenization, these ACRT-driven flows disrupt the persistent undercooling that occurs during THM growth without rotation. An ACRT schedule designed to promote Ekman flows results in lower levels of melt undercooling than a schedule designed to promote Taylor-Görtler flows and greater mixing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stability-based optimization of ACRT for the growth of CZT by the traveling heater method

Numerical modeling is employed to study the effects of several accelerated crucible rotation technique (ACRT) rotation schedules on the growth of cadmium zinc telluride (CZT) by the traveling heater method (THM). In conjunction with these analyses, a measure is developed to estimate the amplitude of disturbances associated with cellular interface growth morphologies that develop over time and across the melt-solid interface. Furthermore, this disturbance amplitude measure is used to assess the ability of ACRT to stabilize the growth interface and thereby minimize inclusion formation processes. Notably, our analyses reveal that interfacial instability is not well correlated with melt undercooling, explaining why classical ACRT approaches to reduce undercooling by mixing may not stabilize growth. Computations show that interface stability is best achieved via rotation schedules that accentuate outward Ekman flows during spin-up and, during spin-down, eliminate Taylor-Görtler flows and minimize inward Ekman flows.

36 MATERIALS SCIENCE↗

Projected network performance for multiple isotopes using next-generation xenon monitoring systems

Since about 2000 (Bowyer et al., 1998), radioxenon monitoring systems have been under development and testing for the verification of the Comprehensive Nuclear Test-Ban Treaty (CTBT). Operation of the systems since then has resulted in development of a next-generation of systems that are nearly ready for operational deployment. By 2010, the need to screen out civilian sources was well known (Auer et al., 2010; Saey, 2009), and isotopic ratio approaches were soon considered (Kalinowski et al., 2014) to identify specific sources. New generation systems are expected to improve the ability to verify the absence of nuclear tests by using isotopic ratios when multiple isotopes are detected. In this work, thousands of releases were simulated to compute the global detection probability of 131m Xe, 133m Xe, 133 Xe, and 135 Xe at 39 noble gas systems in the International Monitoring System (IMS) for both current and next-generation systems. Three release scenarios are defined at 1 h, 1 d, and 10 d past a 1 kt TNT equivalent 235U explosion event. Multiple cases using from one part in a million to the complete release of the xenon isotopic activity are evaluated for each scenario. Coverage maps and global integrals comparing current and next-generation monitoring systems are presented showing that next-generation noble gas systems will create measurable improvements in the IMS. The global detection probability for 133 Xe is shown to be strong in all scenarios, but only modestly improved by next-generation equipment. However, the detection probability for 131m Xe and 133m Xe increased to about 50% in different scenarios, providing a second detectable isotope for many events. As anticipated from shorter sampling intervals, the expected number of detecting samples roughly doubled and the expected number of detecting stations rose by approximately 50% for all release scenarios. Thus, it might be anticipated that future events would consist of multiple 133 Xe detections and one or more second isotope detections. In conclusion, signals of this nature should increase detection confidence, tighten release location estimates, improve rejection of civilian signals, and lessen the impacts from individual systems being offline for maintenance or repair reasons.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Neutron diffraction discriminates between models for the nanoarchitecture of graphene sheets in glassy carbon

Glassy carbon is a chemically inert, biocompatible, disordered material with the graphene sheet as its basic building block. Structural characterisation techniques have so far been unable to provide definitive distinction between the many proposed models for its structure. Computer based simulation methods have made a step forward by predicting structure from interaction potentials, but the results are sensitive to the choice of the potential. Here we use the white neutron beam of the Spallation Neutron Source at Oak Ridge National Laboratory coupled with the SNAP time-of-flight diffractometer to sample the reciprocal space of glassy carbon and calculate accurate radial distribution functions. From the radial distribution information, we determine graphene sheet dimensions, registration, curvature, and assess defect content to rank the proposed structures in terms of their agreement with experimental data. In this work, we find that the recent models derived from empirical potentials give the best agreement.

36 MATERIALS SCIENCE↗

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗

Two-layer cathode architecture for high-energy density and high-power density solid state batteries

Solid state batteries with high-energy density and high-power density require the development of thick and energy dense cathodes. Structured cathode electrodes with a double-layer configuration were enabled using a freeze tape casting technique. A bottom dense layer was utilized to enhance the energy density whereas a top porous layer with vertically aligned walls was utilized to enhance the power density. The porous structure of the power layer was achieved by ice templating this layer on top of the densified energy layer of the cathode. This configuration was found to better utilize the active material of the cathodes and was optimized using numerical simulation and computer modeling. Cells with Li metal anode and LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) at approximately 5 and 20 mg/cm 2 were cycled at 70 °C at different C-rates. Poly(ethylene oxide) (PEO) with lithium bis-trifluoromethanesulfonimide (LiTFSI) was used for the catholyte and the solid-state electrolyte. The structured cathodes exhibited more than double capacity values as well as better Coulombic efficiency compared to non-structured (single-layer) thick cathodes. In conclusion, synchrotron X-ray tomography and scanning electron microscopy were used to characterize the microstructure of the cathodes.

25 ENERGY STORAGE↗

Advanced neutron absorber Ni-Cr-Mo-Gd alloys seawater corrosion mechanism and susceptibility study

Previous studies demonstrated the corrosion susceptibility of advanced neutron absorber (ANA) Ni-Cr-Mo-Gd in a seawater environment but remained inconclusive. In this paper, scanning electron microscopy was employed to identify the corrosion phase in Ni-Cr-Mo-Gd, and computational thermodynamic simulation was used to study phase corrosion potentials. Results showed that the Ni 5 Gd second phase, distributed along FCC grain boundaries in Ni-Cr-Mo-Gd, has lower corrosion potential than the substrate. Furthermore, heat treatment and alloy chemistry adjustment would not ameliorate the ANA corrosion susceptibility. Finally, isolating the Ni 5 Gd from the corrosion medium through advanced manufacturing is a viable way to improve the ANA corrosion resistance.

36 MATERIALS SCIENCE↗

Thermal mechanical assessment of a SiC-SiC-composite clad fuel pin concept in a light water reactor environment

Accident Tolerant Fuels (ATFs) are designed to increase coping time following an accident scenario while preserving or improving current steady state reactor operational performance. A potential ATF concept is SiC-SiC composite claddings. Fuel performance simulations were conducted on a SiC-SiC based cladding concept utilizing a multilayered approach for improved performance. This cladding concept referred to as the Duplex concept is a duplex structure composed of a monolithic SiC layer placed on the outside of a SiC-SiC composite. A liquid metal is added to fuel-cladding gap for improved heat dissipation from the fuel. The monolithic SiC layer is used to improve the coolant corrosion characteristics and protect the SiC-SiC composite layer from exposure to the coolant. The fuel performance code BISON was used to conduct fuel performance simulations on the cladding concepts. Comparisons are made with a current prototypic fuel rod design (UO 2 fuel enclosed in Zircaloy-4 cladding). Representative steady-state cases were considered for normal power and two cycle power histories. Additionally, a PCI ramp case was simulated to analyze potential anticipated operational occurrences. Transient response during a Loss of Coolant Accident and a Reactivity Initiated Accident were also simulated. This computational study demonstrated that for normal operating conditions, the SiC concept cladding performed as well as the baseline for the standard power cases evaluated. The ramping evaluations indicate potential fracturing of the SiC-SiC composite of the composite cladding compared to the Zircaloy-4 cladding due to the temperature gradient and the subsequent differential thermal conductivity degradation and swelling across the composite thickness. In conclusion, the rod fails early at low enthalpy for RIA but survives a LOCA with minimal material loss due to high temperature steam corrosion.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploiting ν-dependence of projected energy correlators in HICs

We extend the recently derived factorization formula for energy-energy correlators to study the analytic structure of general ν-point projected energy correlators in heavy ion collisions. The ν-point projected energy correlators (or, ν-correlators) are an analytically continued family of the integer N-point projected energy correlators, which probe correlations between N final-state particles. By tracking the largest separation (χ) between the N particles, in vacuum, their structure is closely related to the DGLAP splitting functions and exhibits a classical scaling behavior ∼ 1/χ which is modified by resummation through the anomalous dimensions. We show that, in a thermal medium, the ν-correlators display non-trivial angular scaling already at the leading order in perturbation theory. We find that for non-integer values, particularly ν < 1, medium-induced jet function is enhanced compared to ν > 1. This is particularly manifested in the ratios of ν-correlators with respect to the two-point energy correlator, which encodes intrinsic angular information for ν < 1 when compared to large ν values. Moreover, for small-ν values, the ν-correlators appear to saturate at ν = 0.01 . We further confirm our leading-order numerical computations against simulated events from JEWEL for the parton-level production cross-section. Finally, we qualitatively discuss the effect of BFKL resummation for various values of ν.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Metasurface holographic optical traps for ultracold atoms

We propose metasurface holograms as a novel platform to generate optical trap arrays for cold atoms with high quality, efficiency, and thermal stability. We developed design and fabrication methods to create dielectric, phase-only metasurface holograms based on titanium dioxide. We experimentally demonstrated optical trap arrays of various geometries, including periodic and aperiodic configurations with dimensions ranging from 1D to 3D and up to a few hundred trap sites. We characterized the performance of the holographic metasurfaces in terms of the positioning accuracy, size and intensity uniformity of the generated traps, and power handling capability of the dielectric metasurfaces. Our proposed platform has great potential for enabling fundamental studies of quantum many-body physics, and quantum simulation and computation tasks. The compact form factor, passive nature, good power handling capability, and scalability of generating high-quality, large-scale arrays also make the metasurface platform uniquely suitable for realizing field-deployable devices and systems based on cold atoms.

42 ENGINEERING↗

Joint Optimization of Multimodal Transit Frequency and Shared Autonomous Vehicle Fleet Size with Hybrid Metaheuristic and Nonlinear Programming

Shared autonomous vehicles (SAVs) bring competition to traditional transit services but redesigning multimodal transit network can utilize SAVs as feeders to enhance service efficiency and coverage. This paper presents an optimization framework for the joint multimodal transit frequency and SAV fleet size problem, a variant of the transit network frequency setting problem. The objective is to maximize total transit ridership (including SAV-fed trips and subtracting boarding rejections) across multiple time periods under budget constraints, considering endogenous mode choice (transit, point-to-point SAVs, driving) and route selection, while allowing for strategic route removal by setting frequencies to zero. Due to the problem’s non-linear, non-convex nature and the computational challenges of large-scale networks, we develop a hybrid solution approach that combines a metaheuristic approach (particle swarm optimization) with nonlinear programming for local solution refinement. To ensure computational tractability, the framework integrates analytical approximation models for SAV waiting times based on fleet utilization, multimodal network assignment for route choice, and multinomial logit mode choice behavior, bypassing the need for computationally intensive simulations within the main optimization loop. Applied to the Chicago metropolitan area’s multimodal network, our method illustrates a 33.3% increase in transit ridership through optimized transit route frequencies and SAV integration, particularly enhancing off-peak service accessibility and strategically reallocating resources.

Ng, Max↗

Rapid Determination of Uranium Isotopic Abundance from Cotton Swipes: Direct Extraction via a Planar Surface Reader and Coupling to a Microplasma Ionization Source

The collection of solid particulates and liquids from surfaces by the use of cloth swipes is fairly ubiquitous. In such methods, there is a continuous concern regarding the ability to locate and quantitatively sample the analyte species from the material. In this effort, we demonstrate the initial coupling of an Advion Plate Express plate reader to a liquid sampling–atmospheric pressure glow discharge (LS-APGD) microplasma ionization source with an Orbitrap mass spectrometer to perform uranium isotopic analyses of solution residues on cotton swipes. The Plate Express employs a sampling probe head to engage and seal against the swipe surface. Subsequentially, the analyte residues are desorbed and transported within a 2% HNO 3 electrolyte flow to the ionization source. Quantitative recoveries were observed following a single 30 s extraction step, with the absolute mass sampled per extraction being ~100 ng. While the intrasample variability in the analytical responses for triplicate sampling of the same swipe yield ~30% RSD, this lack of precision is offset by the ability to determine isotope ratios for enriched uranium specimens with a precision of better than 10% RSD. Pooled, intersample precision ( n = 9) was found to be <5%RSD across the various sample compositions. Finally, 235 U/ 238 U determinations (ranging from 0.053 to 1.806) were accurate with errors of <10%, absolute. The 234 U- and 236 U-inclusive ratios were determined with similar accuracy in enriched samples. While the driving force for the effort is in the realm of nuclear nonproliferation efforts, the ubiquitous use of cloth swipes across many application areas could benefit from this convenient approach, including the use of versatile, reduced-format mass spectrometer systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Structure-Driven Liquid Microjunction Surface-Sampling Probe Mass Spectrometry

The rhizosphere is the narrow region of soil surrounding the roots of plants that is influenced by root exudates, root secretions, and associated microbial communities. This region is crucial to plant growth and development and plays a critical role in nutrient uptake, disease resistance, and soil transformation. Understanding the function of exogenous compounds in the rhizosphere starts with determining the spatiotemporal distribution of these molecular components. Using liquid microjunction surface-sampling probe mass spectrometry (LMJ-SSP-MS) and microfluidic devices with attached microporous membranes enables in situ, nondisruptive, and nondestructive spatiotemporal measurement of exogenous compounds from plant roots. However, long imaging times (>2 h) can negatively affect plant heath and limit temporal studies. Here, we present a novel strategy to optimize the number and location of sampling sites on these microporous membrane-covered microfluidic devices. In conclusion, this novel, “structure-driven” sampling workflow takes into consideration the channel structure of the microfluidic device to maximize sampling from the channels and minimize acquisition time (~4× less time in some cases while providing similar chemical image accuracy), thus reducing stress on plants during in situ LMJ-SSP-MS analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗