SEARCH · Search NASA

Results for “integrated gradients”

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.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴ speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously previously impossible with conventional tools while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴× speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously—previously impossible with conventional tools—while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Integrating biomechanics in evolutionary studies, with examples from the amphidromous goby model system

The functional capacities of animals are a primary factor determining survival in nature. In this context, understanding the biomechanical performance of animals can provide insight into diverse aspects of their biology, ranging from ecological distributions across habitat gradients to the evolutionary diversification of lineages. To survive and reproduce in the face of environmental pressures, animals must perform a wide range of tasks, some of which entail tradeoffs between competing demands. Moreover, the demands encountered by animals can change through ontogeny as they grow, sexually mature or migrate across environmental gradients. To understand how mechanisms that underlie functional performance contribute to survival and diversification across challenging and variable habitats, we have pursued diverse studies of the comparative biomechanics of amphidromous goby fishes across functional requirements ranging from prey capture and fast-start swimming to adhesion and waterfall climbing. The pan-tropical distribution of these fishes has provided opportunities for repeated testing of evolutionary hypotheses. By synthesizing data from the lab and field, across approaches spanning high-speed kinematics, selection trials, suction pressure recordings, mechanical property testing, muscle fiber-type measurements and physical modeling of bioinspired designs, we have clarified how multiple axes of variation in biomechanical performance associate with the ecological and evolutionary diversity of these fishes. Here our studies of how these fishes meet both common and extreme functional demands add new, complementary perspectives to frameworks developed from other systems, and illustrate how integrating knowledge of the mechanical underpinnings of diverse aspects of performance can give critical insights into ecological and evolutionary questions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Progress on the reevaluation and validation of the n+233U neutron cross sections

The set of 233 U resonance parameters of the ENDF/B-VIII.0 nuclear data library was adopted from the previous ENDF/B-VII.1 evaluation using the external levels to update the thermal values. Adoption of IAEA 2017 thermal standards ( σ f = 533.0 ± 2.2 b, σ c = 44.9 ± 0.9 b, and ν ‾ tot = 2.487 ± 0.011 ) and of the IAEA-recommended thermal-neutron induced prompt fission neutron spectrum (PFNS) with average PFNS energy of 2.030 ± 0.013 MeV requires a re-evaluation of 233 U neutron cross sections in the resolved resonance region. A newly produced evaluation is being tested on benchmarks carefully selected from the Handbook of International Criticality Safety Benchmark Experiments (ICSBEP) which are highly sensitive to 233 U data. An important goal of this work was to eliminate the strong negative gradient of the calculated effective multiplication factors with respect to the epithermal fission fraction observed in the validation of the ENDF/B-VIII.0 library for those assemblies. A significant improvement in integral performance of critical 233 U solutions is observed for the newly proposed evaluation. Further work addressing the fast neutron region is needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Rapid data acquisition and machine learning-assisted composition design of functionally graded alloys via wire arc additive manufacturing

Abstract The lack of high-quality datasets in materials science hinders artificial intelligence (AI)-driven alloy design. To address this challenge, wire arc additive manufacturing (WAAM) was employed to fabricate graded alloys, generating extensive data for machine learning (ML)-assisted property prediction. ML models were developed using high-throughput experiments, computational models, and genetic algorithm to optimize feature selection, successfully predicting hardness and porosity. The ML model demonstrated its efficacy by designing a gradient alloy with enhanced properties. However, scaling up revealed uncertainties in tensile property and porosity due to differences in size and thermal conditions between the designed alloy build and the gradient print used to construct the ML model. This underscores the need for uncertainty quantification and process optimization in WAAM-driven alloy design. Our work advances AI-integrated additive manufacturing, offering a rapid approach to exploring process–structure–property relationships and accelerating materials development.

Wang, Xin↗

DOE OSTI2024

NMPC for Mode-Switching Operation of Reversible Solid Oxide Cell Systems

Solid oxide cells (SOCs) are a promising dual-mode technology that generates hydrogen through high-temperature water electrolysis and generates power through a fuel cell reaction that consumes hydrogen. Reversible operation of SOCs requires a transition between these two modes for hydrogen production setpoints as the demand and price of electricity fluctuate. Moreover, a well-functioning control system is important to avoid cell degradation during mode-switching operation. In this work, we apply nonlinear model predictive control (NMPC) to an SOC module and supporting equipment and compare NMPC performance to classical proportional integral (PI) control strategies, while ramping between the modes of hydrogen and power production. While both control methods provide similar performance in many metrics, NMPC significantly reduces cell thermal gradients and curvatures (mixed spatial-temporal partial derivatives) during mode switching. A dynamic process flowsheet of the reversible SOC system was developed in the open-source, equation-based IDAES modeling framework. Our IDAES dynamic simulation results show that NMPC can ramp the SOC system between hydrogen and power production targets within short mode-switching times. Moreover, NMPC can comply with operating limits in the SOC system more effectively than PI, and only NMPC can directly enforce user-specified limits for mixed spatial-temporal partial derivatives of temperature. This allows for management of the trade-off be-tween operating efficiency and cell degradation, which is dependent on these temperature curvatures.

Li, Mingrui↗

DOE OSTI2026

Climate-Driven Divergence in Biophysical and Economic Impacts of Agrivoltaics

Increasing global demands for food and energy necessitate innovative land-use solutions. Agrivoltaics, colocating solar photovoltaics with agriculture, shows promise, but its widespread adoption faces complex biophysical and economic trade-offs in a changing climate. Here, we develop an integrated biophysical-economic modeling framework to quantify how agrivoltaics affect biophysical and economic impacts across the Midwestern United States under both current and project climate conditions. We find strong regional divergences driven by climate gradients. In the humid eastern Midwest, solar panel shading limits photosynthesis, leading to reduced yields (maize -24%; soybean -16%) and lower farmers' profitability (maize -16%; soybean -2%) compared to conventional agriculture. Conversely, in the semiarid western region, shading alleviates heat and water stress, moderating yield reductions for maize (-12%) and even boosting soybean yields (+6%), resulting in improved economic returns (-6% for maize; +9% for soybean), for a scenario with 33% photovoltaic ground coverage ratio. Although agrivoltaics generate substantial electrical energy across all regions, high upfront installation costs challenge solar developers compared to standalone solar photovoltaics. However, our analysis identifies “win-win” opportunities where soybean-based agrivoltaics in the semiarid region produce economic benefits for both farmers and solar developers, highlighting the necessity for region-specific designs tailored to local climate conditions. Critically, future climate projections indicate eastward expansion of semiarid conditions, broadening areas where agrivoltaics can mitigate crop yield penalties (even boosting yield) and improve overall profitability, especially under high-emission scenarios. The results provide a mechanistic and economically integrated understanding essential for developing evidence-based and region-specific strategies to scale agrivoltaics in a changing climate.

14 SOLAR ENERGY↗

DOE OSTI2025

Accelerating iterative ptychography with an integrated neural network

Electron ptychography is a powerful and versatile tool for high-resolution and dose-efficient imaging. Iterative reconstruction algorithms are powerful but also computationally expensive due to their relative complexity and the many hyperparameters that must be optimised. Gradient descent-based iterative ptychography is a popular method, but it may converge slowly when reconstructing low spatial frequencies. Here, in this work, we present a method for accelerating a gradient descent-based iterative reconstruction algorithm by training a neural network (NN) that is applied in the reconstruction loop. The NN works in Fourier space and selectively boosts low spatial frequencies, thus enabling faster convergence in a manner similar to accelerated gradient descent algorithms. We discuss the difficulties that arise when incorporating a NN into an iterative reconstruction algorithm and show how they can be overcome with iterative training. We apply our method to simulated and experimental data of gold nanoparticles on amorphous carbon and show that we can significantly speed up ptychographic reconstruction of the nanoparticles.

4DSTEM↗

DOE OSTI2024

Data for: Time to anoxia: Observations and predictions of oxygen drawdown following coastal flood events

This dataset contains data and code associated with the paper Patel et al. "Time to anoxia: Observations and predictions of oxygen drawdown following coastal flood events".The coastal terrestrial-aquatic interface (TAI) is a highly dynamic system characterized by strong physical, chemical, and biological gradients. In particular, shifting soil redox conditions and consumption of terminal electron acceptors, due in part to dynamic hydrologic conditions, is a strong driver of carbon availability and transformations across TAIs. However, while redox dynamics are well described, our ability to quantitatively forecast rates of oxic to anoxic shifts in soils with different characteristics and inundation regimes is limited. We integrated field measurements, laboratory incubations, and model simulations to improve mechanistic understanding of oxygen consumption dynamics in coastal soils. Continuous in situ monitoring unexpectedly revealed that flooding caused temporary spikes in subsurface dissolved oxygen followed by rapid consumption in the wetlands. To further investigate these mechanisms in a controlled setting, we performed laboratory incubations using surface and subsurface soils from a TAI gradient (defined here as upland forest to transitional forest to wetland) in Western Lake Erie to measure oxygen consumption rates in TAI soils during flood events. In our experiments, wetland soils reached anoxia the fastest, in ∼ 9 h on average, whereas upland soils turned anoxic in ∼ 18 h. Subsurface upland soils did not turn anoxic even after two weeks of saturation in the lab, and their oxygen consumption patterns suggested carbon and/or nutrient limitation. These results are consistent with in-situ groundwater redox and oxygen measurements in the field, where wetland soils exhibited the highest rates of oxygen consumption along the TAI. Model simulations of oxygen consumption suggested that oxygen consumption had stronger abiotic controls in wetland soils but stronger biotic controls in upland soils, providing a useful framework for future incubation experiments. Microbial activity is a strong driver of oxygen consumption in TAI soils, although it is constrained by the availability of dissolved carbon in subsurface soils.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/COMPASS-DOE/cmps-time_to_anoxia (DOI: 10.5281/zenodo.10815374).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Effect of photon counting shot noise on total internal reflection microscopy

Total internal reflection microscopy (TIRM) measures changes in the distance between a colloidal particle and a transparent substrate by measuring the scattering intensity of the particle illuminated by an evanescent wave. From the distribution of the recorded separation distances, the height-dependent effective potential φ(z) between the colloidal particle and the substrate can be measured. In this work, we show that spatial resolution with which TIRM can measure φ(z) is limited by the photon counting statistics of the scattered laser light. Here, we develop a model to evaluate the effect of photon counting statistics on different potential profiles using Brownian dynamics simulations and experiments. Our results show that the effect of photon counting statistics depends on spatial gradients ∂φ/∂z of the potential, with the result that sharp features tend to be significantly blurred. We further establish the critical role of photon counting statistics and the intensity integration time τ in TIRM measurements, which is a trade-off between narrowing the width of the photon counting distribution and capturing the instantaneous position of the probe particle.

47 OTHER INSTRUMENTATION↗

DOE OSTI2024

Hydra: Computer Vision for Online Data Quality Monitoring

Hydra is a system utilizing computer vision for near real-time data quality monitoring. Currently operational across all of Jefferson Lab’s experimental halls, it reduces the workload of shift takers by autonomously monitoring diagnostic plots during experiments. Hydra uses "off-the-shelf" supervised learning technologies and is supported by a comprehensive MySQL database. To simplify access, web apps have been developed to facilitate both labeling and monitoring of Hydra’s inferences. Hydra can connect with the alarm system and incorporates complete historical tracking, enabling it to identify issues that shift takers could miss. When issues are detected, a natural first question is: "Why does Hydra think there is a problem?" To answer, Hydra employs Gradient-weighted Class Activation Maps (GradCAM) to identify regions of the image that are important for the specific classification. This interpretive layer enhances transparency and trustworthiness, which is essential for integration with experiment workflows and operation. The Hydra system, results, and sociological considerations for deployment will be discussed.

Jeske, Torri↗

DOE OSTI2025

Deep Reinforcement Learning Based Control of Wind Turbines for Fast Frequency Response

In order to fulfill vital auxiliary grid services, such as load regulation, spin and non-spin reserve provision, and frequency support during emergencies, there is often a requirement for certain wind farms to operate in de-loaded modes. Leveraging the swift response capabilities of wind farms, this study demonstrates that reserving power in de-loaded modes can significantly enhance power grid stability and reliability during system contingencies. Controlling wind farms optimally for frequency support is intricate due to the nonlinearity of models and controllers and the complexity of wind farm interactions with power systems. Here, to address this challenge, this paper introduces a novel approach that integrates wind turbines into reinforcement learning-based solutions for frequency response. This innovative methodology utilizes the state-of-the-art reinforcement learning algorithm known as the surrogate-gradient-based evolutionary strategy. The proposed learning-based algorithm provides continuous control of wind farm output to rapidly stabilize system frequency and prevent unnecessary trips of under-frequency load shedding relays. To facilitate efficient training, parallel computing techniques are employed. The proposed methodology is evaluated on a modified IEEE-39 bus system, and simulation results reveal its efficacy in reliably supporting power system frequency and preventing the need for unnecessary load shedding.

Gao, Wei [Argonne National Laboratory (ANL), Argon↗

DOE OSTI2022

Heat conduction in an irregular magnetic field. Part 2. Heat transport as a measure of the effective non-integrable volume

Given the large anisotropy of transport processes in magnetized plasmas, the magnetic field structure can strongly impact heat diffusion: magnetic surfaces and cantori form barriers to transport while chaotic layers and island structures can degrade confinement. When a small but non-zero amount of perpendicular diffusion is included, the structure of the magnetic field becomes less important, allowing pressure gradients to be supported across chaotic regions and island chains. We introduce a metric for the effective volume over which the local parallel diffusion dominates based on the solution to the anisotropic heat diffusion equation. To validate this metric, we consider model fields with a single island chain and a strongly chaotic layer for which analytic predictions of the relative parallel and perpendicular transport can be made. We also analyse critically chaotic fields produced from different sets of perturbations, highlighting the impact of the mode number spectrum on the heat transport. Our results indicate that this metric coincides with the effective volume of non-integrability in the limit $\kappa _{\perp } \rightarrow 0$ , where $\kappa_{\perp}$ is the perpendicular diffusion coefficient. We propose that this metric be used to assess the impact of non-integrability on the heat transport in stellarator equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Binary optimal control by trust-region steepest descent

Abstract We present a trust-region steepest descent method for dynamic optimal control problems with binary-valued integrable control functions. Our method interprets the control function as an indicator function of a measurable set and makes set-valued adjustments derived from the sublevel sets of a topological gradient function. By combining this type of update with a trust-region framework, we are able to show by theoretical argument that our method achieves asymptotic stationarity despite possible discretization errors and truncation errors during step determination. To demonstrate the practical applicability of our method, we solve two optimal control problems constrained by ordinary and partial differential equations, respectively, and one topological optimization problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A real-time distributed solid oxide electrolysis cell (SOEC) model for cyber-physical simulation

System integration and dynamic operability between SOEC and balance-of-plant (BoP) components are major technical challenges before realizing rapid load following of SOEC systems. Cyber-physical simulation (CPS) is a leading-edge digital engineering approach and is regarded as the next step beyond Digital Twins. CPS approach can be used to research SOEC system integration and develop dynamic controls prior to actual pilot testing without using a real SOEC. To seamlessly couple with BoP hardware and access non-observable operational parameters (e.g., local temperature gradient) during transients, a distributed one-dimensional (1D) real-time SOEC model was developed. Its real-time execution was demonstrated for 20 to 640 nodes at the fixed time step of 5 ms. A higher excess air ratio enabled smaller local temperature gradients on SOEC solid materials and faster transients upon current density step change from 0.15 to 0.55 A cm -2 . During the transients, the magnitude of the peak temperature gradient nearly doubled in 10 s from -3.5 to -5.9 °C cm -1 . This represents a significant operating risk that can impact the dynamic operability of SOEC systems. In addition, the local temperature gradient was found to change directions on all nodes in SOEC solid materials, with the greatest impact on the upstream nodes. The SOEC model was also tested at the thermal neutral voltage using actual process air flow parameters as variable model inputs. Variable process air temperatures were found to induce alternating local temperature gradients on SOEC solid materials. These are new operational mechanisms for SOEC degradation relevant for load following operational modes yet distinct from previous reports. To mitigate these unfavorable features, the SOEC can be operated at voltages that are slightly (±20 mV) deviated from the thermal neutral voltage. Here, the corresponding net thermal energy change was less than 1.6% of the electric power consumption. This 1D real-time SOEC model established the basis of cyber-physical simulation of SOEC hybrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE↗

DOE OSTI2021

Suspended Membrane Evaporators Integrating Environmental and Solar Evaporation for Oily Wastewater Purification

Solar-driven evaporation is promising in oily wastewater treatment, in particular for emulsions, but conventional evaporators suffer from pore blocking by residual oil or contamination by volatile oil compounds in the condensed water. In the current research, we develop a suspended membrane evaporator integrating solar evaporation with oil-in-water emulsion separation. The heating and evaporating interface is separated from the rejecting interface to avoid oil escape and improve heat management. A temperature gradient forms on the membrane surface that can promote evaporation performance by combining both solar and environmental evaporation. Such an evaporator achieves a maximum evaporation rate of 1.645 kg/(m 2 ·h) as well as an apparent evaporation efficiency of 111.9%. Moreover, the superhydrophilic and superoleophobic membrane shows excellent oil repellence and emulsion rejection, which can achieve an oil removal efficiency above 98.8% in oil-in-water emulsion separation, and high evaporation rate recovery in cycling tests. Furthermore, a scaled-up membrane evaporator array produces ~8 kg/(m 2 ·d) of clean water from oily wastewater in outdoor experiments, further demonstrating the strong purification performance of this evaporator in oily wastewater treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗