Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized additive models”

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 127 records · Page 7

Computational fluid dynamic analysis of hybrid rocket combustor flowfields

Computational fluid dynamic analyses of the Navier-Stokes equations coupled with solid-phase pyrolysis, gas-phase combustion, turbulence and radiation are performed to study hybrid rocket combustor flowfields. The computational study is closely co-ordinated with a companion experimental program using a planar slab burner configuration with HTPB as fuel and gaseous oxygen. Computational predictions agree reasonably well with measurement data of fuel regression rates and surface temperatures. Additionally, most of the parametric trends predicted by the model are in general agreement with experimental trends. The computational model is applied to extend the results from the lab-scale to a full-scale axisymmetric configuration. The numerical predictions indicate that the full-scale configuration burns at a slower rate than the lab-scale combustor under identical specific flow rate conditions. The results demonstrate that detailed CFD analyses can play a useful role in the design of hybrid combustors.

Venkateswaran, S.↗

Executable Architecture Research at Old Dominion University

Executable Architectures allow the evaluation of system architectures not only regarding their static, but also their dynamic behavior. However, the systems engineering community do not agree on a common formal specification of executable architectures. To close this gap and identify necessary elements of an executable architecture, a modeling language, and a modeling formalism is topic of ongoing PhD research. In addition, systems are generally defined and applied in an operational context to provide capabilities and enable missions. To maximize the benefits of executable architectures, a second PhD effort introduces the idea of creating an executable context in addition to the executable architecture. The results move the validation of architectures from the current information domain into the knowledge domain and improve the reliability of such validation efforts. The paper presents research and results of both doctoral research efforts and puts them into a common context of state-of-the-art of systems engineering methods supporting more agility.

Tolk, Andreas↗

Boson Star from Repulsive Light Scalars and Gravitational Waves

We study properties of boson stars consisting of ultra-light scalar dark matter with repulsive self-interactions. We investigate the origin of the maximum mass of spherically symmetric stable stars which emerges only when solving the full equations of motion in curved space-time, but not when solving the approximated Schrödinger-Newton equations. When the repulsion is weak, the backreaction of the curvature on the scalars acts as an additional source of attraction and can overcome the repulsion, resulting in a maximum star mass and compactness. We also point out that the potential in a UV completed particle physics model of light scalar dark matter is generally more complicated than the widely used φ 4 interaction. Additional interactions beyond φ 4 in the potential can dramatically change the properties of boson stars as well as modify the prospect of LIGO gravitational wave detection for binary mergers of boson stars.

Djuna Croon↗

A Time-Regularized, Multiple Gravity-Assist Low-Thrust, Bounded-Impulse Model for Trajectory Optimization

The multiple gravity assist low-thrust (MGALT) trajectory model combines the medium-fidelity Sims-Flanagan bounded-impulse transcription with a patched-conics flyby model and is an important tool for preliminary trajectory design. While this model features fast state propagation via Keplers equation and provides a pleasingly accurate estimation of the total mass budget for the eventual flight suitable integrated trajectory it does suffer from one major drawback, namely its temporal spacing of the control nodes. We introduce a variant of the MGALT transcription that utilizes the generalized anomaly from the universal formulation of Keplers equation as a decision variable in addition to the trajectory phase propagation time. This results in two improvements over the traditional model. The first is that the maneuver locations are equally spaced in generalized anomaly about the orbit rather than time. The second is that the Kepler propagator now has the generalized anomaly as its independent variable instead of time and thus becomes an iteration-free propagation method. The new algorithm is outlined, including the impact that this has on the computation of Jacobian entries for numerical optimization, and a motivating application problem is presented that illustrates the improvements that this model has over the traditional MGALT transcription.

Optimization↗

A Convolutional Neural Network for Enhancement of Multi-Scale Localization in Granular Metallic Representative Unit Cells

A convolutional neural network was used to enhance the localization of strain and stress for a generalized method of cells model of a metallic microstructure. Enhanced shear strains, measured in terms of the linear regression coefficients as a function of ground truth strains, were improved from inaccurate and uncorrelated (slope=0.003, Rsq=0.000) to accurate and well correlated (slope=0.890, Rsq=0.882) relative to ground truth (slope=1.0, Rsq=1.0). In applying the convolutional neural network, a convolutional stride of 1.0 (padding=’same’) was only modestly effective while strides of 2 or 3 were more effective yet at higher cost. Additional convolutional layers were generally more expensive than additional dense layers, often with limited benefit. The accuracy of enhanced localized shear strains and stress is expected to yield benefits for damage progression models, especially in the context of hierarchical multi-scale methods where the generalized method of cells is applied at the intermediate scale.

Machine Learning↗

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗

Vibration characteristics of Z-ring-stiffened 60 deg conical shell models of a planetary entry spacecraft

An experimental investigation of the vibration characteristics of a 60 deg conical shell model of a planetary entry vehicle is described and the results presented. Model configurations include the shell with or without one or two Z-ring stiffeners and with or without a simulated payload. Tests were conducted with the model clamped at the small diameter and with the model suspended at the simulated payload. Additionally, calculated results obtained from application of several analytical procedures reported in the literature are presented together with comparisons between experimental and calculated frequencies and meridional mode shapes. Generally, very good frequency agreement between experimental and calculated results was obtained for all model configurations. For small values of circumferential mode number, however, the frequency agreement decreased as the number of ring stiffeners increased. Overall agreement between experimental and calculated mode shapes was generally good. The calculated modes usually showed much larger curvatures in the vicinity of the rings than were observed in the experimentally measured mode shapes. Dual resonances associated with modal preference were noted for the shell without Z-ring stiffeners, whereas the addition of stiffeners produced resonances for which the model responded in two or more modes over different sections of the shell length.

Naumann, E. C.↗

A Time-Regularized Multiple Gravity-Assist Low-Thrust Bounded-Impulse Model for Trajectory Optimization

The multiple gravity assist low-thrust (MGALT) trajectory model combines the medium-fidelity Sims-Flanagan bounded-impulse transcription with a patched-conics flyby model and is an important tool for preliminary trajectory design. While this model features fast state propagation via Kepler's equation and provides a pleasingly accurate estimation of the total mass budget for the eventual flight-suitable integrated trajectory it does suffer from one major drawback, namely its temporal spacing of the control nodes. We introduce a variant of the MGALT transcription that utilizes the generalized anomaly from the universal formulation of Kepler's equation as a decision variable in addition to the trajectory phase propagation time. This results in two improvements over the traditional model. The first is that the maneuver locations are equally spaced in generalized anomaly about the orbit rather than time. The second is that the Kepler propagator now has the generalized anomaly as its independent variable instead of time and thus becomes an iteration-free propagation method. The new algorithm is outlined, including the impact that this has on the computation of Jacobian entries for numerical optimization, and a motivating application problem is presented that illustrates the improvements that this model has over the traditional MGALT transcription.

Optimization↗

Diurnal Forcing of Planetary Atmospheres

Much progress has been made on calculations of the Martian seasonal water cycle using the Mars Climate Model developed for this purpose. Two papers, documenting the model and the water transport results obtained with it have been published in the Journal of Geophysical Research - Planets. An additional paper describing results related to the evolution of the seasonal water cycle as a result of orbital changes was published in Advances in Space Research. Since that time, further studies have concentrated on the consequences of the soil adsorption required to match the observed water cycle and its relation to the stability of ground ice and other potential water reservoirs. Earth-related studies have concentrated on incorporating an efficient and realistic microphysical model into the Ames Stratospheric General Circulation Model used to simulate the spread of the ML Pinatubo and other volcanic clouds in the stratosphere. In addition, visualizations of the simulations are being incorporated into a video describing the UARS mission. A paper describing the new stratospheric aerosol microphysics package (and its consequences for volcanic cloud evolution) will be submitted in the near future. The paper will discuss the relative importance of condensation and coagulation to early particle growth and the separation of the cloud by sedimentation of the larger particles. A more general paper which highlights the observation that particle number densities did not increase dramatically after the ML Pinatubo eruption is planned. Simulations of atmospheric transport will be extended to include studies of terrestrial tropospheric tracers using the Fifth-Generation Penn State/NCAR Mesoscale Model.

Houben, Howard C.↗

Bayesian Optimization of Catalysis with In-Context Learning

Large language models (LLMs) can perform accurate classification with zero or few examples through in-context learning (ICL), allowing the model to observe query-relevant examples at inference time and eliminating the need for additional weight updates to generalize beyond its original training data. We extend this capability to regression with uncertainty estimation using frozen LLMs (e.g., GPT-4o, Gemini), enabling Bayesian optimization (BO) in natural language without explicit model training or feature engineering. We apply this to materials discovery by representing materials as synthesis and testing procedures for use in natural language prompts. This Bayesian, design-first approach prioritizes optimization toward target material properties before detailed characterization, in contrast to conventional experimental workflows that often emphasize characterization of suboptimal materials. On benchmarks like aqueous solubility and oxidative coupling of methane (OCM), BO-ICL matches or outperforms Gaussian processes. In live experiments on the reverse water–gas shift (RWGS) reaction, BO-ICL identifies multimetallic catalysts that approach equilibrium CO yield within 6 and 10 iterations from a pool of 3,700 and 360,000 candidates, respectively. Our method redefines materials representation and accelerates discovery, with broad applications across catalysis, materials science, and AI.

Calibration↗

The extratropical 40-day oscillation in the UCLA general circulation model. Part 1: Atmospheric angular momentum

Variations in atmospheric angular momentum (AAM) are examined in a three-year simulation of the large-scale atmosphere with perpetual January forcing. The simulation is performed with a version of the University of California at Los Angeles (UCLA) general circulation model that contains no tropical Madden-Julian Oscillation (MJO). In addition, the results of three shorter experiments with no topography are analyzed. The three-year standard topography run contains no significant intraseasonal AAM periodicity in the tropics, consistent with the lack of the MJO, but produces a robust, 42-day AAM oscillation in the Northern Hemisphere (NH) extratropics. The model tropics undergoes a barotropic, zonally symmetric oscillation, driven by an exchange of mass with the NH extratropics. No intraseasonal periodicity is found in the average tropical latent heating field, indicating that the model oscillation is dynamically rather than thermodynamically driven. The no-mountain runs fail to produce an intraseasonal AAM oscillation, consistent with a topographic origin for the NH extratropical oscillation in the standard model. The spatial patterns of the oscillation in the 500-mb height field, and the relationship of the extratropical oscillation to intraseasonal variations in the tropics, will be discussed in Part 2 of this study.

Marcus, S. L.↗

The Effect of Enhanced Diabatic Heating on Stratospheric Circulation

The objective of this research focuses on the stratospheric dynamical response to the increase in aerosol loading and subsequent enhanced diabatic heating resulting from the eruption of Mt. Pinatubo. The Langley research Center three dimensional general circulation model and modifications made to that model for this study are described (addition of hydrogen fluoride tracer and diabatic heating enhancement). Unperturbed hydrogen fluoride distribution is compared to the hydrogen fluoride distribution measured by HALOE. A comparison of control and perturbed model runs is presented.

Kleb, Mary M.↗

Seasonal-to-Interannual Variability and Land Surface Processes

Atmospheric chaos severely limits the predictability of precipitation on subseasonal to interannual timescales. Hope for accurate long-term precipitation forecasts lies with simulating atmospheric response to components of the Earth system, such as the ocean, that can be predicted beyond a couple of weeks. Indeed, seasonal forecasts centers now rely heavily on forecasts of ocean circulation. Soil moisture, another slow component of the Earth system, is relatively ignored by the operational seasonal forecasting community. It is starting, however, to garner more attention. Soil moisture anomalies can persist for months. Because these anomalies can have a strong impact on evaporation and other surface energy fluxes, and because the atmosphere may respond consistently to anomalies in the surface fluxes, an accurate soil moisture initialization in a forecast system has the potential to provide additional forecast skill. This potential has motivated a number of atmospheric general circulation model (AGCM) studies of soil moisture and its contribution to variability in the climate system. Some of these studies even suggest that in continental midlatitudes during summer, oceanic impacts on precipitation are quite small relative to soil moisture impacts. The model results, though, are strongly model-dependent, with some models showing large impacts and others showing almost none at all. A validation of the model results with observations thus naturally suggests itself, but this is exceedingly difficult. The necessary contemporaneous soil moisture, evaporation, and precipitation measurements at the large scale are virtually non-existent, and even if they did exist, showing statistically that soil moisture affects rainfall would be difficult because the other direction of causality - wherein rainfall affects soil moisture - is unquestionably active and is almost certainly dominant. Nevertheless, joint analyses of observations and AGCM results do reveal some suggestions of land-atmosphere feedback in the observational record, suggestions that soil moisture can affect precipitation over seasonal timescales and across certain large continental areas. The strength of this observed feedback in nature is not large but is still significant enough to be potentially useful, e.g., for forecasts. This talk will address all of these issues. It will begin with a brief overview of land surface modeling in atmospheric models but will then focus on recent research - using both observations and models - into the impact of land surface processes on variability in the climate system.

Koster, Randal↗

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.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

Some Recent Developments in Numerical Modelling at ECMWF

A new atmospheric model was introduced into operational forecasting at the European Center for Medium Range Weather Forecasts, on 21 April 1983. The principal differences between this model and the Center's first operational model were in the adiabatic formulation, which, in the new model, includes use of a spectral representation in the horizontal, a more general vertical coordinate, and a modified, more-efficient, time-stepping scheme. In addition, new programming techniques and standards were adopted to facilitate both the model's use as a research tool and its adaptation to make full use of the features of the recently-acquired CRAY X-MP computer. A number of revisions were also made in detailed aspects of the formulation of the parameterization schemes. The operational change to this new model was accompanied by a second important change, namely the use of a higher envelope orography in the lower boundary conditions of the model. Problems in the operational performance of the model and development preparatory to the future operational implementation of a high resolution version of the model on the CRAY X-MP are discussed.

Simmons, A. J.↗

Effects of Boattail Area Contouring and Simulated Turbojet Exhaust on the Loading and Fuselage-tail Component Drag of a Twin-engine Fighter-type Airplane Model

An investigation of a twin-engine fighter-type airplane model has been conducted in the Langley 16-foot transonic tunnel to determine the effect on drag of a fuselage volume addition incorporating streamline contouring and more extensive boattailing of the engine shrouds. The effect of hot exhausts from the turbojet engines was simulated with hydrogen peroxide gas generators using scaled nonafterburning engine nozzles. Afterbody pressure distributions, base drag coefficients, and forces on the fuselage-tail configurations are presented at Mach numbers from 0.80 to 1.05 angles of attack of 0 degree and 4 degrees for jet pressure ratios from 1 to 7. The effect of jet operation on both the basic and modified models was generally to decrease base pressures but to increase most other afterbody pressures and, therefore, to result in an overall decrease in fuselage-tail component drag. The addition of volume to the basic model reduced the base drag coefficient by 0.0010 with the jets off and 0.0018 at a typical cruise operating condition of a jet pressure ratio of 3, a Mach number of 0.85, and an angle of attack of 4 degrees. The overall jet-off reduction in fuselage-tail component drag due to the volume addition was a maximum of 0.0040 at a Mach number of 0.90 for an angle of attack of 4 degrees.

Foss, Willard E , Jr↗