Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial accelerator”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Application of electron paramagnetic resonance imaging to the characterization of the Ultem(R) exposed to 1 MeV electrons. Correlation of radical density data to tiger code calculations

A major long-term goal of the Materials Division at the NASA Langley Research Center is the characterization of new high-performance materials that have potential applications in the aircraft industry, and in space. The materials used for space applications are often subjected to a harsh and potentially damaging radiation environment. The present study constitutes the application of a novel technique to obtain reliable data for ascertaining the molecular basis for the resilience and durability of materials that have been exposed to simulated space radiations. The radiations of greatest concern are energetic electrons and protons, as well as galactic cosmic rays. Presently, the effects of such radiation on matter are not understood in their entirety. It is clear however, that electron radiation causes ionization and homolytic bond rupture, resulting in the formation of paramagnetic spin centers in the polymer matrices of the structural materials. Since the detection and structure elucidation of paramagnetic species are most readily accomplished using Electron Paramagnetic Resonance (EPR) Spectroscopy, the NASA LaRC EPR system was brought back on-line during the 1991 ASEE term. The subsequent 1992 ASEE term was devoted to the adaptation of the EPR core system to meet the requirements for EPR Imaging (EPRI), which provides detailed information on the spatial distribution of paramagnetic species in bulk media. The present (1994) ASEE term was devoted to the calibration of this EPR Imaging system, as well as to the application of this technology to study the effects of electron irradiation on Ultem(exp R), a high performance polymer which is a candidate for applications in aerospace. The Ultem was exposed to a dose of 2.4 x 10(exp 9) Rads (1-MeV energy/electron) at the LaRC electron accelerator facility. Subsequently, the exposed specimens were stored in liquid nitrogen, until immediately prior to analyses by EPRI. The intensity and dimensions of the EPR Images that were generated for the irradiated specimens showed that the electrons penetrated the material to a depth of approximately 0.125 inch. These data show a very high degree of correlation to the energy deposition profile as predicted by the Tiger Code, a Monte Carlo code that provides guidelines for the transport of electrons in matter. Subsequent efforts will focus on delineating the transport properties of energetic protons in Ultem(R).

Suleman, Naushadalli K.↗
NASA NTRS2005

Precipitation Modeling over the Greenland and Antarctic Ice Sheets and the Relationship to the Surface Mass Balance and Climate

Atmospheric numerical simulation and dynamic retrieval method with atmospheric numerical analyses are used to assess the spatial and temporal variability of Antarctic precipitation for the last two decades. First, the Polar MM5 has been run over Antarctica to study the Antarctic precipitation. With a horizontal resolution of 60km, the Polar MM5 has been run for the period of July 1996 through June 1999 in a series of short-term forecasts from initial and boundary conditions provided by the ECMWF operational analyses. In comparison with climatological maps, the major features of the spatial distribution of Antarctic precipitation are well captured by the Polar MM5. Drift snow effects on redistribution of surface accumulation over Antarctica are also assessed with surface wind fields from Polar MM5 in this study. There are complex divergence and convergence patterns of drift snow transport over Antarctica, especially along the coast. It is found that areas with large drift snow transport convergence and divergence are located around escarpment areas where there is large katabatic wind acceleration. In addition, areas with large snow transport divergence are generally accompanied by areas with large snow transport convergence nearby, indicating that drift snow transport is of local importance for the redistribution of the snowfall

Bromwich, David H.↗
DOE OSTI2026

Multi-resolution enhancement for full-spectrum neural representations

Scientific data acquisition continues to outpace storage and analysis capabilities, making voxel-basedrepresentations increasingly intractable. Implicit neural representations (INRs) offer a promising solutionby encoding signals through coordinate-based neural networks, serving as surrogates of data, withcomputational and storage requirements scaling with network complexity rather than data dimensionality.However, smaller INRs struggle to faithfully represent multiscale structures, high-frequency informationand fine textures that constitute a large proportion of scientific measurements. We propose WIEN-INR, atheoretically guided hierarchical INR framework that distributes modelling across resolution scales andenables improved representation capacity through a novel enhancement network to recover subtle details.This multiscale architecture allows smaller networks to retain the full spatial-frequency content of thesignal as well as preserve training efficiency and lower storage cost. Evaluated on distinct raw experimentalmeasurements across scales and complexities, WIEN-INR represents a practical step towards a broaderadoption of neural representations in scientific workflows, delivering compact, robust and high-fidelityrepresentations.

Ni, Yuan [SLAC National Accelerator Laboratory (SL↗
DOE OSTI2025

A Printed Microscopic Universal Gradient Interface for Super Stretchable Strain‐Insensitive Bioelectronics

Abstract Stretchable electronics capable of conforming to nonplanar and dynamic human body surfaces are central for creating implantable and on‐skin devices for high‐fidelity monitoring of diverse physiological signals. While various strategies have been developed to produce stretchable devices, the signals collected from such devices are often highly sensitive to local strain, resulting in inevitable convolution with surface strain‐induced motion artifacts that are difficult to distinguish from intrinsic physiological signals. Here all‐printed super stretchable strain‐insensitive bioelectronics using a unique universal gradient interface (UGI) are reported to bridge the gap between soft biomaterials and stiff electronic materials. Leveraging a versatile aerosol‐based multi‐materials printing technique that allows precise spatial control over the local stiffnesses with submicron resolution, the UGI enables strain‐insensitive electronic devices with negligible resistivity changes under a 180% uniaxial stretch ratio. Various stretchable devices are directly printed on the UGI for on‐skin health monitoring with high signal quality and near‐perfect immunity to motion artifacts, including semiconductor‐based photodetectors for sensing blood oxygen saturation levels and metal‐based temperature sensors. The concept in this work will significantly simplify the fabrication and accelerate the development of a broad range of wearable and implantable bioelectronics for real‐time health monitoring and personalized therapeutics.

Song, Kaidong [Department of Aerospace and Mechani↗
DOE OSTI2023

Cosmological consequences of first-order general-relativistic viscous fluid dynamics

In this work, we investigate the out-of-equilibrium dynamics of viscous fluids in a spatially flat Friedmann-LemaîtreRobertson-Walker cosmology using the most general causal and stable viscous energy-momentum tensor defined at first order in spacetime derivatives. In this new framework a pressureless viscous fluid having equilibrium energy density ρ can evolve to an asymptotic future solution in which the Hubble parameter approaches a constant while ρ → 0, even in the absence of a cosmological constant (i.e., Λ = 0). Thus, while viscous effects in this model drive an accelerated expansion of the universe, the equilibrium energy density itself vanishes, leaving behind only the acceleration. This behavior emerges as a consequence of causality in first-order theories of relativistic fluid dynamics and it is fully consistent with Einstein’s equations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Embracing fine-root system complexity in terrestrial ecosystem modeling

Projecting the dynamics and functioning of the biosphere requires a holistic consideration of whole-ecosystem processes. However, biases toward leaf, canopy, and soil modeling since the 1970s have constantly left fine-root systems being rudimentarily treated. As accelerated empirical advances in the last two decades establish clearly functional differentiation conferred by the hierarchical structure of fine-root orders and associations with mycorrhizal fungi, a need emerges to embrace this complexity to bridge the data-model gap in still extremely uncertain models. In this work, we propose a three-pool structure comprising transport and absorptive fine roots with mycorrhizal fungi (TAM) to model vertically resolved fine-root systems across organizational and spatial–temporal scales. Emerging from a conceptual shift away from arbitrary homogenization, TAM builds upon theoretical and empirical foundations as an effective and efficient approximation that balances realism and simplicity. A proof-of-concept demonstration of TAM in a big-leaf model both conservatively and radically shows robust impacts of differentiation within fine-root systems on simulating carbon cycling in temperate forests. Theoretical and quantitative support warrants exploiting its rich potentials across ecosystems and models to confront uncertainties and challenges for a predictive understanding of the biosphere. Echoing a broad trend of embracing ecological complexity in integrative ecosystem modeling, TAM may offer a consistent framework where modelers and empiricists can work together toward this grand goal.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗
DOE OSTI2025

FracML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage

Poster on “FRACML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. The accurate characterization of subsurface fracture networks is essential for the secure operation of carbon capture, utilization, and storage (CCUS) projects. A thorough understanding of the spatial distribution of subsurface faults and fractures is crucial for predicting CO2 plume evolution and minimizing risks such as potential leakage into overlying formations or induced seismicity. In this context, robust fracture network quantification plays a pivotal role in reservoir management, providing the data necessary to fine-tune operational parameters, and ensure the environmental and economic viability of CCUS projects. As part of the U.S. Department of Energy’s SMART (Science-informed Machine Learning for Accelerating Real-time Decisions in Subsurface Applications) initiative, we focused on the development and application of a machine learning-based tool (FRACML) designed to quantify and map fracture networks using real-world (non-synthetic) data from an active CO2 injection site. Our objective is to demonstrate the utility of this tool in improving operational efficiency and safety across CCUS sites.

artifical intelligence / machine learning (AI/ML)↗
NASA NTRS1994

Ground and satellite observations of postdawn aurorae near the time of a sudden storm commencement

Meridian scanning photometer measurements taken in the magnetic postdawn sector at Longyearbyen, Svalbard, between 0300 and 0630 UT on December 29, 1981, are analyzed in conjunction with particle and field data retrieved during two near passes of the Dynamics Explorer 2 (DE 2) satellite. The interval included a sudden storm commencement (SSC) at 0455 UT. Pre-SSC optical and particle measurements showed a system of arcs that are spaced at approximately 1.1 deg intervals in magnetic latitude, embedded within the region 1 current system and span the convection reversal. The softer particle precipitation appears to have a source near the flanks of the magnetotail while the harder, more equatorward precipitation originates closer to Earth. During the SSC period the entire sky brightened, with enhanced 630.0-nm emissions extending from the northern horizon to south of magnetic zenith; intense but spatially separated 557.7-nm emissions dominated the southern horizon. DE 2 detected more than an order of magnitude increase and near isotropization of ring current electron fluxes, enhanced precipitation from the plasma sheet and significantly decreases auroral zone convection. Region 1/region 2 currents remained, with wavelike structures superposed. A dual timescale response to the SSC is consistent with ground and satellite measurements. On few minute travel timescales for hydromagnetic waves to pass through the system, magnetospheric particles accelerate and precipitate to increase the ionospheric conductivity. Global, field-aligned currents change more slowly. To maintain similar field-aligned currents with higher ionospheric conductances requires reduced electric fields. After 0520 UT the optical emissions settled into stable, but latitudinally separated bands of 630.0- and 557.7-nm emissions characteristic of cleft and plasma sheet precipitation, respectively.

Egeland, A.↗
NASA NTRS1999

Experimental Studies of the Interaction Between a Parallel Shear Flow and a Directionally-Solidifying Front

It has long been recognized that flow in the melt can have a profound influence on the dynamics of a solidifying interface and hence the quality of the solid material. In particular, flow affects the heat and mass transfer, and causes spatial and temporal variations in the flow and melt composition. This results in a crystal with nonuniform physical properties. Flow can be generated by buoyancy, expansion or contraction upon phase change, and thermo-soluto capillary effects. In general, these flows can not be avoided and can have an adverse effect on the stability of the crystal structures. This motivates crystal growth experiments in a microgravity environment, where buoyancy-driven convection is significantly suppressed. However, transient accelerations (g-jitter) caused by the acceleration of the spacecraft can affect the melt, while convection generated from the effects other than buoyancy remain important. Rather than bemoan the presence of convection as a source of interfacial instability, Hurle in the 1960s suggested that flow in the melt, either forced or natural convection, might be used to stabilize the interface. Delves considered the imposition of both a parabolic velocity profile and a Blasius boundary layer flow over the interface. He concluded that fast stirring could stabilize the interface to perturbations whose wave vector is in the direction of the fluid velocity. Forth and Wheeler considered the effect of the asymptotic suction boundary layer profile. They showed that the effect of the shear flow was to generate travelling waves parallel to the flow with a speed proportional to the Reynolds number. There have been few quantitative, experimental works reporting on the coupling effect of fluid flow and morphological instabilities. Huang studied plane Couette flow over cells and dendrites. It was found that this flow could greatly enhance the planar stability and even induce the cell-planar transition. A rotating impeller was buried inside the sample cell, driven by an outside rotating magnet, in order to generate the flow. However, it appears that this was not a well-controlled flow and may also have been unsteady. In the present experimental study, we want to study how a forced parallel shear flow in a Hele-Shaw cell interacts with the directionally solidifying crystal interface. The comparison of experimental data show that the parallel shear flow in a Hele-Shaw cell has a strong stabilizing effect on the planar interface by damping the existing initial perturbations. The flow also shows a stabilizing effect on the cellular interface by slightly reducing the exponential growth rate of cells. The left-right symmetry of cells is broken by the flow with cells tilting toward the incoming flow direction. The tilting angle increases with the velocity ratio. The experimental results are explained through the parallel flow effect on lateral solute transport. The phenomenon of cells tilting against the flow is consistent with the numerical result of Dantzig and Chao.

Zhang, Meng↗
NASA NTRS2002

Progress In Plasma Accelerator Development for Dynamic Formation of Plasma Liners

An experimental plasma accelerator for magnetic target fusion (MTF) applications under development at the NASA Marshall Space Flight Center is described. The accelerator is a coaxial pulsed plasma thruster (Figure 1). It has been tested experimentally and plasma jet velocities of approx.50 km/sec have been obtained. The plasma jet has been photographed with 10-ns exposure times to reveal a stable and repeatable plasma structure (Figure 2). Data for velocity profile information has been obtained using light pipes and magnetic probes embedded in the gun walls to record the plasma and current transit respectively at various barrel locations. Preliminary spatially resolved spectral data and magnetic field probe data are also presented. A high speed triggering system has been developed and tested as a means of reducing the gun "jitter". This jitter is being characterized and future work for second generation "ultra-low jitter" gun development is being identified.

Thio, Y. C. Francis↗
DOE OSTI2020

A High Precision Measurement of the Proton Charge Radius at JLab

The elastic electron-proton (e-p) scattering and the spectroscopy of hydrogen atoms are the two traditional methods to determine the proton charge radius (r p ). In 2010, a new method using muonic hydrogen (µH) 1 spectroscopy reported a r p result that was nearly ten times more precise but significantly smaller than the values from the compilation of all previous r p measurements. This discrepancy is often referred to as the "proton charge radius puzzle". In order to investigate the puzzle, the PRad experiment (E12-11-106) was first proposed in 2011 and performed in 2016 in Hall B at the Thomas Jefferson National Accelerator Facility, with both 1.1 and 2.2 GeV electron beams. The experiment measured the e-p elastic scattering cross sections in an unprecedented low values of momentum transfer squared region (Q 2 = 2.1× 10 -4 -0.06 (GeV/c) 2 ), with a sub-percent precision. The PRad experiment utilized a calorimetric method that was magnetic-spectrometer-free. Its detector setup included a large acceptance and high resolution calorimeter (HyCal), and two large-area, high spatial-resolution Gas Electron Multiplier (GEM) detectors. To have a better control over the systematic uncertainties, the absolute e-p elastic scattering cross section was normalized to that of the well-known Møller scattering process, which was measured simultaneously during the experiment. For each beam energy, all data with different Q 2 were collected simultaneously with the same detector setup, therefore sharing the same integrated luminosity. The windowless H 2 gas-flow target utilized in the experiment largely removed a typical background source, the target cell windows. The proton charge radius was determined as r p = 0.831±0.007 stat. ±0.012 p fm, which is smaller than the average r p from previous e-p elastic scattering experiments, but in agreement with the µH spectroscopic results within the experimental uncertainties.

Source record↗
DOE OSTI2020

Runaway electron beam dynamics at low plasma density in DIII-D: energy distribution, current profile, and internal instability

Parameters of the post-disruption runaway electron (RE) beam in the collisionless background plasma achieved after deuterium injection are investigated in DIII-D. The spatially resolved RE energy distribution function is measured for the first time during the RE plateau stage by inverting hard X-ray bremsstrahlung spectra. It has maximum energy up to 20 MeV and a non-monotonous feature at 5{6 MeV observed only in the core of the beam supporting the possibility of kinetic instabilities. The RE current profile is reconstructed for the first time using the spatially resolved RE energy distribution. It is found to be more peaked than the pre-disruption plasma current, with higher internal inductance, suggesting preferential formation of REs in the core plasma or potentially a radially inward motion of the REs. The accessed relatively low-current (180 kA) RE beam is found stable, likely due to its elevated q profile. From this base stable equilibrium, instability is accessed by ramping the solenoid current. Under a large applied accelerating loop voltage an internal kink mode is observed. The kink mode leads to sawtooth-like relaxation of the RE current profile, but drives no RE loss. Under conditions of a large decelerating voltage, chirping Alfvénic instabilities driven by REs during the RE plateau stage are observed for the first time. The low-frequency (0.1{7 MHz) Alfvénic instabilities, supposedly compressional Alfvén eigenmodes, are found to be correlated with a RE loss signal. Finally, this provides novel plasma conditions to study these instabilities and investigate opportunities to utilize them for RE control. Bifurcated electron cyclotron emission (ECE) spectrum is measured during the collisionless RE plateau, with a break point at ≈100 GHz suggesting resonant absorption of the ECE at low frequencies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Distributed Optimization for Nonrigid Nano-Tomography

Resolution level and reconstruction quality in nano-computed tomography (nano-CT) are in part limited by the stability of microscopes, because the magnitude of mechanical vibrations during scanning becomes comparable to the imaging resolution, and the ability of the samples to resist radiation induced deformations during data acquisition. In such cases, there is no incentive in recovering the sample state at different time steps like in time-resolved reconstruction methods, but instead the goal is to retrieve a single reconstruction at the highest possible spatial resolution and without any imaging artifacts. Here we propose a distributed optimization solver for tomographic imaging of samples at the nanoscale. Our approach solves the tomography problem jointly with projection data alignment, nonrigid sample deformation correction, and regularization. Projection data consistency is regulated by dense optical flow estimated by Farneback's algorithm, leading to sharp sample reconstructions with less artifacts. Synthetic data tests show robustness of the method to Poisson and low-frequency background noise. We accelerated the solver on multi-GPU systems and validated the method on three nano-imaging experimental data sets.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

Remote Sensing of Solar-Induced Chlorophyll Fluorescence (SIF) in Vegetation: 50 Years of Progress

Remote sensing of solar-induced chlorophyll fluorescence (SIF) is a rapidly advancing front in terrestrial vegetation science, with emerging capability in space-based methodologies and diverse application prospects. Although remote sensing of SIF – especially from space – is seen as a contemporary new specialty for terrestrial plants, it is founded upon a multi-decadal history of research, applications, and sensor developments in active and passive sensing of chlorophyll fluorescence. Current technical capabilities allow SIF to be measured across a range of biological, spatial, and temporal scales. As an optical signal, SIF may be assessed remotely using highly-resolved spectral sensors and state-of-the-art algorithms to distinguish the emission from reflected and/or scattered ambient light. Because the red to far-red SIF emission is detectable non-invasively, it may be sampled repeatedly to acquire spatio-temporally explicit information about photosynthetic light responses and steady-state behaviour in vegetation. Progress in this field is accelerating with innovative sensor developments, retrieval methods, and modelling advances. This review distills the historical and current developments spanning the last several decades. It highlights SIF heritage and complementarity within the broader field of fluorescence science, the maturation of physiological and radiative transfer modelling, SIF signal retrieval strategies, techniques for field and airborne sensing, advances in satellite-based systems, and applications of these capabilities in evaluation of photosynthesis and stress effects. Progress, challenges, and future directions are considered for this unique avenue of remote sensing.

Review↗
NASA NTRS1990

Spatial distribution of H2O(+) in Comet P/Halley

Comet Halley spectroscopic and imaging data near the time of the Vega 1 encounter in March 1986 have been analyzed for H2O(+) distribution and near-nucleus ion physics, using a data-reduction procedure encompassing flat fielding, sky removal, continuum subtraction, contamination correction, absolute flux calibration, and conversion to ion-column densities. A 25,000 km-radius collision zone within which collisions dominate the solar wind's sweeping effect is verified. A production rate of H2O(+) ions is calculated on the basis of the average velocity and column density integral across the tail. An analysis of the tailward ion flow indicates constant acceleration beyond 130,000 km from the nucleus.

Disanti, Michael A.↗
DOE OSTI2024

A Study of the Transition to Turbulence in a Bed of 67 Spherical Pebbles

Packed beds are commonly found in many engineering systems and have been widely studied for decades. A relatively new packed bed system is the Pebble Bed Reactor, a type of generation-IV nuclear reactor. Unlike many of the packed beds encountered in chemical and process engineering applications, Pebble Bed Reactors are larger and operate at significantly higher Reynolds numbers. As a result of these differences, there is a very limited amount of information on the detailed flow physics that exist in these complex geometries. This work seeks to contribute to a growing database of flow data for Pebble Bed Reactor systems by performing Direct Numerical Simulations of the flow in an experimental bed of 67 pebbles for a range of conditions. Simulations are performed at a Prandtl number of 0.66 and Reynolds numbers from 300–600. These Reynolds numbers are chosen to gain additional knowledge on the spatial development of turbulence in these systems. Analysis of the Turbulent Kinetic Energy, turbulence anisotropy, and Turbulent Heat Flux is performed. Results demonstrate significant development of the TKE across the tested range of Reynolds numbers. Examination of both the TKE and THF reveal that development first occurs near the center of the bed and propagates radially as the flow moves further into the bed. Notable regions of negative production of turbulent kinetic energy are observed in regions where flow accelerates around pebble contact points. Furthermore, these regions are found to coincide with regions of 1-component turbulence.Kindly check and confirm, all authors email id is correctly identified.These are correct

Direct numberical simulation↗