Search NASA⌕ Search

SEARCH · Search NASA

Results for “split injection”

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 55 records · Page 3

Transient Relativistic Plasma Grating to Tailor High-Power Laser Fields, Wakefield Plasma Waves, and Electron Injection

In this work, we show the first experiment of a transverse laser interference for electron injection into the laser plasma accelerators. Simulations show such an injection is different from previous methods, as electrons are trapped into later acceleration buckets other than the leading ones. With optimal plasma tapering, the dephasing limit of such unprecedented electron beams could be potentially increased by an order of magnitude. In simulations, the interference drives a relativistic plasma grating, which triggers the splitting of relativistic-intensity laser pulses and wakefield. Consequently, spatially dual electron beams are accelerated, as also confirmed by the experiment.

43 PARTICLE ACCELERATORS↗

Chemo‐mechanical Alterations Induced From CO 2 Injection in Carbonate‐Cemented Sandstone: An Experimental Study at 71 °C and 29 MPa

Abstract Carbon capture, utilization, and storage may lead to mechanical degradation of the subsurface reservoir from fluid‐rock interaction, which could lead to wellbore instability or reservoir compaction. To better understand potential relationship between mechanical degradation with various carbonate cement textures and compositions in sandstone reservoirs, six flow‐through experiments were conducted. Formation water (TDS = 5,390 mg/L) enriched with CO 2 flowed through two types of Pennsylvanian Morrow B Sandstone: an ankerite‐siderite‐cemented sandstone (disseminated cement texture) and a calcite‐cemented sandstone (poikilotopic cement texture). The experiments produced little change in permeability in the ankerite‐siderite‐cemented sandstone, but permeability increased up to more than 1 order of magnitude in the calcite‐cemented sandstone. Ultrasonic measurements and cylinder‐splitting tests (also known as Brazilian tests) suggested negligible mechanical degradation of the ankerite‐siderite‐cemented sandstone. Variable changes, with significant mechanical degradation in the static moduli, were observed in the calcite‐cemented sandstone. Thus, dissolution of the disseminated ankerite‐siderite cement (0.28–0.30%) had minimal impact on modifying the flow network and the mechanical integrity of the sandstone, whereas dissolution of the poikilotopic calcite cement (0.89–1.13%, quantified with fluid chemistry and visualized with X‐ray microcomputed tomography) impacted the mechanical strength of the sandstone by disconnecting framework grains. With the high water‐to‐rock mass ratios (7.3–8.2) and number of pore volumes (147–675) employed in these experiments, potential risks are most relevant to regions near injection wells. Ultimately, the chemo‐mechanical effects induced by CO 2 injection are strongly influenced by the cement texture and composition and the burial history of the reservoir rock.

Wu, Z.↗

Epitaxial heterointerfacial electron bridge synchronizes oxygen evolution activity and stability on a layered double hydroxide surface

Scalable green hydrogen production via electrocatalytic water splitting is largely restricted by the insufficient activity and stability of oxygen evolution reaction (OER) catalysts at the anode. As a class of the most active OER catalysts in alkaline electrolyzers, the application of layered double hydroxides (LDHs) remains a main challenge owing to the unstable lattice oxygen dissolution due to the dominant lattice oxygen-involving OER mechanism during long-term operation. Herein, we found that using an epitaxial hetero-interfacing nickel hydroxide (namely Ni(OH) 2 ) as an electron bridge between an active FeCo LDH and Ni foam support to form an LDH*/NFO catalyst, the electronic storage capacity around the Fermi level (-0.5 to +0.5 eV, e-D FE ) sharply increases from 0.93 per cell to 1.51 per cell. Subsequently, we demonstrate that this high e-D FE enables ceaseless and fast power injection into the kinetic process of intermediate species conversion and inhibits lattice oxygen dissolution in the active FeCo LDH. Consequently, it demonstrated a low OER overpotential of 246 mV at a current density of 100 mA cm -2 and ultrahigh stability for up to 3500 hours with an ultraslow overpotential increase rate of 9.4 × 10 -3 mV h -1 . Therefore, we developed an epitaxial hetero-interfacial electron bridging strategy to synchronize the activity and stability of available catalysts for scalable green hydrogen production via electrocatalytic water splitting.

25 ENERGY STORAGE↗

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time, momentum, and energy resolved pump-probe tunneling spectroscopy of two-dimensional electron systems

Abstract Real-time probing of electrons can uncover intricate relaxation mechanisms and many-body interactions in strongly correlated materials. Here, we introduce time, momentum, and energy resolved pump-probe tunneling spectroscopy (Tr-MERTS). The method allows the injection of electrons at a particular energy and observation of their subsequent decay in energy-momentum space. Using Tr-MERTS, we visualize electronic decay processes, with lifetimes from tens of nanoseconds to tens of microseconds, in Landau levels formed in a GaAs quantum well. Although most observed features agree with simple energy-relaxation, we discovered a splitting in the nonequilibrium energy spectrum in the vicinity of a ferromagnetic state. An exact diagonalization study suggests that the splitting arises from a maximally spin-polarized state with higher energy than a conventional equilibrium skyrmion. Furthermore, we observe time-dependent relaxation of the splitting, which we attribute to single-flipped spins forming skyrmions. These results establish Tr-MERTS as a powerful tool for studying the properties of a 2DES beyond equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mapping the combustion modes of a dual-fuel compression ignition engine

Compression-ignition (CI) engines can produce higher thermal efficiency (TE) and thus lower carbon dioxide (CO 2 ) emissions than spark-ignition (SI) engines. Unfortunately, the overall fuel economy of CI engine vehicles is limited by their emissions of nitrogen oxides (NO x ) and soot, which must be mitigated with costly, resource- and energy-intensive aftertreatment. NO x and soot could also be mitigated by adding premixed gasoline to complement the conventional, non-premixed direct injection (DI) of diesel fuel in CI engines. Several such “dual-fuel” combustion modes have been introduced in recent years, but these modes are usually studied individually at discrete conditions. This paper introduces a mapping system for dual-fuel CI modes that links together several previously studied modes across a continuous two-dimensional diagram. This system includes the conventional diesel combustion (CDC) and conventional dual-fuel (CDF) modes; the well-explored advanced combustion modes of HCCI, RCCI, PCCI, and PPCI; and a previously discovered but relatively unexplored combustion mode that is herein titled “Piston-split Dual-Fuel Combustion” or PDFC. Tests show that dual-fuel CI engines can simultaneously increase TE and lower NO x and/or soot emissions at high loads through the use of Partial HCCI (PHCCI). At low loads, PHCCI is not possible, but either PDFC or RCCI can be used to further improve NO x and/or soot emissions, albeit at slightly lower TE. These results lead to a “partial dual-fuel” multi-mode strategy of PHCCI at high loads and CDC at low loads, linked together by PDFC. Drive cycle simulations show that this strategy, when tuned to balance NO x and soot reductions, can reduce engine-out CO 2 emissions by about 1% while reducing NO x and soot by about 20% each with respect to CDC. This increases emissions of unburnt hydrocarbons (UHC), still in a treatable range (2.0 g/kWh) but five times as high as CDC, requiring changes in aftertreatment strategy.

Engineering↗

Identification of Recombination Losses in CdSe/CdTe Solar Cells from Spectroscopic and Microscopic Time-Resolved Photoluminescence

Due to the lowest-cost and best reliability, CdTe solar cells are the leading thin-film photovoltaic technology. Increasing open-circuit voltage by reducing recombination represents the most promising path toward further improvements. Analysis is needed to identify limitations that cause efficiency losses. To achieve this goal for Cu-doped CdSe/CdTe solar cells, time-resolved spectroscopy and microscopy are developed and applied. Recombination lifetimes and radiative efficiency identify that defect-mediated recombination is the dominant voltage loss mechanism. When carrier lifetimes are averaged over many crystalline grains, they increase from 180 to 430 ns when Al2O3 is applied to the back contact. The quasi-Fermi-level splitting correspondingly increases from 880–905 to 906–931 mV, indicating a pathway to overcome the long-standing 900 mV voltage limitation. However, the dominant recombination losses are attributed to the absorber bulk. From microscopic carrier lifetime measurements, it is identified that space charge fields due to charged grain boundaries (GBs) lead to recombination in the CdTe absorber bulk. At high injection, GB space charge fields are screened, but that occurs above 1 Sun excitation conditions. Alloying with selenium in the near-interface CdSeTe absorber region reduces GB losses and is identified as one of the factors leading to high radiative and power conversion efficiency.

14 SOLAR ENERGY↗

Active Reservoir Management Engineering and Evaluation Summary – Brine Extraction and Storage Test Phase II Final Report – Volumes 1-3

The Energy & Environmental Research Center (EERC) conducted a multiyear program over two phases to perform a field study of the potential efficacy of formation water extraction during carbon dioxide (CO 2 ) injection for CO 2 storage as a means of improving efficiency and/or effectiveness. Phase I was a 2-year study on potential methodologies for a field study of pressure management, culminating in site selection and experimental and site design. Phase II was a 7-year effort that involved constructing site facilities, conducting experiments, and evaluating collected field data. The experimental portion was split into two collaborative components, one for the active reservoir management (ARM) experiment and a second component to host and demonstrate potential water treatment technologies appropriate for high-saline brine at the project site. The results of the overall project are expected to benefit future saline storage projects through validation of techniques that may reduce stress on sealing formations, increase storage capacity, and reduce area of review. The following report focuses on the Phase II ARM experimental implementation and results. A second companion report summarizes the results of the high-saline brine demonstrations.

54 ENVIRONMENTAL SCIENCES↗

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

Colossal anisotropic absorption of spin currents induced by chirality

The chiral induced spin selectivity (CISS) effect, in which the structural chirality of a material determines the preference for the transmission of electrons with one spin orientation over that of the other, is emerging as a design principle for creating next-generation spintronic devices. CISS implies that the spin preference of chiral structures persists upon injection of pure spin currents and can act as a spin analyzer without the need for a ferromagnet. Here, we report an anomalous spin current absorption in chiral metal oxides that manifests a colossal anisotropic nonlocal Gilbert damping with a maximum-to-minimum ratio of up to 1000%. A twofold symmetry of the damping is shown to result from differential spin transmission and backscattering that arise from chirality-induced spin splitting along the chiral axis. These studies reveal the rich interplay of chirality and spin dynamics and identify how chiral materials can be implemented to direct the transport of spin current.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Charge Transfer Dynamics in Dye-Sensitized Photocatalysts Using Metal Complex Sensitizers with Long-Wavelength Visible Light Absorption Based on Singlet–Triplet Excitation

An Os(II) polypyridyl complex was applied as a photosensitizer in dye-sensitized photocatalyst systems based on Pt-intercalated HCa 2 Nb 3 O 10 and Pt-loaded TiO 2 . The Os(II) complex exhibits a spin-forbidden but partially allowed triplet metal-to-ligand charge transfer ( 3 MLCT) transition, enabling broad visible light absorption up to 800 nm, which surpasses that of conventional Ru(II)-based dyes. Despite its shorter excited-state lifetime compared to Ru(II) complexes, efficient electron injection from the excited Os(II) dye into the semiconductor was confirmed. Under visible-light irradiation, the Os(II)- sensitized photocatalysts showed higher H 2 evolution activity than the Ru(II)-sensitized photocatalysts when sodium ascorbate was used as an electron donor, demonstrating effective utilization of long-wavelength visible light. In contrast, negligible H 2 evolution was observed when NaI was employed as a redox mediator for Z-scheme water splitting. Transient absorption spectroscopy revealed that the lack of activity stemmed from inefficient electron transfer from I − to oxidized Os(II). These findings highlight the importance of selecting appropriate redox mediators to fully exploit long-wavelength dyes for overall water splitting under visible light.

artificial photosynthesis↗

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗

Assessment of Potential Future Demands for Hydrogen in the United States

H2@Scale is a U.S. Department of Energy (DOE) initiative that brings together stakeholders to advance the affordable production, transport, storage, and utilization of hydrogen (H2) as an energy carrier to increase revenue opportunities in multiple energy sectors. The focus of the current work is to characterize the growth potential of diverse hydrogen industries in the United States, given research and development (R&D) advancements in hydrogen technologies. Current and emerging hydrogen production technologies utilize diverse energy sources, including natural gas (NG) reforming, as well as renewable and nuclear power for low-temperature and high-temperature water splitting. The produced hydrogen also enables emerging domestic industries that value conventional and renewable hydrogen as an energy carrier for intermediate and end use. The success of H2@Scale (Figure ES.1) depends not only on hydrogen demand from growing existing markets such as petroleum refining and ammonia (NH 3 ) production, but also on the development of new markets such as metals refining, synthetic fuel (synfuel) and chemical production, biofuels, light-duty (LD) and heavy-duty (HD) hydrogen fuel cell (FC) electric vehicles (FCEVs), and injection into NG pipelines, all of which can significantly increase hydrogen demand relative to current levels of approximately10 million metric tons [MMT] annually. This study focused primarily on five demand sectors — synfuels, upgrading of oil/biomass, NH 3 /fertilizer, metals refining, and hydrogen vehicles (transportation) — along with gas infrastructure. For each sector, a hydrogen demand potential was quantified, along with a “threshold price.” The hydrogen demand potential reflects a practical amount of hydrogen that could be used in that sector, barring economic considerations. The threshold price reflects the price at which the consumer would utilize hydrogen in lieu of an alternative that could supply the same performance.

08 HYDROGEN↗

Visible-light H 2 evolution using dye-sensitized TiO 2 : effects of physicochemical properties of TiO 2 on excited carrier dynamics and activity

Dye-sensitized photocatalysts have emerged as promising materials for solar-driven water splitting due to their ability to utilize visible light, in contrast to conventional wide-band-gap semiconductors. However, the relationship between semiconductor properties and charge carrier dynamics remains insufficiently understood. In this study, we investigated Pt/TiO 2 systems sensitized with a visible-light-absorbing Ru(II) polypyridyl complex (RuP), focusing on how the crystal phase and specific surface area of TiO 2 influence excited carrier dynamics and H 2 evolution activity. To isolate the effects of TiO 2 properties, Pt and RuP loadings were standardized across samples. Emission lifetime analysis showed similarly efficient electron injection from RuP to TiO 2 in all cases, suggesting that injection efficiency does not account for observed differences in activity. Transient absorption measurements revealed that back electron transfer (BET) rates depended strongly on the TiO 2 phase, with anatase and P25 exhibiting slower BET and higher activity for H 2 evolution than rutile. The highest apparent quantum yield for H 2 evolution was 12.0% at 450 nm. Among anatase samples, larger surface areas correlated with higher activity, while smaller-area samples exhibited slower BET rates but still low H 2 evolution activity, implying a role for RuP dye–dye interactions in performance loss. This was further supported by improvements in H 2 evolution activity by lowering RuP loading or adding co-adsorbents. Overall, these results demonstrate that both BET suppression and control over RuP dye aggregation are essential for designing efficient dye-sensitized photocatalytic systems.

Harada, Kakeru [Institute of Science Tokyo (Japan)↗

Integrating the Design of Desalination Technologies into Produced Water Network Optimization

The oil and gas energy sector uses billions of gallons of water for hydraulic fracturing each year to extract oil and gas. The water injected into the ground for fracturing along with naturally occurring formation water from oil wells surfaces back in the form of produced water. Produced water can contain high concentrations of total dissolved solids and is unfit for reuse outside the oil and gas industry without desalination. In semi-arid shale plays, produced water desalination for beneficial reuse could play a crucial role in alleviating water shortages and addressing extreme drought conditions. In this poster we co-optimize the design and operation of desalination technologies along with operational decisions across produced water networks. A multi-period produced water network model with simplified split-fraction-based desalination nodes is developed. Rigorous steady-state desalination mathematical models based on mechanical vapor recompression are developed and embedded at the desalination sites in the network model. An optimal common design across all periods is ensured using global capacity constraints. The solution approach is demonstrated for multi-period planning problems on networks from the PARETO open-source library. Model formulation and challenges associated with scalability are discussed.

Naik, Sakshi↗

Toward Emittance Measurements at 11.7 GHz Short-Pulse High-Gradient RF Gun

A short pulse high gradient RF gun has been recently tested at Argonne Wakefield Accelerator (AWA) facility. The carried-out test showed that the 1,5-cell gun was able to inject 3 MeV, up to 100 pC bunches at room tem-perature being fed by 9 ns up to 300 MW 11.7 GHz puls-es. The cathode field was as high as about 400 MV/m. So high field is aimed to mitigate repealing Coulomb forces substantially. In accordance with simulations the emit-tance could be as low as less than 0.2 mcm. To obtain so low emittance in the experiment, the gun is assumed to be equipped with a downstream linac to be fed from the same power extractor as the gun itself. Here we report design of the RF power distribution system splitting RF power among the gun and the linac, results of low-power tests, and emittance measurement plans for upcoming new experiment at AWA.

43 PARTICLE ACCELERATORS↗

Integrating the Design of Desalination Technologies into Produced Water Network Optimization

The oil and gas energy sector uses billions of gallons of water for hydraulic fracturing each year to extract oil and gas. The water injected into the ground for fracturing along with naturally occur-ring formation water from the oil wells surfaces back in the form of produced water. Produced water can contain high concentrations of total dissolved solids and is unfit for reuse outside the oil and gas industry without desalination. In semi-arid shale plays, produced water desalination for beneficial reuse could play a crucial role in alleviating water shortages and addressing extreme drought conditions. In this paper we co-optimize the design and operation of desalination technologies along with operational decisions across produced water networks. A multi-period produced water network model with simplified split-fraction-based desalination nodes is developed. Rigorous steady-state desalination mathematical models based on mechanical vapor recompression are developed and embedded at the desalination sites in the network model. An optimal common design is ensured across all periods using global capacity constraints. The solution approach is demonstrated for multi-period planning problems on networks from the PARETO open-source library. Model formulation and challenges associated with scalability are discussed.

Naik, Sakshi↗