Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “response surface methods”

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

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗
DOE OSTI2022

Update of the Nuclear Criticality Slide Rule Calculations: Plutonium systems – Delayed Fission Gamma

IRSN (France), LLNL (USA) and ORNL (USA) began a long-term collaboration effort in 2015 to update the nuclear criticality Slide Rule for the emergency response to a nuclear criticality accident. The Slide Rule permits the estimation of neutron and gamma dose rates and integrated doses based upon estimated fission yields, as a function of distance from the fission source, and time after criticality accidents for different critical systems. This paper presents results from the fourth phase of the current update of the Slide Rule project, in which delayed fission-product gamma (DFG) dose rates of unreflected plutonium critical systems were compared by several modern 3D radiation transport codes (MCNP, COG, SCALE), using updated flux-to-dose conversion factors. Dose rates are calculated for fissile material at five moderation ratios (H/Pu), at 1 m above the ground as a function of distance (between 30 cm and 1,200 m) from the external surface of the source to the center of the detector, and for periods between 1 s and 1,000 min after the critical instantaneous event. Further efforts have been devoted to the determination of the delayed gamma source, by comparing the time-dependent energy spectra obtained from several methods. Overall, DFG dose rates calculated by each participant led to consistent results. Extra effort is under way to identify the cause of the remaining differences, by comparing precisely the gamma source, and particularly nuclides inventories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measuring He i Stark Line Shapes in the Laboratory to Examine Differences in Photometric and Spectroscopic DB White Dwarf Masses

Abstract Accurate helium White Dwarf (DB) masses are critical for understanding the star’s evolution. DB masses derived from the spectroscopic and photometric methods are inconsistent. Photometric masses agree better with currently accepted DB evolutionary theories and are mostly consistent across a large range of surface temperatures. Spectroscopic masses rely on untested He i Stark line-shape and Van der Waals broadening predictions, show unexpected surface temperature trends, and are thus viewed as less reliable. To test this conclusion, we present in this paper detailed He i Stark line-shape measurements at conditions relevant to DB atmospheres ( T electron ≈12,000–17,000 K, n electron ≈ 10 17 cm −3 ). We use X-rays from Sandia National Laboratories’ Z -machine to create a uniform ≈120 mm long hydrogen–helium mixture plasma. Van der Waals broadening is negligible at our experimental conditions, allowing us to measure He i Stark profiles only. H β , which has been well-studied in our platform and elsewhere, serves as the n e diagnostic. We find that He i Stark broadening models used in DB analyses are accurate within errors at tested conditions. It therefore seems unlikely that line-shape models are solely responsible for the observed spectroscopic mass trends. Our results should motivate the WD community to further scrutinize the validity of other spectroscopic and photometric input parameters, like atmospheric structure assumptions and convection corrections. These parameters can significantly change the derived DB mass. Identifying potential weaknesses in any input parameters could further our understanding of DBs, help elucidate their evolutionary origins, and strengthen confidence in both spectroscopic and photometric masses.

Schaeuble, M. -A. (ORCID:0000000188695203)↗
DOE OSTI2022

Role of Ocean and Atmosphere Variability in Scale-Dependent Thermodynamic Air-Sea Interactions

This study investigates the influence of oceanic and atmospheric processes in extratropical thermodynamic air-sea interactions resolved by satellite observations (OBS) and by two climate model simulations run with eddy-resolving high-resolution (HR) and eddy-parameterized low-resolution (LR) ocean components. Here, in this work, spectral methods are used to characterize the sea surface temperature (SST) and turbulent heat flux (THF) variability and co-variability over scales between 50 and 10,000 km and 60 days to 80 years in the Pacific Ocean. The relative roles of the ocean and atmosphere are interpreted using a stochastic upper-ocean temperature evolution model forced by noise terms representing intrinsic variability in each medium, defined using climate model data to produce realistic rather than white spectral power density distributions. The analysis of all datasets shows that the atmosphere dominates the SST and THF variability over zonal wavelengths larger than ~2,000–2,500 km. In HR and OBS, ocean processes dominate the variability of both quantities at scales smaller than the atmospheric first internal Rossby radius of deformation (R 1 , ~600–2,000 km) due to a substantial ocean forcing coinciding with a weaker atmospheric modulation of THF (and consequently of SST) than at larger scales. The ocean forcing also induces oscillations in SST and THF with periods ranging from intraseasonal to multidecadal, reflecting a red spectrum response to ocean forcing similar to that driven by atmospheric forcing. Such features are virtually absent in LR due to a weaker ocean forcing relative to HR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Method for producing a crystalline film

A method is provided for coating a surface of a material with a film of porous coordination polymer. A first substrate having a first surface to be coated is positioned in a processing chamber such that the first surface is placed in an opposing relationship to a second surface. The second surface may be provided by a wall of the processing chamber, or in some cases the second surface may be provided by a second substrate to be coated. The first substrate is held such that a gap exists between the first and second surfaces, and the gap is filled with at least one reaction mixture comprising reagents sufficient to form the crystalline film on at least the first surface. A thin gap (e.g., less than 2 mm) between the first and second surfaces is effective for producing a high quality film having a thickness less than 100 μm. Confining the volume of the reaction mixture to a thin layer adjacent the substrate surface significantly reduces problems with sedimentation and concentration control. The size, shape, or average thickness of the gap may be adjusted during formation of the film in response to feedback from at least one film growth monitor.

Britt, David K.↗
DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗
DOE OSTI2023

Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas

Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.

58 GEOSCIENCES↗
DOE OSTI2020

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the Emergence of Human Influence on Surface Air Temperature Changes Over India

Human activities in terms of greenhouse gas (GHG) emissions and aerosols resulting from the combustion of fossil fuels have been shown to have affected the temperature of the Earth on global and continental scales. The surface air temperature (TAS) over India has also been observed to be increasing over the last 100 years. Understanding the underlying causes of regional climate change over India can help in developing appropriate mitigation and adaptation strategies. Differentiating signals of externally forced climate changes from the noise of natural internal variability generally becomes more difficult as spatial scale reduces. Therefore detecting and attributing the influence of external forcings such as GHGs and aerosols is harder at local and regional scales. In this study, we applied a detection and attribution (D&A) method to study annual and seasonal mean TAS over the Indian region. We found that the observed warming over India from 1906 to 2005 cannot be explained by natural climate variability alone. We found that the warming is largely driven by the increase in GHGs, and partially offset by regional anthropogenic emissions of aerosols. These results were confirmed for the shorter 1956–2005 period, but results were sensitive to the choice of observational data set. Furthermore, the changes cannot be explained by internal climate variability or natural external forcings alone, but are compatible with the responses to combined anthropogenic GHG and aerosol forcings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Apparatus and method for producing a crystalline film on a substrate surface

An apparatus and method is provided for coating a surface of a material with a film of porous coordination polymer. A first substrate having a first surface to be coated is positioned in a processing chamber such that the first surface is placed in a substantially opposing relationship to a second surface. In some embodiments, the second surface is provided by a wall of the processing chamber, and in other embodiments the second surface is provided by a second substrate to be coated. The first substrate is held such that a gap exists between the first and second surfaces, and the gap is filled with at least one reaction mixture comprising reagents sufficient to form the crystalline film on at least the first surface. A thin gap (e.g., having a thickness less than 2 mm) between the first and second surfaces is effective for producing a high quality film having a thickness less than 100 μm. Confining the volume of the reaction mixture to a thin layer adjacent the substrate surface significantly reduces problems with sedimentation and concentration control. In some embodiments, the size, shape, or average thickness of the gap is adjusted during formation of the film in response to feedback from at least one film growth monitor.

Britt, David K.↗
DOE OSTI2023

Lime slurry treatment of soils developing on abandoned coal mine spoil: Linking contaminant transport from the micrometer to pedon-scale

Historical and abandoned coal mine spoil continues to generate acid- and metal(loid)-rich porewaters and represents a geographically large and diffuse non-point source of contamination to local watersheds. A potentially inexpensive approach to treat these materials and soils developing on them is through the application of lime slurries, to neutralize acidity and encourage the (co)precipitation of metal(loids) with Fe(III)-(oxy)hydroxides, and potentially other metal-oxides and/or Ca-bearing phases. Here, the efficacy of this approach was evaluated through parallel field application and laboratory-based flow-through column experiments. The field site is Huff Run sub-watershed 25 located in Tuscarawas County, Ohio and was chosen in part because it was previously classified as one of the most highly AMD-impacted sub-watersheds in the region. Two locations with historical spoil were chosen and suction lysimeters were installed at 25 and 75 cm depth to monitor porewater composition on two sides at the base of each pile. Half of each slope received seven lime slurry treatments from June through October of 2017. A suite of aqueous (ICP-OES, IC, and TOC-L) and solid phase geochemical and mineralogical approaches (quantitative SEM-EDS and synchrotron μ-XRF) were used to determine how composition, texture, morphology, and spatial distribution of mineral coatings differ in pre- and post-lime treated soils, and how that impacts the distribution and transport of trace metal(loid)s. Mine spoil porewater at site 1 was slightly less alkaline (pH ranging from 7.04 to 7.37) than at site 2 (ranging from 7.55 to 7.71), and average electrical conductivity values at site 1 (316–405 μS cm —1 ) were slightly lower than at site 2 (358–464 μS cm —1 ), although differences between the sites were not significant. Porewater pH and electrical conductivity in all lysimeters decreased over the course of the field season but there was no obvious response to lime treatment at either site or any depth. At site 1, both treatment and depth were significant factors affecting Ca, K, Ni, SO 4 2— , and DOC concentrations while only treatment effects were significant for dissolved Al and Cu (p < 0.05). For all soils, there were no trends in metal concentration observed over time although DOC and SO 4 2— decreased over the field season. Pedon-scale changes in metal porewater concentrations in response to treatment were linked to micrometer-scale changes in mineral surface coatings; specifically, higher concentrations of Ca, Fe, Mn, and Zn were observed in the coatings and no changes were observed in Fe redox speciation, whereas total S decreased likely due to oxidation of S in coal fragments. In contrast to the field experiment, the column experiments exhibited a much greater response in effluent composition with respect to lime treatment. The untreated columns had approximately an order of magnitude more H + leached over the course of the experiment (p < 0.001) and resulted in greater Ca, Al, Cu, and DIC leached and less Mn, Zn, and SO 4 2— . Soils treated with the lime slurry in the column experiments exhibited larger and thicker secondary Fe-coatings, including the addition of Fe-sulfates. Despite clear trends in the laboratory-based column experiment where the lime-to-soil ratio was higher, the effects were either muted or undetected in the field pilot project, suggesting that a higher application rate of lime in the field is needed to achieve a similar effect. This work provides evidence that a less alkaline lime slurry could be a practical and inexpensive method of treating coal mine spoil-impacted soils and represents an important step in linking laboratory-based remediation studies to implemented field-based studies.

58 GEOSCIENCES↗
DOE OSTI2023

Cybersecurity Platform and Certification Framework Development for Extreme Fast Charging (XFC)-Integrated Charging Ecosystem (Final Project Report)

This report summarizes a pioneering effort in Electric Vehicle charging infrastructure ecosystem cybersecurity requirements, assessment methodologies, functional verification, as well as embodiment of the key technologies in the form of hardware and software tools being made available to the public. EPRI led a team of experts, as well as a stakeholder coalition encompassing all key actors in the EV charging infrastructure ecosystem that includes eXtreme Fast Charging (XFC) equipment (defined as 200kW or above). EV charging infrastructure in the United States is a patchwork of networks that have continued to grow organically and have been designed to serve the charging needs of the EV owners, who are their customers. In doing so, each network provider, as well as their connected entities such as the cloud Electric Vehicle Service Providers or EVSPs, utility back office, utility AMI networks, payment networks, as well as Original Equipment Manufacturer (EV manufacturer) telematics networks, have designed systems that may work well individually, but no single entity is responsible for the entire ecosystem to be secure in terms of data exchange. Furthermore, there is no uniformity in how each actor has implemented the cybersecurity requirements since no system-wide cybersecurity requirements existed prior to this project. The final project report describes the technical approach guided by the EV charging infrastructure cybersecurity working group, convened specifically for this project. The technical approach included definition of requirements at the ecosystem level, treated as a ‘system of systems’, and then passed down to individual systems (EVSE, EV, cloud EVSP, utility, and the payment networks), followed by developing the cybersecurity risk and vulnerability assessment methods, that were later applied to real-world cyber-physical systems at EPRI, ANL, and NREL laboratories, to validate both the process and the results. Finally, in a spotlight over the most vulnerable equipment, which is the EV charge station (AC or DC), the team developed a multi-layer cybersecurity implementation in the embedded domain embodied by the open-source Secure Network Interface Card (SNIC) demonstrating the various ways in which the infrastructure can be secured protecting against the identified attack surfaces. Finally, the entire process of EV infrastructure cybersecurity assessment was encapsulated in the Electric Vehicle Charging Cybersecurity Management (EVC2M) online GUI-based tool, that is expected to be released to the public. The report presents the objectives, the technical approach, the key results, and recommendations for future work.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Total metals, carbon, nitrogen & anion concentration data; Slate River & East River floodplains, Crested Butte, CO; May 2022-October 2022

This data package includes processed and undiluted measurements for metal, total carbon, total nitrogen, and anion concentrations from pore water (groundwater) and surface water samples from the Slate River and East River floodplains of Crested Butte, CO, focus field sites for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and October of 2022. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS (metals):Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA. The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solutions (SPEX Certiprep, Metuchen, NJ). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace elements in water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. Lastly, a suitable internal standard (usually Rh, In, Ga or Ge) is added using on-line addition into the sample line and mixing tee.Analysis by Shimadzu TOC-L (TOC/TN):The TOC-L system is a combustion technique where liquid samples are injected and combusted into CO2 for carbon detection by non-dispersive infrared (NDIR) and NO for detection by chemiluminescence. A calibration curve using five standard solutions between 0.1 and 7 ppm for carbon and 0.05 and 3.5 ppm for nitrogen is made for each type of measurement with a linearity >0.99. All samples, standards, and QC’s are prepared in 24mL scintillation vials that have been baked for 4hrs at 475 Cº and made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.Analysis by Ion Chromatography (Anions):The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with sodium carbonate eluent. A calibration curve using five standard solutions between 5 and 250 umol/L is made with a linearity >0.99. Standards and QC’s are prepared in 15mL polypropylene conical tubes, pipetted along with the samples into 1.5mL polypropylene vials. Dilutions are made using RO water (18.2mΩ). QC’s include a calibration blank check (CCB), continuing calibration check (CCC), and a certified reference material check (CRM). All QC’s are within ±10% error and are run before and after each batch of samples. Samples are diluted and rerun if any measurement concentrations are above the highest standard.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗
DOE OSTI2021

Estimating the Acute Health Impacts of Fire-Originated PM 2.5 Exposure During the 2017 California Wildfires: Sensitivity to Choices of Inputs

Exposure to wildfire smoke increases the risk of respiratory and cardiovascular hospital admissions. Health impact assessments, used to inform decision-making processes, characterize the health impacts of environmental exposures by combining preexisting epidemiological concentration–response functions (CRFs) with estimates of exposure. These two key inputs influence the magnitude and uncertainty of the health impacts estimated, but for wildfire-related impact assessments the extent of their impact is largely unknown. We first estimated the number of respiratory, cardiovascular, and asthma hospital admissions attributable to fire-originated PM 2.5 exposure in central California during the October 2017 wildfires, using Monte Carlo simulations to quantify uncertainty with respect to the exposure and epidemiological inputs. We next conducted sensitivity analyses, comparing four estimates of fire-originated PM 2.5 and two CRFs, wildfire and nonwildfire specific, to understand their impact on the estimation of excess admissions and sources of uncertainty. We estimate the fires accounted for an excess 240 (95% CI: 114, 404) respiratory, 68 (95% CI: -10, 159) cardiovascular, and 45 (95% CI: 18, 81) asthma hospital admissions, with 56% of admissions occurring in the Bay Area. Although differences between impact assessment methods are not statistically significant, the admissions estimates' magnitude is particularly sensitive to the CRF specified while the uncertainty is most sensitive to estimates of fire-originated PM 2.5 . Not accounting for the exposure surface's uncertainty leads to an underestimation of the uncertainty of the health impacts estimated. Employing context-specific CRFs and using accurate exposure estimates that combine multiple data sets generates more certain estimates of the acute health impacts of wildfires.

58 GEOSCIENCES↗
DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

Isotope ratio-based quantification of carbon assimilation highlights the role of plastidial isoprenoid precursor availability in photosynthesis

Abstract Background We report a method to estimate carbon assimilation based on isotope ratio-mass spectrometry (IRMS) of 13 CO 2 labeled plant tissue. Photosynthetic carbon assimilation is the principal experimental observable which integrates important aspects of primary plant metabolism. It is traditionally measured through gas exchange. Despite its centrality in plant research, gas exchange performs poorly with rosette growth habits typical of Arabidopsis thaliana , mutant lines with limited biomass, and accounts poorly for leaf shading. Results IRMS-based carbon assimilation values from plants labeled at different light intensities were compared to those obtained by gas exchange, and the two methods yielded similar values. Using this method, we observed a strong correlation between 13 C content and labeling time (R 2 = 0.999) for 158 wild-type plants labeled for 6 to 42 min. Plants cultivated under different light regimes showed a linear response with respect to carbon assimilation, varying from 7.38 nmol 13 C mg −1 leaf tissue min −1 at 80 PAR to 19.27 nmol 13 C mg −1 leaf tissue min −1 at 500 PAR. We applied this method to examine the link between inhibition of the 2 C -methyl- d -erythritol-4-phosphate (MEP) pathway and suppression of photosynthesis. A significant decrease in carbon assimilation was observed when metabolic activity in the MEP pathway was compromised by mutation or herbicides targeting the MEP pathway. Mutants affected in MEP pathway genes 1-DEOXY- d -XYLULOSE 5-PHOSPHATE SYNTHASE (DXS) or 1-HYDROXY-2-METHYL-2-(E)-BUTENYL 4-DIPHOSPHATE SYNTHASE (HDS) showed assimilation rates 36% and 61% lower than wild type. Similarly, wild type plants treated with the MEP pathway inhibitors clomazone or fosmidomycin showed reductions of 52% and 43%, respectively, while inhibition of the analogous mevalonic acid pathway, which supplies the same isoprenoid intermediates in the cytosol, did not, suggesting inhibition of photosynthesis was specific to disruption of the MEP pathway. Conclusions This method provides an alternative to gas exchange that offers several advantages: resilience to differences in leaf overlap, measurements based on tissue mass rather than leaf surface area, and compatibility with mutant Arabidopsis lines which are not amenable to gas exchange measurements due to low biomass and limited leaf surface area. It is suitable for screening large numbers of replicates simultaneously as well as post-hoc analysis of previously labeled plant tissue and is complementary to downstream detection of isotopic label in targeted metabolite pools.

2C‑methyl‑d‑erythritol 4‑phosphate pathway↗
DOE OSTI2024

Covalent integration of polymers and porous organic frameworks

Covalent integration of polymers and porous organic frameworks (POFs), including metal-organic frameworks (MOFs), covalent organic frameworks (COFs) and hydrogen-bonded organic frameworks (HOFs), represent a promising strategy for overcoming the existing limitations of traditional porous materials. This integration allows for the combination of the advantages of polymers, i.e., flexibility, processability and chemical versatility etc., and the superiority of POFs, like the structural integrity, tunable porosity and the high surface area, creating a type of hybrid materials. These resulting polymer-POF hybrid materials exhibit enhanced mechanical strength, chemical stability and functional diversity, thus opening up new opportunities for applications across a large variety of fields, such as gas separation, catalysis, biomedical applications, environmental remediation and energy storage. In this review, an overview of synthetic routes and strategies on how to covalently integrate different polymers with various POFs is discussed, especially with a particular focus on methods like polymerization within, on and among POF structures. To investigate the unique properties and functions of these resultant hybrid materials, the characterization techniques, including nuclear magnetic resonance spectroscopy (NMR), Fourier transform infrared spectroscopy (FTIR), X-ray diffraction (XRD), thermogravimetric analysis (TGA), transmission electron microscopy (TEM) and scanning electron microscopy (SEM), gas adsorption analysis (BET) and computational modeling and machine learning, are also presented. The ability of polymer-POFs to manipulate the pore environments at the molecular level affords these materials a wide range of applications, providing a versatile platform for future advancements in material science. Looking forward, to fully realize the potential of these hybrid materials, the authors highlight the scalability, green synthesis methods, and potential for stimuli-responsive polymer-POF materials as critical areas for future research.

Hossain, Md Amjad↗