Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scaling Studies”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tropical root responses to global changes: A synthesis

Tropical ecosystems face escalating global change. These shifts can disrupt tropical forests' carbon (C) balance and impact root dynamics. Since roots perform essential functions such as resource acquisition and tissue protection, root responses can inform about the strategies and vulnerabilities of ecosystems facing present and future global changes. However, root trait dynamics are poorly understood, especially in tropical ecosystems. We analyzed existing research on tropical root responses to key global change drivers: warming, drought, flooding, cyclones, nitrogen (N) deposition, elevated (e) CO 2 , and fires. Based on tree species- and community-level literature, we obtained 266 root trait observations from 93 studies across 24 tropical countries. We found differences in the proportion of root responsiveness to global change among different global change drivers but not among root categories. In particular, we observed that tropical root systems responded to warming and eCO 2 by increasing root biomass in species-scale studies. Drought increased the root: shoot ratio with no change in root biomass, indicating a decline in aboveground biomass. Despite N deposition being the most studied global change driver, it had some of the most variable effects on root characteristics, with few predictable responses. Episodic disturbances such as cyclones, fires, and flooding consistently resulted in a change in root trait expressions, with cyclones and fires increasing root production, potentially due to shifts in plant community and nutrient inputs, while flooding changed plant regulatory metabolisms due to low oxygen conditions. The data available to date clearly show that tropical forest root characteristics and dynamics are responding to global change, although in ways that are not always predictable. This synthesis indicates the need for replicated studies across root characteristics at species and community scales under different global change factors.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Parallel Adjective High-Order CFD Simulations Characterizing SOFIA Cavity Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge-Kutta, and spatially fth-order accurate WENO- 5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

CFD↗
NASA NTRS2015

Parallel Adaptive High-Order CFD Simulations Characterizing SOFIA Cavitiy Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A tempo- rally fourth-order accurate Runge-Kutta, and a spatially fth-order accurate WENO-5Z scheme were used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

Parallell↗
DOE OSTI2025

Massively parallel phase-field simulations targeting exascale

The interface thickness in the phase-field (PF) method limits its simulation scales. Consequently, large-scale PF simulations become prohibitively expensive for resolving the extremely fine microstructures that typically form during rapid solidification processing. This challenge is significant in predicting microstructure evolution in metal additive manufacturing and has been identified by the United States Department of Energy’s Exascale Computing Project. Here, to address this, we develop a multi-GPU and MPI-based massively parallel simulation code, utilizing state-of-the-art algorithms, software, and libraries, for large-scale three-dimensional (3D) PF simulations. We report the first GPU-parallel PF simulations on Frontier (currently the second TOP500 exascale cluster) and Summit machines, taking dendritic growth as an example problem. We evaluate the parallel performance of our implementation using scaling studies with more than 24 000 GPUs (among the largest known computations to date) and the acceleration performance using large-scale simulations of dendritic growth in 3D. Finally, massively parallel GPUs in these supercomputers enabled the first coupled multiscale simulations of laser melting and subsequent dendritic solidification on the scale of a full melt-pool, demonstrating the feasibility of performing PF simulations with a point total over 2 billion grid points within an acceptable time.

Exascale↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
NASA NTRS1980

Noise suppression due to annulus shaping of a conventional coaxial nozzle

Previous studies have shown that increasing the annulus width of a conventional coaxial nozzle with constant bypass velocity will lower the noise level. In the present model-scale study, the annulus was shaped by an eccentric mounting of the annular nozzle with respect to the conical core nozzle. Acoustic measurements were made in the flyover plane below the widest portion of the annulus and at 90 deg and 180 deg from this point. The model-scale spectra are scaled up to engine size (1.07 m diameter) and the perceived noise levels for the eccentric and concentric coaxial nozzles are compared over a limited range of operating conditions. The implications of the acoustic benefits derived from the eccentric nozzle to practical applications are discussed.

Von Glahn, U.↗
DOE OSTI2020

The Role of Pentaerythritol Tetranitrate (PETN) Aging in Determining Detonator Firing Characteristics

Pentaerythritol tetranitrate (PETN) is one of the most commonly used explosives for both military and mining applications. PETN is used in EBW (exploding bridge wire) and EFI (exploding foil initiator) detonator devices. Due the widespread usage of PETN, studies focusing on how the explosive changes with age, primarily its function in detonators has received particular interest. Aging studies have identified several factors that affect PETN detonator function over time, most notably the specific surface area of the explosive powder. Though small-scale studies have been performed on specific detonator systems with isolated batches of material, what is missing is a large statistically-significant study focusing on PETN aging characteristics - such as powder coarsening - and their influence on detonator performance. Herein we report the first large statistically viable PETN aging study focusing on four batches of PETN powder from the same stock of PETN using two stabilizers, polysaccharide and TriPEON. PETN was aged both as a free-flowing powder as well as in modified EBW RP-2 detonators, for one month at 75 °C. The PETN powder has been analyzed chemically using Fisher specific surface area analysis, SEM imaging and light scattering particle size analysis. Additionally the detonators were tested for performance through voltage-sweep threshold testing and function time measurements. Findings from the study indicate that aging at 75 °C for one month significantly changes the specific surface area and particle size of unstabilized PETN, leading to increases in detonator function time, but not detonator threshold values. Powders stabilized with TriPEON displayed less significant increases in function time, while powder stabilized with polysaccharide exhibited no aging effects, despite the high temperature aging.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1999

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗
DOE OSTI2022

Rotating quantum turbulence in the unitary Fermi gas

Quantized vortices carry the angular momentum in rotating superfluids, and are key to the phenomenon of quantum turbulence. Advances in ultracold-atom technology enable quantum turbulence to be studied in regimes with both experimental and theoretical control, unlike the original contexts of superfluid helium experiments. While much work has been performed with bosonic systems, detailed studies of fermionic quantum turbulence are nascent, despite wide applicability to other contexts such as rotating neutron stars. In this paper, we present a large-scale study of quantum turbulence in rotating fermionic superfluids using an accurate time-dependent density functional theory called the superfluid local density approximation. We identify two different modes of turbulent decay in the dynamical equilibration of a rotating fermionic superfluid, and contrast these results with a computationally simpler description provided by the Gross-Pitaevskii equation, which we find can qualitatively reproduce these decay mechanisms if dissipation is explicitly included. These results demonstrate that dissipation mechanisms intrinsic to fermionic superfluids play a key role in differentiating fermionic from bosonic turbulence, which manifests by enhanced damping of Kelvin waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Facility requirements for helicopter noise research

Development of future helicopter noise prediction methods requires an accurate experimental data base for correlation studies or assessment of theoretical prediction methods. The aerodynamic and acoustic performance criteria for conducting experimental studies of helicopter noise using existing test facilities were defined. Requirements for ground-based facilities were described in addition to the limitations associated with full scale and model scale studies. Flight testing methods were also evaluated briefly and the restrictions associated with these approaches were cited. Finally, a general evaluation of ground-based and flight testing methods is given. Based on this presentation, future investigators can select the experimental approach best suited for generating a desired data base. Also, facility improvements need to extend the state of the art in helicopter noise experimental studies so that problems can be identified.

Schlinker, R. H.↗
DOE OSTI2025

Toucan: A performance portable, scalable implementation of the DECA algorithm

In the field of additive manufacturing (AM), cellular automata (CA) is extensively used to simulate microstructural evolution during solidification. However, while traditional CA approaches are relatively fast, they still require a substantial number of time steps, are limited to moderate volumes, and are relatively difficult to improve through parallelism due to the highly localized nature of the solidification front. Here, to address these issues of time to solution and load balancing, we introduce Toucan, a parallel, performance-portable, and scalable code written in C++ with the Kokkos library that leverages the discrete event inspired cellular automata (DECA) algorithm to perform parallel-in-time (PinT) grain growth simulations. Toucan effectively mitigates load balancing issues by distributing the computational workload more evenly across processors, enhancing scalability and efficiency. We conduct both strong and weak scaling studies on up to 64 GPUs on the Frontier supercomputer, demonstrating that Toucan significantly outperforms the current state-of-the-art, time-stepped CA code, ExaCA, on both single and multi-GPU simulations. Even in AM-specific weak scaling scenarios, Toucan maintains near-ideal scaling, in contrast to the linear increase observed with ExaCA due to the moving laser raster pattern. This study highlights Toucan’s potential to transform microstructural simulations in AM by radically improving both efficiency and scalability over existing methods.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Aerosol observations for climate studies

The Stratospheric Aerosol Measurement II and Stratospheric Aerosol and Gas Experiment satellite systems have provided to date more than five years and almost three years, respectively, of data on atmospheric aerosol profiles on a near-global scale. Studies with these unique data sets are developing a global aerosol climatology for the first time and have shown the existence and quantification of polar stratospheric clouds and tropical stratospheric cirrus. In addition, a tropospheric cirrus climatology is evolving. Since these two experiments were launched, a series of large volcanic eruptions have occurred which have greatly impacted the stratospheric aerosol loading. The aerosol layer produced by the eruption of El Chichon, for example, increased the 30 mb temperatures in the northern tropics by as much as 4 C for six months after the eruption. This paper describes in detail, from a climate perspective, the evolving aerosol and cloud climatologies as a function of space and time, and shows the stratospheric dynamics of volcanic injections and their enhancements on stratospheric optical depth and mass loading.

Mccormick, M. P.↗
DOE OSTI2024

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2020

MicroRNAs as diagnostic and prognostic biomarkers of age-related macular degeneration: advances and limitations

A main cause of vision loss in the elderly is age-related macular degeneration (AMD). Among the cellular, biochemical, and molecular changes linked to this disease, inflammation and angiogenesis appear as being crucial in AMD pathogenesis and progression. There are two forms of the disease: dry AMD, accounting for 80–90% of cases, and wet AMD. The disease usually begins as dry AMD associated with retinal pigment epithelium and photoreceptor degeneration, whereas wet AMD is associated with choroidal neovascularization resulting in severe vision impairment. The new vessels are largely malformed, leading to blood and fluid leakage within the disrupted tissue, which provokes inflammation and scar formation and results in retinal damage and detachment. MicroRNAs are dysregulated in AMD and may facilitate the early detection of the disease and monitoring disease progression. Two recent reviews of microRNAs in AMD had indicated weaknesses or limitations in four earlier investigations. Studies in the last three years have shown considerable progress in overcoming some of these concerns and identifying specific microRNAs as biomarkers for AMD. Further large-scale studies are warranted using appropriate statistical methods to take into account gender and age disparity in the study populations and confounding factors such as smoking status.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1984

Representing Multidimensional Systems Using Visual Displays

Techniques employed to represent multiattribute information in an integrated, object display are reviewed and discussed. A study is proposed to investigate the effects of system parameters such as intervariable correlation on the choice of an optimal display format. The results of a psychophysical scaling study of five different displays are presented.

Casey, E. J.↗
DOE OSTI2022

Microfluidic In-Situ Spectrophotometric Approaches to Tackle Actinides Analysis in Multiple Oxidation States

The study and development of present and future processes for the treatment/recycling of spent nuclear fuels require many steps, from design in the laboratory to setting up on an industrial scale. In all of these steps, analysis and instrumentation are key points. For scientific reasons (small-scale studies, control of phenomena, etc.) but also with regard to minimizing costs, risks, and waste, such developments are increasingly carried out on milli- or microfluidic devices. The logic is the same for the chemical analyses associated with their follow-up and interpretation. Due to this, over the last few years, opto–microfluidic analysis devices adapted to the monitoring of different processes (dissolution, liquid–liquid extraction, precipitation, etc.) have been increasingly designed and developed. In this work, we prove that photonic lab-on-a-chip (PhLoC) technology is fully suitable for all actinides concentration monitoring along the plutonium uranium refining extraction (plutonium, uranium, reduction, extraction, or Purex) process. Several PhLoC microfluidic platforms were specifically designed and used in different nuclear research and development (R&D) laboratories, to tackle actinides analysis in multiple oxidation states even in mixtures. The detection limits reached (tens of µmol·L −1 ) are fully compliant with on-line process monitoring, whereas a range of analyzable concentrations of three orders of magnitude can be covered with less than 150 µL of analyte. Finally, this work confirms the possibility and the potential of coupling Raman and ultraviolet–visible (UV–Vis) spectroscopies at the microfluidic scale, opening the perspective of measuring very complex mixtures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Hot electron preheat in hydrodynamically scaled direct-drive inertial confinement fusion implosions on the NIF and OMEGA

Hot electron preheat has been quantified in warm, directly driven inertial confinement fusion implosions on OMEGA and the National Ignition Facility (NIF), to support hydrodynamic scaling studies. These CH-shell experiments were designed to be hydrodynamically equivalent, spanning a factor of 40 in laser energy and a factor of 3.4 in spatial and temporal scales, while preserving the incident laser intensity of 10 15 W/cm 2 . Experiments with similarly low levels of beam smoothing on OMEGA and NIF show a similar fraction (~0.2%) of laser energy deposited as hot electron preheat in the unablated shell on both OMEGA and NIF and similar preheat per mass (~2 kJ/mg), despite the NIF experiments generating a factor of three more hot electrons (~1.5% of laser energy) than on OMEGA (~0.5% of laser energy). This is plausibly explained by more absorption of hot electron energy in the ablated CH plasma on NIF due to larger areal density, as well as a smaller solid angle of the imploding shell as viewed from the hot electron generating region due to the hot electrons being produced at a larger standoff distance in lower-density regions by stimulated Raman scattering, in contrast to in higher-density regions by two-plasmon decay on OMEGA. Finally, the results indicate that for warm implosions at intensities of around 10 15 W/cm 2 , hydrodynamic equivalence is not violated by hot electron preheat, though for cryogenic implosions, the reduced attenuation of hot electrons in deuterium–tritium plasma will have to be considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗