Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Improving energy yield in photovoltaic modules with photonic structures

When installed outdoors, Si solar cells typically operate 20 – 30 K above ambient conditions. These elevated operating temperatures lead to diminished performance, reducing efficiency by ~0.4%/K for crystalline Si cells and decreasing module lifetime. Since much of the elevated temperature of the modules derives from parasitic absorption of sub-bandgap sunlight by the contacts, encapsulants, or other materials, reflectors that remove the sub-bandgap radiation from the module before it is absorbed promise to reduce the operating temperature. However, there are many competing strategies to reject this light: selective reflectors could be integrated into the outer module glass, at the textured Si-encapsulant interface, or on the rear. Furthermore, many of these approaches are complex and expensive to fabricate. The issue of heating can be even more critical in bifacial modules: high fractions of rear irradiance can lead to higher operating temperatures. Prior to this project, relatively little was known about the optical strategies for temperature reduction in bifacial modules. The goal of this project was to quantitatively assess optical strategies for improving the energy yield via a combination of improved anti-reflection and sub-bandgap light rejection. We studied two complementary concepts: the performance of ideal structures to determine the upper limits to performance, and low-complexity structures made from standard materials as cost-effective implementations. We also studied the performance of these reflectors in different types of modules, including Al BSF, PERC, and bifacial. Here we build from our prior accomplishments in developing mirror optimization techniques, and make use of simulation methods that accurately capture the optical properties of the modules throughout the visible and near-infrared spectrum. We connect these optical properties to energy yield under realistic, outdoor conditions. We also experimentally fabricated mirrors, incorporated them into mini-modules, and performed outdoor testing of performance at NREL. We found definitively that spectrally-selective mirrors on the outside of the module cover glass offer the best performance, both under ideal conditions and for realistic, low- complexity mirrors. Mirrors at the Si cell/encapsulant interface are limited by the multiple reflections needed to remove sub-bandgap light. Mirrors at the cell surface and mirrors on the rear also cannot remove all the sub-bandgap light, since parasitic absorption also occurs in the cover glass, encapsulant, and front contacts, depending on the module configuration. Also, cell texturing places higher demands on any cell-surface mirror because multiple interactions with the mirror are required to reject sub bandgap light which compounds losses. For bifacial modules, mirrors on the back do provide additional cooling and energy yield improvements, and the mirrors can have the same structure. In the monofacial case up to 3.3K of cooling is possible for mirrors on the front, compared to 2.2K at the cell surface and 1.2K at the rear. For bifacial modules, up to 2.4K of cooling is possible when mirrors are included at both the front and back. Our realistic designs utilize 4 – 6 layers of common photovoltaic materials such as SiO 2 , TiO 2 , and SiN x , making them more cost effective. Finally, we used our models to study thermal management in bifacial modules with different back lamination materials, finding that glass-glass laminated modules operated hotter than equivalent glass-polymer backsheet modules, due solely to the increased absorption of light incident from the rear of the module. Improving the energy yield by 2%, which is achievable with these realistic designs, would lead to LCOE reductions of $0.03/kWh by 2030. These benefits make photovoltaic panels more accessible to the public and may reduce associated maintenance costs.

14 SOLAR ENERGY↗
DOE OSTI2022

Electrically controlling and optically observing the membrane potential of supported lipid bilayers

Supported lipid bilayers are a well-developed model system for the study of membranes and their associated proteins, such as membrane channels, enzymes, and receptors. These versatile model membranes can be made from various components, ranging from simple synthetic phospholipids to complex mixtures of constituents, mimicking the cell membrane with its relevant physiochemical and molecular phenomena. In addition, the high stability of supported lipid bilayers allows for their study via a wide array of experimental probes. In this work, we describe a platform for supported lipid bilayers that is accessible both electrically and optically, and demonstrate direct optical observation of the transmembrane potential of supported lipid bilayers. We show that the polarization of the supported membrane can be electrically controlled and optically probed using voltage-sensitive dyes. Membrane polarization dynamics is understood through electrochemical impedance spectroscopy and the analysis of an equivalent electrical circuit model. In addition, we describe the effect of the conducting electrode layer on the fluorescence of the optical probe through metal-induced energy transfer, and show that while this energy transfer has an adverse effect on the voltage sensitivity of the fluorescent probe, its strong distance dependency allows for axial localization of fluorescent emitters with ultrahigh accuracy. We conclude with a discussion on possible applications of this platform for the study of voltage-dependent membrane proteins and other processes in membrane biology and surface science.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing

Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

FIU Project 2: Environmental Remediation Science & Technology [Slides]

FIU’s research under this project involves conducting basic and applied science to fill knowledge gaps and validate potential remediation technologies for contaminated soil and groundwater and the assessment of the fate and transport of contaminants in the environment. The aim of FIU’s research is to reduce the potential for contaminant mobility or toxicity in the surface and subsurface through the development and application of state-of-the-art scientific and environmental remediation technologies at the Hanford Site, Savannah River Site (SRS), and the Waste Isolation Pilot Plant (WIPP), which is the Nation’s only mined geologic repository for permanent disposal of transuranic waste. FIU collaborates with scientists from Pacific Northwest National Laboratory (PNNL), Savannah River National Laboratory (SRNL), Savannah River Ecology Laboratory (SREL), Los Alamos National Laboratory (LANL) and the DOE Carlsbad Field Office (CBFO) in order to plan and execute research that is synergistic with the work being conducted at the sites, and that supports the resolution of critical science and engineering needs which leads to a better understanding of the long-term behavior of subsurface contaminants. The knowledge gained through this research will be used to transform experimental and modeling innovations into practical applications deployed at the sites to support EM’s primary goal of expediting the closure of major contaminated soil and groundwater sites and waste units. Collaborative relationships between FIU and the national laboratories have provided large benefits over the years to FIU, the national laboratories, the DOE complex, and the DOE EM mission. By working closely with the national laboratories, FIU’s research is not only closely aligned with the cleanup mission priorities at the DOE sites, but complements and supports ongoing work at the national laboratories for screening of new remedial technologies. This coordination and leveraging of research efforts results in time- and cost-savings, and will accelerate progress of the DOE EM environmental restoration mission.

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

Combined ambient vibration and surface displacement measurements for improved progressive failure monitoring at a toppling rock slab in Utah, USA

Seismic resonance and surface displacement measurements can be implemented in tandem to improve landslide characterization and progressive failure monitoring. Crack aperture data are frequently used in rock slope stability monitoring and often exhibit recognizable trends prior to failure, such as accelerated crack opening. Alternatively, ambient resonance data offer multiple parameters including modal frequencies, damping, and polarization that can be monitored alongside crack aperture and may respond differently to environmental forcings and complex failure evolution. We analyzed data from continuous ambient vibration monitoring and concomitant crack aperture measurements at the Courthouse Mesa instability, a large toppling sandstone slab in Utah, USA. Three years of data revealed crack aperture increases of 2–4 mm/year with no clearly detectable irreversible changes in modal parameters, including frequency. Annually, frequency and displacement varied by 29% and 19% of the mean, respectively, with average and maximum daily frequency fluctuations of 6.5% and 16%, respectively. These reversible cyclic changes were primarily temperature-driven, but annually, frequency was in-phase with temperature whereas crack aperture lagged temperature changes by ~37 days. Polarization and damping also varied seasonally but were less strongly correlated with temperature. Conceptual 3D finite element modeling demonstrated consistent frequency decreases associated with crack propagation but variable changes in crack aperture measured at a single point; i.e., crack propagation did not always result in increased crack opening but always generated a resonance frequency decrease. Taken together, our data suggest a possible thermal wedging-ratcheting mechanism at the Courthouse Mesa instability, where annual thermoelastic crack closure is impeded by debris infill but the absence of downward crack propagation during the monitoring period is evidenced by no permanent resonance frequency changes. Our study demonstrates that combined seismic resonance and crack aperture data provide an improved description of rock slope instability behavior, supporting refined characterization and monitoring of changes accompanying progressive failure.

58 GEOSCIENCES↗
DOE OSTI2023

Boundary layer moisture variability at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic observatory during marine conditions

Abstract. Boundary layer moisture variability at the Eastern North Atlantic (ENA) site during marine conditions is examined at monthly and daily timescales using 5 years of ground-based observations and output from the European Center for Medium range Weather Forecast (ECMWF) reanalysis model. The annual cycle of the mixed-layer total water budgets is presented to estimate the relative contribution of large-scale advection, local moisture tendency, entrainment, and precipitation to balance the moistening due to surface latent heat flux on monthly timescales. When marine conditions prevail, advection of colder and dry air from the north acts as an important moisture sink (∼ 50 % of the overall budget) during fall and winter driving the seasonality of the budget. Entrainment and precipitation contribute to the drying of the boundary layer (∼ 25 % and ∼ 15 % respectively), and the local change in moisture contributes to a small residual part. On a daily temporal scale, moist and dry mesoscale columns of vapor (∼ 10 km) are analyzed during 10 selected days of precipitating stratocumulus clouds. Adjacent moist and dry columns present distinct mesoscale features that are strongly correlated with clouds and precipitation. Dry columns adjacent to moist columns have more frequent and stronger downdrafts immediately below the cloud base. Moist columns have more frequent updrafts, stronger cloud-top cooling, and higher liquid water path and precipitation compared to the dry columns. This study highlights the complex interaction between large-scale and local processes controlling the boundary layer moisture and the importance of spatial distribution of vapor to support convection and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of spatial resolution on WRF v3.8.1 simulated meteorology over the central Himalaya

The sensitive ecosystem of the central Himalayan (CH) region, which is experiencing enhanced stress from anthropogenic forcing, requires adequate atmospheric observations and an improved representation of the Himalaya in the models. However, the accuracy of atmospheric models remains limited in this region due to highly complex mountainous topography. This article delineates the effects of spatial resolution on the modeled meteorology and dynamics over the CH by utilizing the Weather Research and Forecasting (WRF) model extensively evaluated against the Ganges Valley Aerosol Experiment (GVAX) observations during the summer monsoon. The WRF simulation is performed over a domain (d01) encompassing northern India at 15 km x 15 km resolution and two nests (d02 at 5 km x 5 km and d03 at 1 km x 1 km) centered over the CH, with boundary conditions from the respective parent domains. WRF simulations reveal higher variability in meteorology, e.g., relative humidity (RH = 70.3%-96.1%) and wind speed (WS = 1.1-4.2m s(-1)), compared to the ERA-Interim reanalysis (RH D 80.0%-85.0%, WS D 1.2-2.3m s(-1)) over northern India owing to the higher resolution. WRF-simulated temporal evolution of meteorological variables is found to agree with balloon-borne measurements, with stronger correlations aloft (r = 0.44-0.92) than those in the lower troposphere (r = 0.18-0.48). The model overestimates temperature (warm bias by 2.8 degrees C) and underestimates RH (dry bias by 6.4 %) at the surface in d01. Model results show a significant improvement in d03 (P = 827.6 hPa, T = 19.8 degrees C, RH = 92.3 %), closer to the GVAX observations (P = 801.4 hPa, T = 19.5 degrees C, RH = 94.7 %). Interpolating the output from the coarser domains (d01, d02) to the altitude of the station reduces the biases in pressure and temperature; however, it suppresses the diurnal variations, highlighting the importance of well-resolved terrain (d03). Temporal variations in near-surface P, T, and RH are also reproduced by WRF in d03 to an extent (r>0:5). A sensitivity simulation incorporating the feedback from the nested domain demonstrates the improvement in simulated P, T, and RH over the CH. Our study shows that the WRF model setup at finer spatial resolution can significantly reduce the biases in simulated meteorology, and such an improved representation of the CH can be adopted through domain feedback into regional-scale simulations. Interestingly, WRF simulates a dominant easterly wind component at 1 km x 1 km resolution (d03), which is missing in the coarse simulations; however, the frequency of southeasterlies remains underestimated. The model simulation implementing a highresolution (3 s) topography input (SRTM) improved the prediction of wind directions; nevertheless, further improvements are required to better reproduce the observed local-scale dynamics over the CH.

Singh, Jaydeep↗
DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

FY25 MOOSE Usability Improvements: 3D Meshing Capabilities, Initiation of Geometry Support for Monte Carlo Tools, and Enhancement of MOOSE/Workbench User Input Interactions

Usability improvements have been made to MOOSE and Workbench in FY25 to enhance usability and user workflows. Assorted enhancement have been made to MOOSE’s intrinsic meshing capabilities in order to enable more flexible and complex meshing of nuclear reactor systems, in particular for 3D applications. Mesh generators have been added to perform operations such as batch mesh generation, surface mesh generation, and creation of 3D transition layers. These mesh generation capabilities make it much easier to generate high quality non-extruded 3D meshes. Additionally, work to integrate Monte Carlo reactor physics simulations into MOOSE-based multi-physics workflows has reached another milestone with the implementation of the Constructive Solid Geometry (CSG) base framework. This framework lays the foundation for mesh generators to offer the user a generic CSG output option (as opposed to a finite element mesh). To support users, workshop on the MOOSE Reactor Module was delivered which featured hands-on examples using the NEAMS Workbench on INL’s High Performance Computing system. Recent updates to the NEAMS Workbench, WASP, and the MOOSE language server have introduced several improvements aimed at making MOOSE-based simulation setup and input management faster, more accurate, and easier to use. Key capabilities that have been added include multi-tab-stop autocompletion, visual input diagnostics, developer-directed data visualizations, upgraded ParaView integration, and Workspace-level file tracking. Together, these changes make it easier for users to build, validate, and manage complex MOOSE-based simulation models — especially those involving reusable components, included files, and datasets. The improvements are designed to save time, reduce input errors, and help users get to a successful simulation run faster, with more confidence in the results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Formation and dissociation of shear-induced high-energy dislocations: insight from molecular dynamics simulations

Solid-phase processing (SPP) allows one to create complex microstructures, not achievable via thermal processing alone. The resulting structures exhibit a rich palette of defects, both thermal and non-thermal, including defect substructures, such as dislocation networks. It is essential to understand the mechanisms of deformation and defect structure formation to guide SPP towards achieving desired microstructures and material properties. In this study, large-scale molecular dynamics simulations are used to investigate the effects of inhomogeneous strain distribution, that mimics deformation conditions of tribological tests, on the evolution of defects under severe shear deformation in polycrystalline Al. Analysis of defect nucleation and reaction pathways reveals that strong geometric constraints suppress the nucleation and slide of low energy dislocation 1/2$\langle110\rangle${111} but promote the nucleation and slide of high energy dislocations, such as [$1\bar1$0] (001) and 1/2 [$1\bar1$$\bar2$]($1\bar1$1). Further, a rough contact surface, characteristic to tribological tests, imposes an inhomogeneous stress field leading to inhomogeneous defect substructures due to location-dependent activation of slip systems. The results suggest that high-energy dislocations can dominate the evolution of grain structures in highly constrained environments, which should be considered in modeling plastic deformation and grain refinement during SPP.

42 ENGINEERING↗
DOE OSTI2024

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Coastal-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) Science Plan

Understanding the mechanisms governing the urban atmospheric environment is critical for informing urban populations regarding the impacts of climate change and associated mitigation and adaptation measures. Earth system (climate and weather) models have not yet been adapted to provide accurate predictions of climate and weather variability within cities, nor do they provide well-tested representations of the impacts of urban systems on the atmospheric environment. These limitations are largely due to limited field data available for testing and development of these models. We will deploy the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s first Mobile Facility (AMF1) to the mid-Atlantic region surrounding the city of Baltimore for the Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE). This deployment will create a four-node regional atmospheric observatory network including Baltimore and its three primary surrounding environments – rural, urban, and bay. CoURAGE investigators will study the interactions among the Earth’s surface, the atmospheric boundary layer, aerosols and atmospheric composition, clouds, radiation, and precipitation at each site, and examine how the spatial gradients across the region interact to create the climate conditions in Baltimore. This study will determine the degree to which Baltimore’s atmospheric environment depends on interactive feedbacks in the atmospheric system and the degree to which conditions in Baltimore depend on the surrounding environment. Some topics of interest include how urban land management exacerbates heat waves, the impact of regional mesoscale winds (nocturnal jet, bay breeze) on urban air pollution and cloud cover, and the impact of the urban heat island and aerosol production on heavy precipitation events. Understanding this integrated coast-urban-rural system quantitatively and with good accuracy and precision is critical to informing climate adaptation and mitigation efforts in the city of Baltimore. The understanding gained should be applicable to many similar coastal, mid-latitude urban centers. Another important objective of CoURAGE is to improve the representation of the climate of coastal cities in Earth systems models (ESMs). CoURAGE investigators will use the observations to test current ESMs, identify weaknesses and work towards improved simulations of this complex environment. The ARM core facility will be deployed in the city of Baltimore, complementing the Baltimore Social-Environmental Collaborative (BSEC), a DOE urban integrated field laboratory (UIFL). Ancillary sites will be deployed to rural Maryland northwest of Baltimore, and to the southern end of Kent Island within Chesapeake Bay. The fourth node will be a long-term atmospheric observatory operated in Beltsville, Maryland by Howard University and the Maryland Department of the Environment. Measurements will be conducted for one year, starting in December of 2024. There will be two intensive operational periods (IOPs), one in summer and one in winter, when the ancillary sites will be enhanced with additional balloon launches, tethered balloon system (TBS) operation, and added atmospheric composition measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Microkinetic simulation and fitting of the temperature programmed reaction of methanol on CeO 2 (111): H 2 and H 2 O + V production

The kinetics and mechanism for temperature programmed reaction following adsorption of an adsorbate can be better understood by simulation and fitting with comparison to experiment. A case study is presented here for the chemistry of adsorption of methanol on a CeO 2 (111) surface followed by heating. The gas products observed are CH 3 OH, CH 2 O, H 2 , H 2 O, CO, CO 2 . At low temperatures (< 500 K), there is formation of H 2 and H 2 O, where the H 2 O formation is accompanied by lattice oxygen vacancy (V) formation and is thus important in determining the selectivity towards different products. Microkinetic modeling was performed using a recently published method for fitting to gain mechanistic knowledge of the H 2 and H 2 O + V formation at < 500 K. In the kinetic models used here, most of the H 2 and H 2 O + V V formation can be explained by a mechanism in which a metastable state of hydrogen on the surface (H*) acts as an intermediate. Two possibilities were investigated for the source of the metastable H* intermediate: H* from CH bond breaking of methoxies, or promotion of H+ to H* via electron transfer from ionic methoxies absorbed in oxygen vacancies (CH 3 O - /V). From this study, we consider the latter to be more likely at < 500 K. For the H 2 O formation, it was found to be critical that H 2 O cannot dissociate directly on oxygen vacancies. Finally, catalytic chemistry was observed in simulations, including catalytic formation of oxygen vacancies. Various features of the experimental results were reproduced, including methoxies being the major carbon containing species on the surface at < 500 K.

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

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗