Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “operator complexity”

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 1,081 records · Page 60

Microtructure and Low-Cycle Fatigue Behavior of Additively Manufactured 316L Stainless Steel at 300°C

Metal additive manufacturing (AM) or 3D printing has the potential to transform the nuclear industry by producing high quality components faster and cheaper, thus enhancing the operating performance of current plants and advanced reactors. Two AM 316L tubes - intended to act as surrogates for complex components where nuclear equipment vendors are more likely to consider AM technologies - were printed using a Renishaw AM400 Laser Powder Bed Fusion (L-PBF) system. The porosity of the as-built material was found to be small, 0.06%. In prior research, the fatigue and corrosion fatigue crack growth rate (CGR) response of the AM specimens in the as-built condition was found to be similar to that expected for conventional alloys, and extremely resistant to stress corrosion cracking (SCC). The present research demonstrates that the low-cycle fatigue behavior of the AM material in air is similar to that of wrought stainless steel at LWR temperature. The data suggest that low porosity levels do not have a significant effect on the fatigue performance of AM materials. Taken together, all data generated to date demonstrate that the use of AM alloys in a nuclear environment is plausible. Additional research on environmental fatigue, crack initiation, and SCC needs to be conducted to provide the performance information needed to qualify AM materials for LWR applications.

36 MATERIALS SCIENCE↗

LLNS Request for Variance From 10 CFR 851 for Construction of LLNL Emergency Operations Center

Lawrence Livermore National Security, LLC (LLNS) is submitting the following variance application requesting the use of the State of California’s Occupational Safety and Health Administration (Cal/OSHA) Regulations in lieu of Title 10 Code of Federal Regulations Part 851 (10 CFR 851, the Rule), Worker Safety and Health Program, Subpart B (Program Requirements), Subpart C (Specific Program Requirements), and Appendix A (Worker Safety and Health Functional Areas), solely with respect to the construction of the new Lawrence Livermore National Laboratory (LLNL) Emergency Operations Center (EOC) facility. The request is in support of the National Nuclear Security Administration (NNSA) “pilot” to streamline the delivery of commercial-like line item construction projects under $50 million. This pilot is aligned with the approved NA-50 safety approach for the projects in the pilot. This pilot is to test and demonstrate innovation leading to a reduced cost of construction for building commercial type structures at Department of Energy (DOE)/NNSA facilities. The variance proposes to permit local/regional construction businesses contracted to Management & Operation contractor operated facilities to bid, operate, execute, and comply with safety and health standards they utilize in their business outside of the DOE Complex. This will permit compliance to standards that are as safe and healthful as those in 10 CFR 851 without confusion, additional cost, and inefficiencies associated with the burden of compliance to an unfamiliar regulatory spectrum.

99 GENERAL AND MISCELLANEOUS↗

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗

BOPTEST as a Platform for Building Controls and Grid-Interactive Buildings Workforce Training

Building automation and controls are becoming increasingly complex with the emergence of Grid Integrated Efficient Buildings (GEBs) as well as new highly efficient sequences of operation and data-driven control schemes. However, there remains a significant gap in hands-on training opportunities for building operators and technicians to gain practical experience with advanced control systems in a low-risk environment. This paper presents BOPTEST (Building Optimization Performance Test) as a suitable platform for workforce training in building controls and GEB technologies. BOPTEST provides a suite of standardized building simulation test cases with a REST API, real-time control interfaces through BACnet, semantic models connecting users to building data, and built-in calculation of control metrics and performance indicators. The platform enables trainees to interact with virtual buildings using industry-standard protocols while learning how to implement and innovate control strategies. The training platform is designed to offer a structured and interactive learning experience for building engineers, helping them effectively develop, learn, and retain skills in fault identification, troubleshooting, and correction. The workflow is divided into three main phases: 1) Setup, 2) Exercise, and 3) Review, each comprising specific activities performed by either the instructor or the student. Initial pilot training sessions have yielded positive feedback from instructors and participants and demonstrates that BOPTEST effectively fills an industry need for a low-risk training resource via simulation of real building control systems, allowing trainees to gain practical experience before working in the field. The platform's ability to provide immediate performance feedback while maintaining familiar industry interfaces makes it particularly suitable for workforce development programs. This work provides a replicable model for leveraging building simulation in control education and training.

Paul, Lazlo↗

AI-Based Optimal Design and Controls Can Greatly Reduce Carbon Emissions and Enhance Resilience in Residential Communities in Cold Climates

Net-zero energy residential communities are crucial for achieving decarbonization goals, but the high-penetration photovoltaic (PV) in those communities is posing challenges to the distribution grid. Traditional design and operation of net-zero communities rely on rule-of-thumb methods and may not work in complex scenarios. Artificial intelligence and machine learning methods can optimally size PV for net-zero energy, identify user preferences and usage patterns, and fully unlock the potential of distributed energy resources to address distribution grid issues.

artificial intelligence↗

SaltStone Wastewater Cement Study Using Isothermal Calorimetry, Standard Concrete Characterization Techniques, and CemGEMS - 25169

Cementitious reagents are used to solidify/stabilize aqueous radioactive, hazardous, and mixed salt solutions, and sludges to meet low-level radioactive waste (LLW) and Resource Conservation and Recovery Act (RCRA) requirements for disposal at Department of Energy (DOE). It is flexible enough to solidify radioactive wastewater saturated in complex species that include but are not limited to radioactive isotopes from the bombardment of neutrons in reactor operation, corrosion products from metallic components, and a variety of soluble organic compounds. [1,2]. Waste form testing typically includes processing or fresh properties, cured properties, compressive strength and hydraulic properties, porosity, density, saturated and unsaturated moisture transport, and leachability of contaminants in the waste form pore solution. Properties are collected over a relatively limited time, typically 28 to 365 days [3]. In addition, changes in the waste form as the result of time and changing conditions are important for concrete engineers to predict overall performance of the forms and potential release of contaminants in the disposal process via unintended filtration into the environment [4]. These predictions are determined/calculated characterizing young waste forms (relative to the standard age of concrete) and are based on transport through soluble ions in pore solutions. Characterization methods include X-ray, SEM, and isothermal calorimetry among other methods used to define the composition, amorphous vs. crystalline nature of the components, and the energetic formation mechanisms for multi-phase mineral systems. [5–7] Isothermal calorimetry is a well standardized technique for cements and concretes and can be used to predict the timing and nature of the hydration reactions.[8] The technique can measure long term energetic

Bustamante, Michael E. [Savannah River National La↗

Exact chiral symmetry with quantum signal processing

We give a quantum signal processing (QSP) algorithm for the overlap fermion Hamiltonian which preserves the Ginsparg-Wilson relation up to a controllable error $ε_e$. Quantum simulations of Dirac fermions with exact chiral symmetry are thus nearly free: applying the overlap Hamiltonian costs only a factor logarithmic in $ε_e$ more than the Wilson-Dirac Hamiltonian. Comparing to domain-wall fermions, a mild overhead is found in circuit complexity while reducing qubit costs. We show how QSP effectively constructs an extra dimension when simulating the overlap operator, illustrating that the scaling of quantum algorithms reflects the deeper physics of overlap fermions arising at the boundary of domain-wall fermions.

Lamm, Henry [Fermilab] (ORCID:0000000330330791)↗

Electrostatic dust remediation for future exploration of the Moon

Dust accumulation is one of the critical issues that must be mitigated on in-situ lunar explorations because an in-situ probe is exposed to small dust particles, which are easily attached to it, during its operations. The Lunar Dust Science Definition Team is organized by the Jet Propulsion Lab/California Institute of Technology through NASA’s Biological and Physical Sciences Division to define key science questions and assess dust remediation techniques. Here, we assess three electrostatic remediation technology concepts: electrostatic dust shield; surface electrostatically collecting dust, later called attractive surface; and electron beam — plasma jet inducing electrostatic dust lofting from a surface. We qualitatively investigate their maturity by defining six operational factors: Time and location; Amount of dust removal; Contamination of target surfaces; Operation duration; Installation; and Safety. In addition to these techniques, we discuss a supporting system that loads dust particles onto a test article to examine dust removal efficiency. Here, the results show that further development increases the maturity of all the technologies. While laboratory and theoretical demonstrations reported whether each technology robustly work on the Moon, which hosts a complex, heterogeneous dust environment, we find that it is still uncertain if this is the case because none has been tested in the lunar environment. Particularly, operation duration and safety are critical to be addressed further on both laboratory and spaceflight scales.

42 ENGINEERING↗

Managing the techno-economic impacts of partial string failure in multistring energy storage systems

The role of energy storage systems (ESSs) is becoming increasingly important for today’s electric power systems. Unavailability of an ESS assigned to critical grid services may cause unwanted disruption of those services and hence, may have a significant techno-economic impact. Like any physical equipment, an ESS is vulnerable to various types of faults. Failure of one or more strings in a multistring ESS does not have to be the cause of shutting down the entire ESS. It can still operate with a partial number of strings and continue providing critical services to the grid, if there are no reliability or safety issues and is acceptable under applicable standards. However, it is important to make sure that the control strategies are adaptable to the changes in ESS capacity caused by failed strings. Also, depending on the previous operation and type of failure, the reallocation of duty cycle burden among available strings could be non-uniform. These complexities suggest that mitigation of the impact of partial failure in multistring ESSs is not trivial and needs careful consideration. This is the topic of this paper. The proposed work investigates the impact of partial failure of a large multistring ESS on the assigned service and develops strategies to adjust the ESS control duty cycles for reducing such impacts. In doing so, the paper proposes a novel two-stage framework that plans for the multistring failure using robust optimization theory and then adjusts in real-time using a rule-based algorithm, based on the real-time information on the power availability of the string. Illustrations provided in this work are based on frequency regulation use-case, which is a common application for many utility-scale ESSs. A 750 kilowatt (kW)/1500 kilowatt-hour (kWh) 3-string ESS is used for the demonstration in this work and the efficacy of the proposed method is demonstrated and compared against methods that do not incorporate string failure in their strategy. In particular, here we show that revenue loss of 93% can be incurred when partial string failure is not included in operation and planning. These losses in revenue are reduced by 60% with the proposed method.

25 ENERGY STORAGE↗

Cryo-EM captures the coordination of asymmetric electron transfer through a di-copper site in DPOR

Enzymes that catalyze long-range electron transfer (ET) reactions often function as higher order complexes that possess two structurally symmetrical halves. The functional advantages for such an architecture remain a mystery. Using cryoelectron microscopy we capture snapshots of the nitrogenase-like dark-operative protochlorophyllide oxidoreductase (DPOR) during substrate binding and turnover. DPOR catalyzes reduction of the C17 = C18 double bond in protochlorophyllide during the dark chlorophyll biosynthetic pathway. DPOR is composed of electron donor (L-protein) and acceptor (NB-protein) component proteins that transiently form a complex in the presence of ATP to facilitate ET. NB-protein is an α 2 β 2 heterotetramer with two structurally identical halves. However, our structures reveal that NB-protein becomes functionally asymmetric upon substrate binding. Asymmetry results in allosteric inhibition of L-protein engagement and ET in one half. Residues that form a conduit for ET are aligned in one half while misaligned in the other. An ATP hydrolysis-coupled conformational switch is triggered once ET is accomplished in one half. These structural changes are then relayed to the other half through a di-nuclear copper center at the tetrameric interface of the NB-protein and leads to activation of ET and substrate reduction. These findings provide a mechanistic blueprint for regulation of long-range electron transfer reactions.

Computational biophysics↗

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗

Compact Multicell Superconducting Crab Cavity for ILC

We propose a novel design of a crab cavity (CC) for the ILC project with low parasitic High Order Modes (HOM) losses and preserving the beam emittance, which is critical for operation with high beam current intensity. Multiple electrodes immersed in the hollow waveguide form a trapped-mode resonator. The transverse components of the electromagnetic field of the trapped dipole mode induce a transverse kick and efficiently deflect charged particles passing through the cavity. We present a scalable design of a superconducting Quasi-waveguide Multicell Resonator (QMiR) seamlessly connected with a beam vacuum chamber. The cavity is completely open at both ends, which significantly reduces the maximum loaded quality factor of the higher order modes (HOM), avoids complex HOM couplers, and thus simplifies the mechanical design of the cavity. The same port is used to feed RF power to the operating mode and to extract the same order modes (SOM). Finally, we estimate the expected cryogenic losses, HOM impedance limits, RF input power required, and frequency tuning for a QMiR cavity designed to operate at 2.6 GHz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comments on the quantum field theory of the Coulomb gas formalism

The holomorphic Coulomb gas formalism, as developed by Feigin-Fuchs, Dotsenko-Fateev and Felder, is a set of rules for computing minimal model observables using free field techniques. We attempt to derive and clarify these rules using standard techniques of quantum field theory. We begin with a careful examination of the timelike linear dilaton. Although the background charge of the model breaks the scalar field’s continuous shift symmetry, the exponential of the action remains invariant under a discrete shift because the background charge is imaginary. Gauging this symmetry makes the dilaton compact and introduces winding modes into the spectrum. One of these winding operators corresponds to the anti-holomorphic completion of the BRST current first introduced by Felder, and the full left/right cohomology of this BRST charge isolates the irreducible representations of the Virasoro algebra within the degenerate Fock space of the linear dilaton. The “supertrace” in the BRST complex reproduces the minimal model partition function and exhibits delicate cancellations between states with both momentum and winding. The model at the radius $R=\sqrt{pp^{\prime }}$ has two marginal operators corresponding to the Dotsenko-Fateev “screening charges”. Deforming by them, we obtain a model that might be called a “BRST quotiented compact timelike Liouville theory”. The Hamiltonian of the zero-mode quantum mechanics of this model is not Hermitian, but it is PT -symmetric and exactly solvable. Its eigenfunctions have support on an infinite number of plane waves, suggesting an infinite reduction in the number of independent states in the full quantum field theory. Applying conformal perturbation theory to the exponential interactions reproduces the Coulomb gas calculations of minimal model correlation functions. In contrast to spacelike Liouville, these “resonance correlators” are finite because the zero mode is compact. We comment on subtleties regarding the reflection operator identification, as well as naive violations of truncation in correlators with multiple reflection operators inserted. This work is part of an attempt to understand the relationship between the JT model of two dimen- sional gravity and the worldsheet description of the (2 , p ) minimal string as suggested by Seiberg and Stanford.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Certifying the Performance of Fixative Technologies under Open Air Demolition Activities for D and D - 20096

Department of Energy (DOE) facilities undergoing deactivation and decommissioning (D and D) activities struggle with safety concerns of residual radioactive contamination leftover after gross decontamination efforts have concluded. The positive effects of implementing fixative technologies for decommissioning and maintenance is known across the DOE complex. However, there is no standardized metric in place that quantifies the performance of these technologies such as incombustible fixative platforms that immobilize the contamination in a solid polymer material under normal operating conditions or when exposed to stressors. The primary focus of this research is to provide the empirical data necessary to properly characterize the effects of implementing fixative technologies on mitigating the release of contamination. In addition, this effort has been used to identify and establish uniform testing protocols to quantify operational parameters of fixating platforms and properly credit such systems under a variety of operational conditions and stressors that can arise during open air demolition efforts (e.g. impact, water, thermal, etc.). The established protocols will be provided to the ASTM International E10.03 Subcommittee on Radiological Protection for Decontamination and Decommissioning of Nuclear Facilities and Components for the potential of formal standardization. Attributes to be addressed during testing are material compromise from impact stress, thermal stress for incombustible fixative materials, and immobilization factors (how well an immobilizing platform is capable of retaining fixated material during thermal and impact stressors). In developing these test methods, a greater understanding of the material's behavior under anticipated and unanticipated events, such as decommissioning activities and contingency events as outlined in the Basis of Interim Operations (BIO), can be successfully characterized. There are several key components of the experimental methodology that is essential in empirically certifying fixative technologies. A surrogate contaminant with a unique signature has been utilized with the implementation of a uniform contamination process of test coupons that is quantifiable and replicable for direct fixative comparisons. A modular test chamber has also been used to apply the various stressors of interest. Collection of released contamination that includes airborne and resettled contaminant particles and analyzed using mass spectrometry. This research directly supports DOE complex-wide concerns in terms of final disposition of nuclear facilities and the considerations made during their time between ceasing facility operation and final disposition. The results of this effort aims to provide a standardized road map and logical decision flow to better facilitate fixative technologies into D and D activities required to achieve the desired facility end state. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Demand response event simulator and risk-aware bidding tool for industrial customers

Incentive Based Demand Response (IBDR) program participation delivers financial benefits to the consumers and resiliency benefits to the electricity grid. Effectively participating in these programs as an industrial consumer requires bidding strategies that balance financial risk with operational constraints. Existing bidding tools tend not to fully incorporate stochastic IBDR event modeling, program specific baseline and payment/penalty calculations, or demand reduction process control schemes that account for the cascading impacts of shutdown in complex facilities. Here, this work presents an IBDR event simulator and risk-aware bidding framework tool integrating three key components: a flexible, parameterized demand response event generator that rigorously accounts for program structures and stochasticity, a demand response operational simulation model that generates explicit control strategies for load reduction, and a Monte Carlo simulator to evaluate financial risk for varied capacity bids. A case study at a wastewater treatment plant participating in PG&E's Capacity Bidding Program demonstrates the framework's utility. In the peak capacity price month of August, optimal bidding by the wastewater treatment plant nets a mean IBDR benefit of $101,000 (67% of the August electricity bill) with 0.4% probability of a financial loss. This framework enables industrial operators to make informed bidding decisions, negotiate better program terms with demand response load aggregators, and analyze energy flexibility investments at their facilities. Ultimately, this work reduces participation barriers in IBDR programs and supports the broader goal of enhancing grid reliability and renewable energy integration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Optimal Transfer Operators in Algebraic Two-Level Methods for Nonsymmetric and Indefinite Problems

Consider an algebraic two-level method applied to the 𝑛-dimensional linear system 𝐴⁢𝒙 = 𝒃 using fine-space preconditioner (i.e., “relaxation” or “smoother”) 𝑀, with 𝑀 ≈ 𝐴, restriction and interpolation 𝑅 and 𝑃, and algebraic coarse-space operator 𝐴 𝑐 : = 𝑅 ∗ ⁢𝐴⁢𝑃. Then, what are the best possible transfer operators 𝑅 and 𝑃 of a given dimension 𝑛 𝑐 < 𝑛? Brannick et al. [12] showed that when 𝐴 and 𝑀 are Hermitian positive definite (HPD), the optimal interpolation is such that its range contains the 𝑛 𝑐 smallest generalized eigenvectors of the matrix pencil (𝐴, 𝑀). Recently, in Ali et al. [5] we generalized this framework to the non-HPD setting, by considering both right (interpolation) and left (restriction) generalized eigenvectors of (𝐴, 𝑀) and defining corresponding nonsymmetric transfer operators {𝑅#, 𝑃#}. Tight convergence bounds for {𝑅#, 𝑃#} are derived in spectral radius, as well as a proof of pseudo-optimality. Note, {𝑅#, 𝑃#} are typically complex valued, which is not practical for real-valued problems. Here, in this work, we build on [5], first characterizing all inner products in which the coarse-space correction defined by {𝑅#, 𝑃#} is orthogonal. We then develop tight two-level convergence bounds in these norms, and prove that the underlying transfer operators {𝑅#, 𝑃#} are genuinely optimal. As a special case, our theory both recovers and extends the HPD results from [12]. Finally, we show how to construct optimal, real-valued transfer operators in the case of that 𝐴 and 𝑀 are real valued, but are not HPD. Numerical examples arising from a discretized advection-reaction equation, wave-equation, and Stokes equations are used to verify and illustrate the theory.

97 MATHEMATICS AND COMPUTING↗

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗