Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Contact Multigraph Routing: Overview and Implementation

In Delay Tolerant Networking (DTN), the standard routing algorithm used to navigate time-varying networks has been Contact Graph Routing (CGR). In CGR, a globally distributed list of contacts, periods during which two DTN nodes may communicate, is used to construct a contact graph, in which contacts are vertices. A version of Dijkstra’s algorithm can then be used to find paths through this model of the timevarying network. However, since contact graphs may be large compared to the network, potentially growing with the square of the number of network nodes and linearly with the time interval represented, the resulting algorithm does not scale well with the size of the network or time. Any improvement to the routing algorithm will bring significant returns to scale. In a previous paper, we briefly introduced an alternative to the contact graph model for routing. This alternative model is based on a multigraph (a graph in which there may be multiple edges between a pair of vertices) where vertices represent network nodes instead of contacts. A version of Dijkstra’s algorithm in these multigraph models reduces the time needed to perform the same routing computations done in the existing CGR algorithm. Moreover, a modified version of Yen’s algorithm for multigraphs is included. Our variation of CGR, which we call Contact Multigraph Routing (CMR), provides an in-line replacement for the previously used pathfinding algorithms. This paper describes an implementation created based on the CMR approach, and experimental comparisons to traditional CGR are given. In addition, we explore some additional modifications to the routing pipeline traditionally assumed in CGR. These modifications range from the theoretical to the practical in terms of size and scope. We step forward our understanding of sheaftheoretic networking and describe how to model the routing pipeline using sheaves. We detail some enhanced route selection criteria that addresses some of the added complexity of DTNbased systems. We also include a future works section on future improvements and implementations that would be of service to the broader DTN community.

contact graph routing↗
NASA NTRS1960

Wind-Tunnel Investigation of a Small-Scale Model of an Aerial Vehicle Supported by Tilting Ducted Fans

A wind-tunnel investigation has been made to study the static longitudinal and lateral stability characteristics of a simplified aerial vehicle supported by ducted fans that tilt relative to the airframe. The ducts were in a triangular arrangement with one duct in front and two at the rear in order to minimize the influence of the downwash of the front duct on the rear ducts. The results of the investigation were compared with those of a similar investigation for a tandem two-duct arrangement in which the ducts were fixed (rather than tiltable) relative to the airframe, since the three-duct configuration had been devised in an attempt to avoid some of the deficiencies of the tandem fixed-duct configuration. The results of the investigation indicated that the tilting-duct arrangement had less noseup pitching moment for a given forward speed than the tandem fixed-duct arrangement. The model had less angle-of-attack instability than the tandem fixed-duct arrangement. The model was directionally unstable but had a positive dihedral effect throughout the test speed range.

Smith, Charles C., Jr.↗
DOE OSTI2023

Multi-output multilevel best linear unbiased estimators via semidefinite programming

Multifidelity forward uncertainty quantification (UQ) problems often involve multiple quantities of interest and heterogeneous models (e.g., different grids, equations, dimensions, physics, surrogate and reduced-order models). While computational efficiency is key in this context, multi-output strategies in multilevel/multifidelity methods are either sub-optimal or non-existent. In this paper we extend multilevel best linear unbiased estimators (MLBLUE) to multi-output forward UQ problems and we present new semidefinite programming formulations for their optimal setup. Not only do these formulations yield the optimal number of samples required, but also the optimal selection of low-fidelity models to use. While existing MLBLUE approaches are single-output only and require a non-trivial nonlinear optimization procedure, the new multi-output formulations can be solved reliably and efficiently. Here, we demonstrate the efficacy of the new methods and formulations in practical UQ problems with model heterogeneity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Abrasion-ablation model for neutron production in heavy ion collisions

In intermediate energy nucleus-nucleus collisions, neutron production at forward angles is observed to occur with a Gaussian shape that is centered near the beam energy and extends to energies well above that of the beam. This paper presents an abrasion-ablation model for making quantitative predictions of the neutron spectrum. To describe neutrons produced from the abrasion step of the reaction where the projectile and target overlap, we use the Glauber model and include effects of final-state interactions. We then use the prefragment mass distribution from abrasion with a statistical evaporation model to estimate the neutron spectrum resulting from ablation. Measurements of neutron production from Ne and Nb beams are compared with calculations, and good agreement is found.

NASA Discipline Radiation Health↗
NASA NTRSDate not supplied

Transfer Function Models and Sensitivity Analysis

In many situations, real or induced flaws such as tight cracks with known morphology cannot be manufactured in part configuration specimens or in real parts. Typically, fatigue cracks are manufactured in simple geometry specimens such as flat plates, dog-bone shaped flat or round specimens. If a nondestructive evaluation (NDE) technique is required to provide a reliably detectable target flaw size denoted as α_(90/95) for induced flaws in real part, then the direct method for qualifying the NDE procedure is to use the appropriate induced flaw specimens and perform probability of detection analysis using these flaws. This can be described as direct POD demonstration testing, which may follow guidelines of MIL-HDBK-1823. This paper considers a case, where induced flaws are not available in part configuration specimens. Therefore, a direct POD demonstration study cannot be undertaken. In such situation, general practice for NDE procedure qualification is to use artificial flaws in simple geometry and part configuration specimens, and induced flaws in the chosen simple geometry specimens. Signal response data is taken on all sets of artificial and induced flaws. NDE procedure testing on induced flaw in simple geometry specimen is called NDE demonstration testing here. A transfer function NDE procedure qualification method for forward case calculates predicted induced flaw size for demonstration using a chosen target flaw size. Another transfer function method for inverse case, calculates the target flaw size using a given demonstration flaw size. The transfer function analysis assumes relationship of artificial flaw signal responses in real parts and simple geometry specimens; and induced flaw responses in simple geometry specimens to induced flaws in real parts. The signal response transfer relationships model should be defined before transfer function models can be devised. Assuming that the signal response transfer relationships model is valid, forward and inverse case transfer function methods have been devised. Because of lack of signal response data from induced flaws in real part, 90/95% POD/confidence (P/C) cannot be demonstrated directly. However, the transfer function method may be assessed using simulation to evaluate whether the resulting target flaw size or demonstration flaw size provides adequate confidence to the assumed signal response transfer relationships model. Therefore, the transfer function approach is a risk assessment approach. Both the signal response transfer relationships model and the transfer function model are important in managing risk in results provided by the transfer function NDE technique qualification or assessment. The signal response transfer relationships model needs to be validated with empirical data and then transfer function model needs to be validated for desired P/C on case-by-case basis.

Nondestructive evaluation↗
NASA NTRS1993

Heat transfer in rotating serpentine passages with selected model orientation for smooth or skewed trip walls

Experiments were conducted to determine the effects of model orientation as well as buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. Turbine blades have internal coolant passage surfaces at the leading and trailing edges of the airfoil with surfaces at angles which are as large as +/- 50 to 60 degrees to the axis of rotation. Most of the previously-presented, multiple-passage, rotating heat transfer experiments have focused on radial passages aligned with the axis of rotation. Results from serpentine passages with orientations 0 and 45 degrees to the axis of rotation which simulate the coolant passages for the mid chord and trailing edge regions of the rotating airfoil are compared. The experiments were conducted with rotation in both directions to simulate serpentine coolant passages with the rearward flow of coolant or with the forward flow of coolant. The experiments were conducted for passages with smooth surfaces and with 45 degree trips adjacent to airfoil surfaces for the radial portion of the serpentine passages. At a typical flow condition, the heat transfer on the leading surfaces for flow outward in the first passage with smooth walls was twice as much for the model at 45 degrees compared to the model at 0 degrees. However, the differences for the other passages and with trips were less. In addition, the effects of buoyancy and Coriolis forces on heat transfer in the rotating passage were decreased with the model at 45 degrees, compared to the results at 0 degrees. The heat transfer in the turn regions and immediately downstream of the turns in the second passage with flow inward and in the third passage with flow outward was also a function of model orientation with differences as large as 40 to 50 percent occurring between the model orientations with forward flow and rearward flow of coolant.

Johnson, B. V.↗
NASA NTRS2014

Optimization of an Active Twist Rotor Blade Planform for Improved Active Response and Forward Flight Performance

A study was conducted to identify the optimum blade tip planform for a model-scale active twist rotor. The analysis identified blade tip design traits which simultaneously reduce rotor power of an unactuated rotor while leveraging aeromechanical couplings to tailor the active response of the blade. Optimizing the blade tip planform for minimum rotor power in forward flight provided a 5 percent improvement in performance compared to a rectangular blade tip, but reduced the vibration control authority of active twist actuation by 75 percent. Optimizing for maximum blade twist response increased the vibration control authority by 50 percent compared to the rectangular blade tip, with little effect on performance. Combined response and power optimization resulted in a blade tip design which provided similar vibration control authority to the rectangular blade tip, but with a 3.4 percent improvement in rotor performance in forward flight.

Sekula, Martin K↗
NASA NTRS1987

An investigation of the forward scattering theorem

The calculation of an EM wave's extinction loss during propagation within an inhomogeneous medium, as in active and passive remote sensing modeling, can be undertaken either through the summation of the scattering and absorption losses or through the use of the forward scattering theorem. Attention is presently given to the similarities and differences of these two approaches as a function of dielectric properties of a spherical scatterer and the incident frequency. Scattering loss is obtainable by integrating the magnitude-squared of the scattered field over a spherical surface surrounding the scatterer; the scattered field and the field within the scatterer are computed according to Mie theory.

Karam, M. A.↗
NASA NTRS2009

Cold-Flow Testing of a Proposed Integrated Center-Body Diffuser/Steam Blocker Concept for Plum Brook Station's B-2 Test Facility

The center-body diffuser (CBD) steam blocker (SB) system is a concept that incorporates a set of secondary drive nozzles into the envelope of a CBD, such that both nozzle systems (i.e., the rocket engine and the steam blocking nozzles) utilize the same supersonic diffuser, and will operate either singularly or concurrently. In this manner, the SB performs as an exhaust system stage when the rocket engine is not operating, and virtually eliminates discharge flow on rocket engine shutdown. A 2.25-percent scale model of a proposed SB integrated into a diffuser for the Plum Brook B-2 facility was constructed and cold-flow tested for the purpose of evaluating performance characteristics of various design options. These specific design options addressed secondary drive nozzle design (method of steam injection), secondary drive nozzle location relative to CBD throat, and center-body throat length to diameter (L/D) ratios. The objective of the test program is to identify the desired configuration to carry forward should the next phase of design proceed. The tested scale model can provide data for various pressure ratios; however, its design is based on a proposed B-2 spray chamber (SC) operating pressure of 4.0 psia and a steam supply pressure of 165 psia. Evaluation of the test data acquired during these tests indicate that either the discrete axial or annular nozzle configuration integrated into a CBD, with an annular throat length of 1.5 L/D at the nominal injection position, would be suitable to carry forward from the SB's perspective. Selection between these two then becomes more a function of constructability and implementation than performance. L/D also has some flexibility, and final L/D selection can be a function of constructability issues within a limited range.

Edwards, Daryl A.↗
NASA NTRS1996

Are Planetary Regolith Particles Back Scattering? Response to a Paper by M. Mishchenko

In a recent paper Mishchenko asserts that soil particles are strongly forward scattering, whereas particles on the surfaces of objects in the solar system have been inferred to be back scattering. Mishchenko suggests that this apparent discrepancy is an artifact caused by using an approximate light scattering model to analyse the data, and that planetary regolith particles are actually strong forward scatterers. The purpose of the present paper is to point out the errors in Mishchenko's paper and to show from both theoretical arguments and experimental data that inhomogencous composite particles which are large compared to the wavelength of visible light, such as rock fragments and agglutinates, can be strongly back scattering and are the fundamental scatterers in media composed of them. Such particles appear to be abundant in planetary regoliths and can account for the back scattering character of the surfaces of many bodies in the solar system. If the range of phase angles covered by a data set is insufficient, serious errors in retrieving the particle scattering properties can result whether an exact or approximate scattering model is used. However, if the data set includes both large and small phase angles, approximate regolith scattering models can correctly retrieve the sign of the particle scattering asymmetry.

Hapke, Bruce↗
NASA NTRS1979

Total aircraft flight-control system - Balanced open- and closed-loop control with dynamic trim maps

The availability of the airborne digital computer has made possible a Total Aircraft Flight Control System (TAFCOS) that uses virtually the complete nonlinear propulsive and aerodynamic data for the aircraft to construct dynamic trim maps that represent an inversion of the aircraft model. The trim maps, in series with the aircraft, provide essentially a linear feed-forward path. Basically, open-loop trajectory control is employed with only a small perturbation feedback signal required to compensate for inaccuracy in the aircraft model and for external disturbances. Simulation results for application to an automatic carrier-landing system are presented. Flight-test results for a STOL aircraft operating automatically over a major portion of its flight regime are presented. The concept promises a more rapid and straightforward design from aerodynamic principles, particularly for highly nonlinear configurations, and requires substantially less digital computer capacity than conventional automatic flight-control system designs.

Smith, G. A.↗
NASA NTRS1993

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗
DOE OSTI2021

Initialized Earth system prediction from subseasonal to decadal timescales

Initialized Earth system predictions are made by starting a numerical prediction model in a state as consistent as possible to observations, and running it forward in time for up to ten years. Skillful predictions at time slices from subseasonal to seasonal (S2S), seasonal to interannual (S2I) and seasonal to decadal (S2D) offer information useful for various stakeholders, from agriculture to water resource management, and human and infrastructure safety. In this Review, we examine the processes influencing predictability, and discuss estimates of skill across S2S, S2I and S2D timescales. There are encouraging signs that skillful predictions can be made: at S2S timescales, there has been some skill in predicting the Madden-Julian Oscillation and North Atlantic Oscillation; at S2I in predicting the El Niño-Southern Oscillation; and at S2D, in predicting variability in North Atlantic sea surface temperatures. However, challenges remain, and future work must prioritise reducing model error, more effectively communicating forecasts to users, and increasing process and mechanistic understanding that could increase predictive skill and, in turn, confidence. As numerical models progress towards Earth system models, initialized predictions are expanding to include prediction of sea-ice, air pollution, terrestrial and ocean biochemistry which can bring clear benefit to society and various stakeholders.

climate prediction↗
NASA NTRS1995

A guide to the use of the pressure disk rotor model as implemented in INS3D-UP

This is a guide for the use of the pressure disk rotor model that has been placed in the incompressible Navier-Stokes code INS3D-UP. The pressure disk rotor model approximates a helicopter rotor or propeller in a time averaged manner and is intended to simulate the effect of a rotor in forward flight on the fuselage or the effect of a propeller on other aerodynamic components. The model uses a modified actuator disk that allows the pressure jump across the disk to vary with radius and azimuth. The cyclic and collective blade pitch angles needed to achieve a specified thrust coefficient and zero moment about the hub are predicted. The method has been validated with experimentally measured mean induced inflow velocities as well as surface pressures on a generic fuselage. Overset grids, sometimes referred to as Chimera grids, are used to simplify the grid generation process. The pressure disk model is applied to a cylindrical grid which is embedded in the grid or grids used for the rest of the configuration. This document will outline the development of the method, and present input and results for a sample case.

Chaffin, Mark S.↗
NASA NTRS1997

Magnetospheric Plasmas: A Direct Measurement of the Ionospheric Source

In a paper of more than a decade ago, the estimated the strength of the ionospheric source and its ability to supply the different plasma regions of the Earth's magnetosphere was reported. The launch of the POLAR spacecraft with the Thermal Ion Dynamics Experiment (TIDE) and the active control of spacecraft potential made possible for the first time the direct measurement of low energy ions moving from the ionosphere into the lobes of the magnetotail. This paper presents data taken during the period of operation of the Plasma Source Instrument (PSI) which maintains the spacecraft potential at about 1.8V positive with respect to the ambient plasma. These data give an in-situ observation of the strength and flow direction of the ionospheric plasmas as it moves outward into the lobes of the tail. A particle trajectory model has been used to trace particles back to the location of the ionospheric source and forward to the entry point of the ions to the plasma sheet. A comparison of measured data with the predictions of the earlier modeling effort will be discussed.

Giles, B. L.↗
NASA NTRS2005

Adjoint Sensitivity Analysis of Orbital Mechanics: Application to Computations of Observables' Partials with Respect to Harmonics of the Planetary Gravity Fields

An approach is presented to the inversion of gravity fields based on evaluation of partials of observables with respect to gravity harmonics using the solution of adjoint problem of orbital dynamics of the spacecraft. Corresponding adjoint operator is derived directly from the linear operator of the linearized forward problem of orbital dynamics. The resulting adjoint problem is similar to the forward problem and can be solved by the same methods. For given highest degree N of gravity harmonics desired, this method involves integration of N adjoint solutions as compared to integration of N2 partials of the forward solution with respect to gravity harmonics in the conventional approach. Thus, for higher resolution gravity models, this approach becomes increasingly more effective in terms of computer resources as compared to the approach based on the solution of the forward problem of orbital dynamics.

gravity fields↗
NASA NTRSDate not supplied

Gigatraj: An Atmospheric Trajectory Model

Atmospheric trajectory models have a long history of success in tracking air motions in the lower stratosphere and upper troposphere over periods of up to a few days. Parcels have been traced backwards from observations to identify whatever phenomena (strong convection, volcanic eruptions, rocket launches, etc.) put their signature on them. Parcels have also been initialized at a known event and traced forward to examine their subsequent physical and chemical evolution. We describe a new trajectory model, "gigatraj," that aims to increase exibility by (a) making it straightforward to use new meteorological data sources, including those not based on regular latitude-longitude grids; (b) enabling a run-time choice of vertical coordinate system for kinematic and/or quasi-isentropic calculations; (c) allowing for the output of arbitrary meteorological products, selectable by the user and interpolated to the parcels' locations and times. The model can be run in a serial or parallel processing environment, so that large numbers of parcels can be traced in a reasonable time. Information is presented on model accuracy and performance. The former is demonstrated by runs using both test-pattern winds (comparing expected paths with actual output) and real-world winds (comparing forward and backward runs to characterize how well parcels retrace their paths). Sample cases are also shown, including a reverse domain lling (RDF) calculation illustrating a tropopause fold event. Model output can be displayed using the new Visualization And Lagrangian dynamics Immersive eXtended Reality (VALIXR) system, and an example will be shown. In addition, we describe work to incorporate a version of gigatraj into the Goddard Earth Observing System (GEOS) of NASA's Global Modeling and Assimilation O ce (GMAO) at Goddard Space Flight Center. This enables trajectory calculations to be performed within the running GEOS model at the latter's native time resolution, instead of using the winds from every few hours. It also provides access to all of GEOS's internal variables as they are calculated. This module may be useful, for example, for tracking rapid chemical changes in a Lagrangian framework.

dynamics↗
DOE OSTI2023

A geospatial risk analysis graphical user interface for identifying hazardous chemical emission sources

Background: Performing back trajectory and forward trajectory using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) is a reliable approach for assessing particle transport after release among mid-field atmospheric models. HYSPLIT has an externally facing online interface that allows non-expert users to run the model trajectories without requiring extensive training or programming. However, the existing HYSPLIT interface is limited if simulations have a large amount of meteorological data and timesteps that are not coincident. The objective of this study is to design and develop a more robust tool to rapidly evaluate hazard transport conditions and to perform risk analysis, while still maintaining an intuitive and user-friendly interface. Methods: HYSPLIT calculates forward and backward trajectories of particles based on wind speed, wind direction, and the corresponding location, timestamp, and Pasquill stability classes of the regions of the atmosphere in terms of the wind speed, the amount of solar radiation, and the fractional cloud cover. The computed particle transport trajectories, combined with the online Proton Transfer Reaction-Mass Spectrometry (PTR-MS) data (https://figshare.com/articles/dataset/ARL_Data_from_PROS_station_at_Hanford_site/19993964), can be used to identify and quantify the sources and affected area of the hazardous chemicals’ emission using the potential source distribution function (PSDF). PSDF is an improved statistical function based on the well-known potential source contribution function (PSCF) in establishing the air pollutant source and receptor relationship. Performing this analysis requires a range of meteorological and pollutant concentration measurements to be statistically meaningful. The existing HYSPLIT graphical user interface (GUI) does not easily permit computations of trajectories of a dataset of meteorological data in high temporal frequency. To improve the performance of HYSPLIT computations from a large dataset and enhance risk analysis of the accidental release of material at risk, a geospatial risk analysis tool (GRAT-GUI) is created to allow large data sets to be processed instantaneously and to provide ease of visualization. Results: The GRAT-GUI is a native desktop-based application and can be run in any Windows 10 system without any internet access requirements, thus providing a secure way to process large meteorological datasets even on a standalone computer. GRAT-GUI has features to import, integrate, and convert meteorological data with various formats for hazardous chemical emission source identification and risk analysis as a self-explanatory user interface. The tool is available at https://figshare.com/articles/software/GRAT/19426742.

97 MATHEMATICS AND COMPUTING↗