Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Profiling”

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,387 records · Page 77

The importance of round-robin validation when assessing machine-learning-based vertical extrapolation of wind speeds

The extrapolation of wind speeds measured at a meteorological mast to wind turbine rotor heights is a key component in a bankable wind farm energy assessment and a significant source of uncertainty. Industry-standard methods for extrapolation include the power-law and logarithmic profiles. The emergence of machine-learning applications in wind energy has led to several studies demonstrating substantial improvements in vertical extrapolation accuracy in machine-learning methods over these conventional power-law and logarithmic profile methods. In all cases, these studies assess relative model performance at a measurement site where, critically, the machine-learning algorithm requires knowledge of the rotor-height wind speeds in order to train the model. This prior knowledge provides fundamental advantages to the site-specific machine-learning model over the power-law and log profiles, which, by contrast, are not highly tuned to rotor-height measurements but rather can generalize to any site. Furthermore, there is no practical benefit in applying a machine-learning model at a site where winds at the heights relevant for wind energy production are known; rather, its performance at nearby locations (i.e., across a wind farm site) without rotor-height measurements is of most practical interest. To more fairly and practically compare machine-learning-based extrapolation to standard approaches, we implemented a round-robin extrapolation model comparison, in which a random-forest machine-learning model is trained and evaluated at different sites and then compared against the power-law and logarithmic profiles. We consider 20 months of lidar and sonic anemometer data collected at four sites between 50 and 100 km apart in the central United States. We find that the random forest outperforms the standard extrapolation approaches, especially when incorporating surface measurements as inputs to include the influence of atmospheric stability. When compared at a single site (the traditional comparison approach), the machine-learning improvement in mean absolute error was 28 % and 23 % over the power-law and logarithmic profiles, respectively. Using the round-robin approach proposed here, this improvement drops to 20 % and 14 %, respectively. These latter values better represent practical model performance, and we conclude that round-robin validation should be the standard for machine-learning-based wind speed extrapolation methods.

17 WIND ENERGY↗

Wind and Temperature Consensus at Horn Point, HU-Beltsville, Piney Run (Maryland) in support of CoURAGE

The Maryland Department of the Environment (MDE) operates a ground-based atmospheric profiling network consisting of collocated radar wind profilers (RWP) and radio acoustic sounding systems (RASS) as part of its Ambient Air Monitoring Program. This network provides continuous observations of wind and temperature structure in the lower troposphere to support air quality forecasting, regulatory analysis, and atmospheric research. The network currently includes three fixed sites across Maryland: Horn Point (HP, lower eastern shore) [38.587525°,-76.141006°], Howard University-Beltsville (HUB, central Maryland) [39.055277°, -76.878632°], and Piney Run (PR, western Maryland) [39.705950°, -79.012000°] The network is designed to capture regional variability in atmospheric transport and boundary-layer processes. These systems measure vertical profiles of horizontal wind speed and direction using Doppler radar techniques, with observations typically spanning from ~100 m above ground level up to approximately 2.5–4 km. Measurements are derived from the Doppler shift of backscattered electromagnetic signals, enabling retrieval of wind vectors at multiple altitudes with high temporal resolution (e.g., 30-minute averages reported every 6 minutes). Each radar wind profiler is paired with a Radio Acoustic Sounding System (RASS) to provide profiles of virtual temperature in the lower atmosphere (~100–200 m AGL) by measuring the propagation speed of acoustic waves. Together, the RWP/RASS system yields a coupled data set of thermodynamic and kinematic atmospheric structure, including additional parameters such as vertical velocity, radial velocity, signal-to-noise ratio, and spectral width for advanced analysis. There are two types of files for each station: wind data (files with a "w" prefix) and virtual temperature RASS data (files with a "t" prefix). The wind data files are in the format wYYDDD.cns, where YY is the 2-digit year and DDD is the day of the year. The RASS virtual temperature data files are in the format tYYDDD.cns. Each record has the following header structure: Line 1 : Station Name RASS files Line 2 : RASS rev DeTect_2.0, WINDS files Line 2 : WINDS rev ATI 5.1 Line 3 : N latitude, W longitude, and site elevation (m) Line 4 : Date and begin time of consensus: yy mm dd hh mn ss plus # minutes to add to get UTC Line 5 : Consensus averaging time (minutes); number of beams; number of range gates Line 6 : Number of records required to make consensus (num) total number of records (tot) and the consensus window size (m/s) in the format: num:tot (window) RASS files Line 7 : no. of coded cells, no. of spec, pulse width (ns), and inter-pulse period (µs), WINDS files Line 7 : No. of coded cells, no. of spectra, pulse width (ns), and inter-pulse period (µs), each with a pair of values: first value is for oblique beams, second for vertical RASS files Line 8 : Full scale Doppler value (m/s) Delay to first gate (ns) Number of gates Spacing of gates (ns), WINDS files Line 8 : Full scale Doppler velocity (m/s), oblique and vertical Vertical correction applied to oblique beams? (0 = no, 1 = yes) Delay to first gate (ns), oblique and vertical Number of gates, oblique and vertical Spacing of gates (ns), oblique and vertical Line 9 : Azimuth and elevation (9s indicate vertical beam not used) RASS files Line 10, values : HT = Height above ground (km), T = Uncorrected virtual temperature consensus (deg C), Tc = Corrected virtual temperature consensus (deg C), W = Vertical wind consensus (9s indicate vertical beam not used, w-component, positive upward, m/s), CNT = Number of records that made consensus (for the 3 values in same order), SNR = Average signal to noise ratio (dB) of records in consensus (same order) WINDS files Line 10, values : HT = Height above ground (km), SPD = Wind speed (m/s), DIR = Wind direction (deg E of N from N), RAD = Radial velocities for each beam (m/s) in order given in azimuth and elevation line (positive toward radar; 9s indicate vertical beam not used, CNT = Number of records that made consensus, SNR = Average signal to noise ratio (dB) of records in consensus

{"wind speed and direction",temperature}↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Solar PV

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) research platform. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence data centers and other variable loads. This dataset entry describes the behavior of a 1.25-MW proton exchange membrane MC250 electrolyzer system, manufactured by Nel Hydrogen , [1] when fed historical data generated by the 430-kW, fixed-axis solar photovoltaic (PV) array located at NLR’s Flatirons Campus. (While the electrolyzer balance of plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack.) Solar PV power output data for the 2020 calendar year were categorized on a daily basis by total energy generation and standard deviation. Each day was then ranked by these metrics, and the 25th, 50th, and 100th percentiles were selected. The 75th percentile day did not exhibit sufficient variability to make for a valuable experiment. A similar process was used for the related historical wind dataset . [2] The historical days in 2020 that represented these percentiles are Dec. 19, March 29, and May 4, respectively. The entire solar day’s power profile was then fed through the MC250 electrolyzer. Due to its length, the 100th percentile day experiment was split into two parts, and the final 3 hours of the solar day were not captured. These final 3 hours contained no spikes or dips of interest and simply represented a slow decay of input solar power. Also, a single timestamp (13:13:47 on Jan. 14, 2026) was lost in the hydrogen system supervisory control and data acquisition. Finally, during the 25th percentile experiment (solar day Dec. 19, 2020) data recording was lost from 11:00:13 to 11:14:45. The roughly 15 minutes of the solar profile were rerun at the end of the experiment and spliced into this time slot during post-processing. The electrolysis system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operation of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The historical solar profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. For more details on the statistical analysis process, see the slide deck “Public Reference Data for Megawatt-Scale Hydrogen Electrolysis: NLR Historical Solar PV Analysis and Profile Generation” accessible with this data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single solar PV electrolysis experiment and is formatted as: {technology}_{percentile}_{scaling factor} For instance, “solarPV-430kW_25_2x.zip” reports the experiment using the 25th percentile solar data from the historical 2020 solar PV dataset, scaled to 200%. Scaling factors were applied to the generated solar PV power output files to more closely match the 1.25-MW capacity of the electrolyzer. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and solar power input. A PDF file detailing the historical solar data statistical analysis used to generate the solar profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all experiments combined into one dataset labeled "combined_solarPV_experiments.csv". [1] nelhydrogen.com/product/mc-series-electrolyser . [2] data.nlr.gov/submissions/316 .

08 HYDROGEN↗

Radial Wall Jet Flow Over Sigmoidal Surfaces

Radial wall jet flows across flat smooth surfaces have been studied for decades. These studies show that the radial velocity of these jets decays inversely with distance from the nozzle with modest contribution from friction (Poreh, et al., 1967; Rajaratnam, 1976). However, the extent to which flat surface results apply to curved surfaces remains unclear. In this article we explore the influence of settled particle bed slope on these velocity profiles. We model the step change in thickness as a sigmoidal curve of variable steepness and use conservation of momentum to evaluate the velocity profile. We show that the velocity profile attenuates because of curvature. Next, we evaluate the influence of surface curvature on radial wall jet flow. Jet flows over particle beds often introduce curvature in the particle bed profile but the influence of the developed curvature on the velocity profile has not been explored. Here we develop a solution for steady fixed beds based on conservation of momentum. We find that surface curvature has a significant influence on the velocity decay coefficients and celerity of radial wall jets, provided there is velocity slip in the vicinity of the particle bed interface, which is strictly true for particle surfaces. Conservation of momentum predicts conditions where the forward momentum of the flow is directed completely upward. The solution identifies two new dimensionless groups that determine whether a curved surface is sufficient to block radial flow and force flow vertically.

Surface curvature, jet flow, multiphase flow↗

Developing Near Optimal Control Sequences for Chiller Plants with Water-side Economizers: A Case Study in a Warm and Marine Climate

Various advanced control sequences for chiller plants with water-side economizers (WSE) have been proposed in literature, but the optimization of those controls is limited. It is possible to maximize energy savings by developing near-optimal control sequences, which are dependent on several factors such as the load profile. To address these gaps, we first identify an advanced control sequence and three key control parameters for chiller plants with WSE. Next, optimizations are performed to minimize energy consumption for seven combinations of control parameters. A chiller plant with WSE system in a warm and marine climate is studied and two load profiles are considered. The system and controls are modeled using the Modelica Buildings library. The results show optimizing the selected control parameters can reduce energy consumption by up to 11% depending on the load profile. Specifically, optimizing the cooling tower efficiency threshold in the condenser water reset control can significantly reduce energy savings for the variable load profile by efficiently shifting the load from the cooling tower to the chiller. This paper provides practical guidance for developing near-optimal control sequences for chiller plant with WSE systems considering impacts such as the load profile.

chiller plant↗

Precision Kinematic Sunyaev--Zel'dovich Measurements Across Halo Mass and Redshift with DESI DR2 and ACT DR6: Part I. Luminous Red Galaxies

We present the most precise measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect around luminous red galaxies to date, detecting the signal at $18σ$ significance in both harmonic and configuration space. Our analysis cross-correlates 2.4 million spectroscopic LRGs from the Dark Energy Spectroscopic Instrument (DESI) DR2 sample with Data Release 6 (DR6) of the Atacama Cosmology Telescope (ACT). We develop a novel harmonic-space cross-correlation approach using momentum-weighted kSZ templates, yielding nearly uncorrelated bandpowers within a framework consistent with other large-scale structure analyses. By incorporating the LRG halo occupation distribution (HOD) and its uncertainty, we convert measured galaxy gas profiles into halo gas profiles and provide generalized Navarro-Frenk-White (GNFW) fitting profiles, providing empirical targets for tuning feedback efficiency in hydrodynamical simulations and for baryonic modeling in large-scale structure analyses. We find strong evidence that gas profiles do not trace dark matter, providing direct evidence for gas redistribution beyond gravitational collapse. Comparing to hydrodynamical simulations, our measurements favor feedback efficiencies exceeding those in the Battaglia profile, suggesting more efficient gas ejection in group-scale halos than previously predicted. Splitting by redshift, we detect the kSZ signal at SNR $\approx 5$--$10$ in each of four bins and find amplitude evolution consistent with the expected decline in mean halo mass at fixed comoving number density. Splitting by stellar mass, we study the scaling of kSZ amplitude with galaxy properties. Together with BGS and ELG measurements in Paper II, these results span $0.1 \lesssim z \lesssim 1.6$ across three galaxy populations, demonstrating the potential of spectroscopic kSZ to map circumgalactic gas and constrain baryonic feedback.

Qu, F. J. [KIPAC, Menlo Park; Stanford U., HEPL; S↗

Quantifying the Influence of Size, Shape, and Density of Microplastics on Their Transport Modes: A Modeling Approach

Microplastics (MPs) pose significant risks to marine ecosystems and human health, necessitating accurate predictions of their distributions in aquatic environments for effective risk mitigation. However, understanding MP transport dynamics is challenging because of the inadequate representation of MP characteristics such as size, shape, and density in numerical models. Further, the accuracy of the MP vertical profiles in existing models has not been thoroughly validated. Thus, we developed an MP transport model within the Finite Volume Community Ocean Model framework (FVCOM-MP) by integrating MP characteristics. We validated FVCOM-MP against experimental and analytical data, focusing on various MP transport modes and transitions. FVCOM-MP successfully replicates MP profiles in different transport modes, including the bedload, surface-load, suspended-load, and mixed-load modes. Additionally, we introduce phase diagrams for classifying MP transport modes based on particle characteristics, enhancing our understanding of MP dynamics in aquatic systems. The transport modes for a number of real-world MP particles, including fishing line, plastic bag/bottle fragments, synthetic fibers, tire wear particles, polyvinyl chloride and expanded polystyrene pellets, were analyzed with these phase diagrams.

Microplastic transport, Settling velocity, Rising ↗

Impact of vertical and seasonal variation in leaf traits on simulating soybean canopy photosynthesis via 1D and 3D modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax 25 ), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax 25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax 25 profile achieved comparable performance to the scheme with the measured Vcmax 25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R 2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R 2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

Oxide dissolution and oxygen diffusion scenarios in niobium and implications on the Bean–Livingston barrier in superconducting cavities

We generalize a native Nb 2 O 5 dissolution model [G. Ciovati, Appl. Phys. Lett. 89, 022507 (2006)] to sequential overlayer dissolutions, multilayer dissolution, and realistic temperature profiles, which may be applicable to other materials. The model is applied to secondary ion mass spectrometry depth profile measurements for varying temperature profiles and two-step oxide dissolution in Nb and found to agree well. In the context of the Meissner screening response due to impurity profiles on the length scale of the London penetration depth, the shallow diffusion of O impurities results in a substantial decrease in the peak supercurrent density near the surface. In this framework, oxide dissolution and oxygen diffusion can account for a rise in peak supportable magnetic field in SRF cavities with baking time and a suppression after the optimal baking time is reached, in good agreement with peak-field baking temperatures and times as well as recent quench field measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Behavioral and Population Data-Driven Distribution System Load Modeling

Distribution system residential load modeling and analysis for different geographic areas within a utility or an independent system operator territory are critical for enabling small-scale, aggregated distributed energy resources to participate in grid services under Federal Energy Regulatory Commission Order No. 2222 [1]. In this study, we develop a methodology of modeling residential load profiles in different geographic areas with a focus on human behavior impact. First, we construct a behavior-based load profile model leveraging state-of-the-art appliance models. We simulate human activity and occupancy using Markov chain Monte Carlo methods calibrated with the American Time Use Survey data set. Second, we link our model with cleaned Current Population Survey data from the U.S. Census Bureau. Finally, we populate two sets of 500 households using California and Texas census data, respectively, to perform an initial analysis of the load in different geographic areas with various group features (e.g., different income levels). To distinguish the effect of population behavior differences on aggregated load, we simulate load profiles for both sets assuming fixed physical household parameters and weather data. Analysis shows that average daily load profiles vary significantly by income and income dependency varies by locality.

American Time Use Survey↗

Throughput Estimation of Data Transport Networks From Digital Twin Measurements

Digital twins of networked infrastructures, known as Virtual Infrastructure Twins (VITs), are increasingly used for software development, pre-deployment testing, and design space exploration. While VITs avoid the costs and potential disruptions associated with experiments on operational networks, their throughput measurements are typically not sufficiently accurate for performance profiling of wide-area networks that they emulate. Here, machine learning (ML) methods are developed to transform these inaccurate VIT network throughput measurements to closely match in peak and overall profile of those from a physical testbed or production network. First, a micro kernel network reflecting a physical network is utilized to collect one-time measurements on a host to support this ML transformation. Then, a generic multi-modal ML method is developed to learn a map that transforms measurements from subsequent VITs on the same host to match past, current and follow-on testbed and cloud networks. ML generalization equations are derived to establish its correctness and probabilistically guarantee its generalization accuracy. Experimental results are presented for a variety of VIT hosts with target testbed and cloud networks; they include a case study of a four-site science ecosystem wherein inaccurate convex VIT measurement profiles are transformed into accurate concave profiles of target networks.

97 MATHEMATICS AND COMPUTING↗

The Impact of Nonzenith Elevation Angles on Ground-Based Infrared Thermodynamic Retrievals

Observing thermodynamic profiles within the planetary boundary layer is essential to understanding and predicting atmospheric phenomena because of the significant exchange of sensible and latent heat between the land and atmosphere within that layer. The Atmospheric Emitted Radiance Interferometer (AERI) is a ground-based infrared spectrometer used to obtain the vertical profiles of temperature and water vapor mixing ratio. Most AERIs are only capable of zenith views, although the Marine AERI (M-AERI) has a design that allows it to view various elevation angles. In this study, we quantify the improvement in the information content and accuracy of the retrieved profiles when nonzenith angles are included, as is common with microwave radiometer profilers. The impacts of the additional scan angles are quantified through both a synthetic study and with M-AERI observations from the ARM Cloud Aerosol Precipitation Experiment (ACAPEX) campaign. Furthermore, the simulation study shows that low elevation angles contain more information content for temperature whereas high elevation angles have more information content for water vapor. Outside of very humid environments, the addition of low elevation angles also results in lower root-mean-square errors when compared with high angles for both temperature and water vapor mixing ratio, although this is primarily a result of averaging multiple observations together to reduce instrument noise. Real-world results from the ACAPEX dataset indicate similar results as were found for the simulation study, although not all predicted benefits are realized because of the small sample size and observational uncertainties.

54 ENVIRONMENTAL SCIENCES↗

End-Use Savings Shapes: Public Dataset Release for Residential Round 1 [Slides]

The End-Use Load Profiles project created a public database of 900,000 individual building end-use load profiles. Load profiles were modeled to represent the U.S. building stock as it was in 2018, as nearly as possible based on the best available data. The End-Use Savings Shapes follow-on project adds measure impact profiles for energy efficiency and electrification packages to the public dataset. This presentation details the public dataset release on September 20, 2022.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Porewater Geochemical Assessment of Seismic Indications for Gas Hydrate Presence and Absence: Mahia Slope, East of New Zealand’s North Island

We compare sediment vertical methane flux off the Mahia Peninsula, on the Hikurangi Margin, east of New Zealand’s North Island, with a combination of geochemical, multichannel seismic and sub-bottom profiler data. Stable carbon isotope data provided an overview of methane contributions to shallow sediment carbon pools. Methane varied considerably in concentration and vertical flux across stations in close proximities. At two Mahia transects, methane profiles correlated well with integrated seismic and TOPAS data for predicting vertical methane migration rates from deep to shallow sediment. However, at our “control site”, where no seismic blanking or indications of vertical gas migration were observed, geochemical data were similar to the two Mahia transect lines. This apparent mismatch between seismic and geochemistry data suggests a potential to underestimate gas hydrate volumes based on standard seismic data interpretations. To accurately assess global gas hydrate deposits, multiple approaches for initial assessment, e.g., seismic data interpretation, heatflow profiling and controlled-source electromagnetics, should be compared to geochemical sediment and porewater profiles. A more thorough data matrix will provide better accuracy in gas hydrate volume for modeling climate change and potential available energy content.

03 NATURAL GAS↗

Application of radiochromic gel dosimetry to commissioning of a megavoltage research linear accelerator for small–field animal irradiation studies

To develop and implement an efficient and accurate commissioning procedure for small-field static beam animal irradiation studies on an MV research linear accelerator (Linatron-M9) using radiochromic gel dosimetry. The research linear accelerator (Linatron-M9) is a 9 MV linac with a static fixed collimator opening of 5.08 cm diameter. Lead collimators were manually placed to create smaller fields of 2 × 2 cm 2 , 1 × 1 cm 2 , and 0.5 × 0.5 cm 2 . Relative dosimetry measurements were performed, including profiles, percent depth dose (PDD) curves, beam divergence, and relative output factors using various dosimetry tools, including a small volume ionization chamber (A14), GAFCHROMIC™ EBT3 film, and Clearview gel dosimeters. The gel dosimeter was used to provide a 3D volumetric reference of the irradiated fields. The Linatron profiles and relative output factors were extracted at a reference depth of 2 cm with the output factor measured relative to the 2 × 2 cm 2 reference field. Absolute dosimetry was performed using A14 ionization chamber measurements, which were verified using a national standards laboratory remote dosimetry service. Absolute dosimetry measurements were confirmed within 1.4% (k = 2, 95% confidence = 5%). The relative output factor of the small fields measured with films and gels agreed with a maximum relative percent error difference between the two methods of 1.1 % for the 1 × 1 cm 2 field and 4.3 % for the 0.5 × 0.5 cm 2 field. These relative errors were primarily due to the variability in the collimator positioning. The measured beam profiles demonstrated excellent agreement for beam size (measured as FWHM), within approximately 0.8 mm (or less). Film measurements were more accurate in the penumbra region due to the film’s finer resolution compared with the gel dosimeter. Following the van Dyk criteria, the PDD values of the film and gel measurements agree within 11% in the buildup region starting from 0.5 cm depth and within 2.6 % beyond maximum dose and into the fall-off region for depths up to 5 cm. The 2D beam profile isodose lines agree within 0.5 mm in all regions for the 0.5 × 0.5 cm 2 and the 1 × 1 cm 2 fields and within 1 mm for the larger field of 2 × 2 cm 2 . The 2D PDD curves agree within approximately 2% of the maximum in the typical therapy region (1–4 cm) for the 1 × 1 cm 2 and 2 × 2 cm 2 and within 5% for the 0.5 × 0.5 cm 2 field. This work provides a commissioning process to measure the beam characteristics of a fixed beam MV accelerator with detailed dosimetric evaluation for its implementation in megavoltage small animal irradiation studies. Radiochromic gel dosimeters are efficient small-field relative dosimetry tools providing 3D dose measurements allowing for full representation of dose, dosimeter misalignment corrections and high reproducibility with low inter-dosimeter variability. Altogether, radiochromic gels are valuable for fast, full relative dosimetry commissioning in comparison to films for application in high-energy small-field animal irradiation studies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗

A Study of Intermittent Turbulence in Stable Arctic Boundary Layers

Abstract Polar boundary layers are difficult to model due to the existence of intermittent turbulence within stable layers. Here we present a case study evaluation of coherent structures in a stable boundary layer observed during a series of flights with an uncrewed aircraft system (DataHawk2) on 19 October 2016 at Oliktok Point, Alaska as part of the ERASMUS (Evaluation of Routine Atmospheric Sounding Measurements using Unmanned Systems) field campaign. During a sequence of five flights over a nine-hour period, 57 profiles of atmospheric properties (0–400 m a.g.l) were collected. Turbulence was identified using derived Richardson Number, temperature structure function parameter, and turbulence kinetic energy dissipation. Throughout all flights on this strongly stable day, intermittent turbulence was observed. These turbulent layers showed well-mixed potential temperature profiles embedded within otherwise stable potential temperature profiles; often resulting in a characteristic staircase pattern. Turbulent layers ranged from 1 to 30 m deep, with most individual layers being 1–2 m deep. Vertical propagation velocities of layers in the lower atmosphere were on the order of a few cm s −1 , typical of non-convective environments. In different regions of the profile, turbulence was driven by a different balance of buoyancy and shear forces, with turbulence in the near surface environment driven by strong shear forces overcoming strong resistance to buoyancy, while turbulence in elevated layers characterized by weaker shear forces overcoming weaker resistance to buoyancy. We discuss the potential of such datasets for improving subgrid parameterizations of small-scale turbulence embedded within stable boundary layers.

54 ENVIRONMENTAL SCIENCES↗

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗