Search NASA⌕ Search

SEARCH · Search NASA

Results for “Design Methods Modeling and Validation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

An optimization method for flexural bearing design for high-stroke high-frequency applications

A method for optimization of a high-stroke, high-frequency flexural bearing was completed. The main goal was to design a flexural bearing system as the energy restoration component of a one kW free piston linear engine alternator (LEA) operating at 90 Hz and maximum targeted displacement of 22 mm. Defined electricity output of 1 kW, with more than 30% efficiency, required engine operation under high speed and high stroke conditions. Under such challenging stroke and frequency requirements, a method of design was required to select between a wide range of parameters in flexural bearing design. A parametric CAD model was prepared and used in the optimization package of a finite element analysis (FEA) software. Four main parameters of flexural bearing geometry were identified and included the outer diameter, thickness, number of arms, and sweep angle of spiral cuts. Other design parameters including spiral gap width, shape factor and spiral start points angle were shown to be less important compared to these four main parameters in primary design and were examined separately. A neural network algorithm was used to investigate the interactions and effects of design parameters. The sensitivity of each parameter on output criteria such as maximum stress and natural frequency was examined and design charts developed. Methods for further modification of the spiral curve at end points were introduced to reduce the maximum stress level by up to 20%. Several steel alloys and a titanium-based alloy were identified as proper material candidates for high cycle applications and the final designed spring was manufactured and tested to validate the FEA results of axial stiffness and strain level throughout the flexure arms. The percent of the flexure’s moving mass was estimated to be 29% with FEA analysis and compared well with the experimental value of 31%. The stress/strain results of the FEA analysis matched with the experimental results – within ±5%.

42 ENGINEERING↗

Electromagnetic wave propagation based upon spectral-element methodology in dispersive and attenuating media

SUMMARY We build on mathematical equivalences between Maxwell’s wave equations for an electromagnetic medium and elastic seismic wave equations. This allows us to readily model Maxwell’s wave propagation in the spectral-element codes SPECFEM2D and SPECFEM3D, written for acoustic, viscoelastic and poroelastic seismic wave propagation, providing the ability to handle complex geometries, inherent to finite-element methods and retaining the strength of exponential convergence and accuracy due to the use of high-degree polynomials to interpolate field functions on the elements, characteristic to spectral-element methods (SEMs). Attenuation and dispersion processes related to the frequency dependence of dielectric permittivity and conductivity are also included using a Zener model, similar to shear attenuation in viscoelastic media or viscous diffusion in poroelastic media, and a Kelvin–Voigt model, respectively. Ability to account for anisotropic media is also discussed. Here, we limit ourselves to certain dielectric permittivity tensor geometries, in order to conserve a diagonal mass matrix after discretization of the system of equations. Doing so, simulation of Maxwell’s wave equations in the radar frequency range based on SEM can be solved using explicit time integration schemes well suited for parallel computation. We validate our formulation with analytical solutions. In 2-D, our implementation allows for the modelling of both a transverse magnetic (TM) mode, suitable for surface based reflection ground penetration radar type of applications, and a transverse electric (TE) mode more suitable for crosshole and vertical radar profiling setups. Two 2-D examples are designed to demonstrated the use of the TM and TE modes. A 3-D example is also presented, which allows for the full TEM solution, different antenna orientations, and out-of-plane variations in material properties.

58 GEOSCIENCES↗

Device-scale CFD study for mass transfer coefficient and effective mass transfer area in packed column

Carbon dioxide (CO 2 ) is a significant contributor of greenhouse gases emissions responsible for global warming. The fossil-fueled power plant is one of the major sources of CO 2 emission worldwide and needs to be mitigated. Post-combustion carbon capture by chemical absorption can be a promising technology to reduce CO 2 emission from power plants. The absorption is carried out via countercurrent flow in the packed columns. Although monoethanolamine (MEA) is a popular solvent, it has also limitations such as high volatility of solvent, corrosivity that leads to product degradation and damage to apparatuses, low CO 2 loading capacity, etc. Therefore, alternate solvents are being explored as prospective materials for carbon capture. A water-lean solvent, CO 2 binding organic liquids (CO2BOL), is currently being developed as a prospective material in carbon capture at the Pacific Northwest National Laboratory. Since the development of the material is at the primitive stage, it needs extensive investigation before industrial deployments. Hydrodynamics and mass transfer phenomena in the packed column using CO 2 BOL is unknown. Computational fluid dynamics (CFD) can be a useful tool to understand these phenomena in the packed column. Before employing the extensive CFD investigations, one needs to carefully validate CFD predictions for accurate flow predictions. This report demonstrates the validation of multiphase flow simulations against the experimental results for both structured and random packed columns. The reactive multiphase flow simulations using volume of fluid method were conducted at a wide range of liquid loads and contacts angles. The effective area is one of the critical factors dictating interphase mass transfer due to chemical reactions. Accordingly, CFD-predicted effective areas were compared with the experimental results of Tsai et al. [1] for Mellapak 250.Y structured packing. Both results matched well at different liquid loads in the presence of fast chemical kinetics using caustic solvents. Further, reactive multiphase flow simulations were also conducted in pall ring random packing to compute effective areas in the presence of fast chemical reactions. The CFD-predicted effective area compared well with correlation developed by Song et al. [2]. With successful validations, the current simulation framework can be further used for prediction of hydrodynamics and mass transfer phenomena using CO 2 BOL solvents in the packed column. Subsequently, it would be helpful in design optimization of a packed column. In addition to validation of the CFD model, this report also explains the mechanism for enhanced liquid side mass transfer for a structured packed column as compared to the wetted wall column.

20 FOSSIL-FUELED POWER PLANTS↗

Agglomerate Sizing in Aluminized Propellants Using Digital Inline Holography and Traditional Diagnostics

Aluminized ammonium perchlorate composite propellants can form large molten agglomerated particles that may result in poor combustion performance, slag accumulation, and increased two-phase flow losses. Quantifying agglomerate size distributions are needed to gain an understanding of agglomeration dynamics and ultimately design new propellants for improved performance. Due to complexities of the reacting multiphase environment, agglomerate size diagnostics are difficult and measurement accuracies are poorly understood. To address this, the current paper compares three agglomerate sizing techniques applied to two propellant formulations. Particle collection on a quench plate and backlit videography are two relatively common techniques, whereas digital inline holography is an emerging alternative for three-dimensional measurements. Atmospheric pressure combustion results show that all three techniques are able to capture the qualitative trends; however, significant differences exist in the quantitative size distributions and mean diameters. For digital inline holography, methods are proposed that combine temporally resolved high-speed recording with lower-speed but higher spatial resolution measurements to correct for size–velocity correlation biases while extending the measurable size dynamic range. Finally, the results from this work provide new guidance for improved agglomerate size measurements along with statistically resolved datasets for validation of agglomerate models.

42 ENGINEERING↗

Development of procedures for calculating stiffness and damping properties of elastomers in engineering applications. Part 1: Verification of basic methods

The primary aim of this investigation was verification of basic methods which are to be used in cataloging elastomer dynamic properties (stiffness and damping) in terms of viscoelastic model constants. These constants may then be used to predict dynamic properties for general elastomer shapes and operating conditions, thereby permitting optimum application of elastomers as energy absorption and/or energy storage devices in the control of vibrations in a broad variety of applications. The efforts reported involved: (1) literature search; (2) the design, fabrication and use of a test rig for obtaining elastomer dynamic test data over a wide range of frequencies, amplitudes, and preloads; and (3) the reduction of the test data, by means of a selected three-element elastomer model and specialized curve fitting techniques, to material properties. Material constants thus obtained have been used to calculate stiffness and damping for comparison with measured test data. These comparisons are excellent for a number of test conditions and only fair to poor for others. The results confirm the validity of the basic approach of the overall program and the mechanics of the cataloging procedure, and at the same time suggest areas in which refinements should be made.

Chiang, T.↗

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.

Gee, Ken↗

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗

IAEA Photonuclear Data Library 2019

Here, photo-induced reaction cross section data are of importance for a variety of current or emerging applications, such as radiation shielding design and radiation transport analyses, calculations of absorbed dose in the human body during radiotherapy, physics and technology of fission reactors (influence of photo-reactions on neutron balance) and fusion reactors (plasma diagnostics and shielding), activation analyses, safeguards and inspection technologies, nuclear waste transmutation, medical isotope production and astrophysical applications. To address these data needs the IAEA Photonuclear Data library was produced in 1999, containing evaluated photo-induced cross sections and neutron spectra for 164 nuclides which were deemed relevant for the applications. Since the release of the IAEA Photonuclear Data Library however, new experimental data as well as new methods to assess the reliability of experimental cross sections have become available. Theoretical models and input parameters used to evaluate photo-induced reactions have improved significantly over the years. In addition, new measurements of partial photoneutron cross sections using mono-energetic photon beams and advanced neutron detection systems have been performed allowing for the validation of the evaluations and assessments of the experimental data. Furthermore, technological advances have led to the construction of new and more powerful gamma-beam facilities, therefore new data needs are emerging. We report our coordinated efforts to address these data needs and present the results of the new up-to-date evaluations included in the new updated IAEA Photonuclear Data Library consisting of 219 nuclides. The new library includes 188 new evaluations produced by the CRP evaluators, and one evaluation taken from the JENDL/PD-2016 library, while 20 evaluations were retained from the previous 1999 IAEA Photonuclear Data Library. In most of the cases, the photon energy goes up to 200 MeV. A total of 55 nuclides are new in this library reflecting the progress in measurements but also the developing data needs. In this paper we discuss the new assessment method and make recommendations to the user community in cases where the experimental data are discrepant and the assessments disagree. In addition, in the absence of experimental data, we present model predictions for photo-induced reaction cross section on nuclides of potential interest to medical radioisotope production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Autonomous sputter synthesis of thin film nitrides with composition controlled by Bayesian optimization of optical plasma emission

Autonomous experimentation has emerged as an efficient approach to accelerate the pace of material discovery. Although instruments for autonomous synthesis have become popular in molecular and polymer science, solution processing of hybrid materials, and nanoparticles, examples of autonomous tools for physical vapor deposition are scarce yet important for the semiconductor industry. Here, we report the design and implementation of an autonomous workflow for sputter deposition of thin films with controlled composition, leveraging a highly automated sputtering reactor custom-controlled by Python, optical emission spectroscopy (OES), and a Bayesian optimization algorithm. We modeled film composition, measured by x-ray fluorescence, as a linear function of plasma emission lines monitored during co-sputtering from elemental Zn and Ti targets in an N 2 and Ar atmosphere. A Bayesian control algorithm, informed by OES, navigates the space of sputtering power to fabricate films with user-defined compositions by minimizing the absolute error between desired and measured optical emission signals. We validated our approach by autonomously fabricating Zn x Ti 1-x N y films that deviate from the targeted cation composition by a relative ±3.5%, even for 15 nm thin films, demonstrating that the proposed approach can reliably synthesize thin films with a specific composition and minimal human interference. Moreover, the proposed method can be extended to more difficult synthesis experiments where plasma intensity lines depend non-linearly on pressure, or the elemental sticking coefficients strongly depend on the substrate temperature.

36 MATERIALS SCIENCE↗

Characterization of the Sulfur-Saturated Melt Versions of the HLW HAlG Glasses

This report provides the results from the chemical analyses of a series of sulfur-saturated melt versions of the High-Level Waste High-Aluminum Glass study glasses, a series of simulated nuclear waste glasses designed and fabricated at Pacific Northwest National Laboratory. These data will be used in the development, validation, and implementation of enhanced property/composition models for waste glass vitrification at Hanford. Chemical analyses were performed on a representative sample of each of the sulfur-saturated melt versions of the glasses to allow for comparisons with targeted compositions as well as the measured compositions of the quenched glasses. The relative differences between the targeted and measured concentrations of B 2 O 3 , Cr 2 O 3 , K 2 O, Li 2 O, Na 2 O, and P 2 O 5 for several of the glasses were greater than ±10%. The relative differences between the targeted and measured concentrations of Al 2 O 3 and ZrO 2 for one of the glasses were greater than ±10%. As expected, the measured concentrations of SO 3 in most of the glasses were higher than targeted due to the use of the sulfur saturation method in fabricating these glasses. The wash solutions contained mainly sodium, sulfur, and sulfate ions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dirty dishes or dirty laundry? Comparing two methods for quantifying American consumers' preferences for load management in a smart home

One challenge of transitioning to renewable energy is that household electricity use and renewable generation are often misaligned. Smart home energy management systems hold promise for shifting usage to match generation, but these systems need to be designed with the occupants’ preferences in mind. The purpose of the present research is to compare two approaches for collecting and modeling consumers’ load management preferences, both of which are amenable to use in a home energy management system. Specifically, we examine the performance of Simple Multi-Attribute Rating Technique Exploiting Ranks (SMARTER) and Analytic Hierarchy Process (AHP) in quantifying consumers’ preferences regarding air temperature (air conditioning and heating), water heating, dishwashing, clothes washing and drying, monetary costs, environmental impacts, and comfort/convenience. Two studies are presented: Study 1 examines the SMARTER approach, and Study 2 focuses on the AHP approach. In both studies, online surveys (N SMARTER = 956 and N AHP = 1023) were conducted to elicit preferences from participants across the United States. The preferences modeled by both approaches were validated based on (a) their ability to predict participants’ choices in a Discrete Choice Experiment and (b) their convergence with previous research on load-shifting behavior. The validation procedure suggests that the SMARTER approach is superior in modeling consumers’ preferences for load management. Overall, this research lays the groundwork for designing a smart home interface capable of collecting occupants’ preferences and using those preferences to deliver improved occupant comfort, lower operating costs, reduced environmental impact, and more significant demand response than exists today.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Calcined Waste Storage at the Idaho Nuclear Technology and Engineering Center

This comprehensive report provides definitive information regarding the volume, mass, and composition (chemical and radioactivity) of the calcined waste stored at the Idaho Nuclear Technology and Engineering Center on the Idaho National Laboratory Site. The calcine composition data are required for current regulatory compliance and safe operation of the INTEC’s six Calcined Solids Storage Facilities (CSSFs), as well as future retrieval, treatment, and shipping the calcine to an off-Site facility (such as a geologic repository). This report also contains a description of the calcine storage bins that compose each of the CSSFs, emphasizing the features that may be of concern to future calcine retrieval. The CSSFs were designed by different architectural engineering firms and built at different times. Each CSSF has a unique design, reflecting varying design criteria and lessons learned from historical CSSF operation. The varying CSSF designs will affect future calcine retrieval processes and equipment. This revision (5) of the report presents refinements of the databases used to calculate the volume, mass, chemical content, and radioactivity of calcined waste produced and stored within the CSSFs. This revision corrects minor errors and inconsistencies noted in an extensive validation and verification effort of the prior version of the Microsoft Excel workbooks upon which the calcine composition described in this report is based. Historical calcine samples are insufficient in number and scope of analysis to fully characterize the entire inventory of calcine in the CSSFs. Sample data exist for all of the liquid wastes that were calcined. This report provides calcine composition data based on liquid waste sample analyses, the volume of liquid waste calcined, calciner operating data, and CSSF operating data using several large Microsoft Excel (Microsoft 2010) databases and spreadsheets. The Excel spreadsheets are collectively called the Historical Processing Model (HPM). The calcine composition determined by this method compares favorably with historical calcine sample data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗

Radiative Heat Transfer Capability Implemented in OpenNCC for Conjugate Heat Transfer Applications

Thermal efficiency of gas turbine engine increases as the temperature and pressure at the combustor increases. Consequently, the materials used inside a combustor must survive an increasingly challenging environment. For this reason, accurate assessment of heat transfer is crucial for combustor design. While all three modes of heat transfer are present inside a combustor, the focus of this paper is the thermal radiation. Radiative heat transfer in a gas turbine combustors are particularly interesting from three reasons. Firstly, the radiative heat loss from the combustion region may affect the emission performance. Secondly, the cooling air will protect the liner from convection but not necessary from radiation. Finally, it is less frequently incorporated in CFD analysis than other forms of heat transfer. In this work, radiative heat transfer using discrete ordinate method has been incorporated in OpenNCC (a publicly releasable version of the National Combustion Code) developed at NASA Glenn Research Center. Aside from massively parallel computation capability using MPI and the ability to utilize unstructured mesh, the current implementation includes two types of spectral models, namely, the weighted some of gray gas model and the full spectrum correlated k-distribution model. After presenting the theory and the strategy of implementation, results of validation cases for gray gas and spectral models will be presented. While the implementation of the radiation solver is intended for gas turbine application, the radiation solver can run independently from the convection/combustion solver and the same theory can be applied to other application.

OpenNCC↗

Field Demonstration of The Krauklis Seismic Wave In A Novel Mva Method For Geologic CO 2 Storage

The Energy & Environmental Research Center (EERC) conducted a feasibility study of a new, low-impact geophysics monitoring method designed to incrementally track the saturation front of CO 2 injected into a reservoir. The method leverages a new way of transmitting energy from the surface to the reservoir and employs a new subsurface signal called the Krauklis wave (K-wave) and other guided waves for seismic monitoring. The EERC partnered with Seismos, Inc., the company that holds the exclusive technology license for this method, to advance development of K-wave technology, and Denbury Onshore, LLC, to field-test the system. The system was field-tested at Denbury’s Bell Creek oil field, currently undergoing CO 2 -based enhanced oil recovery (EOR) in southeastern Montana. The K-wave system uses a wellhead-mounted source to propagate energy down the wellbore in the form of a tube wave and through the well perforations into the reservoir. Part of the energy transmitted through the perforations into the reservoir converts to slow moving K-waves and/or other guided waves that are expected to propagate through the formation to neighboring wells. The propagated signal would induce tube waves in the neighboring wells that return to the surface to be recorded with a wellhead-mounted receiver. The propagation time of the signal is directly affected by the presence of CO 2 between the wells, so the difference in transit time of the returning signal between pre- and postinjection data sets could be used to help define the lateral extent of a CO 2 saturation front. Objectives of the study were to demonstrate, validate, and evaluate the K-wave method’s ability to monitor the morphology of injected CO 2 and estimate saturation distribution. Results from the demonstration were to be compared to conventional 4-D surface seismic methods, dynamic reservoir simulations, and CO 2 breakthrough data to validate the technology and provide useful information that may result in improvements. Prior to mobilizing for the full-scale K-wave baseline survey, a field test was conducted to better understand wellhead hardware configurations and test the response from the Bell Creek reservoir using three wells. Useful operational information resulted but signal from the first-generation gas release/gas injection source was not recorded in neighboring wells. After the initial field test, numerical modeling of the K-wave method was initiated using parameters specific to the Bell Creek Field. Modeling results indicated that the fundamental concepts behind the K-wave method exhibited problems and revealed a significantly greater attenuation of the signal at the wellbore perforations than expected, totaling on the order of 10-8 from source well to receiver well over a 400-m well spacing. Modeling also showed that greater source frequency bandwidth was needed to produce guided waves within a 30-foot-thick reservoir layer such as that at Bell Creek and that reservoir layer seismic velocity would need to be slower than the velocities of the bounding layers in order to act as a wave guide opposite to what exists at Bell Creek. Modeling results prompted a second field test with a stronger, second-generation displacement source developed in the interim. Concurrent with the second field test, a 3-D seismic survey was acquired to serve as a baseline for project validation and to serve as a monitor data set for the portion of the project area where CO 2 injection had already commenced. The new source was tested using a four-well geometry and failed to produce a detectable signal, which appeared to confirm modeling results and indicated that the Bell Creek Field is not an ideal site for the K-wave test. It was concluded that a source with sufficient strength and bandwidth could not be built to overcome the factors discovered by the modeling within a reasonable amount of time. Several other sites undergoing CO 2 -based EOR were also evaluated; however, none of these sites was amenable because of geology, limitations of source and receiver hardware, or the advanced state of the CO 2 flood. In addition, two other configurations of the source and receiver layout were considered; however, these design modifications either changed the scope of the project or did not meet the goals of the Carbon Storage Program. In consultation with the U.S. Department of Energy (DOE) Project Manager and Technology Manager, a no-go decision was invoked to forego additional expenditures as the technology would not meet current project objectives. Project activities to date have moved the technology forward by developing modeling equations that can be applied elsewhere, advanced the state of wellhead sensors and sources, and spurred improvements to field operations. The 3-D seismic survey in the project area and the dynamic reservoir simulations done as part of this study have provided value for ongoing and future synergistic activities at Bell Creek related to DOE carbon storage goals.

54 ENVIRONMENTAL SCIENCES↗

Energy-based descriptors to rapidly predict hydrogen storage in metal–organic frameworks

The low volumetric density of hydrogen is a critical limitation to its use as a transportation fuel. Filling a fuel tank with nanoporous materials, such as metal–organic frameworks (MOFs), could greatly improve the deliverable capacity of these tanks if appropriate materials could be found. Yet, since MOFs can be made from many combinations of metal nodes, organic linkers, and functional groups, the design space of possible MOFs is enormous. Experimental characterization of thousands of MOFs is infeasible, and even conventional molecular simulations can be prohibitively expensive for large databases. In this work, we have developed a data-driven approach to accelerate materials screening and learn structure–property relationships. We report new descriptors for gas adsorption in MOFs derived from the energetics of MOF–guest interactions. Using the bins of an energy histogram as features, we trained a sparse regression model to predict gas uptake in multiple MOF databases to an accuracy within 3 g L -1 . The interpretable model parameters indicate that a somewhat weak attraction between hydrogen and the framework is ideal for cryogenic storage and release. Our machine learning method is more than three orders of magnitude faster than conventional molecular simulations, enabling rapid exploration of large numbers of MOFs. As a case study, we applied the method to screen a database of more than 50 000 experimental MOF structures. We experimentally validated one of the top candidates identified from the accelerated screening, MFU-4l. This material demonstrated a hydrogen deliverable capacity of 47 g L -1 (54 g L -1 simulated) when operating at storage conditions of 77 K, 100 bar and delivery at 160 K, 5 bar.

36 MATERIALS SCIENCE↗

Residual Stresses Affect Cell Fragment Movement

Predictive modeling tools such as the finite element method can be of tremendous help in assessing the reliability and long-term performance of photovoltaic modules. In order to obtain accurate results, the proper modeling of materials and manufacturing processes are of utmost importance. Module fabrication introduces thermo-mechanical stresses inside the module laminate, which need to be accounted for as residual stresses in finite element simulations. We found that cell fragment movement and crack opening displacements of fractured silicon cells within modules are affected by those residual stresses. Cell cracking remains a challenging topic in assessing the reliability and durability of damaged modules. Hence, accurately quantifying the separation and movement between cell fragments creates the foundation for establishing reliable lifetime and performance assessments of fractured silicon modules. Here, we present a modeling approach that uses upper and lower bounds to accurately account for the residual stresses introduced by the module lamination process. We designed a four-point flexure coupon test of a laminated, fractured silicon strip to validate our numerical results and found good agreement between our modeling methodology and the experimental data. Finally, we discuss the implications of the residual stresses on the normal crack opening and metallization wear-out of fractured silicon cells.

14 SOLAR ENERGY↗