Search NASA⌕ Search

SEARCH · Search NASA

Results for “control volume methods”

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 631 records · Page 35

Transonic Dynamics Tunnel Force and Pressure Data Acquired on the HSR Rigid Semispan Model

This report describes the aerodynamic data acquired on the High Speed Research Rigid Semispan Model (HSR-RSM) during NASA Langley Transonic Dynamics Tunnel (TDT) Test 520 conducted from 18 March to 4 April, 1996. The purpose of this test was to assess the aerodynamic character of a rigid high speed civil transport wing. The wing was fitted with a single trailing edge control surface which was both steadily deflected and oscillated during the test to investigate the response of the aerodynamic data to steady and unsteady control motion. Angle-of-attack and control surface deflection polars at subsonic, transonic and low-supersonic Mach numbers were obtained in the tunnel?s heavy gas configuration. Unsteady pressure and steady loads data were acquired on the wing, while steady pressures were measured on the fuselage. These data were reduced using a variety of methods, programs and computer systems. The reduced data was ultimately compiled onto a CD-ROM volume which was distributed to HSR industry team members in July, 1996. This report documents the methods used to acquire and reduce the data, and provides an assessment of the quality, repeatability, and overall character of the aerodynamic data measured during this test.

Schuster, David M.↗

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗

Advanced Unstructured Grid Generation for Complex Aerodynamic Applications

A new approach for distribution of grid points on the surface and in the volume has been developed and implemented in the NASA unstructured grid generation code VGRID. In addition to the point and line sources of prior work, the new approach utilizes surface and volume sources for automatic curvature-based grid sizing and convenient point distribution in the volume. A new exponential growth function produces smoother and more efficient grids and provides superior control over distribution of grid points in the field. All types of sources support anisotropic grid stretching which not only improves the grid economy but also provides more accurate solutions for certain aerodynamic applications. The new approach does not require a three-dimensional background grid as in the previous methods. Instead, it makes use of an efficient bounding-box auxiliary medium for storing grid parameters defined by surface sources. The new approach is less memory-intensive and more efficient computationally. The grids generated with the new method either eliminate the need for adaptive grid refinement for certain class of problems or provide high quality initial grids that would enhance the performance of many adaptation methods.

Pirzadeh, Shahyar Z.↗

Advanced Unstructured Grid Generation for Complex Aerodynamic Applications

A new approach for distribution of grid points on the surface and in the volume has been developed. In addition to the point and line sources of prior work, the new approach utilizes surface and volume sources for automatic curvature-based grid sizing and convenient point distribution in the volume. A new exponential growth function produces smoother and more efficient grids and provides superior control over distribution of grid points in the field. All types of sources support anisotropic grid stretching which not only improves the grid economy but also provides more accurate solutions for certain aerodynamic applications. The new approach does not require a three-dimensional background grid as in the previous methods. Instead, it makes use of an efficient bounding-box auxiliary medium for storing grid parameters defined by surface sources. The new approach is less memory-intensive and more efficient computationally. The grids generated with the new method either eliminate the need for adaptive grid refinement for certain class of problems or provide high quality initial grids that would enhance the performance of many adaptation methods.

Pirzadeh, Shahyar↗

High surface area carbon and process for its production

Activated carbon materials and methods of producing and using activated carbon materials are provided. In particular, biomass-derived activated carbon materials and processes of producing the activated carbon materials with prespecified surface areas and pore size distributions are provided. Activated carbon materials with preselected high specific surface areas, porosities, sub-nm (<1 nm) pore volumes, and supra-nm (1-5 nm) pore volumes may be achieved by controlling the degree of carbon consumption and metallic potassium intercalation into the carbon lattice during the activation process.

36 MATERIALS SCIENCE↗

High-fidelity retrieval from instantaneous line-of-sight returns of nacelle-mounted lidar including supervised machine learning

Abstract. Wind turbine applications that leverage nacelle-mounted Doppler lidar are hampered by several sources of uncertainty in the lidar measurement, affecting both bias and random errors. Two problems encountered especially for nacelle-mounted lidar are solid interference due to intersection of the line of sight with solid objects behind, within, or in front of the measurement volume and spectral noise due primarily to limited photon capture. These two uncertainties, especially that due to solid interference, can be reduced with high-fidelity retrieval techniques (i.e., including both quality assurance/quality control and subsequent parameter estimation). Our work compares three such techniques, including conventional thresholding, advanced filtering, and a novel application of supervised machine learning with ensemble neural networks, based on their ability to reduce uncertainty introduced by the two observed nonideal spectral features while keeping data availability high. The approach leverages data from a field experiment involving a continuous-wave (CW) SpinnerLidar from the Technical University of Denmark (DTU) that provided scans of a wide range of flows both unwaked and waked by a field turbine. Independent measurements from an adjacent meteorological tower within the sampling volume permit experimental validation of the instantaneous velocity uncertainty remaining after retrieval that stems from solid interference and strong spectral noise, which is a validation that has not been performed previously. All three methods perform similarly for non-interfered returns, but the advanced filtering and machine learning techniques perform better when solid interference is present, which allows them to produce overall standard deviations of error between 0.2 and 0.3 m s−1, or a 1 %–22 % improvement versus the conventional thresholding technique, over the rotor height for the unwaked cases. Between the two improved techniques, the advanced filtering produces 3.5 % higher overall data availability, while the machine learning offers a faster runtime (i.e., ∼ 1 s to evaluate) that is therefore more commensurate with the requirements of real-time turbine control. The retrieval techniques are described in terms of application to CW lidar, though they are also relevant to pulsed lidar. Previous work by the authors (Brown and Herges, 2020) explored a novel attempt to quantify uncertainty in the output of a high-fidelity lidar retrieval technique using simulated lidar returns; this article provides true uncertainty quantification versus independent measurement and does so for three techniques rather than one.

47 OTHER INSTRUMENTATION↗

Characterization of a selective, zero power sensor for distributed sensing of hydrogen in energy applications

The use of hydrogen as a clean and renewable energy source is increasing rapidly for both vehicle and stationary applications. There are safety concerns for the locations in which hydrogen is made, used, and transported (i.e., pipelines and tanker trucks). Sensors are needed to comply with safety regulations and to enable a smooth and safe rollout of hydrogen as an alternative energy. However, hydrogen sensors do not yet exist that have the combined features of small size and low power for easy deployment coupled with high-volume manufacturability and low cost. This is necessary to accommodate the emerging fixed and mobile markets while retaining critical metrological metrics, including measurement range, detection limits, selectivity, fast response, stability, and long lifetime. An amperometric gas sensor for hydrogen (AGS) has been developed using an innovative manufacturing method. The sensor was designed using scalable fabrication strategies based on “Printed Electronics” (PE) methodology which are compatible with large-scale production. Prototype sensors were batch fabricated with multiple individual elements on a substrate compatible in size with a standard 8-inch wafer to enable high-volume, low-cost manufacturing, thereby leveraging PE and semiconductor fabrication infrastructure. This novel AGS was instrumented with control circuitry and evaluated for hydrogen detection in the range of 0–5000 ppmv H2 in air. Specific performance evaluations included assessment of the sensor measurement range, repeatability, and selectivity.

08 HYDROGEN↗

Predicting the 5-Year Risk of Biochemical Relapse After Postprostatectomy Radiation Therapy in ≥PT2, pN0 Patients With a Comprehensive Tumor Control Probability Model

Purpose: To fit the individual biochemical recurrence-free survival (bRFS) data from patients treated with postprostatectomy radiation therapy (RT) with a comprehensive tumor control probability (TCP) model. Methods and Materials: Considering pre-RT prostate-specific antigen (PSA) as a surrogate of the number of clonogens, bRFS may be expressed as a function of dose-per-fraction–dependent radiosensitivity (α{sub eff}), the number of clonogens for pre-RT PSA = 1 ng/mL (C), and the fraction of patients who relapse because of clonogens outside the treated volume (K), assumed to depend (linearly or exponentially) on pre-RT PSA and Gleason score (GS). Data from 894 node-negative, ≥pT2, pN0 hormone-naive patients treated with adjuvant (n=331) or salvage (n=563) intent were available: 5-year bRFS data were fitted grouping patients according to GS (<7:392, =7:383, >7:119). Results: The median follow-up time, pre-RT PSA, and dose were 72 months, 0.25 ng/mL, and 66.6 Gy (range 59.4-77.4 Gy), respectively. The best-fit values were 0.23 to 0.26 Gy{sup −1} and 10{sup 7} for α{sub eff} and C for the model considering a linear dependence between K and PSA. Calibration plots showed good agreement between expected and observed incidences (slope: 0.90-0.93) and moderately high discriminative power (area under the curve [AUC]: 0.68-0.69). Cross-validation showed satisfactory results (average AUCs in the training/validation groups: 0.66-0.70). The resulting dose-effect curves strongly depend on pre-RT PSA and GS. bRFS rapidly decreases with PSA: the maximum obtainable bRFS (defined as 95% of the maximum) declined by about 2.7% and 4.5% for each increment of 0.1 ng/mL for GS <7 and ≥7, respectively. Conclusions: Individual data were fitted by a TCP model, and the resulting best-fit parameters were radiobiologically consistent. The model suggests that relapses frequently result from clonogens outside the irradiated volume, supporting the choice of lymph-node irradiation, systemic therapy, or both for specific subgroups (GS <7: PSA >0.8-1.0 ng/mL; GS ≥7: PSA >0.3 ng/mL). Early RT should be preferred over delayed RT; the detrimental effect of PSA increase can never be fully compensated by increasing the dose, especially for patients with GS ≥7.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗

Diminished forearm vasomotor response to central hypervolemic loading in aerobically fit individuals

The aim of this study was to test the hypothesis that cardiopulmonary baroreflex control of forearm vascular resistance (FVR) during central hypervolemic loading was less sensitive in exercise trained high fit individuals (HF) compared to untrained average fit individuals (AF). Eight AF (age: 24 +/- 1 yr and weight: 78.9 +/- 1.7 kg) and eight HF (22 +/- 1 yr 79.5 +/- 2.4 kg) voluntarily participated in the investigation. Maximal aerobic power (determined on a treadmill), plasma volume and blood volume (Evans blue dilution method) were significantly greater in the HF than AF (60.8 +/- 0.7 vs. 41.2 +/- 1.9 ml.kg-1.min-1, 3.96 +/- 0.17 vs 3.36 +/- 0.08 1, and 6.33 +/- 0.23 vs 5.28 +/- 0.13 1). Baseline heart rate (HR), central venous pressure (CVP), mean arterial pressure (MAP, measured by an intraradial catheter or a Finapres finger cuff), forearm blood flow (FBF, plethysmography), and FVR, calculated from the ratio (MAP-CVP)/FBF, were not different between the HF and the AF. Lower body negative pressure (LBNP, -5, -10, -15, and -20 torr) and passive leg elevation (LE, 50 cm) combined with lower body positive pressure (LBPP, +5, +10, and +20 torr) were utilized to elicit central hypovolemia and hypervolemia, respectively. Range of CVP (from LBNP to LE+LBPP) was similar in the AF (from -3.9 to +1.9 mm Hg) and HF (from -4.0 to +2.2 mm Hg). However, FVR/CVP was significantly less in the HF (-1.8 +/- 0.1 unit.mm Hg-1) than AF (-34 +/- 0.1 unit.mm Hg-1). The FVR decrease in response to increase in CVP was significantly diminished in the HF (-1.46 +/- 0.45 unit.mm Hg-1) compared to the AF (-4.40 +/- 0.97 unit.mm Hg-1), and during LBNP induced unloading the FVR/CVP of the HF (-2.01 +/- 0.49 unit.mm Hg-1) was less (P < 0.08) than the AF (-3.28 +/- 0.69 unit.mm Hg-1). We concluded that the cardiopulmonary baroreceptor mediated FVR reflex response was significantly less sensitive to changes in CVP in individuals who practice exercise training.

NASA Discipline Number 93-10↗

Exposure to polychlorinated biphenyls and the thyroid gland – examining and discussing possible longitudinal health effects in humans

Background: Many previous studies have dealt with the effect of polychlorinated biphenyls (PCBs) on the thyroid gland, but their findings are inconsistent. One problem of these studies has been their use of cross-sectional designs. Objectives: The aim of the current study is to investigate longitudinal effects of PCBs on the thyroid gland, focusing on: morphological changes in thyroid tissue (i.e. thyroid volume), changes in thyroid hormones and in thyroid antibodies. Methods: A total of 122 individuals (M{sub age}=44.7) were examined over a period of four years (t{sup 1} until t{sup 4}). Medical history was collected via interviews, an ultrasound examination was performed and blood samples were taken to determine plasma PCB levels, thyroid stimulating hormone (TSH), free triiodthyronine (fT3), free thyroxine (fT4), thyroid peroxidase antibodies (TPOab), thyreoglobulin antibodies (TGab) and thyroid-stimulating hormone receptor antibodies (TSHRab). Rank correlation coefficients and mixed effect models were performed controlling for age and total lipids. Results: There were negative correlations between higher chlorinated biphenyls and fT3, cross-sectionally as well as longitudinally. We also found an interaction effect of higher-chlorinated PCBs over time for fT4 as well as TSHRab. In case of high exposure, a decrease in fT4 and an increase in TSHRab level were found over time. In regards to the other variables, our findings yielded no clear results in the examined time period. Conclusion: This is the first study to shows a PCB-related effect on fT3, fT4 and TSHRab over a four year period. The data also suggest that morphological and antibody findings remain inconsistent and do not allow for unambiguous interpretation. - Highlights: • This is a longitudinal study which includes data from four cross sections. • Higher-chlorinated biphenyls are negatively correlated with fT3. • There are interactions of time and higher-chlorinated biphenyls to TSHRab and fT4. • There is an interaction of time and dioxin-like polychlorinated biphenyls to fT4. • No consistent effects were found for thyroid volume, TSH, TPOab and TGab.

54 ENVIRONMENTAL SCIENCES↗

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗

NAECON 87; Proceedings of the IEEE National Aerospace and Electronics Conference, Dayton, OH, May 18-22, 1987. Volumes 1, 2, 3, & 4

The present conference discusses topics in VLSI components and their packaging, signal processing, uses of cartographic data, data transmission, advanced avionics architectures, fiber-optics, information control and display, image processing, airborne radar and fire control, navigation, air data, Kalman filtering, power generation and control, spacecraft power structures, aircraft flying qualities, flight management, fault-tolerant computer architectures, actuation technologies, self-repairing flight control system technology, multivariable control, stability and control methods, and AFTI/F-16 flight test reports. Also discussed are the ADA/JOVIAL language and its applications, software acquisition and testing, advanced software concepts, software management, computer graphics and visual systems softwear, ADA in embedded avionics, 16- and 32-bit architectures, voice interaction applications, human/machine systems analysis, human factors and AI, mental workloads and displays, pilot acceleration protection research, communications system technology, space communications, reliability and maintainability, managerial techniques, engineering management, EM compatibility and nuclear hardening, expert systems, AI language/knowledge representation, expert system implementation, machine vision/optical processing, and advanced AI concepts and architectures.

Avionics↗

In-situ Condition Monitoring of Components in Small Modular Reactors Using Process and Electrical Signature Analysis. Final report, volume 1. Development of experimental flow control loop, data analysis and plant monitoring

The research and development under this project was focused on the following three major objectives: Objective 1: Identification of critical in-vessel SMR components for remote monitoring and development of their low-order dynamic models, along with a simulation model of an integral pressurized water reactor (iPWR). Objective 2: Development of an experimental flow control loop with motor-driven valves and pumps, incorporating data acquisition and on-line monitoring interface. Objective 3: Development of stationary and transient signal processing methods for electrical signatures, machinery vibration, and for characterizing process variables for equipment monitoring. This objective includes the development of a data analysis toolbox. The following is a summary of the technical accomplishments under this project: - A detailed literature review of various SMR types and electrical signature analysis of motor-driven systems was completed. A bibliography of literature is provided at the end of this report. Assistance was provided by ORNL in identifying some key references. - A review of literature on pump-motor modeling and digital signal processing methods was performed. - An existing flow control loop was upgraded with new instrumentation, data acquisition hardware and software. The upgrading of the experimental loop included the installation of a new submersible pump driven by a three-phase induction motor. All the sensors were calibrated before full-scale experimental runs were performed. - MATLAB-Simulink model of a three-phase induction motor and pump system was completed. The model was used to simulate normal operation and fault conditions in the motor-pump system, and to identify changes in the electrical signatures. - A simulation model of an integral PWR (iPWR) was updated and the MATLAB-Simulink model was validated for known transients. The pump-motor model was interfaced with the iPWR model for testing the impact of primary flow perturbations (upsets) on plant parameters and the pump electrical signatures. Additionally, the reactor simulation is being used to generate normal operation data and data with instrumentation faults and process anomalies. A frequency controller was interfaced with the motor power supply in order to vary the electrical supply frequency. The experimental flow control loop was used to generate operational data under varying motor performance characteristics. Coolant leakage events were simulated by varying the bypass loop flow rate. The accuracy of motor power calculation was improved by incorporating the power factor, computed from motor current and voltage in each phase of the induction motor.- A variety of experimental runs were made for steady-state and transient pump operating conditions. Process, vibration, and electrical signatures were measured using a submersible pump with variable supply frequency. High correlation was seen between motor current and pump discharge pressure signal; similar high correlation was exhibited between pump motor power and flow rate. Wide-band analysis indicated high coherence (in the frequency domain) between motor current and vibration signals. - Wide-band operational data from a PWR were acquired from AMS Corporation and used to develop time-series models, and to estimate signal spectrum and sensor time constant. All the data were from different pressure transmitters in the system, including primary and secondary loops. These signals were pre-processed using the wavelet transform for filtering both low-frequency and high-frequency bands. This technique of signal pre-processing provides minimum distortion of the data, and results in a more optimal estimation of time constants of plant sensors using time-series modeling techniques.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Integrated Experimental and Modeling Studies of Mineral Carbonation as a Mechanism for Permanent Carbon Sequestration in Mafic/Ultramafic Rocks

A program of laboratory experiments, modeling and fieldwork was carried out at Yale University, University of Maryland, and University of Hawai‘i, under a DOE Award (DE-FE0004375) to study mineral carbonation as a practical method of geologic carbon sequestration. Mineral carbonation, also called carbon mineralization, is the conversion of (fluid) carbon dioxide into (solid) carbonate minerals in rocks, by way of naturally occurring chemical reactions. Mafic and ultramafic rocks, such as volcanic basalt, are natural candidates for carbonation, because the magnesium and iron silicate minerals in these rocks react with brines of dissolved carbon dioxide to form carbonate minerals. By trapping carbon dioxide (CO 2 ) underground as a constituent of solid rock, carbonation of natural basalt formations would be a secure method of sequestering CO 2 captured at power plants in efforts to mitigate climate change. Geochemical laboratory experiments at Yale, carried out in a batch reactor at 200°C and 150 bar (15 MPa), studied carbonation of the olivine mineral forsterite (Mg 2 SiO 4 ) reacting with CO 2 brines in the form of sodium bicarbonate (NaHCO 3 ) solutions. The main carbonation product in these reactions is the carbonate mineral magnesite (MgCO 3 ). A series of 32 runs varied the reaction time, the reactive surface area of olivine grains and powders, the concentration of the reacting fluid, and the starting ratio of fluid to olivine mass. These experiments were the first to study the rate of olivine carbonation under passive conditions approaching equilibrium. The results show that, in a simple batch reaction, olivine carbonation is fastest during the first 24 hours and then slows significantly and even reverses. A natural measure of the extent of carbonation is a quantity called the carbonation fraction, which compares the amount of carbon removed from solution, during a run, to the maximum amount that could have been removed if the olivine initially present had fully dissolved and the cations released had subsequently precipitated in carbonate minerals. The carbonation fractions observed in batch experiments with olivine grains and powders varied significantly, from less than 0.01 (1%) to more than 0.5 (50%). Over time, the carbonation fractions reached an upper limit after about 24 to 72 hours of reaction, then stayed constant or decreased. The peak Final Scientific/Technical Report DE-FE0004275 | Mineral Carbonation | 4 coincided with the appearance of secondary magnesium-bearing silicate minerals, whose formation competes for magnesium ions in solution and can even promote conditions that dissolve magnesite. The highest carbonation fractions resulted from experiments with low ratios of concentrated solution to olivine, during which amorphous silica spheres or meshes formed, instead of secondary silicate minerals. The highest carbonation fractions appear to result from competing effects. Precipitation of silica layers on olivine reduces the reactive surface area and, thus, the rate of olivine dissolution (which ultimately limits the carbonation rate), but these same silica layers can also inhibit the formation of secondary silicate minerals that consume magnesite formed in earlier stages of carbonation. Simulation of these experiments with simple geochemical models using the software program EQ3/6 reproduces the general trends observed—especially the results for the carbonation fraction in short-run experiments. Although further experimentation and better models are needed, this study nevertheless provides a framework for understanding the optimal conditions for sequestering carbon dioxide by reacting CO 2 -bearing fluids with rocks containing olivine minerals. A series of experiments at the Rock Physics Laboratory at the University of Maryland studied the carbonation process during deformation of thermally cracked olivine-rich rock samples (dunite) saturated with CO 2 brines of varying compositions. A goal of these geomechanical experiments was to see if flow and deformation processes, which accompany natural carbonation reactions in underground settings, work to enhance or inhibit the reactions. The experiments involved hydrostatic compaction, followed by deformation at a constant rate of strain. Sample permeability was monitored during the reactions. Comparison of the samples’ volume changes to their axial strains (shortening) during deformation indicates that samples reacted with CO 2 -saturated brines accommodate more axial compaction, before the onset of dilation (a swelling that precedes rock failure), than samples reacted with distilled water. Analyses of the reacted samples with scanning electron microscope (SEM) images indicate, first, that dissolution of olivine occurring in the initial stages of carbonation can provide pathways to fluid flow that sustain the reaction, and, second, that carbonate minerals precipitated along existing fractures in the rocks may serve as asperities, or roughness on a crack’s surface that restricts its closure. Final Scientific/Technical Report DE-FE0004275 | Mineral Carbonation | 5 In a related study undertaken by one of the principal investigators as a spin-off of the main project, a simple model of (magnesite) crystal growth in the pore space of basalts undergoing carbonation was developed. The model suggests that, under a carefully controlled program of CO 2 injection, carbonate mineral growth can harden the rock formation against earthquakes that might otherwise be induced by the injection of large fluid volumes (Yarushina and Bercovici, 2013). The overall conclusion of the research project is that mineral carbonation of underground mafic and ultramafic rock formations is a viable candidate for long-term sequestration of man-made carbon dioxide. No results obtained during the project indicate that the method is inherently intractable in its implementation; moreover, enormous volumes of basalt near Earth’s surface are candidate locations for large-scale injection programs. The geochemical experiments do indicate, however, that there will be significant engineering challenges in maintaining high rates of carbonation, by delaying the onset of chemical conditions that promote formation of secondary silicate minerals and, therefore, slow down, or even reverse, the carbonation process. It remains an open question as to whether carbonation processes can be sustained for many years in an engineered system operating on a large scale—a scale capable of accommodating millions of tons of CO 2 annually. The development of realistic theoretical models that can systematically describe the combined effects of reactive flow, precipitation and geomechanical deformation is a major barrier to further understanding of the practical viability of mineral carbonation as large-scale method of carbon sequestration.

58 GEOSCIENCES↗

Solvent- and catalyst-free mechanochemical synthesis of alkali metal monohydrides

Alkali metal monohydrides, AH (A = Li–Cs) have been synthesized in quantitative yields at room temperature by reactive milling of alkali metals in the presence of hydrogen gas at 200 bar or less. The mechanochemical approach reported here eliminates problems associated with the malleability of alkali metals — especially Li, Na, and K — and promotes effective solid–gas reactions, ensuring their completion. This is achieved by incorporating a certain volume fraction of the corresponding hydride powder as a process control agent, which allows continuous and efficient milling primarily by coating the surface of metal particles, effectively blocking cold welding. Formation of high-purity crystalline monohydrides has been confirmed by powder X-ray diffraction, solid-state NMR spectroscopy, and volumetric analyses of reactively desorbed H 2 from as-milled samples. The proposed synthesis method is scalable and particularly effective for extremely air-sensitive materials, such as alkali and alkaline earth metal hydrides. Furthermore, the technique may also be favorable for production in continuous reactors operating at room temperature, thereby reducing the total processing time, energy consumption and, hence, the cost of production of these hydrides or their derivatives and composites.

36 MATERIALS SCIENCE↗