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 199 records · Page 11

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

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗

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↗

Observed Changes in the Gravity Field of Mars Due to Seasonal Atmospheric Processes

The atmosphere of Mars deposits approximately 30% of its mass at the winter pole as part of its seasonal cycle of CO2 exchange and sublimes it back to the atmosphere in the spring, thus creating an annual hemispheric cycle of mass re-distribution. Using X-band tracking data of the Mars Global Surveyor (MGS) spacecraft, we have detected the signature of changes in the low degree gravity field from March 1999 through August 2000, corresponding to about three-quarters of a Martian year. The observed variations show a general resemblance to predicted variations from general circulation models. Also observed are irregular changes that appear to be due to transient phenomena in the Martian atmosphere such as large dust storms that provide significant heat into the lower atmosphere, even in the polar regions. In addition, we have identified a change in the rate of rotation of Mars over this same period that appears somewhat smaller than anticipated, but in general agreement with general circulation models.

Smith, David E.↗

Developing a Hybrid Electric Vehicle Eco-Cooperative Adaptive Cruise Control System at Signalized Intersections.

This study develops an eco-driving strategy for hybrid electric vehicles (HEVs) in the vicinity of signalized intersections, entitled HEV Eco-Cooperative Adaptive Cruise Control at Intersections (Eco-CACC-I). The proposed system computes real-time, energy-optimized vehicle trajectories using HEV vehicle dynamics and energy consumption models. In the proposed system, a simple HEV energy model is used to compute the instantaneous fuel consumption. This HEV energy model is selected since it is general, transferable, and can be easily used to compute instantaneous energy consumption levels for HEVs without the additional input of vehicle engine data or complicated power control strategies. In addition, a vehicle dynamics model is used to capture the relationship between speed, acceleration level, and tractive/resistance forces on vehicles. The energy-optimum problem is formulated as an optimization problem with constraints, which is solved using a moving-horizon dynamic programming approach. The proposed HEV Eco-CACC-I system was tested to evaluate its performance for various speed limits, roadway grades, and signal timings. Lastly, the proposed HEV controller was implemented in a microscopic traffic simulation software to test its network-wide performance. The test results from an arterial corridor with three signalized intersections demonstrate that the proposed system can effectively reduce stop-and-go traffic in the vicinity of signalized intersections producing savings of 7.4% in energy consumption, 5.8% in traffic delay and 23% vehicle stops, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Added Value to Global Model Projections of Climate Change by Dynamical Downscaling: A Case Study over the Continental U.S. using the GISS-ModelE2 and WRF Models

Dynamical downscaling is being increasingly used for climate change studies, wherein the climates simulated by a coupled atmosphere-ocean general circulation model (AOGCM) for a historical and a future (projected) decade are used to drive a regional climate model (RCM) over a specific area. While previous studies have demonstrated that RCMs can add value to AOGCM-simulated climatologies over different world regions, it is unclear as to whether or not this translates to a better reproduction of the observed climate change therein. We address this issue over the continental U.S. using the GISS-ModelE2 and WRF models, a state-of-the-science AOGCM and RCM, respectively. As configured here, the RCM does not effect holistic improvement in the seasonally and regionally averaged surface air temperature or precipitation for the individual historical decades. Insofar as the climate change between the two decades is concerned, the RCM does improve upon the AOGCM when nudged in the domain proper, but only modestly so. Further, the analysis indicates that there is not a strong relationship between skill in capturing climatological means and skill in capturing climate change. Though additional research would be needed to demonstrate the robustness of this finding in AOGCM/RCM models generally, the evidence indicates that, for climate change studies, the most important factor is the skill of the driving global model itself, suggesting that highest priority should be given to improving the long-range climate skill of AOGCMs.

Racherla, P. N.↗

Detection and Signal Processing for Near-Field Nanoscale Fourier Transform Infrared Spectroscopy

Researchers from a broad spectrum of scientific and engineering disciplines are increasingly using scattering-type near-field infrared spectroscopic techniques to characterize materials non-destructively with nanoscale spatial resolution. However, a sub-optimal understanding of a technique's implementation can complicate data interpretation and act as a barrier to entering the field. Here the key detection and processing steps involved in producing scattering-type near-field nanoscale Fourier transform infrared spectra (nano-FTIR) are outlined. The self-contained mathematical and experimental work derives and explains: i) how normalized complex-valued nano-FTIR spectra are generated, ii) why the real and imaginary components of spectra qualitatively relate to dispersion and absorption respectively, iii) a new and generally valid equation for spectra which can be used as a springboard for additional modeling of the scattering processes, and iv) an algebraic expression that can be used to extract an approximation to the sample's local extinction coefficient from nano-FTIR. The algebraic model for weak oscillators is validated with nano-FTIR and attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectra on samples of polystyrene and Kapton and further provides a pedagogical pathway to cementing some of the technique's key qualitative attributes.

36 MATERIALS SCIENCE↗

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

Gravitational waves from current-carrying cosmic strings

Cosmic strings are predicted by many Standard Model extensions involving the cosmological breaking of a symmetry with nontrivial first homotopy group and represent a potential source of primordial gravitational waves (GWs). Present efforts to model the GW signal from cosmic strings are often based on minimal models, such as, e.g., the Nambu-Goto action that describes cosmic strings as exactly one-dimensional objects without any internal structure. In order to arrive at more realistic predictions, it is therefore necessary to consider nonminimal models that make an attempt at accounting for the microscopic properties of cosmic strings. With this goal in mind, we derive in this paper the GW spectrum emitted by current-carrying cosmic strings (CCCSs), which may form in a variety of cosmological scenarios. Our analysis is based on a generalized version of the velocity-dependent one-scale (VOS) model, which, in addition to the mean velocity and correlation length of the string network, also describes the evolution of a chiral (light-like) current. Here we are able to show, the solutions of the VOS equations imply a temporarily growing fractional cosmic-string energy density, Ω cs . This results in an enhanced GW signal across a broad frequency interval, whose boundaries are determined by the times of generation and decay of cosmic-string currents. Our findings have important implications for GW experiments in the Hz to MHz band and motivate the construction of realistic particle physics models that give rise to large currents on cosmic strings.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of vertical-tail location on the aerodynamic characteristics at subsonic speeds of a close-coupled canard configuration

The effects were studied of various vertical-tail configurations on the longitudinal and lateral directional-stability characteristics of a general research fighter model utilizing wing-body-canard. The study indicates that the addition of the high canard resulted in an increase in total lift at angles of attack above 4 deg with a maximum lift coefficient about twice as large as that for the wing-body configuration. For the wing-body (canard off) configuration, the center-line vertical tail indicates positive vertical-tail effectiveness throughout the test angle-of-attack range; however, for this configuration none of the wing-mounted vertical-tail locations tested resulted in a positive directional-stability increment at the higher angles of attack. For the wing-body-canard configuration several outboard locations of the wing-mounted vertical tails were found.

Huffman, J. K.↗

A multi-temperature TVD algorithm for relaxing hypersonic flows

In this paper, the extension of a multispecies TVD algorithm, second-order accurate for real-gas flows to a multitemperature formulation is described. The convection algorithm is coupled to internal relaxation processes, and the features of the coupling are examined. The first version consists of a three-temperature model, where translational-rotational, vibrational, and electronic energy modes are separately convected. Although several species are present, there is only one vibrational temperature in this model. The second version generalizes to a vibrational temperature for each molecular specie, with additional couplings between species. The algorithms are applied to a generic two-dimensional flow field, and results are compared with experimental observations.

Cambier, Jean-Luc↗

Satellite Geodetic Constraints On Earthquake Processes: Implications of the 1999 Turkish Earthquakes for Fault Mechanics and Seismic Hazards on the San Andreas Fault

Our principal activities during the initial phase of this project include: 1) Continued monitoring of postseismic deformation for the 1999 Izmit and Duzce, Turkey earthquakes from repeated GPS survey measurements and expansion of the Marmara Continuous GPS Network (MAGNET), 2) Establishing three North Anatolian fault crossing profiles (10 sitedprofile) at locations that experienced major surface-fault earthquakes at different times in the past to examine strain accumulation as a function of time in the earthquake cycle (2004), 3) Repeat observations of selected sites in the fault-crossing profiles (2005), 4) Repeat surveys of the Marmara GPS network to continue to monitor postseismic deformation, 5) Refining block models for the Marmara Sea seismic gap area to better understand earthquake hazards in the Greater Istanbul area, 6) Continuing development of models for afterslip and distributed viscoelastic deformation for the earthquake cycle. We are keeping close contact with MIT colleagues (Brad Hager, and Eric Hetland) who are developing models for S. California and for the earthquake cycle in general (Hetland, 2006). In addition, our Turkish partners at the Marmara Research Center have undertaken repeat, micro-gravity measurements at the MAGNET sites and have provided us estimates of gravity change during the period 2003 - 2005.

Reilinger, Robert↗

An Approach for Utilizing Power Flow Modeling for Simulations of Hybrid Electric Propulsion Systems

This paper describes an approach to creating simulations of the electric components for a hybrid electric propulsion system. The proposed modeling technique is based on power/load flow modeling and is designed to provide a modular framework that includes buses, lines, and other electrical components that can be connected together to form the electrical distribution system. The purpose of this paper is to detail an electric distribution system modeling technique and to demonstrate how these models may be integrated with turbomachinery simulations. These general modeling techniques were created to be utilized for system and control design studies. Additionally, steady-state and dynamic performance for a proposed model example is compared with data from a hardware in the loop simulation.

Chapman, Jeffryes W.↗

An Approach for Utilizing Power Flow Modeling for Simulations of Hybrid Electric Propulsion Systems

This paper describes an approach to creating simulations of the electric components for a hybrid electric propulsionsystem. The proposed modeling technique is based on power/load flow modeling and is designed to provide a modularframework that includes buses, lines, and other electrical components that can be connected together to form the electrical distribution system. The purpose of this paper is to detail an electric distribution system modeling technique and to demonstrate how these models may be integrated with turbomachinery simulations. These general modeling techniques were created to be utilized for system and control design studies. Additionally, steady-state and dynamic performance for a proposed model example is compared with data from a hardware in the loop simulation.

Power Flow↗

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗

Assess and Benchmark Magneto-Inertial Fusion (MIF) Scaling. Final Technical Report for the SNL/Rochester ALPHA Follow-on Project

This project was a follow-on to the Sandia National Laboratories (SNL) and the Laboratory for Laser Energetics (LLE) ARPA-E ALPHA project entitled “Demonstrating Fuel Magnetization and Laser Heating Tools for Low-Cost Fusion Energy”. The primary purpose of this follow-on project was to obtain additional data at the OMEGA facility to help better understand how MagLIF, a platform that has already demonstrated the scientific viability of magneto-inertial fusion, scales across a factor of 1000 in driver energy. A secondary aspect of this project was to extend simulations and analysis at SNL to cover a wider magneto-inertial fusion (MIF) parameter space and test scaling of those models across this wide range of input energies and conditions of the target. This work was successful in improving understanding of how key physics elements of MIF scales and improves confidence in setting requirements for fusion gain with larger drivers. The OMEGA experiments at the smaller scale verified the hypothesis that preheating the fuel plays a significant role in introducing wall contaminants that mix into the fuel and significantly degrade fusion performance. This contamination not only impacts target performance but the optimal input conditions for the target. However, analysis at the Z-scale showed that target performance at high preheat levels is limited by the Nernst effect, which advects magnetic flux from the hot spot, reducing magnetic insulation and consequently reduces the temperature of the fuel. The combination of MagLIF experiments at the disparate scales of OMEGA and Z along with a multiscale 3D simulation analysis has led to new insight into the physical mechanisms responsible for limiting target performance and provides important benchmarks to assess target scaling more generally for MIF schemes. Finally, in addition to the MagLIF related work, a semi-analytic model of liner driven Field Reversed Configuration (FRC) was developed that predicts the fusion gain for such systems. This model was also validated with 2D radiation magneto-hydrodynamic simulations and predicts that fusion gains of near unity could be driven by the Z machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification Testing For Solid-Element Material Models 11-19 in DYNA3D/ParaDyn

This technical report documents the creation and implementation of verification tests for solid-element material models 11 through 19 available in DYNA3D/ParaDyn. The verification tests covered all aspects of each material model, except for the Weibull distribution functionality in material models 15 and 19. General test cases were created to verify the elastic and plastic behavior of the material models. Other additional tests were developed to examine the intricacies of each material model. Each test involved the use a kinematic load case and specification of material parameters necessary to activate corresponding features of the material model. The load cases prescribed the full time history of the kinematic motion for the solid elements, and these loads are independent of the material model or element formulation. When possible, closed form solutions were then derived for each verification test in a continuum setting. The DYNA3D simulations for each test were carried out over a unit time interval, t ϵ [0, 1], and the as implemented DYNA3D response was compared to the closed form solutions evaluated at discrete points in time. A relative error measure was determined for each test to justify the proper implementation of the material model. The relative errors comparing the DYNA3D solution to the analytical solution were, in general, on the order of machine precision except where noted. This signifies the proper implementation of solid-element material models 11-19. In the development of these verification tests, six bugs were found and fixed in the source code. Additionally, this work generated eighteen DYNA3D input decks and answer extraction scripts in the DYNA3D/ParaDyn Software Quality Assurance test suite, which are comprised of a total of 285 solid-element tests. Testing for each material model utilizes two input decks and answer extraction scripts, where one focuses on the linear elastic response and the other examines the inelastic and remaining functionalities of the material model. In total, this work added 285 individual verification test problems in the DYNA3D/ParaDyn test suite.

42 ENGINEERING↗

Microstructure Validation of Graph Theory Model-Derived Cooling Rates in the Wire Arc Additive Manufacturing of ER70S-6 Steel

Wire arc additive manufacturing (WAAM) enables high-rate fabrication of large metallic components, but spatial variations in thermal history can lead to microstructural heterogeneity that requires efficient process models to evaluate. This study evaluates whether cooling rates extracted from a graph theory model (GTM)-based thermal simulation are consistent with the microstructural evolution observed in an ER70S-6 WAAM wall. Thermal histories from the model were analyzed at selected build heights, and cooling rates were extracted from the final thermal excursion through the austenite phase field. Microstructures at corresponding locations were characterized using electron backscatter diffraction (EBSD) to quantify grain size distributions, and pearlite interlamellar spacing was used as an additional indicator of cooling behavior. The modeled cooling rates were highest near the substrate and generally decreased with build height, consistent with the observed reduction in the fine grain fraction and the progressive shift in the grain size distribution as build height increased. Pearlite spacing trends also supported the modeled cooling rate variation. These results indicate that GTM-derived thermal histories can be post-processed into metallurgically meaningful cooling rate estimates for WAAM steel builds and linked to dataset specific empirical grain size distribution relationships for process–thermal history–microstructure assessment.

36 MATERIALS SCIENCE↗