SEARCH · Engineering Papers
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.
This Small Business Voucher CRADA will focus on wind turbine blade fatigue analysis for development of predictive life models. The parties will cooperate to conduct structural validation fatigue testing that produces high quality data enabling validation of blade damage models with sensing technologies developed by Sentient.
Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.
To address the significant productivity challenges associated with the qualification and certification (Q&C) tasks of additively manufactured (AM) parts, which have traditionally relied on rigorous post‐build inspection and testing, we propose an integrated framework that combines model‐based qualification and certification (MBQ&C) with autonomous additive manufacturing (AAM). MBQ&C employs high‐fidelity predictive models, developed within the Integrated Computational Materials Engineering (ICME) paradigm, to simulate process–structure–property–performance relationships for assessing a part’s fitness for use. Since predictive models are commonly machine learning (ML)-based or reduced-order surrogates of validated physics models, they run efficiently, enabling timely inference. In parallel, the self-driving AAM utilises ML-based adaptive, closed‐loop control strategies to avoid, mitigate, or repair defects and anomalies during fabrication, thereby increasing the likelihood of producing acceptable parts. A key feature of the combined AAM-MBQ&C framework is that predictive models explicitly incorporate defects or anomalies that persist after the build, using instance-specific data captured via in-situ sensing. This customisation enables a build‐specific assessment of fitness for use, rather than relying on nominal or generic parameters. Such individualised evaluation provides a robust basis for Q&C-related acceptance decisions relating to each build. Additionally, the rapid solution capabilities of ML or reduced-order models enable the determination of a part’s suitability for service shortly after build completion. As the framework matures, it has the potential to substantially reduce reliance on conventional point‐design approaches—such as time‐consuming post‐build computed tomography scanning and costly destructive testing. Thus, the AAM-MBQ&C framework represents a transformative, scalable strategy for quality assurance of AM components, as parts produced within a stable, validated, and certified envelope can be certified with reduced testing. Key benefits include: (1) significant gains in Q&C productivity through efficient, model-centric assessment; (2) performance-based classification of defects into critical and non-critical categories; (3) the ability to predict potential deviations in the performance of parts affected by real-time, adaptive process control interventions relative to those produced under a certified process, and (4) the enabling of virtual Q&C for service environments that are difficult, hazardous, or impractical to access or reproduce experimentally. Collectively, these capabilities strengthen the business case for AM, particularly for high‐consequence and mission‐critical applications. Finally, although this work focuses on powder-based AM, the proposed techniques could be extended to AM processes employing alternative feedstock forms.
A phenomenological load interaction model referred to as the Multi-Parameter Yield Zone model is presented which accounts for crack growth retardation caused by previous overloads, acceleration due to current overloads, and underload effects. In the present model, load interaction effects are calculated utilizing a residual stress intensity concept. Crack growth retardation and acceleration are accounted for by decreasing or increasing, respectively, the effective stress ratio used in a modified Forman's (1967) crack growth equation. Results of spectrum crack growth predictions for 2219-T851 aluminum are presented.
This study aims to investigate and establish a suitable model that can help to estimate aerosol optical depth (AOD) in order to monitor aerosol variations especially during non-retrieval time. The relationship between actual ground measurements (such as air pollution index, visibility, relative humidity, temperature, and pressure) and AOD obtained with a CIMEL sun photometer was determined through a series of statistical procedures to produce an AOD prediction model with reasonable accuracy. The AOD prediction model calibrated for each wavelength has a set of coefficients. The model was validated using a set of statistical tests. The validated model was then employed to calculate AOD at different wavelengths. The results show that the proposed model successfully predicted AOD at each studied wavelength ranging from 340 nm to 1020 nm. To illustrate the application of the model, the aerosol size determined using measure AOD data for Penang was compared with that determined using the model. This was done by examining the curvature in the ln [AOD]-ln [wavelength] plot. Consistency was obtained when it was concluded that Penang was dominated by fine mode aerosol in 2012 and 2013 using both measured and predicted AOD data. These results indicate that the proposed AOD prediction model using routine measurements as input is a promising tool for the regular monitoring of aerosol variation during non-retrieval time.
NorthWest Research Associates (NWRA) has developed an inverse model for inverting landing aircraft vortex data. The data used for the inversion are the time evolution of the lateral transport position and vertical position of both the port and starboard vortices. The inverse model performs iterative forward model runs using various estimates of vortex parameters, vertical crosswind profiles, and vortex circulation as a function of wake age. Forward model predictions of lateral transport and altitude are then compared with the observed data. Differences between the data and model predictions guide the choice of vortex parameter values, crosswind profile and circulation evolution in the next iteration. Iterations are performed until a user-defined criterion is satisfied. Currently, the inverse model is set to stop when the improvement in the rms deviation between the data and model predictions is less than 1 percent for two consecutive iterations. The forward model used in this inverse model is a modified version of the Shear-APA model. A detailed description of this forward model, the inverse model, and its validation are presented in a different report (Lai, Mellman, Robins, and Delisi, 2007). This document is a User's Guide for the Wake Vortex Inverse Model. Section 2 presents an overview of the inverse model program. Execution of the inverse model is described in Section 3. When executing the inverse model, a user is requested to provide the name of an input file which contains the inverse model parameters, the various datasets, and directories needed for the inversion. A detailed description of the list of parameters in the inversion input file is presented in Section 4. A user has an option to save the inversion results of each lidar track in a mat-file (a condensed data file in Matlab format). These saved mat-files can be used for post-inversion analysis. A description of the contents of the saved files is given in Section 5. An example of an inversion input file, with preferred parameters values, is given in Appendix A. An example of the plot generated at a normal completion of the inversion is shown in Appendix B.
Extended (> 50 years) dry storage is being evaluated by the U.S. Department of Energy (DOE) for the disposition of ~ 18 metric tons of aluminum-clad spent nuclear fuel (ASNF). Transition of the current ASNF inventory into dry storage—using the standard DOE canister—necessitates a rigorous, predictive understanding of the long-term physical and chemical factors that may influence the integrity of the proposed storage canister, including radiolytic molecular hydrogen (H 2 ) generation. Current model predictions employ initial radiolytic yields of H 2 , the values of which change as the cladding’s H 2 -precursor inventory is depleted and H 2 itself becomes progressively more involved in radiolytic and surface dissociation processes. Consequently, the absorbed radiation dose that this steady-state H 2 yield corresponds to is essential for the evaluation and improvement of model predictions. Here, we report our findings on the long-term generation of H 2 from the gamma irradiation (≤ 36 MGy) of corroded AA1100 and AA6061 coupons in helium environments at ambient temperature and ~ 50% RH. Our findings show that AA1100 systems reached steady-state by ~ 36 MGy, while higher doses were necessary for AA6061 systems. This discrepancy was attributed to the AA6061 coupons developing a thicker corrosion layer that led to the trapping of H 2 and its precursors, and potentially additional chemistries, ultimately delaying the depletion of H 2 precursors and the system’s “roll over” point. Further, current model predictions—based on previous AA1100 data—do not show steady-state attainment until above 120 MGy, which is not the case for the AA1100 data collected here. Consequently, the new alloy dependent data presented here are important for the continued improvement of predictive computer models for evaluating the feasibility of extended storage of ASNF in helium backfilled canisters.
Newman crack-closure model and the relevant crack growth program were applied to the analysis of crack growth under constant amplitude and aircraft spectrum loading on a number of aluminum alloy materials. The analysis was performed for available test data of 2219-T851, 2024-T3, 2024-T351, 7075-T651, 2324-T39, and 7150-T651 aluminum materials. The results showed that the constraint factor is a significant factor in the method. The determination of the constraint factor is discussed. For constant amplitude loading, satisfactory crack growth lives could be predicted. For the above aluminum specimens, the ratio of predicted to experimental lives, Np/Nt, ranged from 0.74 to 1.36. The mean value of Np/Nt was 0.97. For a specified complex spectrum loading, predicted crack growth lives are not in very good agreement with the test data. Further effort is needed to correctly simulate the transition between plane strain and plane stress conditions, existing near the crack tip.
This paper introduces a novel physics-informed multi-agents constitutive model to propose prediction in quasi-static constitutive behavior of cross-linked elastomer and the loss of mechanical performance during environmental aging. The presented model is used to simulate the effect of single-mechanism chemical aging (i.e. thermal-inducedor hydrolytic aging) on the behavior of the material in this hybrid framework. Those environmental single-mechanism damages change the polymer matrix over time due to massive chain scission, chain formations, and changing the arrangement of molecules in the polymer matrix. Here we propose a data-driven super-constrained machine-learned engine to represent damage in the polymer matrix and capture the changes in material behavior, including its inelastic features such as Mullins effect and permanent set in the course of aging. We have simplified the 3D stress–strain tensor mapping problem into a small number of super-constrained 1D mapping problems by means of a sequential order reduction. An assembly of multiple replicated conditional neural-network learning-agents (L-agents) is trained to systematically simplify the high-dimensional mapping problem into multiple 1D problems, each represented by a different type of agent. Our hybrid framework is designed to capture the effect of deformation history, aging time, and aging temperature. The model is validated with respect to a comprehensive set of experiments specifically designed to benchmark model capabilities and also against available data in the literature. Thermodynamic consistency and frame independency have been verified. Besides acceptable predictive abilities, a significant reduction of computational cost to predict behavior at multiple states of deformation is the most significant feature of this model.
Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water from ab initio molecular dynamics (AIMD) has been a challenging task. The difficulty mainly comes from a lack of accurate and efficient exchange–correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange–correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice for describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. In conclusion, our AIMD simulations indicate that hydroxide diffuses more slowly than hydronium in water, which is consistent with the experimental results.
A model is formulated for estimating the performance and chamber heat transfer in rocket injectors/chambers operating with gaseous H2-O2 propellants. The model quantifies the combustion performance and chamber heat flux for variables such as chamber length, element type, element area ratio, impingement angle, thrust/element, mixture ratio, moment ratio, element spacing, and physical size. Design equations are given and curves are plotted for evaluation of combustion performance in injectors comprised of F-O-F triplet, premix, coaxial and swirl coaxial element types. Curve plots and equations are also included for estimation of the chamber wall heat fluxes generated by these element types.
Modeling results reported here show that precipitating auroral electrons with sufficient energy to be consistent with the Voyager UVS observations produce bremsstrahlung X-rays with sufficient energy and intensity to be detected by the Solar Flare X-ray and Cosmic-Ray-Burst Instrument on board the Ulysses spacecraft. The detection of such bremsstrahlung X-rays at Jupiter would provide strong evidence for the electron-precipitation mechanism, although it does not rule out the possibility of some heavy ion involvement, and thus makes a significant contribution toward solving the mystery of the Jovian aurora.
Models for Orion parachute performance are based on reconstructions of the Capsule Parachute Assembly System (CPAS) drop test campaign and were documented in the CPAS “Model Memo.” Experience with similar Commercial Crew Program (CCP) parachute systems resulted in some updates to the Orion models in preparation for Artemis missions. The reefing cutter dispersion model for the drogues and mains had been overly-conservative by producing wide timing differences within clusters. A higher-fidelity timing model was generated by separating out in-lot variation and temperature effects. The main parachute inflation model had accounted for some correlations between parameters using complicated 2-D geometric bounding, but the results tended to exaggerate individual peak loads from fast (leading) inflations and under-emphasize actual lagging experience. Several flight tests were reconstructed again with an emphasis on matching peak load magnitudes using a search algorithm. A simpler method for generating inflation parameters uses the 3-D correlated reconstructed “samples” with some random “jitter” applied. Dispersed Monte Carlo inputs are then checked against flight test data to evaluate whether they represent reality.
Statistical analysis of test data is currently the preferred method of creating models for uncharacterized side loading in rocket nozzles. However, statistical analysis requires multiple data sets to be reliable; getting these sets requires building and testing multiple rocket motors, which is not feasible given the incredible time and cost required to do so. By using a combination of geometry, simplified finite element analysis (FEA) analogs, and algebra we can create a limited predictive model from only one test data set. Once created, the model requires only the input from one strain gauge located at each yaw and pitch attachment bracket for the actuation system. This process is not computationally intensive, therefore saving an abundance of time and budget over using full 3D FEA models to develop a similar predictive model. To develop the model, conversion factors between the strain at the gauge locations and corresponding actuator load is derived from the data. These conversion factors are used to create a “conversion ellipse” that bounds the constitutive properties of the nozzle as one lumped factor. Using a 2D cross-sectional analog of the nozzle exit plane, we create a function that allows us to determine the azimuthal location of the unexpected load based on the ratio between yaw and pitch strain readings. Combining this with the “conversion ellipse” yields a function that given the strain reading outputs an induced actuator load, regardless of the azimuthal location of the unexpected load. As previously mentioned, this predictive model is limited in its application as there is inaccuracy introduced by the use of data by only one test as well as simplifying assumptions. This approach makes a useful tool to quickly determine whether uncharacterized actuator loading events can damage actuation systems or induce loading that exceeds design limits.
Cholestatic liver injury is frequently associated with drug inhibition of bile salt transporters, such as the bile salt export pump (BSEP). Reliable in silico models to predict BSEP inhibition directly from chemical structures would significantly reduce costs during drug discovery and could help avoid injury to patients. We report our development of classification and regression models for BSEP inhibition with substantially improved performance over previously published models. We assessed the performance effects of different methods of chemical featurization, data set partitioning, and class labeling and identified the methods producing models that generalized best to novel chemical entities.
Alluvial aquifers by nature are complex caused by varied depositional environments. Developing a reliable groundwater model to represent an alluvial aquifer is non-trivial. Also, relying on a single best calibrated model may not be sufficient because of an inadequate choice of model parameter values. To better understand groundwater dynamics and improve model prediction reliability, this study presents a Bayesian multi-model uncertainty quantification (BMMUQ) framework to account for model parameter uncertainty in complex alluvial groundwater modeling. The methodology was applied to the agriculturally intensive Mississippi River alluvial aquifer (MRAA), Northeast Louisiana. An aquifer architecture was first constructed using 7,259 well logs in the MRAA area which covers three fluvial deposits (alluvium, braided-stream terrace, and braided-stream terrace-loess). A 12-layer MODFLOW model was then developed to address the alluvial aquifer complexity and well calibrated through a genetic algorithm. This study quantified model parameter uncertainty in hydraulic conductivity and specific storage of sand facies. Bayesian model averaging (BMA) with the Expectation Maximization (EM) algorithm was adopted to derive posterior model weights and head variances of 50 alternative conceptual groundwater flow models, and thereby obtains BMA ensemble model predictions instead of only relying on the best calibrated conceptual model. Overall, results show that an estimated around 950 million m3 of groundwater storage loss occurs in 2015 with respect to the beginning of 2004, due to high groundwater demand for irrigation in the MRAA area. Explicitly quantifying model uncertainty can produce more reliable groundwater level predictions from BMA ensemble model. The presented groundwater modeling framework improves our understanding of the MRAA and provides a valuable tool to assist agricultural water management.
Brown dwarf and giant planet atmospheric structure and composition has been studied both by forward models and, increasingly so, by retrieval methods. While indisputably informative, retrieval methods are of greatest value when judged in the context of grid model predictions. Meanwhile retrieval models can test the assumptions inherent in the forward modeling procedure. In order to provide a new, systematic survey of brown dwarf atmospheric structure, emergent spectra, and evolution, we have constructed a new grid of brown dwarf model atmospheres. We ultimately aim for our grid to span substantial ranges of atmospheric metallilcity, C/O ratios, cloud properties, atmospheric mixing, and other parameters. Spectra predicted by our modeling grid can be compared to both observations and retrieval results to aid in the interpretation and planning of future telescopic observations. We thus present Sonora, a new generation of substellar atmosphere models, appropriate for application to studies of L, T, and Y-type brown dwarfs and young extrasolar giant planets. The models describe the expected temperature-pressure profile and emergent spectra of an atmosphere in radiative-convective equilibrium for ranges of effective temperatures and gravities encompassing 200 less than or equal to T(sub eff) less than or equal to 2400 K and 2.5 less than or equal to log g less than or equal to 5.5. In our poster we briefly describe our modeling methodology, enumerate various updates since our group's previous models, and present our initial tranche of models for cloudless, solar metallicity, and solar carbon-to-oxygen ratio, chemical equilibrium atmospheres. These models will be available online and will be updated as opacities and cloud modeling methods continue to improve.
Lignin is a random heteropolymer that has been extensively studied as a renewable source of aromatic precursors for high-value chemicals, biofuels, and bioplastics. A key challenge in lignin valorization is the structural and compositional heterogeneity of lignin feedstocks. Solvent-based approaches are commonly used to fractionate lignin to reduce this heterogeneity, but solvent selection can be challenging due to variability in lignin composition. In this work, we developed computational methods to predict good and poor organic solvents as a function of lignin composition. We analyzed 28 different linear pentamer structures, 18 from known libraries and 10 hypothetical polymers, and calculated their activity coefficients in 50 different organic solvents by using the conductor-like screening model for realistic solvents. We used these data to train a regression model that enabled the extensive investigation of the impact of solvent and monolignol compositions on predicted lignin solubility. The exhaustive exploration of solubility trends using model predictions revealed sets of solvents, identified using Kamlet–Taft parameters, that are predicted to promote lignin dissolution regardless of lignin composition. We further identified solvents expected to selectively isolate lignin fractions enriched in certain subunits. Furthermore, these results establish heuristic guidelines for solvent selection that can be used to tailor fractionation processes for lignin feedstocks of distinct composition or to design new processes that isolate fractions with higher proportions of selected subunits.