SEARCH · Search NASA

Results for “Reduced fracture model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Pressures in Tumuli: A Study of Tumuli Formation

Tumuli form via localized inflation in surface lava flows. These domed features have widths of 10-20 m, lengths of 10-150 m, and heights of 1-9 m. The axial fracture exposes a brittle crust overlying a ductilely deformed layer. The total crustal thickness is typically less than lm. Tumuli are observed on both terrestrial and martian lava flow surfaces, and provide insight on the flow formation processes and rates. Past studies have estimated the inflation pressure using a bending model for a circular, thin elastic plate, assuming small deflection (Rossi and Gudmundson, 1996). This formulation results in unrealistic pressures for some tumuli. We thus examine alternative models, including those with different shapes, bending of the ductile crust, large deflection, plastic deformation, and thick plate bending. Using the thickness of the ductile crust in the equations for thin, circular plates reduces most pressures to reasonable values. Alternative plate shapes do not cause a significant reduction in inflation pressure. Although the large deflection equations should be applicable based on the plate thickness to tumuli height ratios, they give even less realistic pressures. Tumuli with unrealistic pressures appear to have exceeded the critical bending moment, and have relatively thick crusts, requiring thick plate bending models.

Hansen, James E.↗

NASA NTRS2012

Analytical Round Robin for Elastic-Plastic Analysis of Surface Cracked Plates: Phase I Results

An analytical round robin for the elastic-plastic analysis of surface cracks in flat plates was conducted with 15 participants. Experimental results from a surface crack tension test in 2219-T8 aluminum plate provided the basis for the inter-laboratory study (ILS). The study proceeded in a blind fashion given that the analysis methodology was not specified to the participants, and key experimental results were withheld. This approach allowed the ILS to serve as a current measure of the state of the art for elastic-plastic fracture mechanics analysis. The analytical results and the associated methodologies were collected for comparison, and sources of variability were studied and isolated. The results of the study revealed that the J-integral analysis methodology using the domain integral method is robust, providing reliable J-integral values without being overly sensitive to modeling details. General modeling choices such as analysis code, model size (mesh density), crack tip meshing, or boundary conditions, were not found to be sources of significant variability. For analyses controlled only by far-field boundary conditions, the greatest source of variability in the J-integral assessment is introduced through the constitutive model. This variability can be substantially reduced by using crack mouth opening displacements to anchor the assessment. Conclusions provide recommendations for analysis standardization.

Wells, D. N.↗

NASA NTRS1978

Fracture and crack growth in orthotropic laminates

An approximate solution is developed for the determination of the interlaminar normal and shear stresses in the vicinity of a crack in a three dimensional composite containing unidirectional linearly elastic fibers in an infinite linearly elastic matrix. In order to reduce the complexity of the formulation, certain assumptions are made as to the physically significant stresses to be retained. These simplifications reduce the partial differential equations of elasticity to differential-difference equations which are tractable using Fourier transform techniques. The potential for damaged or debonded zones to be generated by an embedded crack is discussed, and stress concentration factors for fibers near the crack are given. Detailed comparisons are made between the present solution, the analogous two dimensional problem, and corresponding shear-lag models.

Goree, J. G.↗

NASA NTRS2007

Ultrahigh-Temperature Ceramics

Ultrahigh temperature ceramics (UHTCs) are a class of materials that include the diborides of metals such as hafnium and zirconium. The materials are of interest to NASA for their potential utility as sharp leading edges for hypersonic vehicles. Such an application requires that the materials be capable of operating at temperatures, often in excess of 2,000 C. UHTCs are highly refractory and have high thermal conductivity, an advantage for this application. UHTCs are potentially applicable for other high-temperature processing applications, such as crucibles for molten-metal processing and high-temperature electrodes. UHTCs were first studied in the 1960 s by the U.S. Air Force. NASA s Ames Research Center concentrated on developing materials in the HfB2/SiC family for a leading-edge application. The work focused on developing a process to make uniform monolithic (2-phase) materials, and on the testing and design of these materials. Figure 1 shows arc-jet models made from UHTC materials fabricated at Ames. Figure 2 shows a cone being tested in the arc-jet. Other variations of these materials being investigated elsewhere include zirconium based materials and fiber-reinforced composites. Current UHTC work at Ames covers four broad topics: monoliths, coatings, composites, and processing. The goals include improving the fracture toughness, thermal conductivity and oxidation resistance of monolithic UHTCs and developing oxidation-resistant UHTC coatings for thermal-protection-system substrates through novel coating methods. As part of this effort, researchers are exploring compositions and processing changes that have yielded improvements in properties. Computational materials science and nanotechnology are being explored as approaches to reduce materials development time and improve and tailor properties.

Johnson, Sylvia M.↗

DOE OSTI2020

Development of Low-Cost, Crack-Tolerant Metallization Using Screen Printing

One of the ways to reduce the cost of solar electricity to 3/kWh, thus reaching parity with fossil-fuel-based generation, is to reduce the degradation rate of solar modules and extend their lifetime well beyond 30 years. The extended module lifetime in turn can positively influence the financial model and the bankability of utility-scale PV projects. Today, the highest-riskpriority solar module degradation mechanism is what is known as hot spots, often induced by cell cracks. In order to address this degradation mechanism, we make use of low-cost, multi-walled carbon nanotubes embedded in commercial screen-printable silver pastes. When the carbon nanotubes are properly functionalized and appropriately incorporated into commercial silver pastes, the resulting metal contacts on solar cells, after screen-printing and firing, show exceptional fracture toughness. These composite metal contacts possess increased ductility, electrical gap-bridging capability up to 50 um, and 'self-healing' to regain electrical continuity even after cycles of complete electrical failure under extreme strain.

14 SOLAR ENERGY↗

NASA NTRS2025

Erosion Behavior of Ti-hBN Multifunctional Coatings in A Custom-Made Planetary Test Rig at Extreme Lunar Temperatures

Spacecraft landings and takeoffs on the lunar surface, along with extreme temperature variations between day and night (-196 to 150° C), cause high-velocity dust impacts and erosion, resulting in the premature failure of structures. Ti/2 vol% hBN coatings were deposited using atmospheric (APS) and vacuum plasma spray (VPS) using cryo-milled powder feedstock to protect the structural components. The erosion performance of coatings at extreme lunar temperature regimes (-150 to 150° C) was evaluated in a custom-made planetary erosion test rig (PETR) at low (50 mph) and high impact velocities (250 mph). The mass loss of VPS coatings was reduced by 50% compared to the APS coatings and 40% compared to the Ti6Al4V substrate. The cryogenic temperature induces brittleness in the material, rendering it susceptible to extreme conditions of material loss. The particle impact-deformation behavior was captured using a high-speed camera to study the erosion mechanism. This analysis revealed chipping in substrates and brittle APS coatings, while particles rebounding and embedding were observed in VPS coatings. Energy calculations, aided by particle trajectory tracking from the high-speed camera, have conclusively shown that VPS coatings absorb 5–10% more energy than APS coatings during erosion tests. A modified erosion index was developed incorporating the fracture toughness and temperatures. New erosion models for brittle and ductile target materials are proposed for developing erosion-resistant material systems.

Abhijith Kunneparambil Sukumaran↗

DOE OSTI2024

Planarizing Spalled GaAs(100) Surfaces by MOVPE Growth

III-V photovoltaic devices have demonstrated exceptional performance across various applications, with controlled crystal fracturing, known as controlled spalling, emerging as a promising method to reduce costs by enabling substrate reuse. Spalling GaAs(100) substrates, a commonly used substrate in III-V photovoltaics, results in faceted ridges that must be planarized to grow high-quality photovoltaic devices. Here, in this study, we demonstrate that a GaAs(100) wafer offcut toward [$0\bar{11}$] and spalled toward [$011$] can be efficiently planarized by growing C:GaAs by metal-organic vapor phase epitaxy (MOVPE) on the surface, with up to 95% of the nominally deposited material used to fill the valleys between ridges. We find that reducing the offcut to 2° enhances the planarizing capability of C:GaAs. A surface morphology model indicates that the density of surface dangling bonds significantly influences the growth evolution of undoped GaAs surfaces. In contrast, the model suggests that the effectiveness of C:GaAs as a smoothing layer stems from modifying the atomic surface structure and, consequently, the associated sticking coefficients of the facets, which can alter the evolution of surface morphology. Our findings provide guidelines for the epitaxial planarization of semiconductor surfaces and improve the understanding of MOVPE growth on nonplanar surfaces.

14 SOLAR ENERGY↗

NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

DOE OSTI2025

Property optimized energy absorber for automotive bumpers utilizing multi-material and structural design strategies

This study proposes a novel design for automotive bumper using optimized lattice structures and multi-materials to balance low-speed collision and high-speed pedestrian impact performance. Different blends of 20 % carbon fiber-reinforced acrylonitrile butadiene styrene with thermoplastic polyurethane were used to tailor material properties. The energy absorber features lattice structures with customized mechanical responses, created by varying the incline angle θ from 0 to 180°. We conducted 576 finite element simulations on a half-scale model to optimize energy absorption and stiffness, leading to 66 optimized designs that met both low-speed and high-speed impact criteria. Two sub-scale optimized energy absorbers with different peak forces—both meeting low-speed impact requirements—were 3D printed and validated through drop-weight testing. The one with lower peak stress demonstrated a more compliant response, exhibiting approximately 90 % lower initial peak force and an increase in energy absorption of around 33 % (from 24 J to 32 J). Compared to the baseline triangular lattice, the optimized absorber increased energy absorption by 68 % from (19 J to 32 J) and reduced peak stress by 70 %. It also showed near-complete recovery with minimal fractures, making it suitable for repeated use. This design improves safety while offering a lightweight, durable, and cost-effective bumper system.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Experimental and Computational Characterization of a Modified Sioutas Cascade Impactor for Respirable Radioactive Aerosols

Oak Ridge National Laboratory is collecting and characterizing aerosols released when spent nuclear fuel (SNF) rods are fractured in bending. An aerosol collection system was designed and tested to collect respirable sized (<10 μm aerodynamic diameter [AED]) particulates inside a hot cell facility. The setup is a modified version of the commercially available Sioutas cascade impactor, to which additional stages were added to expand the aerosol collection range from 2.5 to ~15 μm AED. To accommodate the additional stages and specific test conditions, the operating flow rate for aerosol collection was reduced, and testing was conducted by using pressure drop measurements, surrogate dust collection, and particle size characterization. The fluid flow distribution within the cascade and its stages was simulated in STAR-CCM+, and the stage-wise pressure drops obtained using the computational fluid dynamics model were then compared to experimental data. Lagrangian particle simulations were also performed, and stage-wise collection statistics were obtained from the simulation for comparison with the experimental data obtained using SNF-surrogate dust particles. The results provide valuable insights into the stage-wise particle collection characteristics of the modified cascade impactor and can also be used to improve the prediction accuracy of the manufacturer-determined analytical correlations.

aerosol modeling↗

DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗

NASA NTRS1993

Tension fracture of laminates for transport fuselage. Part 2: Large notches

Tests were conducted on over 200 center-crack specimens to evaluate: (a) the tension-fracture performance of candidate materials and laminates for commercial fuselage applications; and (b) the accuracy of several failure criteria in predicting response. Crack lengths of up to 12 inches were considered. Other variables included fiber/matrix combination, layup, lamination manufacturing process, and intraply hybridization. Laminates fabricated using the automated tow-placement process provided significantly higher tension-fracture strengths than nominally identical tape laminates. This confirmed earlier findings for other layups, and possibly relates to a reduced stress concentration resulting from a larger scale of repeatable material inhomogeneity in the tow-placed laminates. Changes in material and layup result in a trade-off between small-notch and large-notch strengths. Toughened resins and 0 deg-dominate layups result in higher small-notch strengths but lower large-notch strengths than brittle resins, 90 deg and 45 deg dominated layups, and intraply S2-glass hybrid material forms. Test results indicate that strength-prediction methods that allow for a reduced order singularity of the crack-tip stress field are more successful at predicting failure over a range of notch sizes than those relying on the classical square-root singularity. The order of singularity required to accurately predict large-notch strength from small-notch data was affected by both material and layup. Measured crack-tip strain distributions were generally higher than those predicted using classical methods. Traditional methods of correcting for finite specimen width were found to be lacking, confirming earlier findings with other specimen geometries. Fracture tests of two stiffened panels, identical except for differing materials, with severed central stiffeners resulted in nearly identical damage progression and failure sequences. Strain-softening laws implemented within finite element models appear attractive to account for load redistribution in configured structure due to damage-induced crack tip softening

Walker, Tom H.↗

NASA NTRS2012

Elastic Plastic Fracture Analysis of an Aluminum COPV Liner

Onboard any space-launch vehicle, composite over-wrapped pressure vessels (COPVs) may be utilized by propulsion or environmental control systems. The failure of a COPV has the potential to be catastrophic, resulting in the loss of vehicle, crew or mission. The latest COPV designs have reduced the wall-thickness of the metallic liner to the point where the material strains plastically during operation. At this time, the only method to determine the damage tolerance lifetime (safe-life) of a plastically responding metallic liner is through full-scale COPV testing. Conducting tests costs substantially more and can be far more time consuming than performing an analysis. As a result of this cost, there is a need to establish a qualifying process through the use of a crack growth analysis tool. This paper will discuss fracture analyses of plastically responding metallic liners in COPVs. Uni-axial strain tests have been completed on laboratory specimens to collect elastic-plastic crack growth data. This data has been modeled with the crack growth analysis tool, NASGRO 6.20 to predict the response of laboratory specimens and subsequently the complexity of a COPV.

Forth, Scott↗

DOE OSTI2023

Arctic Aerosol Sources and Mixing States Field Campaign Report

The Arctic is warming at a faster rate than anywhere else on Earth, with rapidly shrinking sea ice extent transforming the region. Depending on chemical and physical properties, atmospheric aerosols directly scatter and/or absorb radiation, serve as cloud droplet and/or ice crystal nuclei, and/or reduce the reflectiveness of the snow surface, thereby altering the atmospheric energy budget. There is a wide spread in the magnitude of simulated arctic aerosol radiative forcing, and significant differences in aerosol concentration levels and seasonal cycles often exist between models and observations. Increasing local natural and anthropogenic emissions are significant, with uncertain climate impacts due to complex feedbacks. Model evaluations, however, are limited by the dearth of arctic aerosol observations available and an inadequate understanding of arctic aerosol processes. The majority of previous arctic aerosol observations have been made through intensive spring/summer field campaigns, with few intensive measurement studies focused on the fall-winter transition, a period when freeze-up is occurring later and thinning sea ice is resulting in wintertime ice fracturing. Aerosol monitoring at arctic coastal stations has provided knowledge of long-term seasonal trends in aerosol concentrations. The completed ARM field campaign addresses observational and knowledge gaps through detailed aerosol size and chemical composition measurements during the fall-winter transition in the coastal Arctic and through the entire annual cycle in the central Arctic. The observations are improving our understanding of the sources and processes controlling the aerosol population in the rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2006

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗

DOE OSTI2021

The homogenous alternative to biomineralization: Zn- and Mn-rich materials enable sharp organismal “tools” that reduce force requirements

We measured hardness, modulus of elasticity, and, for the first time, loss tangent, energy of fracture, abrasion resistance, and impact resistance of zinc- and manganese-enriched materials from fangs, stings and other “tools” of an ant, spider, scorpion and nereid worm. The mechanical properties of the Zn- and Mn-materials tended to cluster together between plain and biomineralized “tool” materials, with the hardness reaching, and most abrasion resistance values exceeding, those of calcified salmon teeth and crab claws. Atom probe tomography indicated that Zn was distributed homogeneously on a nanometer scale and likely bound as individual atoms to more than ¼ of the protein residues in ant mandibular teeth. This homogeneity appears to enable sharper, more precisely sculpted “tools” than materials with biomineral inclusions do, and also eliminates interfaces with the inclusions that could be susceptible to fracture. Based on contact mechanics and simplified models, we hypothesize that, relative to plain materials, the higher elastic modulus, hardness and abrasion resistance minimize temporary or permanent tool blunting, resulting in a roughly 2/3 reduction in the force, energy, and muscle mass required to initiate puncture of stiff materials, and even greater force reductions when the cumulative effects of abrasion are considered. We suggest that the sharpness-related force reductions lead to significant energy savings, and can also enable organisms, especially smaller ones, to puncture, cut, and grasp objects that would not be accessible with plain or biomineralized “tools”.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Produced water sharing: Improved economics and reduced community impact – A Pennsylvania case study

Here, hydraulic fracturing for oil and gas extraction from unconventional reservoirs is water intensive. Between water sourced for fracturing purposes and water present in rock formations, operators often produce a greater volume of water than oil or gas. Historically, this surplus of produced water has mainly been disposed of via deep injection wells. Rising disposal costs, seasonally limited water availability, and concerns over induced seismicity have incentivized produced water recycling practices, where an operator uses produced water for hydraulic fracturing operations. The logistical challenges associated with produced water recycling have also encouraged operators to adopt ad-hoc water exchange practices, in which competing operators will exchange produced water for mutual cost savings. In this paper, we investigate the potential benefit from systematic produced water exchange among operators in Northeastern Pennsylvania. We leverage PARETO, a free and open-source modeling framework for produced water management optimization, to quantify the benefits of water exchange practices. In an example drawn from FracFocus data, we find that the adoption of systematic water sharing could improve produced water recycling rates from 49.2% to 99%, decreasing operating and trucking costs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗