Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity model evaluation”

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 253 records · Page 14

Sound Speed as a Candidate for Internal Temperature Monitoring During Solid Phase Processing of Materials

Internal temperature is an important parameter for understanding the solid phase processing (SPP) of materials and improving feedback control of these processes. Friction stir processes, shear assisted processing and extrusion (ShAPE ™ ), and ultrahigh velocity (UHV) cold spray are examples of such processes. Simple thermal models based on temperatures measured at the boundary of the deforming body cannot capture the dynamic nature of temperature generation within the deforming solid. One possibly way of understanding internal temperature is by measuring sound speed, which varies as a function of temperature for materials and directing ultrasonic wave propagation through selected material volumes of interest. In this project, longitudinal wave speed is estimated for aluminum 6061. Piezoelectric materials bonded to a specimen provided data from which longitudinal wave speed was estimated. A discussion evaluated the linear relationship of wave speed versus temperature and postulated how wave speed measurements could be used to estimate the internal temperature of selected material volumes for feedback during SPP of materials. Instrumentation, data and analysis, and a discussion are reviewed.

36 MATERIALS SCIENCE↗

New capabilities for rapid depletion analysis of pebble-bed reactors in SCALE

The SCALE Leap-In for Cores at Equilibrium (SLICE) method leverages capabilities available in the SCALE modeling and simulation suite to facilitate rapid estimation of equilibrium core inventories for flowing-pebble reactor systems in support of a variety of fuel cycle applications. New capabilities that compliment the SLICE method have been developed for the ORIGAMI interface to ORIGEN in SCALE to facilitate rapid depletion calculations for flowing-pebble systems. New features include a more generalized and flexible means of specifying interpolation dimensions, the ability to evaluate pebbles moving through user-defined “axial zones” in each of the pebble’s multiple passes through the core, and the treatment of differential velocities across radial channels (representing both pebble-to-pebble and wall-to-pebble friction effects). ORIGAMI thus provides an efficient user interface to define a pebble’s path and irradiation histories of its multiple passes through the core to calculate the pebble’s time-dependent inventories, which can be useful in various applications domains such as safeguards, criticality safety, and disposal analyses. The consistency of new ORIGAMI approach was verified against the solution obtained using the SLICE method with ORIGEN-ARP, showing excellent agreement.

Skutnik, Steven E. [Oak Ridge National Laboratory ↗

Joule heating in bad and slow metals

Heat supplied to a metal is absorbed by the electrons and then transferred to the lattice. In conventional metals energy is released to the lattice by phonons emitted from the Lindhard continuum. However in a 'bad' metal, with short mean free path, the low energy Lindhard continuum is destroyed. To describe energy transfer to the lattice in these cases we obtain a general Kubo formula for the energy relaxation rate in terms of the electronic density spectral weight \text{lm} \, G^R_{nn}(\omega_{k},k) lm G n n R ( ω k , k ) evaluated on the phonon dispersion \omega_k ω k . We apply our Kubo formula to the high temperature Hubbard model, using recent data from quantum Monte Carlo and experiments in ultracold atoms to characterize \text{lm} \, G^R_{nn}(\omega_{k},k) lm G n n R ( ω k , k ) . We furthermore use recent data from electron energy-loss spectroscopy to estimate the energy relaxation rate of the cuprate strange metal to a high energy optical phonon. As a second, distinct, application of our formalism we consider 'slow' metals. These are defined to have Fermi velocity less than the sound velocity, so that particle-hole pairs are kinematically unable to emit phonons. We obtain an expression for the energy relaxation rate of a slow metal in terms of the optical conductivity.

36 MATERIALS SCIENCE↗

A numerical investigation of coal particle modeling for the inlet section of a pressurized oxy-combustion burner

Concerns over climate change have led to numerous efforts to develop low-carbon technologies, and pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in the power industry. Designing an effective burner for such technologies is vital for the development of these and other coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the pressurized oxy-combustor burner is close to or even higher than the maximum limit for commercial CFD codes (e.g., ANSYS FLUENT). At these high particle volume fractions, the interactions among particles, fluid flow and wall need to be re-evaluated. The present computational work is a first step in a systematic analysis of the influence of various parameters, like method of particle release, release location, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In the POC process, pulverized coal is burned under elevated pressure in an O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing Ansys FLUENT commercial sodtware using the Reynolds-Averaged Navier-Stokes (RANS) approach. It is revealed that for this pilot-scale, pressurized burner, velocity profiles in the near-wall region exhibit some anomalies. In order to investigate the influence of particle loading in the near-wall region, particle release location was investigated. The numerical simulations incorporate the coupling between the turbulent flow and the particles. A sensitivity investigation of particle release location found that by tuning the release location, the velocity profile can be consistent with the pure gas flow velocity profile. More importantly, the particle releasing location also affects flame stability. Particle size was also found to have a significant impact on particle trajectory, flame stability and temperature. Finally, Large Eddy Simulations (LES) were performed and compared with the results from two-dimensional RANS.

Li, Lei↗

Integrating Experiments and Well Logs to Predict Caney Shale Static Mechanical Properties During Production with Supervised Machine Learning

Caney shale is one of the emerging oil reservoirs in Oklahoma. Understanding the impact of effective stress on its mechanical properties is critical for predicting hydraulic fracture geometry and overall hydrocarbon production. The objective of our study is to evaluate the impact of effective stress on the dynamic Young’s modulus using ultrasonic velocity measurements for Caney shale samples. A triaxial cell was utilized to measure ultrasonic (P-wave and S-wave) velocities for ten downhole Caney shale samples under various effective stresses to indirectly assess the impact of pore pressure change. The dynamic Young’s moduli estimated from these measurements were integrated with available conventional well logs (excluding sonic logs) and triaxial test results from Benge et al. (2021) to predict the static Young’s modulus using Random Forest (RF) and Extreme Gradient Boosting (XGBoost) models. The results showed that the estimated dynamic Young’s moduli from ultrasonic measurements were higher than the corresponding static Young’s modulus of cores from the same vertical well at similar depths. With increasing effective stress, the dynamic Young’s modulus increased for all samples. The estimated dynamic-to-static correction factor tended to be higher in zones of high neutron porosity (PHIN) and low density compared to other zones. Finally, SHapley Additive exPlanations (SHAP) for RF and XGBoost models identified depth, gamma ray (GR), and PHIN as key features for predicting the static Young’s modulus. This study enhances our understanding of the dynamic and static Young’s moduli for the Caney shale interval, as a function of effective stress and conventional well logs. The findings from this study can improve predictions of production throughout the well's lifespan by offering insights into the mechanical property degradation resulting from pore pressure depletion.

Kholy, Sherif M.↗

Performance and Wake Characterization of a Model Hydrokinetic Turbine: The Reference Model 1 (RM1) Dual Rotor Tidal Energy Converter

The mechanical power and wake flow field of a 1:40 scale model of the US Department of Energy’s Reference Model 1 (RM1) dual rotor tidal energy converter are characterized in an open-channel flume to evaluate power performance and wake flow recovery. The NACA-63(4)-24 hydrofoil profile in the original RM1 design is replaced with a NACA-4415 profile to minimize the Reynolds dependency of lift and drag characteristics at the test chord Reynolds number. Precise blade angular position and torque measurements were synchronized with three acoustic Doppler velocimeters (ADV) aligned with each rotor centerline and the midpoint between the rotor axes. Flow conditions for each case were controlled to maintain a hub height velocity, uhub= 1.04 ms−1, a flow Reynolds number, ReD= 4.4 × 105, and a blade chord length Reynolds number, Rec= 3.1 × 105. Performance was measured for a range of tip-speed ratios by varying rotor angular velocity. Peak power coefficients, CP= 0.48 (right rotor) and CP= 0.43 (left rotor), were observed at a tip speed ratio, λ= 5.1. Vertical velocity profiles collected in the wake of each rotor between 1 and 10 rotor diameters are used to estimate the turbulent flow recovery in the wake, as well as the interaction of the counter-rotating rotor wakes. The observed performance characteristics of the dual rotor configuration in the present study are found to be similar to those for single rotor investigations in other studies. Similarities between dual and single rotor far-wake characteristics are also observed.

30 DIRECT ENERGY CONVERSION↗

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗

Understanding aerosol–cloud interactions using a single-column model for a cold-air outbreak case during the ACTIVATE campaign

Marine boundary layer clouds play a critical role in Earth's energy balance. Their microphysical and radiative properties are highly impacted by ambient aerosols and dynamic forcings. In this study, we evaluate the representation of these clouds and related aerosol–cloud interaction processes in the single-column version of the E3SM climate model (SCM) against field measurements collected during the NASA ACTIVATE campaign over the western North Atlantic, as well as intercompare results with high-resolution process level models. We show that E3SM SCM reproduces the macrophysical properties of post-frontal boundary layer clouds in a cold-air outbreak (CAO) case well. However, it generates fewer but larger cloud droplets compared to aircraft measurements. Further sensitivity tests show that the underestimation of both aerosol number concentration and vertical velocity variance contributes to this bias. Aerosol–cloud interactions are examined by perturbing prescribed aerosol properties in E3SM SCM with fixed dynamics. Higher aerosol number concentration or hygroscopicity leads to more numerous but smaller cloud droplets, resulting in a stronger cooling via shortwave cloud forcing. This apparent Twomey effect is consistent with prior climate model studies. The cloud liquid water path shows a weakly positive relation with cloud droplet number concentration due to precipitation suppression. This weak aerosol effect on cloud macrophysics may be attributed to the dominant impact of strong dynamical forcing associated with the CAO. Our findings indicate that the SCM framework is a key tool to bridge the gap between climate models, process level models, and field observations to facilitate process level understanding.

54 ENVIRONMENTAL SCIENCES↗

LARGE PARTICLE SEPARATION FROM NON-NEWTONIAN SLURRIES USING BUMP ARRAYS

Here we evaluate the performance of bump arrays to separate large particles from non-Newtonian slurries with Bingham and Cross rheology. Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit. Although these devices have been used to separate a wide variety of large particles from blood cells to sand, partition of large particles from non-Newtonian fluids remains unexplored. Yet, an important set of modestly concentrated slurries, including Hanford nuclear waste, displays non-Newtonian rheology. Here we evaluate the influence of non-Newtonian rheology on the large-small particle size cutoff in bump arrays using a model that explores the influence of yield stresses, ratios of zero and infinite shear viscosities, and Cross’s exponent under strictly laminar well-developed conditions. Surprisingly, we find that viscosity ratios and Cross’s exponent make no significant difference on the particle cutoffs between large particles that bump and small particles that slalom around the posts from entrance to exit. In contrast, we find that yield stresses do significantly affect the size cutoff. As the yield stress increases, the velocity profiles become more plug like lowering the size cutoff. For nuclear waste separations where removing large particles is a priority, increasing yield stresses is conservative.

slurry, solids removal, deterministic later displa↗

Turbulence in The Marine Boundary Layer and Air Motions Below Stratocumulus Clouds at the ARM Eastern North Atlantic Site

Marine stratocumulus clouds are intimately coupled to the turbulence in the boundary layer and drizzle is known to be ubiquitous within them. Six years of data collected at the Atmospheric Radiation Measurement (ARM)’s Eastern North Atlantic site are utilized to characterize turbulence in the marine boundary layer and air motions below stratocumulus clouds. Profiles of variance of vertical velocity binned by wind direction (wdir) yielded that the boundary layer measurements are affected by the island when the wdir is between 90° and 310° (measured clockwise from North where air is coming from). Data collected during the marine conditions (wdir<90 or wdir>310) showed that the variance of vertical velocity was higher during the winter months than during the summer months due to higher cloudiness, wind speeds, and surface fluxes. Furthermore, during marine conditions the variance of vertical velocity and cloud fraction exhibited a distinct diurnal cycle with higher values during the nighttime than during the daytime. Detailed analysis of 32 cases of drizzling marine stratocumulus clouds showed that for a similar amount of radiative cooling at the cloud top, within the sub-cloud layer 1) drizzle increasingly falls within downdrafts with increasing rain rates, 2) the strength of the downdrafts increases with increasing rain rates, and 3) the correlation between vertical air motion and rain rate is highest in the middle of the sub-cloud layer. The results presented herein have implications for climatological and model evaluation studies conducted at the ENA site, along with efforts of accurately representing drizzle-turbulence interactions in a range of atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Tank 49H solids disturbance analysis

Tank 49H serves as the feed tank for the Salt Waste Processing Facility (SWPF). Transfers into Tank 49H may disturb any solids that have settled to the bottom of the tank, resulting in feed that may exceed the insoluble solids content limit of 1,200 mg/L of the SWPF Waste Acceptance Criteria (WAC). During a transfer into Tank 49H, material that free falls into Tank 49H through a downcomer could potentially disturb any solids on the bottom of the tank and scour or suspend solids from a settled solids layer or turbid region. A previous analysis and report evaluated the potential to disturb solids when transferring into Tank 49H through the B4 riser and recommended a minimum tank level of 120 inches to prevent disturbing any solids in the bottom of the tank. The scope of this task is to perform additional fluid flow analysis to determine whether accounting for disturbed particle settling and particle mixing and dispersion during transfer could allow the minimum liquid level to be reduced below 120 inches. The analysis utilized models from the technical literature to calculate the size and shape of the “plunging jet” as a function of input parameters such as initial velocity, initial jet diameter, elevation of the initial jet, liquid level in the tank, and solids depth. The analysis relied on the M-Star® simulations performed for the previous analysis to provide bounding estimates of the amount of solid particles disturbed and used the MStar® software to calculate the dispersion and mixing of the disturbed solids with other liquid in the tank as the solids are transported to the transfer pump. The analysis showed that with a solid particle size of 5 micron or less, a liquid level of 120 inches should be maintained to prevent significant disturbance of the solid layer at the bottom of Tank 49H.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗

Denoising Seismic Waveforms Using a Wavelet-Transform-Based Machine-Learning Method

Seismic waveform data recorded at stations can be thought of as a superposition of the signal from a source of interest and noise from other sources. Frequency‐based filtering methods for waveform denoising do not result in desired outcomes when the targeted signal and noise occupy similar frequency bands. Recently, denoising techniques based on deep‐learning convolutional neural networks (CNNs), in which a recorded waveform is decomposed into signal and noise components, have led to improved results. These CNN methods, which use short‐time Fourier transform representations of the time series, provide signal and noise masks for the input waveform. These masks are used to create denoised signal and designaled noise waveforms, respectively. However, advancements in the field of image denoising have shown the benefits of incorporating discrete wavelet transforms (DWTs) into CNN architectures to create multilevel wavelet CNN (MWCNN) models. The MWCNN model preserves the details of the input due to the good time–frequency localization of the DWT. In this report we use a data set of over 382,000 constructed seismograms recorded by the University of Utah Seismograph Stations network to compare the performance of CNN and MWCNN‐based denoising models. Evaluation of both models on constructed test data shows that the MWCNN model outperforms the CNN model in the ability to recover the ground‐truth signal component in terms of both waveform similarity and preservation of amplitude information. Model evaluation of real‐world data shows that both the CNN and MWCNN models outperform standard band‐pass filtering (BPF; average improvement in signal‐to‐noise ratio of 9.6 and 19.7 dB, respectively, with respect to BPF). Evaluation of continuous data suggests the MWCNN denoiser can improve both signal detection capabilities and phase arrival time estimates.

58 GEOSCIENCES↗

Enhancing reactant selectivity for Ni/Mg reforming catalysts using silicalite-1 shells: A modeling study

The ability of a zeolite shell to enhance the selective conversion of hydrocarbons through diffusional limitations was investigated using a multi-scale model of a fixed-bed reactor. The impact of shell thickness and molecule/pore size on the catalytic performance of silicalite-1 zeolite encapsulated nickel catalyst pellets for steam reforming of C1-C7 hydrocarbons is reported. A reaction–diffusion model using kinetic expressions established in literature was employed. The model was verified through comparison with reported experimental results for steam reforming data over a temperature range of 748 – 1113 K and pressure of 1 – 10 bar. Comparisons are also made against experimental data for steam reforming in the presence of a zeolite shell. Evaluation of the Weisz-Prater criterion for both the core and encapsulated catalyst confirmed mass transfer limitation induced by the utilization of a zeolite shell. The model was used to suggest an optimal thickness that balances diffusional limitations imposed by the zeolite layer on methane versus that of the heavier hydrocarbons. The optimum thickness varied as a function of hydrocarbon size and shape which determined the diffusion rates. For toluene and heptane, a 50 nm thick shell was sufficient to wholly prevent reaction. Hydrocarbons like propane and butane required a shell 7.5 and 5 μm thick. Increasing the gas-hourly-space-velocity from 10,000 to 60,000 h –1 caused a decrease in the optimum shell thickness. Furthermore, this approach can be modified for application to other mixed hydrocarbon systems to predict optimal catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management facility (EMF) stack monitor location was qualified using the LV-S1 scale model stack as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale models and its sampling locations were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the EMF stack was performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on the LV-S1 scale model test DV values, the corresponding stack flow rates for the EMF stack are as listed in Table S1. Table S1. Effluent Management Facility Stack Qualified Flow Range. Stack Parameter EM-1 Minimum Qualified Stack Flow (scfm) 781 Maximum Qualified Stack Flow (scfm) 53,432 The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were less than 5°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, which were all within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes each combination of ACV fan with the DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Development and Evaluation of Velocity Predictive Optimal Energy Management Strategies in Intelligent and Connected Hybrid Electric Vehicles

In this study, a thorough and definitive evaluation of Predictive Optimal Energy Management Strategy (POEMS) applications in connected vehicles using 10 to 20 s predicted velocity is conducted for a Hybrid Electric Vehicle (HEV). The presented methodology includes synchronous datasets gathered in Fort Collins, Colorado using a test vehicle equipped with sensors to measure ego vehicle position and motion and that of surrounding objects as well as receive Infrastructure to Vehicle (I2V) information. These datasets are utilized to compare the effect of different signal categories on prediction fidelity for different prediction horizons within a POEMS framework. Multiple artificial intelligence (AI) and machine learning (ML) algorithms use the collected data to output future vehicle velocity prediction models. The effects of different combinations of signals and different models on prediction fidelity in various prediction windows are explored. All of these combinations are ultimately addressed where the rubber meets the road: fuel economy (FE) enabled from POEMS. FE optimization is performed using Model Predictive Control (MPC) with a Dynamic Programming (DP) optimizer. FE improvements from MPC control at various prediction time horizons are compared to that of full-cycle DP. All FE results are determined using high-fidelity simulations of an Autonomie 2010 Toyota Prius model. The full-cycle DP POEMS provides the theoretical upper limit on fuel economy (FE) improvement achievable with POEMS but is not currently practical for real-world implementation. Perfect prediction MPC (PP-MPC) represents the upper limit of FE improvement practically achievable with POEMS. Real-Prediction MPC (RP-MPC) can provide nearly equivalent FE improvement when used with high-fidelity predictions. Constant-Velocity MPC (CV-MPC) uses a constant speed prediction and serves as a “null” POEMS. Results showed that RP-MPC, enabled by high-fidelity ego future speed prediction, led to significant FE improvement over baseline nearly matching that of PP-MPC.

33 ADVANCED PROPULSION SYSTEMS↗