Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Layer Profiles”

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

Bridging New Observational Capabilities and Process-Level Simulation: Insights into Aerosol Roles in the Earth System

The spatial distribution of ambient aerosol particles significantly impacts aerosol–radiation–cloud interactions, which contribute to the largest uncertainty in global anthropogenic radiative forcing estimations. However, the atmospheric boundary layer and lower free troposphere have not been adequately sampled in terms of spatiotemporal resolution, hindering a comprehensive characterization of various atmospheric processes and impeding our understanding of the Earth system. To address this research data gap, we have leveraged the development of uncrewed aerial systems (UAS) and advanced measurement techniques to obtain mesoscale spatial data on aerosol microphysical and optical properties around the U.S. Southern Great Plains (SGP) atmospheric observatory. Our study also benefits from state-of-the-art laboratory facilities that include three-dimensional molecular imaging techniques enabled by secondary ion mass spectrometry and nanogram-level chemical composition analysis via micronebulization aerosol mass spectrometry. Through our study, we have developed a framework for observation–modeling integration, enabling an examination of how various assumptions about the organic–inorganic components mixing state, inferred from chemical analysis, affect clouds and radiation in observation-constrained model simulations. By integrating observational constraints (derived from offline chemical analysis of the aerosol surface using collected samples) with in situ UAS observations, we have identified a prominent role of organic-enriched nanometer layers located at the surface of aerosol particles in determining profiles of aerosol optical and hygroscopic properties over the SGP observatory. Furthermore, we have improved the agreement between predicted clouds and ground-based cloud lidar measurements. This UAS–model–laboratory integration exemplifies how these new advanced capabilities can significantly enhance our understanding of aerosol–radiation–cloud interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Realization of a gas puff imaging system on the Wendelstein 7-X stellarator

A system for studying the spatiotemporal dynamics of fluctuations in the boundary of the W7-X plasma using the “Gas-Puff Imaging” (GPI) technique has been designed, constructed, installed, and operated. This GPI system addresses a number of challenges specific to long-pulse superconducting devices, such as W7-X, including the long distance between the plasma and the vacuum vessel wall, the long distance between the plasma and diagnostic ports, the range of last closed flux surface (LCFS) locations for different magnetic configurations in W7-X, and management of heat loads on the system’s plasma-facing components. The system features a pair of “converging–diverging” nozzles for partially collimating the gas puffed locally ≈135 mm radially outboard of the plasma boundary, a pop-up turning mirror for viewing the gas puff emission from the side (which also acts as a shutter for the re-entrant vacuum window), and a high-throughput optical system that collects visible emission resulting from the interaction between the puffed gas and the plasma and directs it along a water-cooled re-entrant tube directly onto the 8 × 16 pixel detector array of the fast camera. The DEGAS 2 neutral code was used to simulate the Hα (656 nm) and HeI (587 nm) line emission expected from well-characterized gas-puffs of H2 and He and excited within typical edge plasma profiles in W7-X, thereby predicting line brightnesses used to reduce the risks associated with system sensitivity and placement of the field of view. Operation of GPI on W7-X shows excellent signal-to-noise ratios (>100 at 2 Mframes/s) over the field of view for minimally perturbing gas puffs. The GPI system provides detailed measurements of the two-dimensional (radial and poloidal) dynamics of plasma fluctuations in the W7-X edge and scrape-off layer and in and around the magnetic islands outside the LCFS that make up the island divertor configuration employed on W7-X.

Instruments & Instrumentation↗
DOE OSTI2020

Exploration of Sn 70 Li 30 alloy as possible material for flowing liquid metal plasma facing components

As an advanced alternative to solid materials, Liquid Metals (LM) may offer more resilient and feasible Plasma Facing Components (PFCs). Particularly, regarding the unavoidable material erosion/degradation produced by particle/heat fluxes in future fusion devices where much longer duty cycles are expected. Furthermore, configurations that propose a flowing LM surface can add the advantage of a continuously fresh and clean layer facing the plasma. Although lithium is the most widely tested option, tin-lithium (SnLi) alloys have been proposed to attempt to combine the positive characteristics of both pure elements and ameliorate the specific issues of lithium. In this work, the potential use of Sn70Li30 alloy in such flowing concepts has been explored by addressing several preliminary and mandatory aspects for its utilization. Key issues such as wettability and compatibility of the alloy with relevant substrates have been studied in a multidisciplinary approach. The data obtained from deposited liquid tin-lithium droplets indicates approximate wetting temperatures of 360 °C, 390 °C and 405 °C for the fresh alloy on 316 stainless steel, molybdenum, and tungsten, respectively. However, the alloy contamination appeared to strongly affect the wetting characteristics of materials, increasing their wetting temperature by ~130 °C in the worst observed cases. Interestingly, in some instances, the instability of the liquid alloy surface was observed in the form of sudden gaseous ejection. The deposited droplets were posteriorly characterized in terms of absolute composition and depth profile by Inductively Coupled PlasmaOptical Emission Spectroscopy (ICP-OES) and Secondary Ion Mass Spectrometry (SIMS-ToF). Additionally, the nature and composition of the boundaries between the substrates and alloy microparticles was investigated by Scanning Electron Microscopy (SEM), Energy Dispersive X-Ray Spectroscopy (EDS), and 3D Laser microscopy. The overall results of this post-mortem characterization revealed that first signs of corrosion induced by both alloy elements (lithium-chromium association and iron-tin intermetallic mixing) were present on 316 stainless steel after short exposures (≤3 h) at temperatures lower than 550 °C. Conversely, molybdenum and tungsten showed good compatibility with the alloy in equivalent conditions. The global implications of these results are finally addressed, focusing on the future perspectives and the more viable scenarios for the eventual utilization of these alloys in flowing liquid metal configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Effect of the integral length scales of turbulent inflows on wind turbine loads

As wind turbines become larger, the fluctuations in the inflow become increasingly influential in the turbine structural loading. These fluctuations are characterized by the integral length scale, which represents the average size of the largest energy-containing turbulent eddies. Current design standards neglect the varying integral length scales that characterize inflows of wind turbines in operation. Using large-eddy simulations, we generate turbulent inflows of varying integral length scales and quantify the loads of the IEA 15-MW reference wind turbine. Results illustrate that the impact of turbulence on rotor and tower loads is up to 10 times greater than the impact of the mean shear profile. Increasing integral length scales from 0.3x to 0.5x the rotor diameter reduces blade root flapwise moments and rotor and tower loads. Increasing integral length scales from 0.5x to 1.7x the rotor diameter increases the rotor aerodynamic thrust force and the blade root flapwise shear loads and decreases the tower base torsional moment and the tilting and yawing rotor aerodynamic moments. Additionally, turbulence intensity has a greater impact on wind turbine loads than integral length scales. In conclusion, findings indicate that design standards should consider varying integral length scales for accurate wind turbine loading characterization in turbulent inflow conditions.

17 WIND ENERGY↗
DOE OSTI2024

PSU-BSEC Doppler Lidar Processed Scans

This data repository includes the processed vertically staring measurements (stare files) and the angled scans (profile files) necessary to calculate horizontal winds retrieved by the Pennsylvania State University (PSU) Doppler Lidar as a part of the Baltimore Social-Environmental Collaborative Urban Integrated Field Laboratory (BSEC UIFL). The PSU Doppler Lidar measures aerosol backscatter intensity (m-1 sr-1), signal-to-noise ratio, and radial velocity (i.e., vertical velocity in the case of stare files, units: m/s) at approximately 1 Hz temporal resolution and 30 m spatial resolution. Within the "stare_data" subdirectory there exist example figures of all retrieved data. For further information, please email Nicholas Prince, nec5299@psu.edu.

Air Quality↗
DOE OSTI2023

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗
DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Offshore low-level jet observations and model representation using lidar buoy data off the California coast

Abstract. Low-level jets (LLJs) occur under a variety of atmospheric conditions and influence the available wind resource for wind energy projects. In 2020, lidar-mounted buoys owned by the US Department of Energy (DOE) were deployed off the California coast in two wind energy lease areas administered by the Bureau of Ocean Energy Management: Humboldt and Morro Bay. The wind profile observations from the lidars and collocated near-surface meteorological stations (4–240 m) provide valuable year-long analyses of offshore LLJ characteristics at heights relevant to wind turbines. At Humboldt, LLJs were associated with flow reversals and north-northeasterly winds, directions that are more aligned with terrain influences than the predominant northerly flow. At Morro Bay, coastal LLJs were observed primarily during northerly flow as opposed to the predominant north-northwesterly flow. LLJs were observed more frequently in colder seasons within the lowest 250 m a.s.l. (above sea level), in contrast with the summertime occurrence of the higher-altitude California coastal jet influenced by the North Pacific High, which typically occurs at heights of 300–400 m. The lidar buoy observations also validate LLJ representation in atmospheric models that estimate potential energy yield of offshore wind farms. The European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) was unsuccessful at identifying all observed LLJs at both buoy locations within the lowest 200 m. An extension of the National Renewable Energy Laboratory (NREL) 20-year wind resource dataset for the Outer Continental Shelf off the coast of California (CA20-Ext) yielded marginally greater captures of observed LLJs using the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary layer (PBL) scheme than the 2023 National Offshore Wind dataset (NOW-23), which uses the Yonsei University (YSU) scheme. However, CA20-Ext also produced the most LLJ false alarms, which are instances when a model identified an LLJ but no LLJ was observed. CA20-Ext and NOW-23 exhibited a tendency to overestimate the duration of LLJ events and underestimate LLJ core heights.

17 WIND ENERGY↗
DOE OSTI2022

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

LES simulations of a vacuum membrane distillation channel with geometric alterations

3D LES simulations were carried out to study the performance of a vacuum membrane distillation module. A wiggly wall profile with/without embedded stiffeners was considered to alleviate axial and radial temperature polarization, the cause of performance decrease in membrane distillation. Results of the flow field show the wiggles create unsteady vortex shedding inducing intense mixing in the channel. Vortex shedding intensity increases as the Reynolds number increases or stiffeners are added. 3D results show a 56% improvement in flux when moving from a flat sheet membrane to a wiggly membrane with stiffeners at a constant Reynolds number, corresponding to an improvement from 11.7 to 41.7 in the Nusselt number, showing that the alterations improved the flux performance by enhancing the heat transfer along the membrane surface and therefore alleviating temperature polarization, a critical bottleneck in membrane distillation systems. Further, a merit criterion was defined based on the Nusselt number and friction factor, and a 40% increase in merit was shown switching from a flat to a wiggly channel, while a 97% merit increase was seen going from flat to wiggly with stiffeners. A 42% enhancement in the flux was also seen moving from a straight channel to a wiggly channel at a higher Reynolds number which highlights the importance of strategically choosing the mass flow rate, as well as inducing flow separation and vortex shedding in the channel to promote mixing and dissipate the thermal boundary layer. The variation in results between the wiggly module with stiffeners and the wiggly module at high Reynolds numbers suggests enhancing mixing structures can be more impactful on flux for this geometry than increasing the flow rate to a turbulent/transitional regime. However, both are preferable for peak system flux performance. Furthermore, a 2D approximation was used to perform simulations on more extended channels to examine the length degradation. The modules with wiggly channels performed at the same flux level with a doubling of the length. On the other hand, the flat sheet modules experienced length degradation by temperature polarization and dropped in flux yield by around 13% of the short-channel value. This work illustrates that modeling the system and understanding how the performance decreases as the membrane surface area increases are critical for a larger module (scaling up from a lab to a prototype module) to maintain high flux performance.

42 ENGINEERING↗
DOE OSTI2023

Near-Cloud Atmospheric Ingredients for Deep Convection Initiation

A lack of routine environmental observations located near deepening cumulus congestus clouds limits verification of important theorized and simulated updraft–environment interaction processes occurring during deep convection initiation (CI). We analyze radiosonde profiles collected during several hundred CI events near a mountain range in central Argentina during the CACTI field campaign. Statistical analyses illustrate environmental conditions supporting radar-observed CI outcomes that span a spectrum of convective cell depths, widths, and durations, as well as events lacking precipitating convection. Tested environmental factors include a large variety of sounding-derived measurements of CAPE, CIN, moisture, terrain-relative winds, vertical shear, and lifted parcel properties, with supplemental model reanalysis of background larger-scale vertical motion. CAPE and CIN metrics do not consistently differentiate CI success from failure. Only a few environmental factors contain consistent monotonic relationships among the spectrum of cloud depths achieved during CI: (i) the depth and strength of background ascent, and (ii) the component of low-level flow oriented parallel to the ridgeline. These metrics suggest that the ability of the surrounding flow to lift parcels to their LFC and terrain-modified flow are consistently relevant processes for CI. Low- to midlevel relative humidity strongly discriminated between CI and non-CI events, likely reflecting entrainment-driven dilution processes. However, we could not confidently conclude that relative humidity similarly discriminated robust from marginal CI events. Circumstantial evidence was found linking cell width, an important cloud property governing the probability of CI, to LCL height, boundary layer depth, depth and magnitude of the CIN layer, and ambient wind shear.

Storm environments↗
DOE OSTI2024

Low-Level Liquid-Bearing Clouds Contribute to Seasonal Lower Atmosphere Stability and Surface Energy Forcing over a High-Mountain Watershed Environment

Measurements of atmospheric structure and surface energy budgets distributed along a high-altitude moun tain watershed environment near Crested Butte, Colorado, from two separate, but coordinated, field campaigns, Surface Atmosphere Integrated field Laboratory (SAIL) and Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH), are analyzed. This study identifies similarities and differences in how clouds influence the radiative budget over one snow-free summer season (2022) and two snow-covered seasons (2021/22; 2022/23) for this alpine location. A relationship between lower-tropospheric stability stratification and longwave radiative flux from the presence or absence of clouds is identified. When low clouds persisted, often with signatures of supercooled liquid in win ter, the lower troposphere experienced weaker stability, while radiatively clear skies that are less likely to be influenced by liquid droplets were associated with appreciably stronger lower-tropospheric stratification. Corresponding surface turbu lent heat fluxes partitioned differently based upon the cloud–stability stratification regime derived from early morning radiosounding profiles. Combined with the differences in the radiative budget largely resulting from dramatic seasonal dif ferences in surface albedo, the lower atmosphere stratification, surface energy budget, and near-surface thermodynamics are shown to be modified by the effective longwave radiative forcing of clouds. The diurnal evolution of thermodynamics and surface energy components varied depending on the early morning stratification state. Thus, the importance of quies cent versus synoptically active large-scale meteorology is hypothesized as a critical forcing for cloud properties and associ ated surface energy budget variations. The physical relationships between clouds, radiation, and stratification can provide a useful suite of metrics for process understanding and to evaluate numerical models in such an undersampled, highly com plex terrain environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

High-Resolution Wind Resource Data Set of the Greater Puerto Rico Region

In February 2022, the U.S. Department of Energy and six national laboratories launched the Puerto Rico Grid Resilience and Transitions to 100% Renewable Energy Study (PR100). PR100 aims to provide a comprehensive analysis of possible pathways for Puerto Rico's energy future, with a goal of 100% renewable energy by 2050. As a part of the renewable energy potential assessment in this project, we developed 20 years (2001-2020) of data using a numerical weather prediction (NWP) model for onshore and offshore wind resource assessment for the Puerto Rico. The research steps in developing the long-term wind resource data sets based on the NWP model were: 1. Model wind resource based on the Weather Research and Forecasting (WRF) model. 2. Develop WRF model configurations for Puerto Rico. 3. Test WRF with 11 different physics parameterizations for planetary boundary layer (PBL). 4. Assess WRF output from the different PBL schemes against observations. 5. Select a final model configuration which can produce the modeled wind speed with sufficient accuracy. 6. Produce 20 years wind resource data sets for Puerto Rico region. In the first stage of our framework for developing wind resource data, we developed a WRF model configuration using two nested domains (9 km and 3 km) to cover Puerto Rico and U.S. Virgin Islands and downscale the ERA5 reanalysis data (0.25 degrees x 0.25 degrees; hourly interval) to a 3-km domain. For the second stage, we implemented one-year simulations focused on using 11 different PBL physics parameterizations to find a combination of WRF physics parameterizations that could provide accurately modeled wind speed for Puerto Rico. We also analyzed the sensitivity of the modeled wind speed to PBL schemes for onshore and offshore locations. The WRF output resulting from the 11 WRF experiments using different PBL parameterizations were evaluated against observations obtained from the National Data Buoy Center (NDBC) as well as at hub height for a location for which measurements were available. A final model setup selected through the validation with observational data was used to produce 20 years of data with 3-km spatial and 5-minute temporal resolution. The WRF model output was post-processed to include wind profiles and basic atmospheric variables in a format that can be easily used for downstream modeling. The 20 years of wind resource data will be made available through NREL and support the estimation of wind energy development costs for the PR100 study.

17 WIND ENERGY↗
DOE OSTI2022

Limitations of reanalysis data for wind power applications

Wind energy resource estimates commonly depend on simulated wind speed profiles generated by reanalysis or weather models due to the lack of long time series measurements with sufficient coverage at relevant heights (roughly 90 m above ground). However, modeled data, including reanalyses, can be noisy and display a wide range of biases and errors, variously attributed to terrain effects, poor coverage of assimilated inputs, and model resolution. Wind generation records, if available at high temporal and geographical resolution, can provide a proxy for wind measurements and allow for evaluation of reanalyses and weather model wind time series. We use a 7-year-long data set of hourly, plant-level generation records from over 100 wind plants across Texas to evaluate two commonly used reanalysis data sets (MERRA2 and ERA5). Additionally, we use 1-year of records (2019) to evaluate an operational, high-resolution regional weather modeling product (HRRR v3). We find that across the region, and across all modeling products, the modeled representation of wind generation (i.e., wind speeds at hub heights passed through a power curve) has relatively small mean errors when aggregated daily, but that accuracy and hourly correlation have a strong diurnal sensitivity. Accuracy and correlation systematically decline through the evening and markedly improve after sunrise. These diurnal patterns persist even in the highest resolution model tested (HRRR v3). We hypothesize the nighttime decline in accuracy is mostly due to poorly represented boundary layer conditions, perhaps related to model representation of stability, while other uncertainties (such as wake effects) play a secondary role.

17 WIND ENERGY↗
DOE OSTI2025

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Texas A&M University Mobile Facility Measurements during TRACER (Field Campaign Report)

One of the main goals of the U.S Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Tracking Aerosol Convection Interactions Experiment (TRACER) near Houston, Texas is to improve understanding of how meteorology and aerosols impact storm dynamical and microphysical processes in deep convection to better constrain and improve their model representation. The Houston area is strongly influenced by sea- and bay-breeze circulations that generate convergence and help to initiate and organize deep convection. To properly isolate and understand the roles of varied meteorological conditions and cloud condensation nuclei (CCN) and ice nucleation particles (INP) distributions in different air masses, co-located thermodynamic, kinematic, and aerosol vertical profile observations are needed. The focus of this campaign was to provide key measurements in air masses both in front of and behind sea/bay breeze fronts moving through the greater Houston area to sample the airmass heterogeneity. The overarching scientific goal of this campaign is to understand how the vertical distributions of both CCN and INP correspond to the inflow layer of deep convection in maritime, background continental, and polluted continental air masses, and how these variations influence deep convection. To fully sample the heterogeneity in both meteorological conditions and aerosols across the sea-breeze front (SBF), Texas A&M University (TAMU) deployed a InterMet 3050A 403 MHz mobile unit launching iMet-4 radiosondes and the new Rapid Onsite Atmospheric Measurement Van (ROAM-V) for aerosol sampling during the TRACER intensive operational period (IOP) from June to September 2022. The suite of instruments deployed on ROAM-V included a condensation particle counter (CPC; GRIMM Model 5.403 CPC), scanning mobility particle sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A differential mobility analyzer), cloud condensation nuclei counter (Droplet Measurement Technologies CCN counter), micropulse lidar (Droplet Measurement Technologies micropulse lidar [miniMPL]), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the Global Positioning System (GPS) on the phone application “My Altitude”. The DRUM data were collected as part of a closely related ARM field campaign and also supported by DOE Atmospheric System Research grant DE-SC0021047. The TAMU team sampled these airmass heterogeneities by strategically choosing deployment sites in a different airmass than the ARM fixed sites in La Porte and Guy, Texas. On days when the sea/bay breeze boundary was pushing inland, the TAMU team would usually sample the airmass on the maritime side of the SBF at a coastal site in Galveston, Texas in the early afternoon (1730-1900 UTC) and then move inland ahead of the SBF to sample the airmass on the continental side during late afternoon (2030-2230). Figure 1 shows the Galveston maritime site and the array of sites for the late afternoon continental measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Revisiting the validity of eddy viscosity models for predicting airflow over water waves

In this study, we revisit the validity of eddy viscosity models for predicting wave-induced airflow disturbances over ocean surface waves. We first derive a turbulence curvilinear model for the phase-averaged Navier–Stokes equations, extending the work of Cao, Deng & Shen (2020 J. Fluid Mech. 901, A27), by incorporating turbulence stress terms previously neglected in the linearised viscous curvilinear model. To verify our formulation, we perform a priori tests by numerically solving the model using mean wind and turbulence stress profiles from large-eddy simulations (LES) of airflow over waves across various wave ages. Results show that including turbulence stress terms improves wave-induced airflow predictions compared with the previous viscous curvilinear model. We further show that using a standard mixing-length eddy viscosity yields inaccurate predictions at certain wave ages, as it fails to capture wave-induced turbulence, which fundamentally differs from mean shear-driven turbulence. The LES data show that accurate representations of wave-induced stresses require a complex-valued eddy viscosity. The maximum magnitude of this eddy viscosity scales as ∼𝑢 𝜏 ⁢𝜁 𝑖𝑛𝑛𝑒𝑟 , where 𝑢 𝜏 is the friction velocity and 𝜁 𝑖𝑛𝑛𝑒𝑟 is the inner-layer thickness, the height at which the eddy-turnover time matches the wave advection time scale. This scaling aligns with the prediction by Belcher & Hunt (1993 J. Fluid Mech. 251, 109–148). Overall, the findings demonstrate that traditional eddy viscosity models are inadequate for capturing wave-induced turbulence. More sophisticated turbulence models are essential for the accurate prediction of airflow disturbances and form drag in wind–wave interaction models.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗