Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Alternating minimization”

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 865 records · Page 48

Thermophysical Property Measurement and Materials Research in the NASA/MSFC Electrostatic Levitator

Containerless processing is an important tool for thermophysical property measurements and materials research. The freedom from a crucible allows processing of liquid materials in a metastable undercooled state, as well as allowing processing of high temperature and highly reactive melts. Electrostatic levitation (ESL) is a containerless method which provides a number of unique advantages, including the decoupling of positioning force from sample heating, the ability to operate in ultra-high vacuum or at moderate gas pressure (approximately 3 atm), and the ability to process non-conducting materials. ESL also has the potential to reduce internal flow velocities below those possible with electromagnetic, acoustic, or aero-acoustic techniques. The ESL facility at NASA's Marshall Space Flight Center (MSFC) is in use for thermophysical property measurements and materials research by a number of different internal and external investigators. The methods for obtaining access to the facility, as well as a summary of current capabilities and some future directions will be discussed. In electrostatic levitation, the acceleration of gravity (or residual acceleration in reduced gravity) is opposed by the action of an applied electric field on a charged sample. This positioning method is applicable to any material which can be electrically charged, whether solid or liquid, conducting, or insulating. Because the position of the sample is unstable, a 3-dimensional active control loop rapidly adjusts the applied field to maintain levitation and minimize motion of the sample. Heating lasers melt the sample, and may be adjusted to maintain specified thermal profiles. Microgravity allows electrostatic levitation to work even more effectively. With the need to cancel less than 1 milli-g, the applied field required is reduced from approximately 10 MV/m to approximately 10 kV/m. Alternatively, a microgravity ESL can position larger samples than is possible on the ground, or it can position samples which maintain their charge poorly. Microgravity also reduces the effects of buoyant convection and sedimentation. In the MSFC ESL, many different classes of experiments and measurements have been demonstrated in collaboration with many investigators. These capabilities include heat capacity measurement, phase diagram determination, and nucleation and growth experiments including nucleation rate, TTT curves, and phase selection. We have demonstrated the capability of triggering recalescence of an undercooled liquid sample without loss of levitation. Currently under development are techniques for measuring the surface tension, viscosity, and density of levitated samples, as well as other thermophysical properties. Expanded levitation capabilities including various gas environments for sample processing, sample rotation control and management of thermal gradients are under development for the MSFC ESL. Containerless processing allows research on materials which may not be processed any other way. Electrostatic levitation extends the range of containerless processing to include many materials and classes of materials which cannot be processed by conventional electromagnetic, acoustic, or aerodynamic levitation techniques.

Rogers, Jan↗

Flexible Cryogenic Temperature and Liquid-Level Probes

Lightweight, flexible probes have been developed for measuring temperatures at multiple locations in tanks that contain possibly pressurized cryogenic fluids. If the fluid in a given tank is subcritical (that is, if it consists of a liquid and its vapor), then in one of two modes of operation, the temperature measurements made by a probe of this type can be used to deduce the approximate level of the liquid. The temperature sensors are silicon diodes located at intervals along a probe. If the probe is to be used to measure a temperature gradient along a given axis in the tank, then the probe must be mounted along that axis. In the non-liquid-level-sensing temperature-measurement mode, a constant small electric current is applied to each diode and the voltage across the diode . a known function of the current and temperature . is measured as an indication of its temperature. For the purpose of this measurement, "small electric current" signifies a current that is not large enough to cause a significant increase in the measured temperature. More specifically, the probe design calls for a current of 10 A, which, in the cryogenic temperature range of interest, generates heat at a rate of only about 0.01 mW per diode. In the liquid-level-sensing mode, one applies a larger current (30 mA) to each diode so as to heat each diode appreciably (with a power of about 36 mW in the temperature range of interest). Because the liquid cools the diode faster than does the vapor, the temperature of the diode is less when diode is immersed in the liquid than when it is above the surface of the liquid. Thus, the temperature (voltage) reading from each diode can be used to determine whether the liquid level is above or below the diode, and one can deduce that the liquid level lies between two adjacent diodes, the lower one of which reads a significantly lower temperature. The aforementioned techniques for measuring temperature and deducing liquid level are not new. What is new here are the designs of the probes and of associated external electronic circuitry. In each probe, the diodes and the lead wires are embedded in a strong, lightweight, flexible polyimide strip. Each probe is constructed as an integral unit that includes a multipin input/output plug or socket for solderless connection of the lead wires to the external circuitry. The polyimide strip includes mounting tabs with holes that can accommodate rivets, screws, or other fasteners. Alternatively, a probe can be mounted by use of an epoxy. A probe can be manufactured to almost any length or width, and the diodes can be embedded at almost any desired locations along and across the polyimide strip. In designing a probe for a specific application, one seeks a compromise between (1) minimizing the number of diodes in order to minimize the complexity of input/output connections and external electronic circuitry while (2) using enough diodes to obtain the required precision. Optionally, to minimize spurious heating of the cryogenic fluid, the external circuitry can be designed to apply power to the probe only during brief measurement intervals. Assuming that the external circuitry is maintained at a steady temperature, a power-on interval of only a few seconds is sufficient to obtain accurate data on temperatures and/or the height of the liquid/vapor interface.

Haberbusch, Mark↗

Cleanser, Detergent, Personal Care Product Pretreatment Evaluation

The purpose of the Cleanser, Detergent, Personal Care Product, and Pretreatment Evaluation & Selection task is to identify the optimal combination of personal hygiene products, crew activities, and pretreatment strategies to provide the crew with sustainable life support practices and a comfortable habitat. Minimal energy, mass, and crew time inputs are desired to recycle wastewater during long duration missions. This document will provide a brief background on the work this past year supporting the ELS Distillation Comparison Test, issues regarding use of the hygiene products originally chosen for the test, methods and results used to select alternative products, and lessons learned from testing.

Adam, Niklas↗

Cleanser, Detergent, Personal Care Product, and Pretreatment Evaluation

The purpose of the Cleanser, Detergent, Personal Care Product, and Pretreatment Evaluation & Selection task is to identify the optimal combination of personal hygiene products, crew activities, and pretreatment strategies to provide the crew with sustainable life support practices and a comfortable habitat. Minimal energy, mass, and crew time inputs are desired to recycle wastewater during long duration missions. This document will provide a brief background on the work this past year supporting the ELS Distillation Comparison Test, issues regarding use of the hygiene products originally chosen for the test, methods and results used to select alternative products, and lessons learned from testing.

Adam, Niklas↗

Flexible Cryogenic Temperature and Liquid-Level Probes

Lightweight, flexible probes have been developed for measuring temperatures at multiple locations in tanks that contain possibly pressurized cryogenic fluids. If the fluid in a given tank is subcritical (that is, if it consists of a liquid and its vapor), then in one of two modes of operation, the temperature measurements made by a probe of this type can be used to deduce the approximate level of the liquid. The temperature sensors are silicon diodes located at intervals along a probe. If the probe is to be used to measure a temperature gradient along a given axis in the tank, then the probe must be mounted along that axis. In the temperature-measurement mode, a constant small electric current is applied to each diode and the voltage across the diode a known function of the current and temperature is measured as an indication of its temperature. For the purpose of this measurement, small electric current signifies a current that is not large enough to cause a significant increase in the measured temperature. More specifically, the probe design calls for a current of 10 A, which, in the cryogenic temperature range of interest, generates heat at a rate of only about 0.01 mW per diode. In the liquid-level-sensing mode, one applies a larger current (30 mA) to each diode so as to heat each diode appreciably (with a power of about 36 mW in the temperature range of interest). Because the liquid cools the diode faster than does the vapor, the temperature of the diode is less when the diode is immersed in the liquid than when it is above the surface of the liquid. Thus, the temperature (voltage) reading from each diode can be used to determine whether the liquid level is above or below the diode, and one can deduce that the liquid level lies between two adjacent diodes, the lower one of which reads a significantly lower temperature. The aforementioned techniques for measuring temperature and deducing liquid level are not new. What is new here are the designs of the probes and of associated external electronic circuitry. In each probe, the diodes and the lead wires are embedded in a strong, lightweight, flexible polyimide strip. Each probe is constructed as an integral unit that includes a multipin input/output plug or socket for solderless connection of the lead wires to the external circuitry. The polyimide strip includes mounting tabs with holes that can accommodate rivets, screws, or other fasteners. Alternatively, a probe can be mounted by use of an epoxy. A probe can be manufactured to almost any length or width, and the diodes can be embedded at almost any desired location along and across the polyimide strip. In designing a probe for a specific application, one seeks a compromise between (1) minimizing the number of diodes in order to minimize the complexity of input/output connections and external electronic circuitry while (2) using enough diodes to obtain the required precision. Optionally, to minimize spurious heating of the cryogenic fluid, the external circuitry can be designed to apply power to the probe only during brief measurement intervals. Assuming that the external circuitry is maintained at a steady temperature, a power-on interval of only a few seconds is sufficient to obtain accurate data on temperatures and/or the height of the liquid/vapor interface.

Haberbusch, Mark↗

Alternative Explanations for Extreme Supersolar Iron Abundances Inferred from the Energy Spectrum of Cygnus X-1

Here we study a 1-200 keV energy spectrum of the black hole binary Cygnus X-1 taken with NuSTAR and Suzaku. This is the first report of a NuSTAR observation of Cyg X-1 in the intermediate state, and the observation was taken during the part of the binary orbit where absorption due to the companion's stellar wind is minimal. The spectrum includes a multi-temperature thermal disk component, a cutoff power-law component, and relativistic and nonrelativistic reflection components. Our initial fits with publicly available constant density reflection models (relxill and reflionx) lead to extremely high iron abundances (greater than 9.96 and 10.6 minus 0.9 (sub) to plus 1.6 (sup) times solar, respectively). Although supersolar iron abundances have been reported previously for Cyg X-1, our measurements are much higher and such variability is almost certainly unphysical. Using a new version of reflionx that we modified to make the electron density a free parameter, we obtain better fits to the spectrum even with solar iron abundances. We report on how the higher density (n (sub e) equal to 3.98 minus 0.25 (sub) to plus 0.12 (sup) times 10 (sup 20) per cubic centimeter) impacts other parameters such as the inner radius and inclination of the disk.

Tomsick, John A.↗

Revisiting a minimally destructive analytic approach for determining electrochemical kinetic parameters: Measuring aluminum corrosion across a wide pH range based on the Butler-Volmer equation

Here, this study revisits the three-point sampling of the simplified Butler-Volmer equation to address the limitations of strong potentiodynamic polarization, which can introduce irreversible damage and uncertainty in corrosion analysis. The method extracts electrochemical kinetic parameters while minimizing polarization effects, evaluates noise sensitivity relative to overpotential, and accounts for errors from signal noise, OCP drift, ohmic resistance, and mass-transfer constraints. Verified against the Tafel extrapolation method for aluminum corrosion across a wide pH range, this low-polarization approach enables accurate evaluations with specific error estimates, offering a robust alternative to linear polarization resistance methods that assume constant Tafel slopes.

36 MATERIALS SCIENCE↗

Trace levels of PtRu on NiMo foam electrodes towards hydrogen evolution reaction

Electrochemical water splitting plays a critical role in high purity hydrogen production. To lower capital cost and energy consumption, efficient catalysts are required to boost the hydrogen evolution reaction (HER), especially in alkaline media. Although PtRu alloy materials are still considered the state-of-the-art catalyst for the HER under alkaline media, large-scale application is hindered by its scarcity and high cost. NiMo alloy catalysts can be potential Pt-free alternatives to drive the sluggish kinetics of the HER under alkaline media, but their HER activity is still not ideal due to the limitation of Volmer step. Herein, we investigated the impact of trace levels of PtRu catalysts on NiMo electrodes towards the HER by incorporating advantages of both PtRu and NiMo catalysts to minimize cost and maintain a promising HER activity. Trace levels of PtRu catalysts from 0.0173 to 0.2648 mg cm −2 were loaded on commercial NiMo foam electrodes by a facile electrodeposition technique utilizing ppm concentration levels of PtRu in electrolyte. The results show that 0.0173 mg cm −2 PtRu loading on NiMo electrode can significantly enhance HER activity compared to the pristine NiMo foam electrode (Tafel slopes drop from ∼170 to ∼60 mV dec −1 ) due to the remarkable decline of energy barriers towards water dissociation and hydrogen adsorption. Meanwhile, the rate determining step (RDS) switches from the Volmer step to the Heyrovsky step when the PtRu loading increases to 0.0675 mg cm −2 . Promising catalytic stability is observed at 100 mA cm −2 over the course of 50 h with a PtRu loading of 0.1198 mg cm −2 . In conclusion, this work demonstrates a potential strategy to decrease the cost and simultaneously maintain superior HER performance by integrating trace levels of PtRu catalyst with NiMo foam electrodes.

Electrochemical water splitting↗

Kinetic energy classification and smoothing for compact B-spline basis sets in quantum Monte Carlo

Quantum Monte Carlo calculations of defect properties of transition metal oxides have become feasible in recent years due to increases in computing power. As the system size has grown, availability of on-node memory has become a limiting factor. Saving memory while minimizing computational cost is now a priority. The main growth in memory demand stems from the B-spline representation of the single particle orbitals, especially for heavier elements such as transition metals where semi-core states are present. Despite the associated memory costs, splines are computationally efficient. In this paper, we explore alternatives to reduce the memory usage of splined orbitals without significantly affecting numerical fidelity or computational efficiency. We make use of the kinetic energy operator to both classify and smooth the occupied set of orbitals prior to splining. By using a partitioning scheme based on the per-orbital kinetic energy distributions, we show that memory savings of about 50% is possible for select transition metal oxide systems. Finally, for production supercells of practical interest, our scheme incurs a performance penalty of less than 5%.

97 MATHEMATICS AND COMPUTING↗

SU-F-T-347: An Absolute Dose-Volume Constraint Based Deterministic Optimization Framework for Multi-Co60 Source Focused Radiotherapy

Purpose: Treatment plan optimization in multi-Co60 source focused radiotherapy with multiple isocenters is challenging, because dose distribution is normalized to maximum dose during optimization and evaluation. The objective functions are traditionally defined based on relative dosimetric distribution. This study presents an alternative absolute dose-volume constraint (ADC) based deterministic optimization framework (ADC-DOF). Methods: The initial isocenters are placed on the eroded target surface. Collimator size is chosen based on the area of 2D contour on corresponding axial slice. The isocenter spacing is determined by adjacent collimator sizes. The weights are optimized by minimizing the deviation from ADCs using the steepest descent technique. An iterative procedure is developed to reduce the number of isocenters, where the isocenter with lowest weight is removed without affecting plan quality. The ADC-DOF is compared with the genetic algorithm (GA) using the same arbitrary shaped target (254cc), with a 15mm margin ring structure representing normal tissues. Results: For ADC-DOF, the ADCs imposed on target and ring are (D100>10Gy, D50,10, 0<12Gy, 15Gy and 20Gy) and (D40<10Gy). The resulting D100, 50, 10, 0 and D40 are (9.9Gy, 12.0Gy, 14.1Gy and 16.2Gy) and (10.2Gy). The objectives of GA are to maximize 50% isodose target coverage (TC) while minimize the dose delivered to the ring structure, which results in 97% TC and 47.2% average dose in ring structure. For ADC-DOF (GA) techniques, 20 out of 38 (10 out of 12) initial isocenters are used in the final plan, and the computation time is 8.7s (412.2s) on an i5 computer. Conclusion: We have developed a new optimization technique using ADC and deterministic optimization. Compared with GA, ADC-DOF uses more isocenters but is faster and more robust, and achieves a better conformity. For future work, we will focus on developing a more effective mechanism for initial isocenter determination.

60 APPLIED LIFE SCIENCES↗

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

Novel low profile antenna candidates for EHF portable terminals

This paper presents three low profile antenna candidates for EHF portable communication terminals. The first structure is a planar, multilayer microstrip antenna utilizing electromagnetic coupling to minimize the complexity and losses associated with the feed network. The second candidate is a medium gain (15 dB) radiating cavity antenna utilizing a thick metallic ground plane. This element is amenable to device integration because a heat sink is incorporated and an area is available for fastening modular RF components. Thirdly, as an alternative to microstrip antennas, dielectric resonator antennas (DRA's) are presented. A broadband (28 percent) DRA is discussed. Experimental measurements are presented for all three antenna candidates.

D. Roscoe↗

Development of COPVS for High pressure, In-Space, Cryogenic Fuel Storage

Polymeric composite overwrapped pressure vessels (COPVs) provide an attractive material system to support developing commercial launch business and alternate fuel ventures. However to be able to design with these materials, the mechanical behavior of the materials must be understood with regards to processing, performance, damage tolerance, and environment. For the storage of cryogenic propellants, it is important to evaluate the materials performance and impact damage resistance at cryogenic temperatures in order to minimize weight and to ensure safety and reliability. To evaluate the ultimate performance, various polymeric COPV's have been statically burst tested at cryogenic conditions before and after exposure to irradiation. Materials selected for these COPVs were based on the measured mechanical properties of candidate resin systems and fibers that were also tested at cryogenic conditions before and after exposure to irradiation. The correlation of COPV burst pressures with the constituent material properties has proven to be a valuable screening method for selection of suitable candidate materials with resistance to material degradation due to exposure to temperature and radiation.

DeLay, Tom↗

High Performance COPVs for In-Space Storage of High Pressure Cryogenic Fuels

Polymeric composite overwrapped pressure vessels (COPVs) provide an attractive material system to support developing commercial launch business and alternate fuel ventures. However to be able to design with these materials, the mechanical behavior of the materials must be understood with regards to processing, performance, damage tolerance, and environment. For the storage of cryogenic propellants, it is important to evaluate the materials performance and impact damage resistance at cryogenic temperatures in order to minimize weight and to ensure safety and reliability. As part of this study, material tests of candidate fiber and resin systems were used as the basis for the selection of the material combinations for evaluation in a COPV at cryogenic conditions. This comprehensive approach has also been expanded to address issues with impact damage tolerance and material degradation due to environmental factors. KEY WORDS: Cryogenic testing, evaluation and applications for pressure vessels, COPVs, tanks, or storage vessels.

Schneider, Judy↗

Path Forward for Groundwater Compliance Strategy at the Tuba City, Arizona, Disposal Site - 19581

The Tuba City, Arizona, Disposal Site is located in the Navajo Nation and near the Hopi Reservation. A uranium mill operated at the site from 1956 to 1966. During the operational period, mill tailings were slurried to an unlined impoundment. Excess slurry water was transferred to unlined evaporation ponds. Groundwater contamination resulting from the infiltration of fluids in the mill tailings impoundment and evaporation ponds was characterized in the 1990's. The site's Groundwater Compliance Action Plan was prepared in 1999, and a groundwater extraction and treatment system was constructed in 2002. Issues with aquifer productivity and treatment plant reliability led to the suspension of operations in 2014. The limitations of the remediation system prompted a realization that the site conceptual model should be reevaluated. An integrated study to better understand the geochemical and hydrologic factors affecting contaminant transport, to support the development of a new groundwater compliance strategy, was initiated in 2017 and is ongoing. In addition to the groundwater science components, the study included the evaluation of enhanced evaporation methods. The integrated evaluation will lead to refinement of the site groundwater flow model and development of a contaminant transport model and subsequent development of remedial alternatives. Remedial alternatives will be evaluated to determine a new groundwater compliance strategy. The objectives of geochemical testing were to confirm and quantify geochemical reactions controlling uranium fate and transport. These reactions are important because they control the distribution of uranium on the solid phase; potential for release from the solid phase; mobility of uranium in the water phase; and potential for uranium redistribution onto the currently uncontaminated solid phase. A key outcome is the change in understanding of the mass of uranium in the subsurface that is potentially soluble. These results contradict the previous basis for remediation: that uranium was almost exclusively present in soluble form, within the boundaries of the contaminant plume. Testing showed that uranium sorption downgradient from the plume was measurable but relatively minimal. The movement of other mill-related constituents (nitrate, sulfate, alkalinity) results in the dissolution of naturally occurring uranium and a higher than background concentration of uranium in groundwater. These constituents are not subject to sorption, and their downgradient concentrations are reduced only by diffusion and dispersion. The geochemical testing provides conceptual and quantitative information on potential plume migration. While the testing revealed the potential availability of uranium contaminant sources, the aquifer buffering capacity, pH, and presence of oxygen in groundwater may, individually or in combination, be controlling factors for plume movement. Evaluation of evaporative treatment alternatives was undertaken to determine one or more sustainable and efficient methods capable of treating 151 liters per minute (L/min) (40 gallons per minute), the annual average flow rate. In addition to establishing the baseline capacity for evaporation using the site's existing pond, prospective methods of enhanced evaporation included preheating water before inflow to the evaporation pond, vertical circulation in the pond, dye addition, and Wind-Aided Intensified Evaporation. Consideration of a 151 L/min evaporative treatment arose from previous analysis of extraction well production rates and contaminant mass removal, which indicated that about 80% contaminant removal, relative to full-scale operation of the site's groundwater treatment plant, could be achieved with a smaller capacity system. The lower-flow evaporative option also represented a potentially more sustainable treatment system in comparison to restarting or reconfiguring the groundwater treatment plant. Aquifer pumping tests were designed to develop or improve estimates of hydraulic properties such as transmissivity, horizontal and vertical hydraulic conductivity, storativity, and storage coefficient. Testing involved pumping from one well at a constant flow rate for up to 5 days while monitoring groundwater drawdown in the pumping well and five surrounding wells and then shutting down the pumping well to observe groundwater recovery. The aquifer pumping test results provide insight regarding the hydraulics of the saturated zone. The saturated zone is in a low-conductivity formation with preferential flow in the horizontal direction. The aquifer exhibits low vertical conductance. The use of AQTESOLV software indicates that the saturated zone is a layered system due to the cross-bedded sand-grain orientation within the Navajo Sandstone. These indicators of limited vertical flow provide a potential justification for the observed depth of the contaminant plume of 15-46 meters (50-150 feet) below ground surface). Given the relatively low vertical conductance, groundwater extraction should be focused on the upper portions of the saturated zone. Aquifer response testing provided an improved understanding of aquifer characteristics. However, the primary objective, to identify optimal groundwater extraction strategies, has not been fully addressed. A list of path forward activities for hydrologic analyses was developed, with an overarching objective of providing input parameters for groundwater flow and contaminant transport modeling. Defensible predictions of contaminant transport, through modeling, are critical to the development of a new site groundwater compliance strategy. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Efficient sparse state preparation via quantum walks

Continuous-time quantum walks (CTQWs) on dynamic graphs, referred to as dynamic CTQWs, are a recently introduced universal model of computation that offers a new paradigm in which to envision quantum algorithms. In this work, we develop an algorithm that converts single-edge and self-loop dynamic CTQWs to the gate model of computation. We use this mapping to introduce an efficient sparse quantum state preparation framework based on dynamic CTQWs. Our approach utilizes combinatorics techniques such as minimal hitting sets, minimum spanning trees, and shortest Hamiltonian paths to reduce the number of controlled gates required to prepare sparse states. We show that our framework encompasses the current state of the art ancilla-free sparse state preparation method by reformulating this method as a CTQW. This CTQW-based framework offers an alternative to the uniformly controlled rotation method used by Qiskit by requiring fewer CX gates when the target state has a polynomial number of non-zero amplitudes.

dynamic continuous time quantum walks↗

Modeling the formation of globular cluster systems in the Virgo cluster

The mass distribution and chemical composition of globular cluster (GC) systems preserve fossil record of the early stages of galaxy formation. The observed distribution of GC colors within massive early-type galaxies in the ACS Virgo Cluster Survey (ACSVCS) reveals a multi-modal shape, which likely corresponds to a multi-modal metallicity distribution. We present a simple model for the formation and disruption of GCs that aims to match the ACSVCS data. This model tests the hypothesis that GCs are formed during major mergers of gas-rich galaxies and inherit the metallicity of their hosts. To trace merger events, we use halo merger trees extracted from a large cosmological N-body simulation. We select 20 halos in the mass range of 2 × 10{sup 12} to 7 × 10{sup 13} M {sub ☉} and match them to 19 Virgo galaxies with K-band luminosity between 3 × 10{sup 10} and 3 × 10{sup 11} L {sub ☉}. To set the [Fe/H] abundances, we use an empirical galaxy mass-metallicity relation. We find that a minimal merger ratio of 1:3 best matches the observed cluster metallicity distribution. A characteristic bimodal shape appears because metal-rich GCs are produced by late mergers between massive halos, while metal-poor GCs are produced by collective merger activities of less massive hosts at early times. The model outcome is robust to alternative prescriptions for cluster formation rate throughout cosmic time, but a gradual evolution of the mass-metallicity relation with redshift appears to be necessary to match the observed cluster metallicities. We also affirm the age-metallicity relation, predicted by an earlier model, in which metal-rich clusters are systematically several billion younger than their metal-poor counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗