Search NASA⌕ Search

SEARCH · Search NASA

Results for “sensitivity analyses”

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 181 records · Page 10

Simultaneous optimization of simulated moving bed adsorption and distillation for 2,3‐butanediol recovery

Abstract A combined simulated moving bed (SMB) and distillation separation scheme is developed to recover 2,3‐butanediol (BDO) from a dilute fermentation broth. The scheme was integrated into a lignocellulosic biorefinery that produces hydrocarbon fuels from corn stover with BDO as an intermediate. BDO recovery is one of the most challenging processes in this biorefinery; and given the high associated energy duties, direct distillation is considered cost‐prohibitive. An alternative separation is SMB adsorption in nanoporous materials, which can reject 90% of the water and reduce subsequent distillation costs. Rigorous models were used to optimize the SMB and distillation simultaneously. The separation can be added to the biorefinery while keeping the projected minimum fuel selling price (MFSP) below $0.66 USD (US dollars) per liter gasoline‐equivalent ($2.50/GGE, gallon gasoline equivalent). Finally, sensitivity analyses were conducted to assess the effects of cost and lifetime of the adsorbent, titer concentration, and BDO purity.

09 BIOMASS FUELS↗

Influence of the double bond position in combustion chemistry of methyl butene isomers: A shock tube and laser absorption study

Shock tube experiments have been carried out on 2–methyl–1–butene (2M1B), 2–methyl–2–butene (2M2B), and 3–methyl–1–butene (3M1B)—the three isomers of methyl butene compound. Here, carbon monoxide (CO) time–histories and ignition delay times are obtained behind reflected shockwaves over the temperature range of 1350–1630 K and pressures of 8.3–10.5 atm with stoichiometric mixtures of 0.075% fuel in O 2 /Ar. Comparative ignition study reveals that 3M1B ignites significantly faster than the other two isomers, while 2M1B dissociates earlier but ignites later than 2M2B. Possible mechanisms for this behavior are discussed with ignition delay time sensitivity and reaction path analysis. In addition, time–resolved CO measurements are compared with three different reaction mechanisms from the literature. Sensitivity analyses have been carried out to identify important reactions that need attention to accurately predict the chemistry of these isomers. Further investigation into the rates of unimolecular fuel decomposition reactions and C 3 H 3 + O 2 = CH 2 CO + HCO reaction are suggested based on the current investigation.

09 BIOMASS FUELS↗

ReactionMechanismSimulator.jl: A modern approach to chemical kinetic mechanism simulation and analysis

Abstract We present ReactionMechanismSimulator.jl (RMS), a modern differentiable software for the simulation and analysis of chemical kinetic mechanisms, including multiphase systems. RMS has already been applied to problems in combustion, pyrolysis, polymers, pharmaceuticals, catalysis, and electrocatalysis. RMS is written in Julia, making it easy to develop and allowing it to take advantage of Julia's extensive numerical computing ecosystem. In addition to its extensive library of optimized analytic Jacobians, RMS can generate and use Jacobians computed using automatic differentiation and symbolically generated analytic Jacobians. RMS is demonstrated to be faster than Cantera and Chemkin in several benchmarks. RMS also implements an extensive set of features for analyzing chemical mechanisms, including a library of easy‐to‐call plotting functions, molecular structure resolved flux diagram generation, crash analysis, traditional sensitivity analysis, transitory sensitivity analysis, and an automatic mechanism analysis toolkit. RMS implements efficient adjoint and parallel forward sensitivity analyses. We also demonstrate the ease of adding new features to RMS.

Johnson, Matthew S.↗

Structural stability and artificial buckling modes in topology optimization

Abstract This paper demonstrates how a strain energy transition approach can be used to remove artificial buckling modes that often occur in stability constrained topology optimization problems. To simulate the structural response, a nonlinear large deformation hyperelastic simulation is performed, wherein the fundamental load path is traversed using Newton’s method and the critical buckling load levels are estimated by an eigenvalue analysis. The goal of the optimization is to minimize displacement, subject to constraints on the lowest critical buckling loads and maximum volume. The topology optimization problem is regularized via the Helmholtz PDE-filter and the method of moving asymptotes is used to update the design. The stability and sensitivity analyses are outlined in detail. The effectiveness of the energy transition scheme is demonstrated in numerical examples.

Dalklint, Anna (ORCID:0000000346195205)↗

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗

Convolutional neural networks for intra-hour solar forecasting based on sky image sequences

Accurate and timely solar forecasts play an increasingly critical role in power systems. Compared to longer forecasting timescales, very short-term solar forecasting has lagged behind in both research and practice. In this paper, we propose deep convolutional neural networks (CNNs) to provide operational intra-hour (10-minute-ahead to 60-minute-ahead) solar forecasts. We develop two CNN structures inspired by a widely-used CNN architecture. The CNNs are tailored to our solar forecasting regression tasks and rely solely on sky image sequences. Case studies based on six years of data (over 150,000 data points) demonstrate that the best CNN model has forecast skill scores of 20%-39% over the naive persistence of cloudiness benchmark, even at these very short timescales. The CNNs also have consistently superior performance when compared to shallow machine learning models with meteorological predictors, where the improvement averages around 7%. The sensitivity analyses show that the sky image length, resolution, and weather conditions have impacts on the deep learning model accuracy. In our intra-hour problem with specific setups, two sky images with a 10-minute 128 x 128 resolution yield the most accurate forecasts. Current limitations, future work, and deployment challenges and solutions are also discussed.

14 SOLAR ENERGY↗

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗

Optimization based process modeling of an anaerobic membrane bioreactor system: Application to swine wastewater

To maintain current levels of consumption in the economy with the dwindling supply of non-renewable material and energy, alternative resource streams more traditionally viewed as waste streams must be considered. Fermentation of high-strength wastewaters is one such pathway that allows for the recovery of energy, nitrogen, phosphorus, and carbon compounds. Anaerobic membrane bioreactors (AnMBRs) are an emerging technology that allow for the digestion of wastewater in a much smaller footprint than traditional anaerobic digesters. Adoption of this technology into industry has been limited by membrane capital and cleaning costs, but these costs may be offset through the recovery of valuable products. To evaluate the viability of AnMBR technology in the context of swine wastewater treatment, an optimization-based process model built upon Anaerobic Digestion Model No. 1 (ADM1) has been developed. Modeling results show that a swine wastewater stream provides potential for net positive energy generation from the AnMBR system in most cases. Sensitivity analyses around important variables were conducted to determine focus areas for future research into AnMBR technology and evaluate the robustness of the model to microbial variables that may change with different microbial communities.

09 BIOMASS FUELS↗

Technoeconomic analysis and life cycle assessment of purification processes for captured CO 2 streams

Captured carbon dioxide (CO 2 ) streams contain impurities that must be removed to meet specifications for safe transport, storage, and utilization. Among these impurities, oxygen poses challenges due to its high reactivity and potential to cause corrosion, motivating stringent purity limits below 10 ppmv. Building on recent experimental demonstrations of catalytic oxygen removal using hydrogen (H 2 ), carbon monoxide (CO), methanol (CH 3 OH), and methane (CH 4 ) as reducing agents, this study presents a technoeconomic (TEA) and life cycle assessment (LCA) of these four catalytic purification pathways. Process flowsheets were developed and simulated in Aspen Plus for CO 2 streams representative of both low-temperature and high-temperature capture processes, with integrated heat recovery and energy optimization. Results showed that total purification costs were dominated by feedstock procurement and electricity consumption. Among the studied reducing agents, the CH 4 -assisted route achieved the lowest purification cost and highest CO 2 recovery. Sensitivity analyses showed that the H 2 route became competitive at H 2 prices below $\$$0.56/kg to $\$$0.84/kg, depending on the CO 2 feed temperature conditions. In conclusion, environmental impacts were primarily driven by indirect CO 2 emissions from raw material production and utility consumption.

CO2 pipeline specifications↗

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗

Ignition delay time and speciation of dibutyl ether at high pressures

Dibutyl ether (DBE) is a promising biofuel due to its high cetane number (~ 100) and high volumetric energy density (31.6 MJ/L). It could either be used directly in compression ignition engines or blended with other conventional or renewable fuels. Oxidation and pyrolysis kinetics of DBE are not well known, particularly at high pressures. In this work, we have experimentally investigated the chemical kinetics of DBE in three domains: (a) ignition delay time measurements in a rapid compression machine over T = 550–650 K, P = 10, 20, 40 bar, φ = 0.5, 1; (b) ignition delay time measurements in a shock tube over T = 900–1300 K, P = 20, 40 bar, φ = 0.5, 1; (c) laser-based carbon monoxide speciation measurements in a shock tube during DBE pyrolysis and oxidation over T = 1100–1400 K, P = 20 bar. Pressure time-histories measured in RCM experiments exhibited unique 3-stage and 4-stage ignition behavior predominantly at fuel-lean conditions. Experimental data were compared with the predictions of two recent chemical kinetic models of DBE. Sensitivity analyses were carried out to identify key reactions which may have caused the discrepancy between experiments and simulations. We found that the rate of decomposition of DBE may need to be revisited to improve the oxidative and pyrolytic predictions of DBE kinetic model.

09 BIOMASS FUELS↗

New insights into fuel blending effects: Intermolecular chemical kinetic interactions affecting autoignition times and intermediate-temperature heat release

Fuel blending effects on chemically-dominated fuel properties, such as gasoline anti-knock quality, are influenced by fundamental chemical kinetic interactions between the blending agent and the base fuel. Historically, quantification of such interactions has focused on changes to the radical pool, including $\dot{O}$H and HO 2 , while intermolecular interactions pertaining to carbonated, non-fuel-specific intermediates are typically overlooked. In this regard, this work aims to derive new insight into intermolecular chemical kinetic interactions that are intrinsic to fuel blending effects via a case study on blends of 0–30% ethanol (by volume) into FGF-LLNL (a multi-component gasoline surrogate for FACE-F research gasoline) using a rapid compression machine at a diluted/stoichiometric fuel loading, compressed pressure of 40 bar and low- to intermediate-temperature regimes that are representative of boosted SI engine operation. Ethanol blending effects on the intermediate temperature heat release (ITHR) of FGF-LLNL are characterized using experimental measurements, where ethanol is found to promote the extent of ITHR and suppress the transition from ITHR to main ignition. Chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model. Sensitivity analyses on ITHR characteristics further corroborate the ethanol blending effects, and highlight the significant dependence of ITHR on both fuel-specific and non-fuel-specific reactions. An approach allowing comprehensive characterization of the complex intermolecular chemical kinetic interactions between constitutes in a fuel blend is then proposed. Application of the approach to FGF-LLNL/E0–E30 reveals that ethanol perturbs the heat release and autoignition characteristics of FGF-LLNL not only by directly changing the $\dot{O}$H and HO 2 radical pools via fuel-specific reactions, but also through intermolecular interactions where participating intermediates can be produced and consumed by various sub-chemistries. Disabling the intermolecular interactions in carbonated species between ethanol and FGF-LLNL sub-chemistries leads to somewhat slower ITHR evolution and lower ignition reactivity. The role of the individual intermolecular interaction is also characterized using the proposed approach. Lastly, implications of intermolecular interactions for future studies that aim to improve model performance are highlighted, where it is found that further investigations on the core C0-C4 chemistry are needed for developing highly accurate chemistry models for complex fuel blends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automatically generated model for light alkene combustion

Light alkenes are common combustion intermediates for a variety of fuels. Therefore, understanding their oxidation and pyrolysis chemistry is key to building detailed mechanisms for heavier fuels. This work was focused on the development and evaluation of a detailed kinetic mechanism suitable for the combustion of light alkenes up to C 4 without the use of tuned parameters, instead, the parameter values come from first principles or direct measurements. The generated mechanism accurately estimates the laminar burning velocity (S u ) and ignition delay time (IDT) of light alkenes available in the literature, which represent fundamental combustion properties at a wide range of conditions. Because each parameter is thought to have a physically realistic value, not tuned to these measurements, the new model could be used as a sub-mechanism in models for other applications. The reaction network was generated with the open-source Reaction Mechanism Generator (RMG) software. Sensitivity analyses were performed under wide ranges of temperatures and pressures, allowing for the identification of the most impactful species and reactions. Based on these, a comprehensive thermochemistry database, including calculations on 550 molecules performed in this work at the CBS-QB3 level of theory, and a kinetic library, including theoretically-derived reaction rates retrieved from the literature, were built and used in the mechanism generation. The developed mechanism was compared against several existing detailed kinetic mechanisms for ethene, propene, 1-butene, 2-butene, and isobutene. Here, the newly generated model is the most accurate among the ones analyzed, in terms of fractional bias and normalized mean square error. Hence, this new model was used to analyze the chemistry of alkene combustion. Key rate coefficients were compared, to identify the cause of deviations between the models and possible areas for further improvements.

09 BIOMASS FUELS↗

A wide range experimental and kinetic modeling study of the oxidation of 2,3-dimethyl-2-butene: Part 1

2,3-Dimethyl-2-butene (TME) is a potential fuel additive with high research octane number (RON) and octane sensitivity (S), which can improve internal combustion engine performance and efficiency. However, the combustion characteristics of TME have not been comprehensively investigated. Thus, it is essential to study the combustion characteristics of TME and construct a detailed chemical kinetic model to describe its combustion. In this paper, two high-pressure shock tubes and a constant-volume reactor are used to measure ignition delay times and laminar flame speeds of TME oxidation. Further, the ignition delay times were measured at equivalence ratios of 0.5, 1.0, and 2.0 in “air”, at pressures of 5 and 10 bar, in the temperature range of 950 – 1500 K. Flame speeds of the TME/ “air” mixtures were measured at atmospheric pressure, at a temperature of 325 K, for equivalence ratios ranging from 0.78 to 1.31. Two detailed kinetic mechanisms were constructed independently using different methodologies; the KAUST TME mechanism was constructed based on NUIGMech1.1, and the MIT TME mechanism was built using the Reaction Mechanism Generator (RMG). Both mechanisms were used to simulate the experimental results using Chemkin Pro. In the present work, reaction flux and sensitivity analyses were performed using the KAUST mechanism to determine the critical reactions controlling TME oxidation at the conditions studied.

33 ADVANCED PROPULSION SYSTEMS↗

Flame dynamics and kinetic coupling of ammonia and dimethyl-ether in non-premixed cool and warm flames at elevated pressure

Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.

42 ENGINEERING↗

Optimization of battery swapping infrastructure for e-commerce drone delivery

Drone delivery is widely-researched to alter the current e-commerce delivery convention for providing short delivery lead times. Yet, flight range of drones, constrained by the available battery technology, sets a milestone toward realizing the sole-drone delivery. To tackle with the flight range limitation, locating automated battery swapping machines (ABSM) have been proposed and a few studies modeled the problem. Using the ABSMs, drones can be loaded with fully-charged batteries along the route to a demand location. Here, we introduce a mixed-integer nonlinear program to model the problem. The objective of the program is to optimally select ABSM locations, determine the delivery-mode choices (drone-only, truck-only, and mixed delivery) of demand locations, find drone delivery routes, and approximate the baseline requirements for the number of drones and batteries needed. The program minimizes the overall delivery system costs including: ABSM, delivery, drone ownership, battery inventory, and service congestion. A cutting-plane method is developed to find exact solutions in finite iterations. Computational experiments showed that the method quickly yields the optimal solution to instances with less than 60 ABSM candidates and 20 demand locations. A case study shows that the optimal drone delivery infrastructure can save almost 20% cost compared to the conventional truck-only delivery. Sensitivity analyses were conducted to reveal the impact of key parameters in the decision-making and found that a decrease in the ABSM and drone costs highly affect the system cost.

25 ENERGY STORAGE↗

Comparative techno-economic analysis of synthetic renewable natural gas production via reactive CO 2 capture and conversion

Reactive CO 2 capture and conversion (RCC) is an emerging carbon management strategy that integrates CO 2 capture and conversion and avoids intermediate CO 2 purification. In this study, we design an RCC process to capture atmospheric CO 2 and react it with renewable hydrogen to produce synthetic renewable natural gas (SRNG), which serves as a carbon-neutral energy source and a chemical form of long-duration renewable energy storage. We assess the technological potential of RCC through process modeling, techno-economic, carbon footprint, and sensitivity analyses. Our findings demonstrate that RCC offers energy savings and comparable cost to separated capture and conversion processes. The cost is dominated by renewable hydrogen and material replacement cost. SRNG produced via RCC is competitive with existing low-carbon natural gas technologies and presents a promising low-cost option for long-duration energy storage. This work highlights the potential for deploying RCC technologies within a circular carbon economy and the scientific and technical challenges that must be overcome for material and technology developers.

03 NATURAL GAS↗

Cost optimization of low-salt-rejection reverse osmosis

Low-salt-rejection reverse osmosis (LSRRO) is an emerging membrane-based desalination technology for concentrating brines with potentially lower energy consumption and cost than thermally driven processes. We assess the technoeconomic performance of LSRRO by using an optimization model that minimizes the levelized cost of water (LCOW). We develop detailed models for the RO module and other equipment using WaterTAP (Water treatment Technoeconomic Assessment Platform), an open-source simulation and optimization software for water treatment. We improve upon previous LSRRO cost optimization studies by optimizing the salt permeability for each stage, which has significant implications for technoeconomic performance. We present the cost-optimal design, operation, and performance for three cases: 35 g/L TDS feed at 70% recovery, 70 g/L at 55%, and 125 g/L at 35%. The cost-optimal LCOW and specific energy consumption (SEC) for these cases are 0.70, 1.89, and 7.41 $/m 3 and 3.9, 8.4, and 30.4 kWh/m 3 , respectively. For each case, we use sensitivity analyses to provide guidance on key factors that impact economic viability. We also present the cost-optimal LCOW and SEC for feed concentrations between 5 and 220 g/L TDS and recoveries between 30 and 90%. So, we compare LSRRO to other high-salinity desalination technologies and find that it may be cost-competitive below feed concentrations of 125 g/L.

54 ENVIRONMENTAL SCIENCES↗