Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model selection”

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,711 records · Page 95

Modeling Subsurface Performance of a Geothermal Reservoir Using Machine Learning

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells—increasing or decreasing the fluid flow rates across the wells—and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. In this study, we describe a new approach combining reservoir modeling and machine learning to produce models that enable such a strategy. Our computational approach allows us, first, to translate sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy, and second, to find optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an “open-source” reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 h, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 s. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs.

15 GEOTHERMAL ENERGY↗

A study of the diffusional behavior of a two-phase metal matrix composite exposed to a high temperature environment

The progress of diffusion-controlled filament-matrix interaction in a metal matrix composite where the filaments and matrix comprise a two-phase binary alloy system was studied by mathematically modeling compositional changes resulting from prolonged elevated temperature exposure. The analysis treats a finite, diffusion-controlled, two-phase moving-interface problem by means of a variable-grid finite-difference technique. The Ni-W system was selected as an example system. Modeling was carried out for the 1000 to 1200 C temperature range for unidirectional composites containing from 6 to 40 volume percent tungsten filaments in a Ni matrix. The results are displayed to show both the change in filament diameter and matrix composition as a function of exposure time. Compositional profiles produced between first and second nearest neighbor filaments were calculated by superposition of finite-difference solutions of the diffusion equations.

Tenney, D. R.↗

A LEAST ABSOLUTE SHRINKAGE AND SELECTION OPERATOR (LASSO) FOR NONLINEAR SYSTEM IDENTIFICATION

Identification of parametric nonlinear models involves estimating unknown parameters and detecting its underlying structure. Structure computation is concerned with selecting a subset of parameters to give a parsimonious description of the system which may afford greater insight into the functionality of the system or a simpler controller design. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of nonlinear systems. The LASSO minimises the residual sum of squares by the addition of a 1 penalty term on the parameter vector of the traditional 2 minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudolinear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. The performance of this LASSO structure detection method was evaluated by using it to estimate the structure of a nonlinear polynomial model. Applicability of the method to more complex systems such as those encountered in aerospace applications was shown by identifying a parsimonious system description of the F/A-18 Active Aeroelastic Wing using flight test data.

Kukreja, Sunil L.↗

Two- and three-nucleon contact interactions and ground-state energies of light- and medium-mass nuclei

Classes of two-nucleon (2N) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet np scattering length and deuteron binding energy at LO, and np and pp scattering data in the laboratory-energy ranges 0–15 MeV at NLO and 0–25 MeV at N3LO. These interactions are regularized by including two Gaussian cutoffs, one for T = 0 and the other for T = 1 channels. The cutoffs are taken to vary in the ranges R 0 = (1.5 –2.3) fm and R 1 = (1.5 –3.0) fm. The 780 (1100) data points up to 15 (25) MeV energy, primarily differential cross sections, are fitted by the NLO (N3LO) models with a χ 2 / datum about 1.7 or less (well below 1.5), when harder cutoff values are adopted. As a first application, we report results for the binding energies of nuclei with mass numbers A = 3 –6 and 16 obtained with selected LO and NLO 2N models both by themselves as well as in combination with a LO three-nucleon (3N) contact interaction. The latter is characterized by a single low-energy constant that is fixed to reproduce the experimental 3 H binding energy. The inclusion of the 3N interaction largely removes the sensitivity to cutoff variations in the few-nucleon systems and leads to predictions for the 3 He and 4 He binding energies that cluster around 7.8 and 30 MeV, respectively. However, in 16 O this cutoff sensitivity remains rather strong. Finally, predictions at LO only are also reported for medium-mass nuclei with A = 40 , 48, and 90.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Metamodeling Time-Temperature Creep Parameters

Abstract There exist many time-temperature parameter (TTP) models for creep rupture prediction of components including the Larson–Miller (LM), Manson–Haferd (MH), Manson–Brown (MB), Orr–Sherby–Dorn (OSD), Manson–Succop (MS), Graham–Walles (GW), Chitty–Duval (CD), Goldhoff–Sherby (GS) models. It remains a challenge to determine which model is “best”, capable of accurate interpolation and physically realistic extrapolation of creep rupture data for a given material. In this study, metamodeling is applied to create a unified TTP metamodel that combines and regresses into twelve TTP models (eight existing and four newly derived). An analysis of the mathematical problems that exist in TTP models is provided. A matlab code is written that can: (1) calibrate the material constants of any of the twelve TTP models (using the metamodel); (2) determine the most suitable stress-parameter function; (3) and report the normalized mean square error (NMSE) of rupture predictions for a given material database. Using the metamodel, and code, a design engineer can make an intelligent selection of the “best” TTP model for creep resistant design. This process is demonstrated using four isotherms of alloy P91 creep rupture data. To assess the influence of material, further validation is performed on alloys Hastelloy X, 304SS, and 316SS. It is determined that the “best” model is dependent on material type and the quality and quantity of available data.

Engineering↗

A numerical evaluation of the ambient air temperature in the Electron-Ion Collider tunnel

The Electron-Ion Collider (EIC) is a next-generation collider-accelerator that may require consistent operating temperature conditions for the beams within the accelerator tunnels to maintain stable operation. Variations in ambient temperature within the tunnel can cause thermal expansion of beampipe and component supports and can negatively affect the tunnel equipment, impacting the stability of the beamline. Modifications will be made to the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) to create the EIC, which necessitates a temperature model that addresses these modifications. To approach this problem, the consistency of temperature changes in different tunnel sections was first evaluated by plotting RHIC tunnel temperature data at various times of the day and year. From this data, a tunnel section was selected and a 2D temperature model was created for RHIC, EIC, and EIC with added cooling configurations. Soil temperature data was analyzed to determine the maximum, average, and mode soil temperatures, which were used as boundary conditions in different temperature scenarios. Computational fluid dynamics modeling was used to create 2D temperature profiles for the configurations. From this model, the predicted temperatures indicate that further analysis is required to validate the boundary conditions and benchmark the current conditions to allow the prediction of the tunnel ambient conditions at EIC. This research can be used as a preliminary model to create an EIC tunnel cooling system that will increase the operational stability of the EIC. As a result of my work this summer, I have become familiar with computational fluid dynamics, including creating fluid dynamic simulations using ANSYS Fluent and related software. I have also learned about the project process required for planning large-scale engineering projects.

43 PARTICLE ACCELERATORS↗

Explicit Dynamic Impact Analysis of the Building 3525 Packages

This report describes the methodology used for the safety basis analysis of the Building (Bldg.) 3525 casks and presents preliminary results for a selected geometry. The specific problem analyzed is a free fall load case onto an unyielding target at various orientations. The various orientations were analyzed to ensure that the packages were subjected to the orientation that causes the most damage. The analysis approach, material modeling, finite element modeling, and results for the selected case are presented, along with a brief discussion of conservatism and uncertainties. The results for the analyzed load cases show that through selected administrative controls, no breach of the container or loss of shielding will occur in the package.

42 ENGINEERING↗

Navigating Transition-Metal Chemical Space: Artificial Intelligence for First-Principles Design

Conspectus The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom. Although first-principles density functional theory (DFT) remains the workhorse of computational chemistry in mechanism deduction and property prediction, it is of limited use here. DFT is both far too computationally costly for widespread exploration of transition-metal chemical space and also prone to inaccuracies that limit its predictive performance for localized d electrons in transition-metal complexes. These challenges starkly contrast with the well-trodden regions of small-organic-molecule chemical space, where the analytical forms of molecular mechanics force fields and semiempirical theories have for decades accelerated the discovery of new molecules, accurate DFT functional performance has been demonstrated, and gold-standard methods from correlated wavefunction theory can predict experimental results to chemical accuracy. The combined promise of transition-metal chemical space exploration and lack of established tools has mandated a distinct approach. In this Account, we outline the path we charted in exploration of transition-metal chemical space starting from the first machine learning (ML) models (i.e., artificial neural network and kernel ridge regression) and representations for the prediction of open-shell transition-metal complex properties. The distinct importance of the immediate coordination environment of the metal center as well as the lack of low-level methods to accurately predict structural properties in this coordination environment first motivated and then benefited from these ML models and representations. Once developed, the recipe for prediction of geometric, spin state, and redox potential properties was straightforwardly extended to a diverse range of other properties, including in catalysis, computational “feasibility”, and the gas separation properties of periodic metal–organic frameworks. Interpretation of selected features most important for model prediction revealed new ways to encapsulate design rules and confirmed that models were robustly mapping essential structure–property relationships. Encountering the special challenge of ensuring that good model performance could generalize to new discovery targets motivated investigation of how to best carry out model uncertainty quantification. Distance-based approaches, whether in model latent space or in carefully engineered feature space, provided intuitive measures of the domain of applicability. With all of these pieces together, ML can be harnessed as an engine to tackle the large-scale exploration of transition-metal chemical space needed to satisfy multiple objectives using efficient global optimization methods. In practical terms, bringing these artificial intelligence tools to bear on the problems of transition-metal chemical space exploration has resulted in ML-model assessments of large, multimillion compound spaces in minutes and validated new design leads in weeks instead of decades.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Software Searches for Better Spacecraft-Navigation Models

ADAPT is a computer program that searches for better mathematical models for spacecraft navigation. The task of tuning trajectory-determination models for interplanetary navigation is complex, requiring an intensive search of multiple dynamical and nondynamical models that yield trajectory solutions with minimal errors. By automating the search, ADAPT eases the task of human analysts and enables them to consider wider ranges of potential solutions. ADAPT uses genetic algorithms to search a range of relevant parameters in a user-selected design space to arrive at values for those parameters that best fit the measured spacecraft-tracking data. The user s guide for ADAPT reviews the theoretical basis of the program and presents two example applications. One example is that of selecting a solar-radiation model for the Mars Pathfinder (MPF) mission using MPF tracking data and an extended Kalman filter from prior spacecraft-navigation software. The second example is of the use of tracking data from the Stardust spacecraft mission combined with a pseudo-epoch-state batch filter and an empirical small-forces model to find improved impulse models for use during Stardust attitude adjustments.

Ely, Todd↗

Selection of Thermal Worst-Case Orbits via Modified Efficient Global Optimization

Efficient Global Optimization (EGO) was used to select orbits with worst-case hot and cold thermal environments for the Stratospheric Aerosol and Gas Experiment (SAGE) III. The SAGE III system thermal model changed substantially since the previous selection of worst-case orbits (which did not use the EGO method), so the selections were revised to ensure the worst cases are being captured. The EGO method consists of first conducting an initial set of parametric runs, generated with a space-filling Design of Experiments (DoE) method, then fitting a surrogate model to the data and searching for points of maximum Expected Improvement (EI) to conduct additional runs. The general EGO method was modified by using a multi-start optimizer to identify multiple new test points at each iteration. This modification facilitates parallel computing and decreases the burden of user interaction when the optimizer code is not integrated with the model. Thermal worst-case orbits for SAGE III were successfully identified and shown by direct comparison to be more severe than those identified in the previous selection. The EGO method is a useful tool for this application and can result in computational savings if the initial Design of Experiments (DoE) is selected appropriately.

Moeller, Timothy M.↗

Biaxial steel plated concrete constitutive models for composite structures: Implementation and validation

The Steel-plated Concrete (SC) technique is an alternative construction technique with faster onsite construction speed, reduced construction time, and increased structural performance. Aiming to predict the force transfer mechanism of SC elements, an innovative constitutive model package is proposed by implementing experimental-based biaxial steel plate concrete models into the nonlinear finite element (FE) model “Membrane Model of SC (MM-CS).” First, the formulation and implementation of the MM-SC is illustrated in detail, including the equilibrium and compatibility equations and the implementation of constitutive models. Next, various experimental data of SC members are selected and simulated using the proposed MM-SC model. In this research, two types of structures, SC panels and framed SC shear walls, and three types of loading conditions, including uniaxial compression, pure shear, and combined axial-flexural-shear tests, are analyzed respectively. Good agreements were obtained between the reported results and the FE simulation results in terms of yield capacity and ultimate capacity, proving the reliability of the implemented analytical constitutive material models and the biaxial membrane model formulations.

42 ENGINEERING↗

A large-scale comparison of Artificial Intelligence and Data Mining (AI&DM) techniques in simulating reservoir releases over the Upper Colorado Region

In recent years, the Artificial Intelligence and Data Mining (AI&DM) models have become popular tools in assisting various aspects of reservoir operation. However, the practical uses are still rarely reported. Comparison experiment of many AI&DM models over a large number of reservoir cases is particularly valuable to help reservoir operators first examine the usefulness and transferability of different AI&DM models, and then identify the most stable and reliable AI&DM model in assist of various decision-making processes. In this study, a total of 12 AI&DM models with different parameterizations and simulation scenarios are comprehensively tested out and compared in simulating the controlled reservoir outflows of 33 reservoir cases over the Upper Colorado Region, United States. Results show that the Random Forecast and the Long-Short-Term-Memory model could consistently derive the best statistical performance than other models under the baseline simulation scenario. The employed AI&DM models could obtain satisfactory statistical interquartile ranges (25–75%) between [0.6–0.9], [0.3–0.8], and [0.2–0.8], for CORR, NSE, and KGE measurements, respectively, and [1.5–6.5], [–15 to 20], and [0.5–8.5] for the normalized RMSE, PBIAS and RSR measurements, respectively. Results also show Multi-Layer Perceptron model and Extreme Gradient Boosting Tree Algorithm produced more stable and superior performance than other models under more complex input scenarios. We also found that the performance of different AI&DM models are closely relevant to the reservoir elevations, sizes, and functionalities. Discussions were made about the sensitivity of AI&DM models’ parameterizations and the key advantages of AI&DM models over the rule-based reservoir models. We further identify that the main advantage of AI&DM models is the flexibility in designing input structures, whereas the rule-based simulation model is rather limited. Future studies were suggested regarding the best way reservoir operators and researchers could use, select, and apply different AI&DM models in simulating reservoir releases under different natural and modeling environments. Finally, this comparison study also serves as a reference and a piece of groundwork for further promoting the practical uses of AI&DM models in assisting reservoir operation.

54 ENVIRONMENTAL SCIENCES↗

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗

Development of Electro-chemical Battery Model for Plug-and-Play Eco-system Library

Energy storage components are fundamental to the concept of an Integrated Energy System (IES). They serve to store surplus energy during low-demand periods for later release when other IES components (i.e., Secondary Energy Source, Balance of Plant, etc.) would otherwise have to operate flexibly. This provision for storage avoids high-amplitude power ramps in these components thereby limiting thermal and mechanical stresses to their internals and providing for extended service life. This report describes a dynamic model that has been developed for an electrochemical battery. The lithium-ion (Li-ion) cell was selected as representative technology. The battery model was developed in the Dymola simulation environment and meets the requirements of the ecosystem plug-and-play library. The model accurately describes the electric dynamic response of a Li-ion battery for an imposed charging/discharging power profile. The corresponding physical limitations related to over-power scenarios, and the impact of the residual state of charge are accounted for in the model. A literature review of the major degradation processes affecting Li-ion batteries was performed. Given the purposes of the CTD-IES project, the progressive fade of the installed capacity, the reduction of the round-trip efficiency, and the limits on the number of charging/discharging cycles are aspects that need to be taken into account in techno-economic analyses. The modeling of these degradation phenomena becomes crucial when predictions over long time horizons (capacity expansion) are made. For each one of these phenomena, a brief description is given, and some figures to be implemented in the HERON optimization algorithm are presented.

25 ENERGY STORAGE↗

FFT-local gravimetric geoid computation

Model computations show that changes of sampling interval introduce only 0.3 cm changes, whereas zero padding provides an improvement of more than 5 cm in the fast Fourier transformation (FFT) generated geoid. For the Global Positioning System (GPS) survey of Franklin County, Ohio, the parameters selected as a result of model computations, allow large reduction in local data requirements while still retaining the cm accuracy when tapering and padding is applied. The results are shown in tables.

Nagy, Dezso↗

Manufacturing of Thermoset Polyimide Composites by Laser Sintering

Selective Laser Sintering (SLS) is an additive manufacturing technique that builds 3D models layer by layer using a laser to selectively melt cross sections in powdered polymeric materials, following sequential slices of the computer-aided design (CAD) model. SLS generally uses thermoplastic polymeric powders such as polyamides. The resultant 3D-printed objects are often weaker in their strength compared to traditionally processed materials, due to their higher porosity. This paper described the process development of using melt-processable imide oligomers terminated with reactive 4-phenylethynylphthalic anhydride (4-PEPA) to conduct laser sintering (LS). The first successful 3D-printing of high temperature RTM370 thermoset polyimide carbon fiber composites were further post-cured to promote additional crosslinking for achieving higher temperature (T g = 370°C) capability. Another novel imide oligomer, RTM385-SLS, formulated with a complex melt viscosity [ƞ*] of ~104-105 poise is also suitable for LS. RTM385-SLS resin powder was mixed with 20-25% of hexagonal boron nitride (h-BN) and subjected to LS to print out “Green” specimens which could be further post-cured to afford a thermally conductive but electrically insulating composites with high T g of 385 °C. The cured composite specimens were then subjected to mechanical testing, thermal conductivity and porosity measurements as well as SEM characterization.

additive Manufacturing↗

Assembly of the Populus Microbiome Is Temporally Dynamic and Determined by Selective and Stochastic Factors

Recent work shows that the plant microbiome, particularly the initial assembly of this microbiome, influences plant health, survival, and fitness. Here, we characterize the initial assembly of the Populus microbiome across ten genotypes belonging to two poplar species in a common garden using 16S rRNA gene and ITS2 region amplicon sequencing of the leaf endosphere, leaf surface, root endosphere, and rhizosphere. We sampled these microbiomes three times throughout the first growing season and found that the composition of the microbiome changed dramatically over time across all plant-associated habitats and host genotypes. For archaea and bacteria, these changes were dominated by strong homogenizing selection (accounting for 29 to 62% of pairwise comparisons). However, fungal assembly was generally characterized by multiple ecological assembly processes (i.e., a mix of weak selective and dispersal processes). Interestingly, genotype, while a significant moderator of microbiome composition, generally explained less variation than sample date across plant-associated habitats. We defined a set of core genera that accounted for, on average, 36% of the microbiome. The relative abundance of this core community was consistent over time. Additionally, using source tracking modeling, we determined that new microbial taxa colonize from both aboveground and belowground sources, and combined with our ecological assembly null models, we found that both selective and dispersal processes explained the differences between exo- (i.e., leaf surface and rhizosphere) and endospheric microbiomes. Taken together, our results suggest that the initial assembly of the Populus microbiome is time-, genotype-, and habitat-dependent and is moderated by both selective and stochastic factors.

16S rRNA↗

Predicting the Activity and Selectivity of Bimetallic Metal Catalysts for Ethanol Reforming using Machine Learning

Machine learning is ideally suited for the pattern detection in large uniform datasets, but consistent experimental datasets on catalyst studies are often small. Here we demonstrate how a combination of machine learning and first-principles calculations can be used to extract knowledge from a relatively small set of experimental data. The approach is based on combining a complex machine-learning model trained on an extensive computational library of transition-state energies with simple linear regression models of experimental catalytic activities and selectivities from the literature. Using the combined model, we identify the key C–C bond scission reactions involved in ethanol reforming and perform a computational screening for ethanol reforming on monolayer bimetallic catalysts with architectures TM-Pt-Pt(111) and Pt-TM-Pt(111) (TM = 3d transition metals). The model also predicts four promising catalyst compositions for future experimental studies. In conclusion, the approach is not limited to ethanol reforming but is of general use for the interpretation of experimental observations as well as for the computational discovery of novel catalytic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗