Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atmospheric model”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dairy Methane Emissions in California's San Joaquin Valley Inferred With Ground‐Based Remote Sensing Observations in the Summer and Winter

Abstract The dairy industry in the San Joaquin Valley (SJV) is one of California’s largest methane (CH 4 ) sources. Reducing dairy emissions is a priority for the state’s climate change plans. Observations of current dairy CH 4 emissions are key to monitoring actions taken toward this goal. To help support this, we present new ground‐based measurements of atmospheric column‐averaged CH 4 mixing ratio (XCH 4 ) gradients across a group of 600 dairies in the central SJV using EM27/SUN solar spectrometers. We used measurements from the 2019 summer and 2020 winter seasons for a top‐down emission inversion based on the WRF‐STILT model. Our top‐down estimates of the region’s dairy emissions range from 90% to 183% of the current CALGEM inventory’s emissions of 277 Gg/yr. In contrast to the strong temperature dependence found by earlier dairy CH 4 emission studies, we also find that our top‐down emissions during the winter measurement days are comparable to the summer measurement days, possibly due to seasonal changes in dairy management practices and meteorological conditions. Furthermore, we find significant interday variability in our measurements and find that our emission estimates overlap with earlier top‐down studies and bottom‐up inventories in this region. Our study demonstrates how analysis of ground‐based remotely sensed CH 4 gradient observations can help improve our understanding of CH 4 sources at scales relevant to mitigation policy. It also reflects the need for long‐term monitoring of CH 4 emissions in the region and at individual facilities to better understand their emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Radiological Safety Analysis Computer (RSAC 7.2) Presentation

This presentation will be for the graduate seminar class for the University of New Mexico. We will discuss the Radiological Safety Analysis Computer (RSAC) version 7.2, a critical tool in radiological consequence analysis within the nuclear safety domain. We will begin with an introduction to radiological consequence analysis and the importance of dose modeling in nuclear safety, referencing key regulatory frameworks such as DOE-STD-3009 and 10 CFR 20. The presentation will provide an overview of RSAC, including its history, primary functions, and the inputs and outputs involved in its calculations. A detailed examination of RSAC's atmospheric dispersion model, based on the Gaussian plume methodology, will be presented, highlighting how meteorological factors influence the results. We will discuss the regulatory alignment of RSAC with DOE and NRC guidelines, emphasizing the assumptions required for compliance and the balance between conservative and realistic inputs. The presentation will also address the limitations and assumptions inherent in RSAC, comparing it with other consequence analysis codes like MACCS, HotSpot, and RASCAL. An example RSAC run will be demonstrated, showcasing a hypothetical accident scenario and interpreting the resulting dose calculations. Finally, we will introduce new ideas to enhance the current RSAC framework. This future work aims to provide multi-industry adaptability and integration with new dispersion models, addressing the evolving needs of DOE and NRC. The presentation will conclude with a discussion on RSAC's ongoing role in nuclear safety and the opportunities for future innovation.

computer program↗
DOE OSTI2026

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Numerical simulation of asteroid geometry variance on airburst threat

For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.

Airburst↗
DOE OSTI2020

Modeling the smoky troposphere of the southeast Atlantic: a comparison to ORACLES airborne observations from September of 2016

The southeast Atlantic is home to well-defined smoke outflow from Africa coinciding vertically with extensive marine boundary-layer cloud decks, both reaching their climatological maxima in spatial extent around September. A framework is put forth for evaluating the performance of a range of global and regional aerosol models against observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) airborne mission in September 2016. The sparse airborne observations are first aggregated into 2o grid boxes and into three vertical layers: the cloud-topped marine boundary layer (MBL), the layer from cloud top to 3 km, and the 3-6 km layer. Aerosol extensive properties simulated for the entire study region for all September suggest that the 2016 ORACLES observations are reasonably representative of the regional monthly average, with systematic deviations of 30% or less. All six models typically place the bottom of the smoke layer at lower altitudes than do the airborne lidar observations by 300-1400 m, whereas model aerosol top heights are within 0-500 m of the observations. All but one of the models that report carbonaceous aerosol masses underestimate the ratio of particulate extinction to the masses, a proxy for mass extinction efficiency, in 3-6 km. Notable findings on individual models include that WRF-CAM5 predicts the mass of black carbon and organic aerosols with minor (~10% or less) biases. GEOS-5 overestimates the carbonaceous particle masses in the MBL by a factor of 3-6. Extinction coefficients in the free troposphere (FT) and above-cloud aerosol optical depth (ACAOD) are 10-30% lower in WRF-CAM5, 30-50% lower in GEOS-5, 10-40% higher in GEOS-Chem, 10-20% higher in EAM-E3SM except for the practically unbiased 3-6 km extinction, and 20-70% lower in the Unified Model, than the airborne in situ, lidar and sunphotometer measurements. ALADIN-Climate also underestimates the ACAOD, by 30%. GEOS-5 and GEOS-Chem predict carbon monoxide in the MBL with small (10% or less) negative biases, despite their overestimates of carbonaceous aerosol masses. Overall, this study highlights a new approach to utilizing airborne aerosol measurements for model diagnosis.

Shinozuka, Yohei↗
DOE OSTI2025

AMR-Wind: A Performance-Portable, High-Fidelity Flow Solver for Wind Farm Simulations

We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.

17 WIND ENERGY↗
DOE OSTI2022

How long can luminous blue variables sleep? A long-term photometric variability and spectral study of the Galactic candidate luminous blue variable MN 112

ABSTRACT Luminous blue variables (LBVs) are massive stars that show strong spectral and photometric variability. The questions of what evolutionary stages they represent and what exactly drives their instability are still open, and thus it is important to understand whether LBVs without significant ongoing activity exist, and for how long such dormant LBVs may ‘sleep’. In this paper we investigate the long-term variability properties of the LBV candidate MN 112, by combining its optical and infrared spectral data covering 12 years with photometric data covering nearly a century, as acquired from both modern time-domain sky surveys and historical photographic plates. We analyse the spectra, derive the physical properties of the star by modelling its atmosphere, and use a new distance estimate from Gaia data release 3 (DR3) to determine the position of MN 112 both inside the Galaxy and in the Hertzsprung–Russell diagram. The distance estimation has almost doubled in comparison with Gaia DR2. Because of this, MN 112 moved to upper part of the diagram, and according to our modelling it lies on an evolutionary track for a star with initial mass $M_*=70~\rm M_\odot$ near the Humphreys–Davidson limit. Given the absence of any significant variability, we conclude that the star is a dormant LBV that has now been inactive for at least a century.

Maryeva, O. V. (ORCID:0000000314424755)↗
DOE OSTI2026

Seismic Observations of the OSIRIS-REx Sample Return Capsule Reentry: Deployment, Signal Characteristics, and Wavefield Phenomenology

Controlled spacecraft reentries from interplanetary trajectories provide rare, well-characterized hypersonic sources for advancing seismoacoustic observation techniques. Here we present seismic observations of the OSIRIS-REx sample return capsule (SRC) reentry on 24 September 2023, recorded by 16 three-component nodal seismometers deployed near Eureka, Nevada, at ground distances of 7–20 km from the capsule trajectory. Air-to-ground coupled signals are detected at all stations, exhibiting impulsive onsets consistent with ballistic shock arrivals from the descending Mach cone. We characterize the seismic wavefield through signal amplitude, period, waveform cross-correlation, and array processing. Signal periods decrease systematically with increasing distance from the trajectory within the airport array, indicating that higher-frequency content becomes more prominent at greater offsets, opposite to expectations from geometric spreading and atmospheric absorption. Seismic array processing identifies frequency-dependent back-azimuth variations whose origin remains unresolved; possible contributing factors include source geometry, scattering by fine-scale layered structure in the stratosphere, and near-surface effects. These observations document a spatially complex seismic wavefield from a well-characterized hypersonic line source and provide constraints for future modeling of atmospheric propagation and air-to-ground coupling.

OSIRIS-REx↗
DOE OSTI2022

The RAM-SCB model and its applications to advance space weather forecasting

Energetic particle fluxes that are part of the Earth’s ring current and radiation belts can intensify significantly during space weather events like geomagnetic storms and could cause severe damage to satellite-based technologies. Understanding the physical processes that control their dynamics and improving our capability for their prediction is thus extremely important. In the context of space weather applications and user needs, this paper provides a brief description of our kinetic ring current-atmosphere interactions model with self-consistent magnetic field (RAM-SCB) and its further extension to implement a self-consistent electric (E) field. Specific examples that demonstrate RAM-SCB capabilities and limitations to reproduce the near-Earth space weather environment are given. In this paper, the current status of RAM-SCB is assessed and plans for its further improvement are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2024

Reduced‐Order Probabilistic Emulation of Physics‐Based Ring Current Models: Application to RAM‐SCB Particle Flux

Abstract In this work, we address the computational challenge of large‐scale physics‐based simulation models for the ring current. Reduced computational cost allows for significantly faster than real‐time forecasting, enhancing our ability to predict and respond to dynamic changes in the ring current, valuable for space weather monitoring and mitigation efforts. Additionally, it can also be used for a comprehensive investigation of the system. Thus, we aim to create an emulator for the Ring current‐Atmosphere interactions Model with Self‐Consistent magnetic field (RAM‐SCB) particle flux that not only improves efficiency but also facilitates forecasting with reliable estimates of prediction uncertainties. The probabilistic emulator is built upon the methodology developed by Licata and Mehta (2023), https://doi.org/10.1029/2022sw003345 . A novel discrete sampling is used to identify 30 simulation periods over 20 years of solar and geomagnetic activity. Focusing on a subset of particle flux, we use Principal Component Analysis for dimensionality reduction and Long Short‐Term Memory (LSTM) neural networks to perform dynamic modeling. Hyperparameter space was explored extensively resulting in about 5% median symmetric accuracy across all data sets for one‐step dynamic prediction. Using a hierarchical ensemble of LSTMs, we have developed a reduced‐order probabilistic emulator (ROPE) tailored for time‐series forecasting of particle flux in the ring current. This ROPE offers accurate predictions of omnidirectional flux at a single energy with no pitch angle information, providing robust predictions on the test set with an error score below 11% and calibration scores under 8% with bias under 2% providing a significant speed up as compared to the full RAM‐SCB run.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Seasonal Patterns in Mobile Colistin Resistance Gene Variants in Wastewater Bioaerosols and Liquid Sludge

To address the growing threat of global antimicrobial resistance, a one-health approach is needed to understand the complex socioecological cycling of antibiotic resistance genes. In this study, a metagenomics approach using DNA shotgun sequencing, metagenome assembly, and antibiotic resistance gene (ARG) annotation was used to examine seasonal patterns in the abundance of mobile colistin resistance (mcr) gene variants in bioaerosols and liquid sludge in three wastewater treatment plants (WWTPs). ARGs represented 0.2–0.8 and 0.1–0.2% of the bioaerosol and liquid sludge metagenomes, respectively, while mcr genes represented 0–0.3 and 0–0.5% of the identified ARGs in bioaerosol and liquid sludge metagenomes. Seven of the ten known mcr variants were detected in wastewater bioaerosol and liquid samples, with mcr-5 and mcr-8 being the most prevalent across all seasons and sites. Furthermore, additional functional and taxonomic annotation of mcr-containing metagenomic contigs showed that mcr genes were often located on contigs with other co-occurring ARGs and mobile genetic elements and may be harbored by opportunistic human pathogens and other bacterial taxa not previously associated with mcr genes. Atmospheric dispersion modeling showed that mcr-containing bioaerosols can be transported kilometers away from the WWTPs, resulting in the possible dissemination of these ARGs into surrounding environments and communities.

Antimicrobial agents↗
DOE OSTI2022

Spatial heterogeneity effects on land surface modeling of water and energy partitioning

Abstract. Understanding the influence of land surface heterogeneity on surface water and energy fluxes is crucial for modeling earth system variability and change. This study investigates the effects of four dominant heterogeneity sources on land surface modeling, including atmospheric forcing (ATM), soil properties (SOIL), land use and land cover (LULC), and topography (TOPO). Our analysis focused on their impacts on the partitioning of precipitation (P) into evapotranspiration (ET) and runoff (R), partitioning of net radiation into sensible heat and latent heat, and corresponding water and energy fluxes. An initial set of 16 experiments were performed over the continental US (CONUS) using the E3SM land model (ELMv1) with different combinations of heterogeneous and homogeneous datasets. The Sobol' total and first-order sensitivity indices were utilized to quantify the relative importance of the four heterogeneity sources. Sobol' total sensitivity index measures the total heterogeneity effects induced by a given heterogeneity source, consisting of the contribution from its own heterogeneity (i.e., the first-order index) and its interactions with other heterogeneity sources. ATM and LULC are the most dominant heterogeneity sources in determining spatial variability of water and energy partitioning, mainly contributed by their own heterogeneity and slightly contributed by their interactions with other heterogeneity sources. Their heterogeneity effects are complementary, both spatially and temporally. The overall impacts of SOIL and TOPO are negligible, except TOPO dominates the spatial variability of R/P across the transitional climate zone between the arid western and humid eastern CONUS. Accounting for more heterogeneity sources improves the simulated spatial variability of water and energy fluxes when compared with ERA5-Land reanalysis dataset. An additional set of 13 experiments identified the most critical components within each heterogeneity source, which are precipitation, temperature, and longwave radiation for ATM, soil texture, and soil color for SOIL and maximum fractional saturated area parameter for TOPO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Identifying Characteristics of Wildfire Towers and Troughs

Wildfire behavior is dictated by the complex interaction of numerous physical phenomena including dynamic ambient and fire-induced winds, heat transfer, aerodynamic drag on the wind by the fuel and combustion. These phenomena create complex feedback effects between the fire and its surroundings. In this study, we aim to study the mechanisms by which buoyant flame dynamics along with vortical motions and instabilities control wildfire propagation. Specifically, this study employs a suite of simulations conducted with the physics-based coupled fire-atmosphere behavior model (FIRETEC). The simulations are initialized with a fire line and the fires are allowed to propagate on a grass bed, where the fuel heights and wind conditions are varied systematically. Flow variables are extracted to identify the characteristics of the alternating counter-rotational vortices, called towers and troughs, that drive convective heat transfer and fire spread. These vortices have previously been observed in wildfires and laboratory fires, and have also been observed to arise spontaneously in FIRETEC due to the fundamental physics incorporated in the model. However, these past observations have been qualitative in nature and no quantitative studies can be found in the literature which connected these coherent structures fundamental to fire behavior with the constitutive flow variables. To that end, a variety of state variables are examined in the context of these coherent structures under various wind profile and grass height conditions. Identification of various correlated signatures and fire-atmosphere feedbacks in simulations provides a hypothesis that can be tested in future observational or experimental efforts, potentially assisting experimental design, and can aid in the interpretation of data from in situ detectors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using STAX data to predict IMS radioxenon concentrations

The noble gas collection and measurement stations in the International Monitoring System (IMS) are heavily influenced by releases from medical isotope production facilities. The ability to reliably model the movement of radioxenon from the points of release to these IMS samplers has improved enough that a routine aspect of the analysis of IMS radioxenon data should be the prediction of the effect of releases from civilian nuclear facilities on the sample concentrations. In this work, predicted concentrations at IMS noble gas systems in Germany and Sweden based on measured releases from Institute for Radioelements (IRE) in Belgium and atmospheric transport modeling for a four-month period are presented and discussed.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Improved constraints on hematite refractive index for estimating climatic effects of dust aerosols

Abstract Uncertainty in desert dust composition poses a big challenge to understanding Earth’s climate across different epochs. Of particular concern is hematite, an iron-oxide mineral dominating the solar absorption by dust particles, for which current estimates of absorption capacity vary by over two orders of magnitude. Here, we show that laboratory measurements of dust composition, absorption, and scattering provide valuable constraints on the absorption potential of hematite, substantially narrowing its range of plausible values. The success of this constraint is supported by results from an atmospheric transport model compared with station-based measurements. Additionally, we identify substantial bias in simulating hematite abundance in dust aerosols with current soil mineralogy descriptions, underscoring the necessity for improved data sources. Encouragingly, the next-generation imaging spectroscopy remote sensing data hold promise for capturing the spatial variability of hematite. These insights have implications for enhancing dust modeling, thus contributing to efforts in climate change mitigation and adaptation.

Environmental Sciences & Ecology↗
DOE OSTI2020

Near-canopy horizontal concentration heterogeneity of semivolatile oxygenated organic compounds and implications for 2-methyltetrols primary emissions

Semivolatile oxygenated organic compounds (SV-OVOCs) are important atmospheric species, in particular for the production and chemistry of atmospheric particulate matter and related impacts on air quality and climate. In this study, SV-OVOCs were collected in the horizontal plane of the roughness layer over the tropical forest in the central Amazon during the wet season of 2018. A sampler mounted to a copter-type, hovering unmanned aerial vehicle was used. Underlying the collection region, a plateau forest transitioned into a slope forest across several hundred meters. The concentrations of pinonic and pinic acids, which are monoterpene oxidation products, had no statistical difference over the two forests. By comparison, across the study period, differences in the concentration of 2-methyltetrols, which are products of isoprene oxidation, ranged from -70% to +480% over the two forests. The chemical lifetime of 2-methyltetrols in the atmosphere is sufficiently long that heterogeneity in the isoprene emission rate from the two forests followed by atmospheric oxidation does not explain the concentration heterogeneity of 2-methyltetrols. Standing waves and local meteorology also cannot account for the heterogeneity. Of the possibilities considered, the most plausible explanation is the direct emission from the forest of 2-methyltetrols produced through biological processes within the plants. Under this explanation, the rate of direct SV-OVOC emissions should be modulated by forest type and related environmental stressors. Direct emissions of SV-OVOCs should be more broadly considered for constraining and improving models of atmospheric composition, transport, and chemistry over tropical forests.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Generation of reactive species in water film dielectric barrier discharges sustained in argon, helium, air, oxygen and nitrogen

Activation of liquids with atmospheric pressure plasmas is being investigated for environmental and biomedical applications. When activating the liquid using gas plasma produced species (as opposed to plasmas sustained in the liquid), a rate limiting step is transport of these species into the liquid. To first order, the efficiency of activating the liquid is improved by increasing the ratio of the surface area of the water in contact with the plasma compared to its volume—often called the surface-to-volume ratio (SVR). Maximizing the SVR then motivates the plasma treatment of thin films of liquids. In this paper, results are discussed from a computational investigation using a global model of atmospheric pressure plasma treatment of thin water films by a dielectric barrier discharge (DBD) sustained in different gases (Ar, He, air, N 2 , O 2 ). The densities of reactive species in the plasma activated water (PAW) are evaluated. The residence time of the water in contact with the plasma is increased by recirculating the PAW in plasma reactor. Longer lived species such as H 2 O 2aq and NO 3 - aq accumulate over time (aq denotes an aqueous species). DBDs sustained in Ar and He are the most efficient at producing H 2 O 2aq , DBDs sustained in argon produces the largest density of NO 3 - aq with the lowest pH, and discharges sustained in O 2 and air produce the highest densities of O 3aq . Finally, comparisons to experiments by others show agreement in the trends in densities in PAW including O 3aq , OH aq , H 2 O 2aq and NO 3 - aq , and highlight the importance of controlling desolvation of species from the activated water.

42 ENGINEERING↗