Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large Scale Additive”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Techno Economic Analysis for Production of Renewable Natural Gas and Value-Added Chemicals from Forest Biomass Residues (CRADA Final Report)

West Biofuels - in collaboration with the University of California San Diego, the University of California Davis, the National Laboratory of the Rockies (NLR), the Colorado School of Mines’ Center for Hydrate Research, Placer County Air Pollution Control District, the Sierra Business Council, and the Southern California Gas Company (SoCalGas)-will demonstrate an innovative pathway to convert forest biomass to renewable gas (RG). The project, Production of Pipeline Grade Renewable Natural Gas and Value-Added Chemicals from Forest Biomass Residues, will utilize an existing pilot-scale gasification system and data and lessons learned from a lab-scale catalytic reactor design, develop, and demonstrate an integrated pilot-scale RG process with a scaled-up catalytic reactor and hydrate-based gas separation process. This robust and simple process uses a commercial proven gasifier, a commercially available catalyst that can tolerate gas contaminants, and a gas separation process that simplifies the downstream processing. Current testing has shown that the process can yield a large amount of methane rich renewable gas in addition to a mixture of propanol, ethanol and other higher alcohols and the relative amounts can be controlled by shifting the process conditions. This process is novel and groundbreaking because the product mixture of RG and valuable alcohols byproducts makes RG production economically feasible at $\$$12 per million British thermal units (MMBtu) or less using forest biomass from high hazard zones.

09 BIOMASS FUELS↗
DOE OSTI2026

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗
DOE OSTI2024

The Dual Nature of Entrainment-Mixing Signatures Revealed through Large-Eddy Simulations of a Convection-Cloud Chamber

Abstract Entrainment of subsaturated air into a cloud can influence its optical and microphysical properties in various ways, depending on the droplet evaporation and turbulent mixing time scales. Previous experiments in the Pi convection-cloud chamber have revealed that, given a fixed entrained air property, the mixing of entrained subsaturated air results in complete evaporation of some cloud droplets, with the rest remaining unchanged. This is a signature of inhomogeneous mixing. While comparing the results of entrainment with varying air properties, the mixing signature appears as if the subsaturated air is well mixed with the cloud to evenly reduce the droplets’ size. In other words, taken together, the experiments appear to have the signature of homogeneous mixing. To explore these results in a greater depth, we conduct large-eddy simulations combined with a bin microphysics scheme. Our results reproduce the similar signatures of inhomogeneous and homogeneous mixing, implying that LES can resolve the inhomogeneous mixing when the grid spacing is smaller than the entrained air parcel. Additionally, we observe that increasing the aerosol injection rate enhances the signature of inhomogeneous mixing, while coarser grid spacing diminishes it. Finally, the change in wall fluxes in response to various entrained air properties confirms that the homogeneous signature seen in the analysis of an ensemble of simulations is the result of various equilibrium states. This further strengthens the suggestion that the homogeneous mixing signature found in aircraft observations near the cloud top may result from combining entrainment events of different intensities, possibly caused by various-sized eddies. Significance Statement Large-eddy simulation and size-resolved microphysics can resolve time scales for turbulent mixing and evaporation and, therefore, are well suited for reproducing, extending, and interpreting the entrainment experiment in the Pi convection-cloud chamber. Our simulation results confirm (i) the inhomogeneous mixing signature for an individual entrainment event and (ii) the appearance of homogeneous mixing in an ensemble of entrainment episodes. Furthermore, we demonstrate that the inhomogeneous mixing signature is more pronounced in a polluted cloud, but coarser grid spacing in simulations may compromise the accuracy of this signature. Last, the homogeneous mixing signature results from various equilibrium states established for different entrainment intensities and adjusted wall fluxes, which are challenging to measure experimentally but can be easily analyzed in the simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Experimental Examination of Additively Manufactured Patterns on Structural Nuclear Materials for Digital Image Correlation Strain Measurements

Abstract Background There are a limited number of commercially available sensors for monitoring the deformation of materials in-situ during harsh environment applications, such as those found in the nuclear and aerospace industries. Such sensing devices, including weldable strain gauges, extensometers, and linear variable differential transformers, can be destructive to material surfaces being investigated and typically require relatively large surface areas to attach (> 10 mm in length). Digital image correlation (DIC) is a viable, non-contact alternative to in-situ strain deformation. However, it often requires implementing artificial patterns using splattering techniques, which are difficult to reproduce. Objective Additive manufacturing capabilities offer consistent patterns using programmable fabrication methods. Methods In this work, a variety of small-scale periodic patterns with different geometries were printed directly on structural nuclear materials (i.e., stainless steel and aluminum tensile specimens) using an aerosol jet printer (AJP). Unlike other additive manufacturing techniques, AJP offers the advantage of materials selection. DIC was used to track and correlate strain to alternative measurement methods during cyclic loading, and tensile tests (up to 1100 µɛ) at room temperature. Results The results confirmed AJP has better control of pattern parameters for small fields of view and facilitate the ability of DIC algorithms to adequately process patterns with periodicity. More specifically, the printed 100 μm spaced dot and 150 μm spaced line patterns provided accurate measurements with a maximum error of less than 2% and 4% on aluminum samples when compared to an extensometer and commercially available strain gauges. Conclusion Our results highlight a new pattern fabrication technique that is form factor friendly for digital image correlation in nuclear applications.

Novich, K. A. (ORCID:0000000204466022)↗
DOE OSTI2025

CO 2 Transport Infrastructure Outlook in the United States

Carbon capture and storage (CCS) represents one of the most important methods to mitigate anthropogenic carbon emissions at a large scale, playing a key role in meeting climate change targets (Bui et al., 2018) and for net-zero CO 2 by 2050 scenarios in the United States (Browning et al., 2023). This technology involves capturing CO 2 emissions from industrial processes, transporting them via pipelines, trucks, rails, or ships, and ultimately storing them in underground geological sites, such as saline aquifers or depleted oil reservoirs. Thus, to encourage carbon reduction initiatives, the U.S. Congress enacted the Bipartisan Budget Act in 2018, reforming the 45Q tax credit to benefit operators storing CO 2 in geologic formations (Jones and Sherlock, 2021). Additionally, the 2022 Inflation Reduction Act further expanded these incentives, providing additional support for CCS initiatives (Hackett and Kuehn, 2023). Although numerous studies describe the importance of optimal CO 2 transportation to support the decision-making of CCS projects aligned with the objective of net-zero emissions by 2050 (Abramson and Christensen, 2021; Chen and Pawar, 2023; Greig and Pascale, 2021), further efforts are required to optimize the transport infrastructure for national-scale CCS deployment. Therefore, in this study, we examine three nationwide scenarios with the SimCCS 3.0 tool (Ma et al., 2022, 2023, 2024) along with a novel geospatial splitting approach developed by Velasco-Lozano et al. (Velasco- Lozano et al., 2024a, 2024b). We present optimized pipeline networks that meet the dynamic evolution of annual capture amounts, describing the required total pipeline lengths at each stage as a function of the pipeline diameters. Thus, the cases presented demonstrate the feasibility of CO 2 pipeline infrastructure for large-scale CCS projects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

SynthEsizing Novel H2 Sensors for Operational Resilience in Pipeline Infrastructure (SENSOR) (CRADA Final Report)

Hydrogen (H₂) is gaining attention as a versatile energy carrier with potential applications across industrial processes, power generation, and transportation. However, its practical deployment, particularly in large-scale distribution systems, faces significant infrastructure challenges. Transporting hydrogen through dedicated pipelines or blending it into existing natural gas networks can lead to serious issues, such as leakage due to the small size of hydrogen molecules and material degradation in pipelines through embrittlement. These technical risks raise safety concerns and could limit the integration of hydrogen into current energy infrastructures. Additionally, using hydrogen-enriched gas mixtures in combustion systems like turbines and engines introduces new performance and compatibility challenges that must be resolved before widespread use becomes feasible.

08 HYDROGEN↗
DOE OSTI2026

Materials Studies of Niobium Thin Films for Quantum Circuit Applications: Progress and Challenges

Niobium (Nb) films have emerged as a crucial material in the development of superconducting qubits, which are key components in quantum computing technology. Here, this review provides a comprehensive examination of Nb films from a materials perspective, focusing on their intrinsic properties, fabrication methods/techniques, and their influence on qubit performance, particularly through surface and interface driven loss mechanisms. We discuss the key material properties that are essential for qubit operation. Various deposition techniques for Nb thin films, such as sputtering, evaporation, molecular beam epitaxy, and atomic layer deposition, are explored, alongside their impact on film quality, uniformity, and qubit performance. Additionally, the influence of surface roughness, thin-film thickness, and substrate materials on quantum coherence is analyzed. Challenges such as defects and material degradation in Nb films are reviewed, along with strategies to mitigate these issues. Finally, we present the latest advancements and future directions in Nb film research, including potential improvements to enhance qubit coherence and scalability for large-scale quantum computing systems. Ultimately, a deeper understanding of surface and interface phenomena is essential for pushing the limits of qubit performance and realizing next-generation quantum technologies.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

High-significance detection of correlation between the unresolved gamma-ray background and the large-scale cosmic structure

Our understanding of the γ-ray sky has improved dramatically in the past decade, however, the unresolved γ-ray background (UGRB) still has a potential wealth of information about the faintest γ-ray sources pervading the Universe. Statistical cross-correlations with tracers of cosmic structure can indirectly identify the populations that most characterize the γ-ray background. In this study, we analyze the angular correlation between the γ-ray background and the matter distribution in the Universe as traced by gravitational lensing, leveraging more than a decade of observations from the Fermi-Large Area Telescope (LAT) and 3 years of data from the Dark Energy Survey (DES). We detect a correlation at signal-to-noise ratio of 8.9. Most of the statistical significance comes from large scales, demonstrating, for the first time, that a substantial portion of the UGRB aligns with the mass clustering of the Universe as traced by weak lensing. Blazars provide a plausible explanation for this signal, especially if those contributing to the correlation reside in halos of large mass (∼ 10 14 M ⊙ ) and account for approximately 30–40% of the UGRB above 10 GeV. Additionally, we observe a preference for a curved γ-ray energy spectrum, with a log-parabolic shape being favored over a power-law. We also discuss the possibility of modifications to the blazar model and the inclusion of additional γ-ray sources, such as star-forming galaxies, misalinged active galactic nuclei, or particle dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Laser system design and critical technologies for the NSF OPAL project

The three-year NSF OPAL project is completing the preliminary design of a 2 x 25-PW all-optical parametric chirped-pulse amplification (OPCPA) system capable of delivering 500-J, 20-fs pulses. The design includes two Alpha beamlines that can be co-timed and pointed into two experimental areas for high-intensity research in several configurations. Additionally, a more flexible 2-PW Beta beam can be used in place of one Alpha beam to allow a wider range of focusing geometries and synchronized electron-beam generation for diverse experimental applications. Here, building on the development of the mid-scale MTW-OPAL prototype system, the NSF OPAL project also focuses on advancing essential technologies, such as large-grating fabrication techniques and actively-cooled disk amplifiers.

42 ENGINEERING↗
DOE OSTI2025

Atomic-scale insights into topotactic transformations in an extra-large-pore zeolite using time-resolved 3D electron diffraction

Understanding the atomic-scale structural dynamics of phase transformations is crucial for developing materials and tailoring their properties. However, many materials are obtained as polycrystalline powders with large unit cells and/or complex structures, making it challenging to investigate detailed structural changes using conventional X-ray diffraction techniques. Here we employ time-resolved three-dimensional electron diffraction to reveal the topotactic reactions and transformations that convert the extra-large-pore silicate zeolite ECNU-45 into ECNU-46. ECNU-45 features three-dimensional interconnecting 24 × 10 × 10-ring channels, while ECNU-46 consists of one-dimensional 24-ring channels connected to 10-ring pockets. ECNU-45 and ECNU-46 are both examples of pure silicate zeolites with pore openings larger than 22-ring. Our findings indicate changes at six distinct tetrahedral silicon sites, involving atom displacement, addition and removal of framework atoms through bond breakage and formation. This work presents the synthesis of zeolites and also provides atomic-level insights into the dynamic processes of topotactic reactions. Our results have implications for advancing materials engineering and understanding complex solid-state reactions at an atomic scale.

Inorganic chemistry↗
DOE OSTI2024

CORRLA-RS

The CORRLA-RS package provides a suite of statistical methods for sampling multidimensional distributions and to conduct sensitivity and correlation analysis of large scale data in the Rust programming language. The software provides a unique solution to multidimensional constrained sampling problems utilizing a combination of parallelized Markov Chain Monte Carlo methods and traditional rejection sampling. The sensitivity and correlation analysis methods are backed by a high performance randomized singular value decomposition implementation which enables datasets larger than the random access memory (RAM) size to be analyzed. Additionally, CORRLA-RS implements the active subspace identification method using a KD-Tree and the randomized singular value decomposition acting in concert.

Gurecky, William [Oak Ridge National Laboratory (O↗
DOE OSTI2024

A thermal-driven graupel generation process to explain dry-season convective vigor over the Amazon

Large-eddy simulations (LESs) are conducted for each day of the intensive observation periods (IOPs) of the Green Ocean Amazon (GoAmazon) field campaign to characterize the updrafts and microphysics within deep convective cores while contrasting those properties between Amazon wet and dry seasons. Mean Doppler velocity (V dop ) values simulated using LESs are compared with 2-year measurements from a radar wind profiler (RWP) as viewed by statistical composites separated according to wet- and dry-season conditions. In the observed RWP and simulated LES V dop composites, we find more intense low-level updraft velocity, vigorous graupel generation, and intense surface rain during the dry periods compared with the wet periods. To investigate coupled updraft–microphysical processes further, single-day golden cases are selected from the wet and dry periods to conduct detailed cumulus thermal tracking analysis. Tracking analysis reveals that simulated dry-season environments generate more droplet-loaded low-level thermals than wet-season environments. This tendency correlates with seasonal contrasts in buoyancy and vertical moisture advection profiles in large-scale forcing. Employing a normalized time series of mean thermal microphysics, the simulated cumulus thermals appear to be the primary generator of cloud droplets. When subsequent thermals penetrate the ice crystal layer, droplets within the thermals interact with entrained ice crystals, which enhances riming in the thermals. This appears to be a production pathway of graupel/hail particles within simulated deep convective cores. In addition, less-diluted dry-case thermals tend to be elevated higher, and graupel grows further during sedimentation after spilling out from thermals. Therefore, greater concentrations of low-level moist thermals likely result in more graupel/hail production and associated dry-season convective vigor.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Real-Time event reconstruction for Nuclear Physics Experiments using Artificial Intelligence

Charged track reconstruction is a critical task in nuclear physics experiments, enabling the identification and analysis of particles produced in high-energy collisions. Machine learning (ML) has emerged as a powerful tool for this purpose, addressing the challenges posed by complex detector geometries, high event multiplicities, and noisy data. Traditional methods rely on pattern recognition algorithms like the Kalman filter, but ML techniques, such as neural networks, graph neural networks (GNNs), and recurrent neural networks (RNNs), offer improved accuracy and scalability. By learning from simulated and real detector data, ML models can identify and classify tracks, predict trajectories, and handle ambiguities caused by overlapping or missing hits. Moreover, ML-based approaches can process data in near-real-time, enhancing the efficiency of experiments at large-scale facilities like the Large Hadron Collider (LHC) and Jefferson Lab (JLAB). As detector technologies and computational resources evolve, ML-driven charged track reconstruction continues to push the boundaries of precision and discovery in nuclear physics. In these proceedings, we highlight advancements in charged track identification leveraging Artificial Intelligence within the CLAS12 detector, achieving a notable enhancement in experimental statistics compared to traditional methods. Additionally, we showcase real-time event reconstruction capabilities, including the inference of charged particle properties, such as momentum, direction, and species identification, at speeds matching data acquisition rates. These innovations enable the extraction of physics observables directly from the experiment in real-time.

Gavalian, Gagik (ORCID:0000000267385457)↗
DOE OSTI2026

Identifying Barriers to Solar and Storage Hybrids: Modeled vs. empirical wholesale market value and net-value for co-located solar + storage projects [Slides]

Large-scale (1MW+) co-located solar and battery storage projects are expanding rapidly in the United States, but their realized contribution to the bulk power system remains poorly understood because public project-level operating data are limited. The Lawrence Berkeley National Laboratory estimates the wholesale market value of 280 operational photovoltaic-plus-storage (PV+S) projects across the seven ISOs/RTOs and 19 additional balancing authorities, representing roughly 95% of the U.S. PV+S fleet in 2024. We model optimized hourly dispatch under energy, capacity, and ancillary-service market opportunities and compare the resulting value with standalone PV value, project-specific levelized cost estimates, and empirical operating or revenue data where available.

14 SOLAR ENERGY↗
DOE OSTI2025

Predicting multi-nodal in-nozzle particle interactions in high-viscosity fluid mediums for acoustophoretic direct-ink writing of line-patterned composites

Patterned functional materials offer improved properties (electrical, thermal, etc.) over their bulk counterparts in many applications, including energy storage, flexible electronics, and sensors. However, manufacturing approaches for patterning materials over large areas with features on the order of hundreds of microns or less are limited. Acoustophoresis, which uses acoustic forces to control particle arrangement in a fluid medium, is a pathway to address this challenge. This process is dependent on particle and fluid properties and enables patterning of a broad range of materials. Herein, a model with experimental validation is presented to demonstrate that acoustophoresis can be combined with direct-ink writing (DIW) to fabricate line patterns over large cm-scale areas. An in-nozzle particle interaction model was developed to investigate the impact of processing conditions on multi-nodal acoustophoretic DIW. The model predicts patterned line widths within a factor of two relative to experimental results for a high viscosity case study. Here, the model was used to investigate the impact of frequency, particle loading, particle radius, and acoustic pressure on line width and patterning time, providing critical feedback regarding the processing conditions suitable for a target application. Model results illustrate that frequency has the greatest impact on line patterns: increasing from 1 to 3 MHz resulted in a greater than 65% reduction in line width and a greater than 85% reduction in patterning time. Additionally, experiments were conducted with an alumina-epoxy ink and a ~21 cm 2 area pattern was rastered in ~5.5 minutes, demonstrating a path towards large-area line-patterned composite fabrication.

25 ENERGY STORAGE↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimization of Desalination Systems with Detailed Water Chemistry through Integration of Reaktoro in WaterTAP

Chemistry predictions are critical for an accurate estimation of performance and costs in desalination process models, which allows for the estimation of the value of new technologies and the viability of treating new water sources. Herein, we present how an implicit function formulation can be used to integrate the chemical modeling package, Reaktoro, into the techno-economic assessment and modeling platform, WaterTAP. This approach resolves the critical issues of integrating large-scale thermodynamic models and databases into equation-oriented process models while allowing more flexibility relative to previously presented surrogate-based methods. We describe how this integration into Pyomo and WaterTAP models is implemented and used through the open-source package Reaktoro-PSE . We first validate this integration approach by performing optimization on a previously presented desalination treatment train with softening and acid addition as the pretreatment steps. Then, to demonstrate the value of this approach, we extend the cost-optimization problem to include the simultaneous addition of lime and soda ash for softening, and HCl and H 2 SO 4 in the acidification steps. Finally, we were able to confirm the previously established results that were obtained by using surrogate models and demonstrate that the implicit function approach enables exploration of different feedwater compositions and a larger number of chemicals and their combinations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗