Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “residual 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

A conceptual basis for the design of damage-tolerant structural systems

Researchers define damage-tolerant structural systems as those systems which not only have adeqate intact strength to withstand initial failure but also adequate residual strength to minimize the possibility of, and hence the consequences of, further failure. The incorporation of damage tolerance cannot be done in total isolation of the function being required of the system and the costs associated with obtaining improved damage tolerance. The approach, therefore, is to formulate multiple-objective, multi-level decision support problems (DSP), the solutions of which represent a compromise between higher costs and higher damage tolerance. Mulitple-objective decision support problems are easily solved in the linear domain. These formulations, however, include both linear and nonlinear constraints and goals, which in the past, have not been considered due to the resulting complexity. Here, researchers: (1) present a complete discussion and description of decision support problems; (2) identify what further research needs to be done in order to obtain information that is required but not known for solving problems using these models; and (3) identify what needs to be done to implement this prototype method in practice.

Mistree, F.↗
NASA NTRS1990

Comparison of two on-orbit attitude sensor alignment methods

Compared here are two methods of on-orbit alignment of vector attitude sensors. The first method uses the angular difference between simultaneous measurements from two or more sensors. These angles are compared to the angular differences between the respective reference positions of the sensed objects. The alignments of the sensors are adjusted to minimize the difference between the two sets of angles. In the second method, the sensor alignment is part of a state vector that includes the attitude. The alignments are adjusted along with the attitude to minimize all observation residuals. It is shown that the latter method can result in much less alignment uncertainty when gyroscopes are used for attitude propagation during the alignment estimation. The additional information for this increased accuracy comes from knowledge of relative attitude obtained from the spacecraft gyroscopes. The theoretical calculations of this difference in accuracy are presented. Also presented are numerical estimates of the alignment uncertainties of the fixed-head star trackers on the Extreme Ultraviolet Explorer spacecraft using both methods.

Krack, Kenneth↗
NASA NTRS2001

On-Orbit Performance of Autonomous Star Trackers

This paper presents the results of a performance study of the autonomous star trackers (ASTs) on the IMAGE and the EO-1 spacecraft. IMAGE is a spinning spacecraft without gyros or redundant precision attitude sensors, so the statistical properties of the AST are estimated simply by comparing the output observed quaternions with a rigid rotator model with constant angular momentum. The initial conditions are determined by a least-squares fit to minimize the AST residuals. An additional fit is used to remove the remaining systematic error and to obtain the inherent sensor noise. Gyro rate data are available for the EO-1 mission, so the AST noise statistics are obtained from the residuals after solving for an epoch attitude and gyro bias also using a least-squares method.

Airapetian, V.↗
NASA NTRS2012

Spectral Line Parameters Including Temperature Dependences of Self- and Air-Broadening in the 2 (left arrow) 0 Band of CO at 2.3 micrometers

Temperature dependences of pressure-broadened half-width and pressure-induced shift coefficients along with accurate positions and intensities have been determined for transitions in the 2<--0 band of C-12 O-16 from analyzing high-resolution and high signal-to-noise spectra recorded with two different Fourier transform spectrometers. A total of 28 spectra, 16 self-broadened and 12 air-broadened, recorded using high- purity (greater than or equal to 99.5% C-12-enriched) CO samples and CO diluted with dry air(research grade) at different temperatures and pressures, were analyzed simultaneously to maximize the accuracy of the retrieved parameters. The sample temperatures ranged from 150 to 298K and the total pressures varied between 5 and 700 Torr. A multispectrum nonlinear least squares spectrum fitting technique was used to adjust the rovibrational constants (G, B, D, etc.) and intensity parameters (including Herman-Wallis coefficients), rather than determining individual line positions and intensities. Self-and air-broadened Lorentz half-width coefficients, their temperature dependence exponents, self- and air-pressure-induced shift coefficients, their temperature dependences, self- and air-line mixing coefficients, their temperature dependences and speed dependence have been retrieved from the analysis. Speed-dependent line shapes with line mixing employing off-diagonal relaxation matrix element formalism were needed to minimize the fit residuals. This study presents a precise and complete set of spectral line parameters that consistently reproduce the spectrum of carbon monoxide over terrestrial atmospheric conditions.

Devi, V. Malathy↗
DOE OSTI2024

Spot pattern welding scanning strategy for sensor embedding and residual stress reduction in laser-foil-printing additive manufacturing

Here, this paper aims to present spot pattern welding (SPW) as a scanning strategy for laser-foil-printing (LFP) additive manufacturing (AM) in place of the previously used continuous pattern welding (CPW) (line-raster scanning). The SPW strategy involves generating a sequence of overlapping spot welds on the metal foil, allowing the laser to form dense and uniform weld beads. This in turn reduces thermal gradients, promotes material consolidation and helps mitigate process-related risks such as thermal cracking, porosity, keyholing and Marangoni effects. 304L stainless steel (SS) feedstock is used to fabricate test specimens using the LFP system. Imaging techniques are used to examine the melt pool dimensions and layer bonding. In addition, the parts are evaluated for residual stresses, mechanical strength and grain size. Compared to CPW, SPW provides a more reliable heating/cooling relationship that is less dependent on part geometry. The overlapping spot welds distribute heat more evenly, minimizing the risk of elevated temperatures during the AM process. In addition, the resulting dense and uniform weld beads contribute to lower residual stresses in the printed part. To the best of the authors’ knowledge, this is the first study to thoroughly investigate SPW as a scanning strategy using the LFP process. In general, SPW presents a promising strategy for securing embedded sensors into LFP parts while minimizing residual stresses.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Neutron Characterization for Additive Manufacturing

Oak Ridge National Laboratory (ORNL) is leveraging decades of experience in neutron characterization of advanced materials together with resources such as the Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR) shown in Fig. 1 to solve challenging problems in additive manufacturing (AM). Additive manufacturing, or three-dimensional (3-D) printing, is a rapidly maturing technology wherein components are built by selectively adding feedstock material at locations specified by a computer model. The majority of these technologies use thermally driven phase change mechanisms to convert the feedstock into functioning material. As the molten material cools and solidifies, the component is subjected to significant thermal gradients, generating significant internal stresses throughout the part (Fig. 2). As layers are added, inherent residual stresses cause warping and distortions that lead to geometrical differences between the final part and the original computer generated design. This effect also limits geometries that can be fabricated using AM, such as thin-walled, high-aspect- ratio, and overhanging structures. Distortion may be minimized by intelligent toolpath planning or strategic placement of support structures, but these approaches are not well understood and often "Edisonian" in nature. Residual stresses can also impact component performance during operation. For example, in a thermally cycled environment such as a high-pressure turbine engine, residual stresses can cause components to distort unpredictably. Different thermal treatments on as-fabricated AM components have been used to minimize residual stress, but components still retain a nonhomogeneous stress state and/or demonstrate a relaxation-derived geometric distortion. Industry, federal laboratory, and university collaboration is needed to address these challenges and enable the U.S. to compete in the global market. Work is currently being conducted on AM technologies at the ORNL Manufacturing Demonstration Facility (MDF) sponsored by the DOE's Advanced Manufacturing Office. The MDF is focusing on R&D of both metal and polymer AM pertaining to in-situ process monitoring and closed-loop controls; implementation of advanced materials in AM technologies; and demonstration, characterization, and optimization of next-generation technologies. ORNL is working directly with industry partners to leverage world-leading facilities in fields such as high performance computing, advanced materials characterization, and neutron sciences to solve fundamental challenges in advanced manufacturing. Specifically, MDF is leveraging two of the world's most advanced neutron facilities, the HFIR and SNS, to characterize additive manufactured components.

Watkins, Thomas↗
DOE OSTI2025

Ruthenium/Carbon Nanocomposites for Efficient Hydrogen Electrocatalysis: Impacts of Halide Residues

Ruthenium has emerged as a promising substitute for platinum toward the hydrogen evolution/oxidation reaction (HER/HOR). Herein, ruthenium/carbon composites are prepared by magnetic induction heating (300 A, 10 s) of RuCl3, RuBr3 or RuI3 loaded on hollow N-doped carbon cages (HNC). The HNC-RuCl3-300A sample consists of Ru nanoparticles (dia. 1.96 nm) and abundant Cl residues. HNC-RuBr3-300A possesses a larger nanoparticle size (≈19.36 nm) and lower content of Br residues. HNC-RuI3-300A contains only bulk-like Ru agglomerates with a minimal amount of I residues, due to reduced Ru-halide bonding interactions. Among these, HNC-RuCl3-300A exhibits the best HER activity in alkaline media, with a low overpotential of only -26 mV to reach 10 mA cm-2, even outperforming Pt/C, and can be used as the cathode catalyst for anion exchange membrane water electrolyzer (along with commercial RuO2 as the anode catalyst), producing 0.5 A cm- 2 at 1.88 V for up to 100 h, a performance markedly better than that with Pt/C. HNC-RuCl3-300A also exhibits the best HOR activity, with a half-wave potential (+18 mV) even lower than that of Pt/C (+35 mV). These activities are ascribed to the combined contributions of small Ru nanoparticles and Ru-to-halide charge transfer that weaken H adsorption.

Yu, Bingzhe↗
NASA NTRS2000

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗
DOE OSTI2021

Residual stresses and microstructure within Allvac 718Plus laser powder bed fusion bars

The residual stresses within Allvac 718Plus bars built using LBPF with different laser powers, speeds and table displacements were measured using diffraction and mechanical methods and modelled with the baseplate attached and then removed. The residual stress profiles within the bars from the top surface down through the bulk were all quite similar, becoming less tensile due to a strain hardening mechanism. Table displacement has the greatest impact on residual stress, decreasing with increasing displacement/powder layer thickness. There was good agreement amongst the modeling and measurements. The microstructures were examined and varied slightly with energy density with higher densities having larger grains and enhanced post solidification diffusion. We found energy density had minimal impact on the residuals stresses within the parameters to produce dense material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Laboratory Evaluation to Increase Effectiveness of Field-Scale Soil Flushing in the Hanford 100 Areas

This laboratory study was initiated to develop an improved technical understanding of the key controlling geochemical and physical processes of mobilizing residual Cr(VI) in Hanford sediments to provide increased effectiveness in soil flushing activities planned for Hanford 100 Areas. This was accomplished by (1) quantifying Cr mass, release rate, and surface phase changes as Cr is leached from different Hanford sediments; (2) maximizing Cr leaching during soil flushing by evaluating different leach solutions and leach strategies; and (3) minimizing leaching of remaining residual Cr after soil flushing ends. To quantify geochemical controls on Cr leaching from sediments, sequential liquid extractions were used to identify aqueous, adsorbed, and solid Cr surface phases (precipitates or Cr incorporated into other phases) before and after water-saturated and unsaturated leach experiments with pH 8 artificial groundwater and with other amendment solutions. Changes in the Cr release rate from sediments were also correlated to changes in Cr surface phases. During unsaturated infiltration, water advection occurs primarily through larger pores, with less movement through smaller grain layers leading to a slower release of Cr trapped in smaller pores next to smaller grains. To quantify these physical controls on Cr leaching from sediments, 10 ft high 1-D infiltration experiments were conducted at differing leach solution application rates, and Cr leaching during application and subsequent residual water flow was quantified. Qualitative 2-D experiments were also conducted to evaluate surfactant addition to increase flow in low-permeability zones. Cr-contaminated Hanford sediments from the 100-H, 100-K, 100-D, and 200 East Areas exhibited a wide variety of leach behavior from rapid release of a small amount of Cr to fast (i.e., minutes) and slow (i.e., hours or more) release of high Cr. Sediments with low labile Cr (< 0.2 µg/g) exhibited fast Cr release but less overall release, and most (>95%) leached from the sediments within a few pore volumes in water-saturated columns. With nearly no sorption of chromate, aqueous and sorbed Cr (as chromate) were released quickly from sediment within the first few pore volumes, resulting in high leach concentrations, fast initial Cr release rates, and changes in pre- and post-leach extractions. In contrast, sediments with high labile Cr (0.2 to 20 µg/g) exhibited a combination of fast and slow Cr release due to multiple Cr surface phases contributing to Cr release from the sediment at different rates with greater overall release. This resulted in elevated Cr effluent concentrations even after dozens of pore volumes were leached. In addition, the Cr release rate decreased with increasing time during leaching. A decrease in labile Cr (aqueous Cr, adsorbed Cr, pH 5 acetate dissolved precipitates) in post-leach extractions compared to pre-leach extractions showed initial aqueous and adsorbed Cr release, followed by dissolution of high-solubility precipitates such as CaCrO 4 (if present), then slow dissolution of calcite that may contain some chromate and possibly BaCrO 4 . The labile Cr in all sediments correlated well with the Cr release rate (r 2 = 0.81). Labile Cr also correlated well with the leached mass in water-saturated columns (r 2 = 0.98), and therefore can be used for prediction of potential removal with flushing applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2007

Asymmetric Bulkheads for Cylindrical Pressure Vessels

Asymmetric bulkheads are proposed for the ends of vertically oriented cylindrical pressure vessels. These bulkheads, which would feature both convex and concave contours, would offer advantages over purely convex, purely concave, and flat bulkheads (see figure). Intended originally to be applied to large tanks that hold propellant liquids for launching spacecraft, the asymmetric-bulkhead concept may also be attractive for terrestrial pressure vessels for which there are requirements to maximize volumetric and mass efficiencies. A description of the relative advantages and disadvantages of prior symmetric bulkhead configurations is prerequisite to understanding the advantages of the proposed asymmetric configuration: In order to obtain adequate strength, flat bulkheads must be made thicker, relative to concave and convex bulkheads; the difference in thickness is such that, other things being equal, pressure vessels with flat bulkheads must be made heavier than ones with concave or convex bulkheads. Convex bulkhead designs increase overall tank lengths, thereby necessitating additional supporting structure for keeping tanks vertical. Concave bulkhead configurations increase tank lengths and detract from volumetric efficiency, even though they do not necessitate additional supporting structure. The shape of a bulkhead affects the proportion of residual fluid in a tank that is, the portion of fluid that unavoidably remains in the tank during outflow and hence cannot be used. In this regard, a flat bulkhead is disadvantageous in two respects: (1) It lacks a single low point for optimum placement of an outlet and (2) a vortex that forms at the outlet during outflow prevents a relatively large amount of fluid from leaving the tank. A concave bulkhead also lacks a single low point for optimum placement of an outlet. Like purely concave and purely convex bulkhead configurations, the proposed asymmetric bulkhead configurations would be more mass-efficient than is the flat bulkhead configuration. In comparison with both purely convex and purely concave configurations, the proposed asymmetric configurations would offer greater volumetric efficiency. Relative to a purely convex bulkhead configuration, the corresponding asymmetric configuration would result in a shorter tank, thus demanding less supporting structure. An asymmetric configuration provides a low point for optimum location of a drain, and the convex shape at the drain location minimizes the amount of residual fluid.

Ford, Donald B.↗
DOE OSTI2024

Biocompatible organosolv fractionation via a novel alkaline lignin-first strategy towards lignocellulose valorization

Simultaneous valorization of both carbohydrate and lignin fractions in lignocellulose remains a great challenge. Herein, a novel lignin-first strategy using triethylene glycol (TEG) under alkaline conditions for effective biomass fractionation producing highly digestible carbohydrates and reactive lignin was developed. Delignification was over 80% and fermentable sugar yields were close to 90% after pretreatment at 90°C. The biocompatibility of TEG allowed direct enzymatic hydrolysis of the solid residue without washing, thus minimizing wastewater generation. Furthermore, the obtained lignin (TEGL) had an uncondensed structure with well-preserved β-O-4 linkage, leading to near-equal aromatic monomer yields compared to cellulolytic enzyme lignin after catalytic-free pyrolysis, demonstrating high valorization potential. Finally, the proposed TEG solvent system is promising for a green and sustainable biorefinery process to achieve the complete utilization of lignocellulose.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

An optical-input Maximum Likelihood Estimation feedback system demonstrated on tokamak horizontal equilibrium control

A readily parallelized Maximum Likelihood Estimation (MLE) algorithm with linear computational complexity is demonstrated in real time using only measurements from an extreme ultraviolet (EUV) diagnostic to control the horizontal position of a tokamak plasma. A set of trial emissivity profiles are parameterized by the control quantity of interest (R m ), and the MLE is identified from the profile which minimizes the signal reconstruction residual. The algorithm depends on an empirically determined likelihood function with exponential form. EUV emission (λ ≈ 15eV-1keV) is captured in a poloidal plane by four 16-channel AXUV diodes mounted at different poloidal angles with radial and angular resolution sufficient to discern plasma equilibrium evolution in HBT-EP. Calculations of the plasma major radius by the system are consistent within diagnostic uncertainty for the majority of the discharge with those of: a weighted average of vertical soft X-ray or EUV chords, magnetic sensors, and an equilibrium reconstruction. The feedback system corrects for a horizontal displacement of the major radius equal to 20% of the plasma minor radius by adjusting the vertical field produced from 40 in-vessel control coils in real time. The MLE calculation is performed on a GPU in a 15 μs cycle, with similar performance in this application to a simple weighted average of vertical chords. Finally, results demonstrate horizontal position control using magnetic actuators and an optical observer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spatial and temporal variability of tritium in precipitation within South Africa and it's bearing on hydrological studies

Tritium, the radioactive isotope of hydrogen, has been used to understand groundwater recharge processes for decades. The current variation of tritium in the atmosphere is largely attributed to stratospheric production and fall out rates as well as global circulation phenomena controlling the hydrological cycle. Global controls on the variability in atmospheric tritium activity are poorly suited to explain local variation and tritium activities in precipitation are often assumed to be uniform over both local and regional catchments and watersheds. This assumption can result in both over and under estimation of modern recharge within an aquifer when using tritium as the recharge proxy. In order to minimize the inherent prediction residuals associated with tritium based recharge investigations, the variability of tritium in precipitation was modelled from 127 spatial precipitation samples taken over a two year period, combined with a 76 precipitation sample group-set taken over a one year period in a single location. Precipitation events were traced backward in time, from the point of collection, using HYSPLIT modelling to ascertain the origins of moisture content as well as the altitudes of moisture origin reached along the particle track. Tritium activities, collected over a one year period in Paarl, range from 0.45 to 4.16 TU and have a mean of 1.59 TU. Spatial storm events in the Western Cape in 2017 and 2018 had a range from 0 to 2.2 and 0.37 to 3.27 TU, respectively, with mean activities of 1.18 (n = 34) and 1.25 TU (n = 32). Both storm events had similar tritium variability (σ = 0.5 n = 35 and 0.48 n = 32). Regional precipitation events had the largest range of tritium activities (0.55–12.2 TU). Although not all tritium activities can be explained by interrogating the water mass origin, this study suggests that approximately 90% of events can be completely or partially attributed to the origin of the water mass. The variability of tritium, both spatially and temporally, was higher than expected, confirming that when uniform tritium inputs are used, the groundwater system would provide inaccurate modern recharge estimates. Higher spatial resolution of tritium variation in precipitation for a particular region will improve our ability to relate tritium activities in groundwater to local precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗
DOE OSTI2023

Vadose Zone Soil Flushing for Chromium Remediation: A Laboratory Investigation to Support Field‐scale Application

Cr(VI) flushing from the vadose zone to the groundwater (with subsequent Cr(VI) removal in groundwater by pump-and-treat system) is a promising remedial technique that has recently been used at field scale. This laboratory study was conducted to provide the technical basis to design a field soil flushing strategy. The objectives were to (1) quantify the relationship between sediment Cr(VI) and Cr(III) mass and release rates and subsequent Cr(VI) leaching; (2) investigate different methodologies to maximize Cr(VI) leaching, and (3) investigate methods to minimize leaching of remaining residual Cr. Characterization of Cr-contaminated sediments (Hanford Site, WA) exhibited Cr(VI) showed that leach rates that were correlated to different Cr surface phases. Sediments with low leachable Cr(VI) (<2 μg/g) leached Cr rapidly, so slow infiltration of water in a single pulse was sufficient to leach most Cr. In contrast, sediments with high Cr (2 to 200 μg/g) released some Cr(VI) quickly but 10 to 50% Cr(VI) slowly (tens to hundreds of hours). Efficient unsaturated leaching of these sediments required a different infiltration strategy that includes: multiple slow leach pulses with time between flushing cycles; the use of a surfactant to increase Cr leaching from low-permeability zones, and the use of a reductant (Na-dithionite or Ca-polysulfide) in the final leach water was highly effective at decreasing residual Cr leaching. This study clearly demonstrated that the methodology of basing laboratory Cr flushing on parameters such as Cr release mass and rates could be used to improve the efficiency of soil flushing at field scale.

Accelerated Cr Flushing↗
DOE OSTI2021

Optimization-based algorithms for nonlinear mechanics and frictional contact

An optimization-based strategy for solving nonlinear mechanics problems is proposed. In contrast to typical nonlinear equation solver algorithms that aim to find zeros in the residual force function, we minimize an energy (or energy-like) function to encourage solutions which are locally stable equilibria. These smooth and potentially non-convex objective functions are minimized using a preconditioned conjugate-gradient trust-region algorithm. Contact is formulated as an inequality constrained minimization problem, and is solved with an augmented Lagrangian algorithm. Friction is included in the approach via a regularized quasi-potential energy, and other dissipative behavior is included through the use of variational constitutive updates. Finally, to accelerate convergence rates for the Lagrange multipliers, we propose a novel multiplier update algorithm utilizing the Fischer-Burmeister function, and demonstrate super-linear solver convergence for some applications.

42 ENGINEERING↗
NASA NTRS1976

Advanced helium regulator for a fluorine propellant system

The space storable propulsion module is an advanced high performance (375 seconds Isp minimum) planetary spacecraft propulsion system with a mission life of 5-10 years. The propellants used are liquid fluorine and amine fuel. This application requires high pressure regulator accuracy to optimize propellant depletion characteristics. An advanced regulator concept was prepared which is compatible with both fuel and oxidizer and which features design concepts such as redundant bellows, all-metallic/ceramic construction, friction-free guidance of moving parts and gas damping. Computer simulation of the propulsion module performance over two mission profiles indicated satisfactory minimization of those propellant residual requirements imposed by regulator performance variables.

Wichmann, H.↗