Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Beam characterization at NSRL for radiobiological experiments—phase 1

An experimental campaign was carried out at the NASA Space Radiation Laboratory to perform an additional, independent dosimetric characterization of the beams of protons, helium and carbon ions for radiobiological experiments. The campaign was undertaken by the request and with the support from the National Cancer Institute, U.S. In this initial phase, the goals were to obtain a first assessment of the dosimetric reproducibility of the beam control system, including analysis of spatial homogeneity and evaluation of ion beam contamination. They should facilitate the design of further experimental campaigns for beam characterization for radiobiological experiments. Measurements included reference dosimetry with comparison of in-house and external ionization chambers and electrometers, lateral-dose profile measurements in air, depth-dose profile in a water tank, evaluation of water equivalent thickness of a HDPE binary range shifter and estimation of impurities of the investigated helium-ion beam. The experiments and results are presented.

43 PARTICLE ACCELERATORS↗
NASA NTRS1993

Low-gravity electrodeposition and growth of polymer thin films with large third-order optical nonlinearities by electrochemical processes for devices: Thiophene-based polymers

It has been proposed that NLO thin film properties may be improved by low-gravity processing. Strong candidates for NLO thin film applications are the polythiophenes. Polymeric thiophenes are attractive materials due to their ease of preparation, stability, and high X(exp 3). A simple and convenient method for preparation of polythiophenes is electrochemical oxidation. We will apply some of our experience and lessons learned in low-gravity metal, metal/cermet electrode position to improve the quality of polythiophene(s) thin films. In low gravity electrode position of Ni at a high rate on an Au substrate often results in the production of an x-ray non diffracting surface. Cobalt metal deposition does not give this result nor does Ni when deposited similarly on a glassy carbon substrate. Co/Ni alloy composition produced during electrode position is strongly dependent upon the amount of convection. Code position of neutral inert cermets with metals is influenced significantly by the presence of gravity and the size of the cermets. Tracks left in the 1-g surfaces by unsuccessful particle occlusion indicate suspension of the large particles is not the only reason for poor volume percentages of the larger particles in the deposits. All size particles are more homogeneously distributed in the deposits in low-gravity electrocodeposition than in 1-g. Low gravity gives larger volume percentages for the larger particles in the deposits, while 1-g gives larger volume percentages for the smaller particles. Intermediate size particles give mixed results. The experimental cells were constructed with flat electrode end plates such that 1-g bench reference electrode positions could be carried out at various orientations with respect to gravity. A series of bench studies using similar designed cells are suggested so that convection modification can be applied to electrochemical thin film preparation. Convection effects can then be coupled with other parameter variations in current, temperature, concentration, solvent, electrolyte, and anode substrates to optimize the properties of NLO thin films before resorting to low-gravity processing.

Penn, Benjamin G.↗
DOE OSTI2020

Heterogeneity effects on nondestructive assay measurements of enrichment in UF 6 cylinders

A set of measurements were performed using multiple mechanically cooled HPGe detectors at six positions around a cylinder uranium hexafluoride to assess if matrix inhomogeneity is detectible and effects the determined apparent uranium enrichment. These measurements were performed on nine 30B cylinders of uranium hexafluoride. Six of the cylinders appeared to be homogeneous and had similar apparent measured uranium enrichment, as calculated by the FRAM software, at all positions. However, three cylinders appeared to have local inhomogeneity based on the observed results. This was manifested as very low apparent enrichments, often <10% of the declared enrichment, on select positions. The side with a low apparent enrichment shows elevated daughter isotope activity (e.g., elevated 1001 keV and 766 keV peak count rate) and slightly reduced 235 U activity (lower 185.7 keV count rate). We hypothesize that these heterogenous cylinders may have experienced asymmetric solar heating, which caused volatile UF 6 to sublime preferentially away from the warmed side. A similar phenomenon was observed at a second facility where dose rate measurements confirm that the “low-enriched side” of cylinders show elevated gamma ray dose rates, likely from the removal of UF 6 attenuation and concentration of daughter-products. Care should be taken during uranium enrichment verification when applying methods that include gamma-rays associated with daughter nuclides to the determination of uranium enrichment on cylinders that are stored in sunlight and have asymmetric dose readings around the cylinder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2015

The Ocean Colour Climate Change Initiative: II. Spatial and Temporal Homogeneity of Satellite Data Retrieval Due to Systematic Effects in Atmospheric Correction Processors

The established procedure to access the quality of atmospheric correction processors and their underlying algorithms is the comparison of satellite data products with related in-situ measurements. Although this approach addresses the accuracy of derived geophysical properties in a straight forward fashion, it is also limited in its ability to catch systematic sensor and processor dependent behaviour of satellite products along the scan-line, which might impair the usefulness of the data in spatial analyses. The Ocean Colour Climate Change Initiative (OC-CCI) aims to create an ocean colour dataset on a global scale to meet the demands of the ecosystem modelling community. The need for products with increasing spatial and temporal resolution that also show as little systematic and random errors as possible, increases. Due to cloud cover, even temporal means can be influenced by along-scanline artefacts if the observations are not balanced and effects cannot be cancelled out mutually. These effects can arise from a multitude of results which are not easily separated, if at all. Among the sources of artefacts, there are some sensor-specific calibration issues which should lead to similar responses in all processors, as well as processor-specific features which correspond with the individual choices in the algorithms. A set of methods is proposed and applied to MERIS data over two regions of interest in the North Atlantic and the South Pacific Gyre. The normalised water leaving reflectance products of four atmospheric correction processors, which have also been evaluated in match-up analysis, is analysed in order to find and interpret systematic effects across track. These results are summed up with a semi-objective ranking and are used as a complement to the match-up analysis in the decision for the best Atmospheric Correction (AC) processor. Although the need for discussion remains concerning the absolutes by which to judge an AC processor, this example demonstrates clearly, that relying on the match-up analysis alone can lead to misjudgement.

Ocean Colour↗
NASA NTRSDate not supplied

Deep Convective Cloud Calibration Sensitivity Studies in Support of Radiometrically Scaling GEO Imagers With VIIRS

The NASA CERES SYN1deg product provides the scientific community regional hourly TOA and surface broadband fluxes and clouds. For consistent geostationary (GEO) derived fluxes and clouds the GEO imagers are radiometrically scaled to the Aqua-MODIS calibration reference. The CERES project utilizes GEO and MODIS or VIIRS analogous channel coincident, collocated, and co-angled radiance pairs as the primary method to inter-calibrate the GEO imagers. Tropical deep convective clouds (DCC) are bright, near Lambertian, top of atmosphere pseudo invariant Earth targets that do not rely on coincident ray-matched radiance pairs to radiometrically scale sensors to a common calibration reference. The DCC invariant target (DCC-IT) methodology collectively analyzes all tropical DCC identified pixel radiances by way of probability density function (PDF) distributions. Perfectly inter-calibrated sensor pairs should reveal nearly identical PDF distributions given the same DCC identification criterion. The PDF median, mean, mode, and inflection point statistics were tested as a function of DCC identification criterion using SNPP-VIIRS and Himawari-8 AHI 0.65μm channel radiances during January 2019. It was found that the PDF inflection point provided inter-calibration factors within 0.25% that were nearly independent of DCC identification criterion. The PDF median provided inter-calibration factors within 0.25% for the coldest BT and most stringent homogeneity factors. The PDF mean and mode statistics were inadequate under any DCC conditions. It is critical for the DCC pixel radiances to be anisotropically corrected. The DCC-IT methodology will also be tested for other visible and SWIR bands.

DCC↗
DOE OSTI2020

Crystal growth and scintillation properties of new ytterbium-activated scintillators Cs 4 CaI 6 :Yb and Cs 4 SrI 6 :Yb

Yb 2+ was investigated as a potential alternative activator for Cs 4 SrI 6 and Cs 4 CaI 6 scintillators for the first time, as opposed to the previously studied Eu2+ activator. Single crystals with nominal Yb 2+ concentrations of 0.5 mol%, 1 mol%, and 3 mol% were grown in Ø7 mm ampoules using the vertical Bridgman method. Luminescence and scintillation properties were evaluated as a function of Yb 2+ concentration and host compound. Here, the 5d → 4f electronic transition of Yb 2+ was observed for both the Sr- and Ca-containing compositions. X-ray induced radioluminescence emission was centered between 449 nm and 463 nm depending on Yb 2+ concentration and host compound. Both the spin-allowed and spin-forbidden transitions were observed in photoluminescence emission spectra and were centered at 446 nm and 476 nm, respectively, for both compounds. The best scintillation performance was achieved with Cs 4 CaI 6 :Yb 1 mol%, which had a 3.5% energy resolution at 662 keV and 43,000 ph/MeV light yield. To our knowledge this is the best energy resolution ever reported for a Yb 2+ -doped scintillator. Additional Cs 4 CaI 6 :Yb 1% crystals were grown in Ø12 mm ampoules to investigate the size dependence of scintillation properties, as well as crystal homogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Influence of Ink Formulation and Drying Conditions on Ionomer Distribution in High-Performance Roll-to-Roll-Coated Gas-Diffusion Electrodes

To enable mass production of fuel cell membrane electrode assemblies (MEAs) catalyst layers production will require continuous roll-to-roll (R2R) coating processes. Gas diffusion electrodes (GDEs) are advantageous for mass production because the catalyst layer can be directely coated on the microporous layer of the gas diffusion media without the need for a decal-transfer process. It is known that the water-to-alcohol ratio in the catalyst ink influences the interactions of the ionomer with the catalyst leading to different distributions of ionomer in spray-coated catalyst layers. It is also known that during drying of colloidal mixtures, like fuel cell inks, factors such as drying rate, particle size, and agglomeration influence how the materials distribute themselves throughout the thickness of the dired film. Thus far there have only been limited studies to understand how process conditions such as ink formulation and drying temperature influence the distribution of ionomer and catalyst coated using scalable methods. This understanding is especially important for GDEs since it is known that having a sufficient amount of ionomer at the catalyst layer-membrane interface is critical for high performance. In this study we have focused on determining how the ratio of water to 1-propanol in the catalyst ink ink and drying temperature influence the distribution of ionomer throughout the thickness of the catalyst layer. Using a combination of Kelvin probe and x-ray photoelectron spectroscopy we show that an ionomer-rich surface is promoted by a higher drying rate and a water-rich catalyst ink. In contrast, a 1-propanol catalyst ink leads to a lower concentration of ionomer on the top surface. Using x-ray computed tomography, we are able to characterize the ionomer distribution throughout the thickness of the layer. We find that, in addition to promoting an ionomer-rich top surface, water-rich inks lead to a more homogenous distribution of ionomer, whereas a 1-propanol-rich ink leads to a more irregular distribution. It is found that MEA performance is improved by selecting conditions and ink formulations that promote ionomer enrichment at the top surface to facilitate a good interface with the membrane. MEAs prepared with a 75 wt% water catalyst ink with a 0.9 I/C have equivalent performance to spray-coated GDEs. Critically, these R2R-coated GDEs do not need an additional ionomer overlayer like the spray-coated GDEs do, reducing the number of processing steps in a manufacturing setting. This work shows that with the appropriate selection of materials, ink formulation, and processing conditions gas-diffusion electrodes are a viable pathway for fuel cell manufacturing.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗
NASA NTRS2001

Portfolios in Stochastic Local Search: Efficiently Computing Most Probable Explanations in Bayesian Networks

Portfolio methods support the combination of different algorithms and heuristics, including stochastic local search (SLS) heuristics, and have been identified as a promising approach to solve computationally hard problems. While successful in experiments, theoretical foundations and analytical results for portfolio-based SLS heuristics are less developed. This article aims to improve the understanding of the role of portfolios of heuristics in SLS. We emphasize the problem of computing most probable explanations (MPEs) in Bayesian networks (BNs). Algorithmically, we discuss a portfolio-based SLS algorithm for MPE computation, Stochastic Greedy Search (SGS). SGS supports the integration of different initialization operators (or initialization heuristics) and different search operators (greedy and noisy heuristics), thereby enabling new analytical and experimental results. Analytically, we introduce a novel Markov chain model tailored to portfolio-based SLS algorithms including SGS, thereby enabling us to analytically form expected hitting time results that explain empirical run time results. For a specific BN, we show the benefit of using a homogenous initialization portfolio. To further illustrate the portfolio approach, we consider novel additive search heuristics for handling determinism in the form of zero entries in conditional probability tables in BNs. Our additive approach adds rather than multiplies probabilities when computing the utility of an explanation. We motivate the additive measure by studying the dramatic impact of zero entries in conditional probability tables on the number of zero-probability explanations, which again complicates the search process. We consider the relationship between MAXSAT and MPE, and show that additive utility (or gain) is a generalization, to the probabilistic setting, of MAXSAT utility (or gain) used in the celebrated GSAT and WalkSAT algorithms and their descendants. Utilizing our Markov chain framework, we show that expected hitting time is a rational function - i.e. a ratio of two polynomials - of the probability of applying an additive search operator. Experimentally, we report on synthetically generated BNs as well as BNs from applications, and compare SGSs performance to that of Hugin, which performs BN inference by compilation to and propagation in clique trees. On synthetic networks, SGS speeds up computation by approximately two orders of magnitude compared to Hugin. In application networks, our approach is highly competitive in Bayesian networks with a high degree of determinism. In addition to showing that stochastic local search can be competitive with clique tree clustering, our empirical results provide an improved understanding of the circumstances under which portfolio-based SLS outperforms clique tree clustering and vice versa.

Mengshoel, Ole J.↗
NASA NTRS1991

Resonant Acoustic Determination of Complex Elastic Moduli

A simple, inexpensive, yet accurate method for measuring the dynamic complex modulus of elasticity is described. Using a 'free-free' bar selectively excited in three independent vibrational modes, the shear modulus is obtained by measuring the frequency of the torsional resonant mode and the Young's modulus is determined from measurement of either the longitudinal or flexural mode. The damping properties are obtained by measuring the quality factor (Q) for each mode. The Q is inversely proportional to the loss tangent. The viscoelastic behavior of the sample can be obtained by tracking a particular resonant mode (and thus a particular modulus) using a phase locked loop (PLL) and by changing the temperature of the sample. The change in the damping properties is obtained by measuring the in-phase amplitude of the PLL which is proportional to the Q of the material. The real and imaginary parts or the complex modulus can be obtained continuously as a function of parameters such as temperature, pressure, or humidity. For homogeneous and isotropic samples only two independent moduli are needed in order to characterize the complete set of elastic constants, thus, values can be obtained for the dynamic Poisson's ratio, bulk modulus, Lame constants, etc.

Brown, David A.↗
DOE OSTI2022

THE IMPACT OF A STATIC-MIXING NOZZLE ON UNIFORMITY IN MATERIAL EXTRUSION LARGE-SCALE ADDITIVE MANUFACTURING

There are many methods of incorporating more than one material in Additive Manufacturing (AM) processes. Oak Ridge National Laboratory has developed a unique solution that enables in-situ material switching by developing a dual-hopper feed system for Cincinnati’s Big Area Additive Manufacturing (BAAM) system. Continuous extrusion during a step-change in material feedstock results in a unique blended material transition region that exhibits a heterogeneous internal morphology. To improve mixing of materials during extrusion, a customized static-mixing nozzle was created for use with the BAAM. Single-bead transitions from Material A to B and B to A were printed with the mixing nozzle at a specified screw speed. Compositional analysis tracked the progression of the material transition as a function of extrudate volume. The resulting transition curves were compared against a standard nozzle configuration. Optical microscopy of cross sections also demonstrated that the static-mixing nozzle promoted a more uniform bead geometry as well as a more homogeneous internal structure throughout the material transition.

Brackett, James↗
NASA NTRS1989

Explicit modeling of composite plates and beams in the dynamics of multibody systems

The state of the art dynamic response analysis of flexible multibody systems is currently restricted to elastic bodies with homogeneous materials. The requirements for high speed operation has made it necessary to use lightweight multi layered composite bodies in robotic systems and space structure applications. Dynamic modeling and analysis of such systems are particularly important since the effects of body flexibility to the performance are likely to be more pronounced. The eight-noded isoperimetric quadrilateral element with independent rotational and displacement degrees of freedom is extended to laminated composite elements. The element includes an arbitrary number of bonded layers, each of which may have a different thickness. The transverse shear deformation which is a predominant factor in the analysis of laminated composite structures is taken into account in developing the stiffness and mass matrices. The corresponding 3-D mode shapes are then incorporated to the multibody system dynamical equations. Floating body reference frames allow the selection of different boundary conditions, and the dynamical equations contain all the nonlinear interactions between the rigid and elastic motion. Example simulations are presented to illustrate the methods proposed.

Amirouche, F. M. L.↗
DOE OSTI2022

Thermal Neutron Scattering Law Evaluation for Zirconium Carbide and Critical Mass Calculations

Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗
DOE OSTI2022

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

2025 Review and Revision of Federal Guidance Report 15

Federal Guidance Report No. 15 (FGR 15), External Exposure to Radionuclides in Air, Water and Soil, published in 2019 and referred to below as FGR 2019, provides age-specific effective dose rate coefficients for reference persons externally exposed to each of 1252 radionuclides homogenously distributed in environmental media. Soon after completion of FGR 2019, the International Commission on Radiological Protection (ICRP) published a similar report (ICRP Publication 144, 2020) addressing the same radionuclides and many of the external exposure scenarios addressed in FGR 2019. In 2023 Argonne National Laboratory (ANL) published a review of dose coefficients in FGR 2019, concluding that “The external effective dose from beta radiation is not appropriately accounted for in FGR 15 for both low- and high-energy beta emitters in air, in soil, and on soil surfaces.” That conclusion was based largely on comparisons of dose coefficients in FGR 2019 with values in ICRP Publication 144 but also on consideration of some unexpected patterns of equivalent dose rate coefficients across tissues and of effective dose rate coefficients across different soil depths, for a selected set of low-energy beta emitters. In response to the ANL report, the Center for Radiation Protection Knowledge (CRPK) at Oak Ridge National Laboratory (ORNL) performed an extensive reexamination of the methods and published values of FGR 2019. CRPK found that coding errors had resulted in inaccurate estimates, primarily overestimates, in many of the dose coefficients tabulated in FGR 2019 but that many of the differences between results in FGR 2019 and ICRP Publication 144 could be traced to differences in methodology. CRPK has corrected all errors in FGR 2019 associated with the coding errors and has taken the opportunity to improve a major portion of the remaining dose coefficients in FGR 2019, primarily through additional Monte Carlo calculations resulting in improved statistics for tissue equivalent dose rate coefficients for exposure to monoenergetic sources. This has eliminated the need for extrapolation of results from high and medium energies to low energies, as was done in FGR 2019. In addition, inconsistencies in the methodology identified in the review, such as computational representation of the adult female, were eliminated. The revised effective dose rate coefficients are consistent with values in ICRP Publication 144 except for differences in values clearly arising from differences in methodology; and differences in values for a relatively small set of very low-energy radionuclides with highly uncertain dose coefficients regardless of methodology.

54 ENVIRONMENTAL SCIENCES↗