Search NASA⌕ Search

SEARCH · Search NASA

Results for “exit time”

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 127 records · Page 7

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗

Disturbance of hibernating bats due to researchers entering caves to conduct hibernacula surveys

Estimating population changes of bats is important for their conservation. Population estimates of hibernating bats are often calculated by researchers entering hibernacula to count bats; however, the disturbance caused by these surveys can cause bats to arouse unnaturally, fly, and lose body mass. We conducted 17 hibernacula surveys in 9 caves from 2013 to 2018 and used acoustic detectors to document cave-exiting bats the night following our surveys. We predicted that cave-exiting flights (i.e., bats flying out and then back into caves) of Townsend’s big-eared bats (Corynorhinus townsendii) and western small-footed myotis (Myotis ciliolabrum) would be higher the night following hibernacula surveys than on nights following no surveys. Those two species, however, did not fly out of caves more than predicted the night following 82% of surveys. Nonetheless, the activity of bats flying out of caves following surveys was related to a disturbance factor (i.e., number of researchers × total time in a cave). We produced a parsimonious model for predicting the probability of Townsend’s big-eared bats flying out of caves as a function of disturbance factor and ambient temperature. That model can be used to help biologists plan for the number of researchers, and the length of time those individuals are in a cave to minimize disturbing bats.

59 BASIC BIOLOGICAL SCIENCES↗

Operation of a Modular 3D-Pixelated Liquid Argon Time-Projection Chamber in a Neutrino Beam

The 2x2 Demonstrator, a prototype for the Deep Underground Neutrino Experiment (DUNE) liquid argon (LAr) Near Detector, was exposed to the Neutrinos from the Main Injector (NuMI) neutrino beam at Fermi National Accelerator Laboratory (Fermilab). This detector is a prototype of a new modular design for a liquid argon time-projection chamber (LArTPC), comprising a two-by-two array of four modules, each further segmented into two optically isolated LArTPCs. The 2x2 Demonstrator features a number of pioneering technologies, including a low-profile resistive field shell to establish drift fields, native 3D ionization pixelated imaging, and a high-coverage dielectric light readout system. The 2.4-tonne active mass detector is flanked upstream and downstream by supplemental solid-scintillator tracking planes, repurposed from the MINERvA experiment, which track ionizing particles exiting the argon volume. The antineutrino beam data collected by the detector over a 4.5 day period in 2024 include over 30,000 neutrino interactions in the LAr active volume—the first neutrino interactions reported by a DUNE detector prototype. During its physics-quality run, the 2x2 Demonstrator operated at a nominal drift field of 500 V/cm and maintained good LAr purity, with a stable electron lifetime of approximately 1.25 ms. This paper describes the detector and supporting systems, summarizes the installation and commissioning, and presents the initial validation of collected NuMI beam and off-beam self-triggers. In addition, it highlights observed interactions in the detector volume, including candidate muon antineutrino events.

47 OTHER INSTRUMENTATION↗

Inside out: Measuring the effect of wood anatomy on the efflux and assimilation of xylem–transported CO 2

Carbon dioxide concentrations (including aqueous CO 2 , carbonic acid, bicarbonate, and carbonate) in woody tissues can be up to 750 times higher than atmospheric CO 2 concentrations, ranging from <1% to >26% vs atmospheric ~0.04% (Fig. 1) (Teskey et al. 2008). CO 2 formed through respiration is generally assumed to diffuse to the atmosphere from tissues adjacent to where it is produced. Here, this CO 2 buildups in the stem due to the diffusional barriers in woody and bark tissues. CO 2 in the stem has three fates: 1) it can be refixed for photosynthesis, 2) it can be used for anaplerotic reactions, or 3) it can exit the plant either adjacent to where it was produced (radial diffusion) or in an area remote from its point of origin (xylem-transported CO 2 ). Not accounting for assimilation of xylem-transported CO 2 may result in underestimating total plant photosynthesis. Alternatively, overlooking the transport of CO 2 away from its point of origin complicates the estimation of respiration in stems, branches, or even leaves. For example, CO 2 efflux from the stem may not solely represent stem respiration, but it may represent CO 2 generated through respiration in areas remote from the point of efflux (Stutz et al., 2017). Additionally, wood anatomy and branching architecture likely influence how much and how far xylem-transported CO 2 travels (Fig. 2). Few studies have combined measurements of both the efflux of xylem-transported CO 2 from trees together and the assimilation of xylem-transported CO 2 . Excitingly in this issue of Plant, Cell & Environment, Salomón et al. (pp. ) demonstrate the importance of xylem-transported CO 2 in plants with different wood anatomies by comparing the amount of xylem-transported CO 2 used for photosynthesis to the amount effluxed.

54 ENVIRONMENTAL SCIENCES↗

Particle Scale Impacts on Deconstruction Energy of Pine Residues

The goal of this Case Study was to quantify the impacts of variable moisture and ash on hammer mill throughput and energy consumption and on generation of fines that are not able to be fed to conversion, as compared to a status quo Base Case system. Also considered was convertible carbon content (minimum carbon specification) and maximum ash content and the delivered feedstock cost impacts of not being able to feed residue not meeting both specifications to the conversion reactor. Laboratory data on the impacts of input particle size and moisture content on the exit particle size were received from FCIC Subtask 5.2: Preprocessing, High Temperature Conversion Preprocessing from their single particle impact population balance modeling study (Tiasha Bhattacharjee, INL). Additional throughput and energy consumption data were obtained from FCIC Subtask 5.2 (Jordan Klinger, INL) for the same grinder with a 6 mm screen in place. These data were utilized to develop the necessary response surface equations to perform throughput analysis using discrete event simulation. Because the ash contents in the separated fines had not been analyzed in the laboratory at the time of the model runs, we chose to assume that the ash distributed proportionally with total mass into the overs and unders in the disk screen following grinding.

energy consumption↗

A reduced-order modeling of a tubular solar reactor for long duration thermochemical energy storage

The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.

14 SOLAR ENERGY↗

Optimizing Ammonia Separation via Reactive Absorption for Sustainable Ammonia Synthesis

Metal halide salts such as magnesium chloride have been shown to be promising candidates for ammonia storage materials for energy storage and agriculture applications due to their ability to incorporate several moles of ammonia per mole of salt. Ammonia exiting a synthesis reactor can be separated from nitrogen and hydrogen by absorption into magnesium chloride. Such an absorption can be more complete and hotter than separation via ammonia condensation, the current standard in the Haber–Bosch process. Here, we discuss the optimal conditions for the cyclic uptake and release of ammonia from the supported magnesium chloride absorbents. An automated system was designed for measuring the nonequilibrium working capacity of the absorbent, as well as the impact of important operating conditions such as absorption and desorption temperature, pressure, and desorption time. Measurements of absorption and desorption kinetics provide insight into the mechanisms involved. The temperatures and pressures during absorption and desorption were designed to use minimal energy input to maximize the uptake and release of ammonia within a reasonable amount of time. In a laboratory-scale bed, absorption has a small unused bed length, so it is largely independent of mass transfer; it is dominated by how fast ammonia is fed to the bed. Yet, desorption is restricted both by the speed of heating the bed and by diffusion out of the absorbent. These measurements provide guidelines for ammonia separations and cycling sorbent materials on a larger scale.

10 SYNTHETIC FUELS↗

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗

Multi-sensor fusion for robust autonomous flight in indoor and outdoor environments with a rotorcraft micro-aerial vehicle (MAV)

The subject matter described herein includes a modular and extensible approach to integrate noisy measurements from multiple heterogeneous sensors that yield either absolute or relative observations at different and varying time intervals, and to provide smooth and globally consistent estimates of position in real time for autonomous flight. We describe the development of the algorithms and software architecture for a new 1.9 kg MAV platform equipped with an IMU, laser scanner, stereo cameras, pressure altimeter, magnetometer, and a GPS receiver, in which the state estimation and control are performed onboard on an Intel NUC 3rd generation i3 processor. We illustrate the robustness of our framework in large-scale, indoor-outdoor autonomous aerial navigation experiments involving traversals of over 440 meters at average speeds of 1.5 m/s with winds around 10 mph while entering and exiting buildings.

42 ENGINEERING↗

Ion Manipulation from Liquid Xe to Vacuum: Ba-Tagging for a nEXO Upgrade and Future 0 νββ Experiments

Neutrinoless double beta decay (0𝜈𝛽𝛽) provides a way to probe physics beyond the Standard Model of particle physics. The upcoming nEXO experiment will search for 0𝜈𝛽𝛽 decay in 136 Xe with a projected half-life sensitivity exceeding 10 28 years at the 90% confidence level using a liquid xenon (LXe) Time Projection Chamber (TPC) filled with 5 tonnes of Xe enriched to ∼90% in the 𝛽𝛽-decaying isotope 136 Xe. In parallel, a potential future upgrade to nEXO is being investigated with the aim to further suppress radioactive backgrounds and to confirm 𝛽𝛽-decay events. This technique, known as Ba-tagging, comprises extracting and identifying the 𝛽𝛽-decay daughter 136 Ba ion. One tagging approach being pursued involves extracting a small volume of LXe in the vicinity of a potential 𝛽𝛽-decay using a capillary tube and facilitating a liquid-to-gas phase transition by heating the capillary exit. The Ba ion is then separated from the accompanying Xe gas using a radio-frequency (RF) carpet and RF funnel, conclusively identifying the ion as 136 Ba via laser-fluorescence spectroscopy and mass spectrometry. Simultaneously, an accelerator-driven Ba ion source is being developed to validate and optimize this technique. The motivation for the project, the development of the different aspects, along with the current status and results, are discussed here.

a-tagging↗

Single-nucleon knockout cross sections for reactions producing resonance states at or beyond the drip line

Single-nucleon knockout cross sections from fast secondary beams of the proton-drip-line nuclei 9 C, 13 O, and 17 Ne on a 9 Be target have been studied with emphasis on the production of resonance states. These states were identified by their invariant mass, and resonances with two-, three-, and five-body exit channels were examined. The measured cross sections for these states were compared with eikonal-model predictions using shell-model or variational Monte Carlo spectroscopic factors. The experimental yields were found to be suppressed relative to the model predictions, especially when a well-bound neutron or proton is removed. This suppression exceeds that found systematically in measured inclusive cross sections to particle-bound final states. In neutron knockout from 9 C and 13 O projectiles, this suppression of the unbound ground-state residuals yield is a factor of two to three times larger than that found in the bound final-state studies. Modifications to the structure of these systems due to coupling of the shell-model configurations to the continuum is expected to contribute to this extra suppression, especially when the final state is a near-threshold resonance. Furthermore, other considerations including the role of nuclear dynamics may be required to explain all the observed trends.

6 ≤ A ≤ 19↗

Overall Cooling Effectiveness of Internally Cooled Additively Manufactured Blades

Power generation gas turbine blades must utilize highly effective cooling schemes to withstand moderate heat loads for long operating times. Prior work indicates that gas turbines in the 5-10 MW range primarily utilize internal cooling to prevent further increase of aerodynamic losses that result from film cooling, but few studies focus solely on internal cooling. In this study, three additively manufactured internal cooling designs were tested in the 50 percent span section of the National Experimental Turbine (NExT) blade: a baseline cooling design with pin-fin arrays and rib turbulators, a double wall design with impingement jets and serpentine channels, and an incremental impingement design. The designs were tested in a high-speed linear cascade at exit Reynolds numbers ranging from 700,000 to 1,400,000 at an exit Mach of 0.8. In addition, freestream turbulence and build angle effects on cooling performance were studied.

gas turbine blades↗

Fuel Stratification Influence on NO x Emission in a Premixed Axial Reacting Jet-in-Crossflow at High Pressure

Three reacting jet-in-crossflow (JiC) methane/air flames were numerically investigated in a lean axially staged combustor at a pressure of five atmospheres. A detailed chemistry Star-CCM+ computational fluid dynamics (CFD) model was used with 53 species considered and the result of turbulence-governed finite-rate modeling was validated with in-house experimental data. An optically accessible test section features three side windows, allowing local flow and flame analysis with particle image velocimetry (PIV) and CH* chemiluminescence as well as pressure, temperature, and species exit measurements. The research objective was to predict and verify NO x formation of the premixed 12.7 mm axial jet. Three headend temperature levels were investigated along with three premixed jets at lean (φ Jet = 0.75), near-stoichiometric (φ Jet = 1.07), and rich (φ Jet = 1.78) axial fuel line equivalence ratio. Based on the matching exit emission concentration, global emission benefits were investigated by adjustment of the fuel stratification. The perfectly premixed methane/air flames of this study were shown to ignite at the lee-side of the jet. For the elevated headend temperature level T = 1800 K, the flame extended beyond the windward jet trajectory and caused high axial NO production. For industry application, a firing temperature of 1920 K was achieved with a NO x optimized fuel split of 25%, combining a lean headend (φ Headend = 0.61) with a rich (φ Jet = 1.78) jet equivalence ratio. As a result, this operating point allowed minimization of the combustor residence time at temperatures above 1700 K as well as combustion in a compact flame at the jet lee-side along the counter rotating vortex pair.

03 NATURAL GAS↗

ORBIT simulations of fast ion power loads on the wall of the Divertor Tokamak Test

Abstract Neutral beam injection in tokamaks produces a population of fast ions, which interact with 3D magnetic fields in a variety of ways, often resulting in energetic particle losses in very short times to the wall. Careful design of neutral beams and active control of error fields helps to keep these losses to a minimum. Nevertheless, past experience in tokamaks in the 1980s and detailed simulations for future machines, such as International Thermonuclear Experimental Reactor, suggest that very localized fast ion losses (‘hot spots’) can be present, even if the overall losses are low. In this paper, we discuss this issue in the Divertor Tokamak Test (DTT) project, and we show that in the standard single-null full-power scenario of the DTT, fast ions produce two hot spots, corresponding approximately to the beam injection and exit toroidal angles: the former being mainly due to prompt losses/passing particles, while the latter is due to trapped ions. However, the maximum power load in these spots is of the order ∼ 63 kW m −2 , below the tolerance of plasma-facing components of the machine.

Physics↗

MACHINE LEARNING-ENABLED PREDICTION OF TRANSIENT INJECTION MAP IN AUTOMOTIVE INJECTORS WITH UNCERTAINTY QUANTIFICATION

Accurate prediction of injection profiles is a critical aspect of linking injector operation with engine performance and emissions. However, highly resolved injector simulations can take one to two weeks of wall-clock time, which is incompatible with engine design cycles with desired turnaround times of less than a day. Hence, it is important to reduce the time-to-solution of the internal flow simulations by several orders of magnitude to make it compatible with engine simulations. This work demonstrates a data-driven approach for tackling the computational overhead of injector simulations, whereby the transient injection profiles are emulated for a side-oriented, single-hole diesel injector using a Bayesian machine-learning framework. First, an interpretable Bayesian learning strategy was employed to understand the effect of design parameters on the total void fraction field. Then, autoencoders are utilized for efficient dimensionality reduction of the flowfields. Gaussian process models are finally used to predict the spatiotemporal void fraction field at the injector exit for unknown operating conditions. The Gaussian process models produce principled uncertainty estimates associated with the emulated flowfields, which provide the engine designer with valuable information of where the data-driven predictions can be trusted in the design space. The Bayesian flowfield predictions are compared with the corresponding predictions from a deep neural network, which has been transfer-learned from static needle simulations from a previous work by the authors. The emulation framework can predict the void fraction field at the exit of the orifice within a few seconds, thus achieving a speed-up factor of up to 38 x 10(6) over the traditional simulation-based approach of generating transient injection maps.

machine learning↗

Effects of frequency and pulse width on electron density, hydrogen peroxide generation, and perfluorooctanoic acid mineralization in a nanosecond pulsed discharge gas-liquid plasma reactor

Plasma electron density and temperature were characterized in a continuous flowing gas-liquid film reactor with argon carrier gas by time-resolved optical emission spectroscopy. The plasma parameters were studied as a function of time for varying pulse widths and frequencies. Pulse frequency was varied between 1 and 10 kHz at 16 kV (input voltage) and 40 ns (pulse width) using an Eagle Harbor Technologies, Inc. (EHT) power supply and 5–100 kHz using an Airity Technologies, LLC (AT) power supply. The pulse width was varied between 40 and 200 ns at 16 kV, 2 kHz with the EHT power supply. Optimal frequencies of 5 and 20 kHz were observed for peak electron density with EHT and AT power supplies, respectively. The peak electron density increased with increasing pulse width between 40 and 200 ns using the EHT power supply. Hydrogen peroxide exiting the reactor in the liquid phase increased with discharge power irrespective of the power supply or pulse parameters. Mineralization of 12.5, 50, and 200 ppm perfluorooctanoic acid (PFOA) dissolved in DI water to fluoride (F - ) correlated to the peak electron density. Glycerol, a liquid-phase hydroxyl radical scavenger, depleted hydrogen peroxide but did not affect PFOA mineralization. CO, a gas-phase hydroxyl radical scavenger, led to a reduction in the formation of F - production, suggesting hydroxyl radicals in the gas-liquid film play a necessary, but not singular, role in mineralization of PFOA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-Fouling, Low Cost Electrolytic Coagulation & Disinfection for Treating Flowback and Produced Water for Reuse

Executive Summary This Final Report is composed of two major sections. The first section presents experimental results obtained in the laboratory during Budget Period 1. The second section presents results from the field test conducted during Budget Period 2. Laboratory Results This research investigated a novel electrochemical process for producing a ferric iron coagulant for use in treating flowback and produced water from hydraulic fracturing and oil production operations. The treatment system improves the effectiveness and lowers the cost of coagulation processes using Fe3+ as the coagulant. The electrolytic coagulant generation (ECG) system uses an electrochemical cell to produce acid and base from oilfield brine solutions. The acid is used to dissolve scrap iron to provide a Fe3+ coagulating agent. The base is used to neutralize the treated water. Compared to conventional electrocoagulation (EC), the main advantage of the ECG system is an order of magnitude lower cost for the source of iron. The second advantage over conventional EC is that it can deliver Fe3+ doses greater than 1 mM, since it is not limited by the amount of dissolved oxygen in the water required to oxidize ferrous to ferric iron. The capital costs for conventional EC and the ECG system are similar, but the operational costs for the ECG system are an order of magnitude lower than conventional EC. The combined costs for iron and electrical energy for treating 1 m3 of FPW with 1 mM Fe3+ is estimated to be $0.87 for conventional EC, and $0.087 for the ECG system. The estimated all-in cost for treating FPW with a 2 mM Fe3+ dose is $0.73/m3 ($0.12/bbl). Field Test Results The field test was conducted at the Paul Foster Central Tank Battery (CTB) in Lea County, New Mexico from November 10, 2022 through December 15, 2022. The feed water to the system was produced water from the Tatanka 1H formation. After approximately two weeks of testing, the initial batch of produced water had been treated and no untreated produced water was available. Thus, after this time, the feed water to the system consisted of previously treated water (i.e., recycled water). The recycled water had nearly all colloidal particles removed, and had a much lower alkalinity due to precipitation of carbonate minerals during the first pass through the system. Although the recycled water was not an ideal test solution due to its low particulate concentrations, its lower alkalinity did allow us to identify the main problem with the treatment system. The main problem with the treatment system was caused by the high alkalinity of the initial feed water (5.4 meq/L) that consumed a significant fraction of the electrochemically generated acid. This resulted in pH values exiting the iron contact tank that were too high to dissolve enough iron to effectively treat the produced water. Tests performed with recycled water with lower alkalinity did not have this issue, and dissolved iron concentrations greater than 20 mM could be achieved. One consistent observation was that effluent water from the iron contact tank was always free of particulates, even when fed with circumneutral solutions. This suggests that there is no need to dissolve high concentrations of iron if all the water to be treated is passed through the scrap iron canister. In this case, dissolved O2 and hypochlorous acid can promote sufficient iron corrosion to provide an effective coagulating agent – even in neutral pH water. This solves the problem resulting from highly alkaline produced water. Modifications to the design of the treatment system were made based on the field test results. These modifications will add minimal additional cost, and were tested in bench-scale laboratory experiments. In short-term testing, the modified treatment process was able to remove colloidal FeS particulates to levels below detection. Long-term, steady state testing will be required to determine whether the modified process is suitable for commercial treatment systems.

54 ENVIRONMENTAL SCIENCES↗

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗