Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional neutrality”

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 109 records · Page 6

The sensitivity of the Fitch wind farm parameterization to a three-dimensional planetary boundary layer scheme

Abstract. Wind plant wake impacts can be estimated with a number of simulation methodologies, each with its own fidelity and sensitivity to model inputs. In turbine-free mesoscale simulations, hub-height wind speeds often significantly vary with the choice of a planetary boundary layer (PBL) scheme. However, the sensitivity of wind plant wakes to a PBL scheme has not been explored because, as of the Weather Research and Forecasting model v4.3.3, wake parameterizations were only compatible with one PBL scheme. We couple the Fitch wind farm parameterization with the new NCAR 3DPBL scheme and compare the resulting wakes to those simulated with a widely used PBL scheme. We simulate a wind plant in pseudo-steady states under idealized stable, neutral, and unstable conditions with matching hub-height wind speeds using two PBL schemes: MYNN and the NCAR 3DPBL. For these idealized scenarios, average hub-height wind speed losses within the plant differ between PBL schemes by between −0.20 and 0.22 m s−1, and correspondingly, capacity factors range between 39.5 %–53.8 %. These simulations suggest that PBL schemes represent a meaningful source of modeled wind resource uncertainty; therefore, we recommend incorporating PBL variability into future wind plant planning sensitivity studies as well as wind forecasting studies.

17 WIND ENERGY↗

Inverse Segmented Motor Drive Using Dual ANPC Inverters for Common-Mode Voltage and Neutral-Point Current Cancellation

This article proposes an inverse segmented motor drive (SgMD) utilizing dual active neutral point clamped (ANPC) inverters. In the proposed configuration, the neutral point current and common-mode (CM) voltage is topologically canceled, achieving zero total neutral point current and CM voltage under ideal conditions. Also, the zero total neutral point current minimizes the neutral point voltage imbalance in ANPC inverters. The mechanisms behind neutral point current and CM voltage cancellation in the proposed inverse SgMD are first introduced. The modifications to the motor windings for implementing the inverse SgMD are explained, showing that a standard motor can be readily adapted for the proposed configuration. A space vector modulation (SVM) scheme tailored for the proposed topology is presented, along with a carrier-based implementation. Simulation results validate that the proposed topology can achieve zero total CM voltage and neutral point current. It is also shown that the proposed inverse SgMD can reduce neutral point voltage fluctuation by about 90% and RMS current stress in the dc-link capacitors by about 43% compared to the conventional SgMD.

Lee, Sangwhee [ORNL] (ORCID:0000000335490057)↗

An expanded kinetics model describing injector region desorption

During the initial stages of electron beam injection, absorbed or entrained material can be released into the nascent beam pulse by thermal desorption or electron beam impact. These processes, and subsequent ionization and dissociation, are examined assuming a prototypical temporal beam pulse. An expanded set of kinetic equations that includes thermal desorption, electron-stimulated ion emission, electron impact dissociation, electron impact ionization, and thermal decomposition of adsorbed water molecules is developed. Estimates of the three most important ionized species H 2 O + , H + and OH + are provided for typical injector conditions. Electron impact dissociative ionization of neutral H 2 O is the dominant source of H + under these conditions. Since the abundance of thermally desorbed neutral H 2 O has an exponential rate, the production of H+ is likewise controlled by this exponential desorption rate. Thus thermal management of the confining walls could suppress proton production in the injector region. Thermal decomposition of H 2 O will not be a contributing physical process in the injector region on time scales less than a μs. The formalism is sufficiently flexible to include entrained hydrogen.

42 ENGINEERING↗

Rejection of low-molecular weight neutral organics is highly sensitive to reverse osmosis system design and operation

A computational model was developed to investigate the significance of system design and operating conditions on the rejection of neutral, low-MW organics by reverse osmosis for potable reuse. Here, the model demonstrated that the decrease in local rejection as net driving pressure decreases is substantially greater for moderately rejected compounds than for highly rejected compounds. At recovery values less than 70%, the local permeate concentration can exceed the pressure vessel feed concentration for moderately rejected compounds. System-level rejection of moderately rejected compounds is likewise substantially more sensitive to operating conditions than highly rejected compounds. The findings highlight a drawback of relying on rejection results from bench-scale testing that operates at low recovery, which invariably has higher rejection than full-scale systems operating at similar pressure. The analysis demonstrates a trade-off in which the low-pressure, high-recovery operation desired for potable reuse systems can be detrimental to the removal of low-MW neutral organics. The removal of low-MW neutral organics can be improved if organics rejection is explicitly evaluated during the design process.

42 ENGINEERING↗

Large eddy simulations of offshore wind turbine wakes for two floating platform types

The growing prospect for large farms of floating offshore wind turbines requires a better understanding of wake effects for floating turbines. In this work, large eddy simulations with an actuator line model are used to study the wake of the NREL 5 MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCwind semi-submersible platforms. The simulations are carried out in the Simulator fOr Wind Farm Applications (SOWFA) coupled with OpenFAST for the platform and turbine motion. The wake location, deficit, and turbulence levels are compared for the two floating platforms and equivalent fixed-turbine cases. The effects of neutral versus stable atmospheric conditions are also compared. Most notably, floating-turbine wakes are deflected upwards compared to fixed-turbine wakes, because of mean platform pitch. The spar wake deflects upwards more than the semi-submersible, while the stable atmosphere increases this vertical deflection compared to the neutral. The time-varying rotor motions do not significantly affect the mid-to-far wake, though the stable atmosphere shows larger fixed-floating differences in horizontal wake fluctuations.

17 WIND ENERGY↗

Observations of the Marine Atmospheric Boundary Layer’s Response to a Solar Eclipse

The atmospheric response to the solar eclipse of 8 April 2024 in North America is investigated with a specific focus on the marine atmospheric boundary layer (MABL). We leverage measurements collected during the Third Wind Forecast Improvement Project (WFIP3), including Doppler lidars, sonic anemometers, and thermodynamic profiler data to investigate the atmospheric response across sites that experienced partial eclipse conditions with nearly 90% obscuration. Using these measurements, we examine eclipse-induced changes in key meteorological parameters, such as temperature, wind speed, and turbulent fluxes. Most previous eclipse studies have been conducted over land, whereas this study provides new observations for both coastal and marine environments, offering additional insight into eclipse-driven variability in the MABL. The findings confirm a notable decrease in downwelling shortwave radiation during the eclipse, which results in rapid cooling of surface air. The temperature reduction ranges from $1.2^\circ \text {C}$ to $1.4^\circ \text {C}$ in coastal regions and from $0.3^\circ \text {C}$ to $0.5^\circ \text {C}$ over the ocean. This analysis suggests that the MABL’s higher thermal inertia compared to coastal regions moderates the temperature decrease during the eclipse. Wind speed exhibits a more complex behavior, as it is influenced by both the MABL and preexisting synoptic conditions. Although a reduction in wind speed is observable up to approximately 140 m above ground level (AGL) at more inland sites, at other locations closer to the coast, this reduction is constrained to the lowest 100 m AGL. Turbulence parameters retrieved from sonic anemometers, such as turbulence kinetic energy, turbulent heat flux, and friction velocity, decrease during the eclipse at coastal sites, accompanied by a brief transition of atmospheric stability from unstable to neutral or weakly stable conditions. For the open-ocean sites, the variability in turbulence statistics and atmospheric stability is minimal during the occurrence of the eclipse.

16 TIDAL AND WAVE POWER↗

High manganese redox variability and manganate predominance in temperate soil profiles as determined by X-ray absorption spectroscopy

Manganese speciation is a key to understanding the fate of contaminants, nutrients, and organic matter in soils. To date, quantification of Mn species in bulk soils has been performed mainly by sequential extraction methods and rarely supported by spectroscopic analysis. In order to obtain quantitative information on the Mn species inventory of soils, we investigated 46 soil horizons (<2-mm fraction, 45.1–2,280 mg/kg Mn) of nine typical Central European soils (Cambisols, Chernozems, Luvisols, Podzol, Stagnosol) by chemical Mn analyses and Mn K-edge X-ray absorption spectroscopy, and related speciation results to major soil properties. Amounts of Mn 2+ , Mn 3+ , and Mn 4+ and the average oxidation state of Mn were evaluated by linear combination fitting (LCF) of X-ray absorption near edge structure (XANES) spectra. Additionally, we used extended X-ray absorption fine structure (EXAFS) spectroscopy to identify and quantify major Mn species. For this, EXAFS spectra of 20 organic and mineral soil samples from five soils (Cambisols, Chernozem, Luvisol) were analyzed by LCF and shell fitting. XANES analyses revealed a high Mn redox variability in organic surface layers, with Mn 2+ being most abundant (≤100%, $\overline{x}$ = 54%), followed by Mn 3+ (≤80%, $\overline{x}$ = 32%) and Mn 4+ (≤55%, $\overline{x}$ = 14%). Mineral soil horizons contained significantly less Mn 2+ (≤56%, $\overline{x}$ = 23%), about equal quantities of Mn 3+ (≤68%, $\overline{x}$ = 31%), and were enriched in Mn 4+ (≤89%, $\overline{x}$ = 46%). EXAFS analyses implied the presence of six major Mn species groups: manganates, organically complexed Mn, Mn(III) oxyhydroxides, silicate-bound Mn, Mn oxides without tunnel- or layer structure, and physisorbed Mn. In litter horizons, Mn was mainly present in organic complexes (58–91%, $\overline{x}$ = 78%) and as physisorbed Mn (≤15%), but individual horizons also comprised manganates, Mn(III) oxyhydroxides, and silicate-bound Mn. Manganates, likely mixtures of phyllomanganates with hexagonal layer symmetry and tectomanganates, dominated in all mineral soil horizons (37–94%, $\overline{x}$ = 67%). Correlation analysis showed that manganates dissolve completely during dithionite-citrate and acid ammonium oxalate extractions, and suggested that Mn 4+ -rich manganates preferentially form under less acidic soil conditions, partly by oxidation of organically complexed Mn(II), and that they are enriched in the soil clay fraction. Mineral soil horizons also contained minor quantities of organically complexed Mn (≤39%, $\overline{x}$ = 11%), silicate-bound Mn (≤30%, $\overline{x}$ = 8%), Mn(III) oxyhydroxides (≤37%, $\overline{x}$ = 7%), Mn oxides without tunnel- or layer structure (≤18%, $\overline{x}$ = 5%), and physisorbed Mn (≤14%, $\overline{x}$ <1%). The detection of Mn(III) oxyhydroxides such as feitknechtite (β-MnOOH) or groutite (α-MnOOH) as well as the spinel hausmannite (Mn 3 O 4 ) in acidic soils is remarkable, since their formation is normally linked to neutral or alkaline pH conditions. Minor contributions of silicate-bound Mn indicate the release of Mn from primary minerals already at early stages of soil formation, and low concentrations of physisorbed Mn suggest that exchangeable Mn is rapidly converted into manganates in oxic soils. In conclusion, the predominance of manganates in mineral soils has far-reaching implications for the functioning of soils and biogeochemical element cycles, as these minerals play an important role in metal binding, plant nutrition, and redox-related processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis and Characterization of Silver-Modified Nanoporous Silica Materials for Enhanced Iodine Removal

In aquatic environments, the presence of iodine species, including radioactive isotopes like 129 I and I 2 , poses significant environmental and health concerns. Iodine can enter water resources from various sources, including nuclear accidents, medical procedures, and natural occurrences. To address this issue, the use of natural occurring nanoporous minerals, such as zeolitic materials, for iodine removal will be explored. This study focuses on the adsorption of iodine by silver-modified zeolites (13X-Ag, 5A-Ag, Chabazite-Ag, and Clinoptilolite-Ag) and evaluates their performance under different conditions. All materials were characterized using scanning electron microscopey (SEM), energy-dispersive X-ray spectroscopy (EDS), powdered X-ray diffraction (P-XRD), Fourier-transform infrared spectrometry (FTIR), and nitrogen adsorption studies. The results indicate that Chabazite-Ag exhibited the highest iodine adsorption capacity, with an impressive 769 mg/g, making it a viable option for iodine removal applications. 13X-Ag and 5A-Ag also demonstrated substantial adsorption capacities of 714 mg/g and 556 mg/g, respectively, though their behavior varied according to different models. In contrast, Clinoptilolite-Ag exhibited strong pH-dependent behavior, rendering it less suitable for neutral to slightly acidic conditions. Furthermore, this study explored the impact of ionic strength on iodine adsorption, revealing that Chabazite-Ag is efficient in low-salinity environments with an iodine adsorption capacity of 51.80 mg/g but less effective in saline conditions. 5A-Ag proved to be a versatile option for various water treatments, maintaining its iodine adsorption capacity across different salinity levels. In contrast, Clinoptilolite-Ag exhibited high sensitivity to ionic competition, virtually losing its iodine adsorption ability at a NaCl concentration of 0.1 M. Kinetic studies indicated that the pseudo-second-order model best describes the adsorption process, suggesting chemisorption mechanisms dominate iodine removal. Chabazite-Ag exhibited the highest initial adsorption rate with a k 2 value of 0.002 mg g -1 h -1 , emphasizing its superior adsorption capabilities. Chabazite and Clinoptilolite, naturally occurring minerals, provide eco-friendly solutions for iodine adsorption. Chabazite superior iodine removal highlights its value in critical applications and its potential for addressing pressing environmental challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Catabolic pathway acquisition by rhizosphere bacteria readily enables growth with a root exudate component but does not affect root colonization

Horizontal gene transfer (HGT) is a fundamental evolutionary process that plays a key role in bacterial evolution. The likelihood of a successful transfer event is expected to depend on the precise balance of costs and benefits resulting from pathway acquisition. Most experimental analyses of HGT have focused on phenotypes that have large fitness benefits under appropriate selective conditions, such as antibiotic resistance. However, many examples of HGT involve phenotypes that are predicted to provide smaller benefits, such as the ability to catabolize additional carbon sources. We have experimentally simulated the consequences of one such HGT event in the laboratory, studying the effects of transferring a pathway for catabolism of the plant-derived aromatic compound salicyl alcohol between rhizosphere isolates from the Pseudomonas genus. We find that pathway acquisition enables rapid catabolism of salicyl alcohol with only minor disruptions to the existing metabolic and regulatory networks of the new host. However, this new catabolic potential does not confer a measurable fitness advantage during competitive growth in the rhizosphere. We conclude that the phenotype of salicyl alcohol catabolism is readily transferable but is selectively neutral under environmentally relevant conditions. We propose that this condition is common and that HGT of many pathways will be self-limiting because the selective benefits are small.

59 BASIC BIOLOGICAL SCIENCES↗

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

High-fidelity wind farm simulation methodology with experimental validation

The complexity and associated uncertainties involved with atmospheric-turbine-wake interactions produce challenges for accurate wind farm predictions of generator power and other important quantities of interest (QoIs), even with state-of-the-art high-fidelity atmospheric and turbine models. A comprehensive computational study was undertaken with consideration of simulation methodology, parameter selection, and mesh refinement on atmospheric, turbine, and wake QoIs to identify capability gaps in the validation process. For neutral atmospheric boundary layer conditions, the massively parallel large eddy simulation (LES) code Nalu-Wind was used to produce high-fidelity computations for experimental validation using high-quality meteorological, turbine, and wake measurement data collected at the Department of Energy/Sandia National Laboratories Scaled Wind Farm Technology (SWiFT) facility located at Texas Tech University’s National Wind Institute. The wake analysis showed the simulated lidar model implemented in Nalu-Wind was successful at capturing wake profile trends observed in the experimental lidar data.

17 WIND ENERGY↗

Insights into determining pore size properties of ultrafiltration membranes

The selectivity of porous membranes is often characterized using solute rejection tests, where membranes are challenged with dilute aqueous solutions of neutral solutes at operating conditions that minimize concentration polarization and fouling. In single solute tests, a membrane is challenged with one molecular weight (MW) solute at a time from low to high MW. Since single solute methods are time-intensive, mixed solute tests have become more common, where a mixture of several MW solutes challenges a membrane at once. However, the presence of large solutes in a mixture increases the rejection of smaller solutes. Furthermore, there are no universally accepted operating conditions or standard methods used by membrane manufacturers or researchers for the experiments, leading to difficulty in pore size and pore characteristic comparisons. In this paper, commercial ultrafiltration membranes were challenged with single and mixed solute polyethylene glycol (PEG) and dextran aqueous solutions. First, rejection values determined using total organic carbon (TOC) and high-performance liquid chromatography (HPLC) from single solute filtration experiments are compared. Differences in rejection curves obtained by the two techniques are attributed to solute polydispersity. Mixed solute filtration experiments with binary mixtures of solutes showcased co-solute interactions, which increase with both the size and weight percent of large solute in the mixture. Mixed solute filtration experiments at varying operating conditions (i.e., stir speed and flux) were conducted to determine operating conditions that mitigate co-solute interactions. Stir speed had a minimal effect on co-solute interactions. In contrast, low flux conditions can help minimize co-solute interactions, leading to pore size distributions that closely resemble results observed in single solute filtration using narrowly dispersed solutes. Additionally, at low flux conditions, the predicted membrane pore size distributions utilizing mixed solute experiments with PEG and dextran were similar.

36 MATERIALS SCIENCE↗

Use of a Multiple Hydride Donor To Achieve an n-Doped Polymer with High Solvent Resistance

We report the ability to insolubilize doped semiconducting polymer layers can help enable the fabrication of efficient multilayer solution-processed electronic and optoelectronic devices. Here, we present a promising approach to simultaneously n-dope and largely insolubilize conjugated polymer films using tetrakis[{4-(1,3-dimethyl-2,3-dihydro-1H-benzo[d]imidazol-2-yl)phenoxy}methyl]methane (tetrakis-O-DMBI-H), which consists of four 2,3-dihydro-1H-benzoimidazole (DMBI-H) n-dopant moieties covalently linked to one another. Doping a thiophene-fused benzodifurandione-based oligo(p-phenylenevinylene)-co-thiophene polymer (TBDOPV-T) with tetrakis-O-DMBI-H results in a highly n-doped film with bulk conductivity of 15 S cm –1 . Optical absorption spectra provide evidence for film retention of ~93% after immersion in o-dichlorobenzene for 5 min. The optical absorption signature of the charge carriers in the n-doped polymer decreases only slightly more than that of the neutral polymer under these conditions, indicating that the exposure to solvent also results in negligible dedoping of the film. Moreover, thermal treatment studies on a tetrakis-O-DMBI-H-doped TBDOPV-T film in contact with another undoped polymer film indicate immobilization of the molecular dopant in TBDOPV-T. This is attributed to the multiple electrostatic interactions between each dopant tetracation and up to four nearby anionic doped polymer segments.

36 MATERIALS SCIENCE↗

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Enhanced pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal high-βp scenario for the National Spherical Torus Experiment Upgrade (NSTX-U) and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (βN > 5), and significant bootstrap fraction (fBS > 0.6) at Ip/BT = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. The EP H-mode was most often observed as a transition following a large edge-localized mode in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with an elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the kinetic ballooning mode and kink-peeling mode, can be responsible for the increased particle transport in EP H-mode.

Battaglia, D. J. (ORCID:0000000188979740)↗

Enhanced Pedestal H-mode at low edge ion collisionality on NSTX

The Enhanced Pedestal (EP) H-mode regime is an attractive wide-pedestal ELM-free high-betap scenario for NSTX-U and next-step devices as it achieves enhanced energy confinement (H98y,2 > 1.5), large normalized pressure (betaN > 5) and significant bootstrap fraction (f_BS > 0.6) at I_p/B_T = 2 MA/T. This regime is realized when the edge ion collisionality becomes sufficiently small that a positive feedback interaction occurs between a reduction in the ion neoclassical energy transport and an increase in the particle transport from pressure-driven edge instabilities. EP H-mode was most often observed as a transition following a large ELM in conditions with low edge neutral recycling. It is hypothesized that the onset of pressure-driven instabilities prior to the full recovery of the neutral density leads to a temporary period with elevated ion temperature gradient that triggers the transition to EP H-mode. Linear CGYRO and M3D-C1 calculations are compared to beam emission spectroscopy (BES) and magnetic spectroscopy in order to describe the evolution of the edge particle transport mechanisms during the ELM recovery and the saturated EP H-mode state. The observations are consistent with the hypothesis that the onset of pressure-driven edge instabilities, such as the KBM and kink-peeling, can be responsible for the increased particle transport in EP H-mode.

edge and boundary physics↗

Tropospheric Temperature

The 2020 annual global lower tropospheric temperature (LTT) tied with 2016 as the highest on record. The annual average LTT was 0.49°–0.72°C above the 1981–2010 average, depending on dataset, and 10%–16% of Earth’s surface experienced record high temperatures (Plate 2.1e). Such expansive and record warmth is notable because it occurred even though the ENSO exhibited neutral or La Niña conditions throughout the year (Fig. 2.7a) and is thus consistent with the background upward trend since 1958.

54 ENVIRONMENTAL SCIENCES↗

A Graphics Processing Unit (GPU) Approach to Large Eddy Simulation (LES) for Transport and Contaminant Dispersion

Recent advances in the development of large eddy simulation (LES) atmospheric models with corresponding atmospheric transport and dispersion (AT&D) modeling capabilities have made it possible to simulate short, time-averaged, single realizations of pollutant dispersion at the spatial and temporal resolution necessary for common atmospheric dispersion needs, such as designing air sampling networks, assessing pollutant sensor system performance, and characterizing the impact of airborne materials on human health. The high computational burden required to form an ensemble of single-realization dispersion solutions using an LES and coupled AT&D model has, until recently, limited its use to a few proof-of-concept studies. An example of an LES model that can meet the temporal and spatial resolution and computational requirements of these applications is the joint outdoor-indoor urban large eddy simulation (JOULES). A key enabling element within JOULES is the computationally efficient graphics processing unit (GPU)-based LES, which is on the order of 150 times faster than if the LES contaminant dispersion simulations were executed on a central processing unit (CPU) computing platform. JOULES is capable of resolving the turbulence components at a suitable scale for both open terrain and urban landscapes, e.g., owing to varying environmental conditions and a diverse building topology. In this paper, we describe the JOULES modeling system, prior efforts to validate the accuracy of its meteorological simulations, and current results from an evaluation that uses ensembles of dispersion solutions for unstable, neutral, and stable static stability conditions in an open terrain environment.

54 ENVIRONMENTAL SCIENCES↗