Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Assessment of 3D Codes for Predicting Liner Attenuation in Flow Ducts

This paper presents comparisons of seven propagation codes for predicting liner attenuation in ducts with flow. The selected codes span the spectrum of methods available (finite element, parabolic approximation, and pseudo-time domain) and are collectively representative of the state-of-art in the liner industry. These codes are included because they have two-dimensional and three-dimensional versions and can be exported to NASA's Columbia Supercomputer. The basic assumptions, governing differential equations, boundary conditions, and numerical methods underlying each code are briefly reviewed and an assessment is performed based on two predefined metrics. The two metrics used in the assessment are the accuracy of the predicted attenuation and the amount of wall clock time to predict the attenuation. The assessment is performed over a range of frequencies, mean flow rates, and grazing flow liner impedances commonly used in the liner industry. The primary conclusions of the study are (1) predicted attenuations are in good agreement for rigid wall ducts, (2) the majority of codes compare well to each other and to approximate results from mode theory for soft wall ducts, (3) most codes compare well to measured data on a statistical basis, (4) only the finite element codes with cubic Hermite polynomials capture extremely large attenuations, and (5) wall clock time increases by an order of magnitude or more are observed for a three-dimensional code relative to the corresponding two-dimensional version of the same code.

Watson, W. R.↗
NASA NTRS2015

NASA Galactic Cosmic Radiation Environment Model: Badhwar-O'Neill (2014)

The Badhwar‐O'Neill (BON) Galactic Cosmic Ray (GCR) flux model is used by NASA to certify microelectronic systems and in the analysis of radiation health risks for human space flight missions. Of special interest to NASA is the kinetic energy region below 4.0 GeV/n due to the fact that exposure from GCR behind shielding (e.g., inside a space vehicle) is heavily influenced by the GCR particles from this energy domain. The BON model numerically solves the Fokker‐Planck differential equation to account for particle transport in the heliosphere due to diffusion, convection, and adiabatic deceleration under the assumption of a spherically symmetric heliosphere. The model utilizes a GCR measurements database from various particle detectors to determine the boundary conditions. By using an updated GCR database and improved model fit parameters, the new BON model (BON14) is significantly improved over the previous BON models for describing the GCR radiation environment of interest to human space flight.

O'Neill, P. M.↗
NASA NTRS2019

Atmospheric Chemistry Modeling and Air Quality Forecasting Using Machine Learning

Atmospheric chemistry models are a central tool to study the impact of chemical constituents on the environment, vegetation and human health. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O1D to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a chemistry model.We have investigated a machine learning approach to emulate the chemistry instead of solving the differential equations numerically. From a one-month simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry.This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, etc. We discuss speed and accuracy of our approach, and highlight some potential future directions for improving it.

Keller, Christoph A.↗
NASA NTRS2019

Atmospheric Chemistry Modeling Using Machine Learning

Atmospheric chemistry models are a central tool to study the impact of chemical constituents on the environment, vegetation and human health. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O1D to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a chemistry model. We have investigated a machine learning approach to emulate the chemistry instead of solving the differential equations numerically. From a one-month simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry. This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, etc. We discuss speed and accuracy of our approach, and highlight some potential future directions for improving it.

Keller, Christoph A.↗
NASA NTRS2019

Machine Learning Application to Atmospheric Chemistry Modeling

Atmospheric chemistry models are a central tool to study the impact of chemical constituents on the environment, vegetation and human health. These models split the atmosphere in a large number of grid-boxes and consider the emission of compounds into these boxes and their subsequent transport, deposition, and chemical processing. The chemistry is represented through a series of simultaneous ordinary differential equations, one for each compound. Given the difference in life-times between the chemical compounds (milli-seconds for O (sup 1) D (Deuterium) to years for CH4) these equations are numerically stiff and solving them consists of a significant fraction of the computational burden of a chemistry model. We have investigated a machine learning approach to emulate the chemistry instead of solving the differential equations numerically. From a one-month simulation of the GEOS-Chem model we have produced a training dataset consisting of the concentration of compounds before and after the differential equations are solved, together with some key physical parameters for every grid-box and time-step. From this dataset we have trained a machine learning algorithm (regression forest) to be able to predict the concentration of the compounds after the integration step based on the concentrations and physical state at the beginning of the time step. We have then included this algorithm back into the GEOS-Chem model, bypassing the need to integrate the chemistry. This machine learning approach shows many of the characteristics of the full simulation and has the potential to be substantially faster. There are a wide range of application for such an approach - generating boundary conditions, for use in air quality forecasts, chemical data assimilation systems, etc. We discuss speed and accuracy of our approach, and highlight some potential future directions for improving it.

Keller, Christoph A.↗
NASA NTRS2019

Real-Time Modeling of Thermal Interactions in Cryogenic Ball Bearings

Based on the integration of classical differential equations of motion of bearing elements, a real-time dynamic simulation of ball bearing performance is coupled with simultaneous modeling of thermal interactions between the bearing elements. The transient heat generations are time averaged over a significantly larger time step before they are used to compute the changes in bearing temperature fields. Such an averaging algorithm eliminates the numerical difficulties associated with simultaneous integration of mechanical thermal differential equations with vastly different time scales. Under stable operation, the step change in temperatures at each thermal time step converges to a steady-state solution. Following is a list of accomplishments in the present investigation: (1) A real-time dynamic performance of both 440C and hybrid ball bearings in a LOX environment is modeled in both transient and steady-state time domains, as the inner race accelerates to a prescribed operating speed. (2) The frictional interactions at the ball to race contacts are modeled by experimentally measured traction coefficients as a function of slide-to-roll ratios. (3) Bearing heat generation, as defined by the heat transferred to circulating LOX, is validated against experimental values determined by measured LOX temperatures under prescribed pressures and flow rates. The predicted heat generation is in good agreement with the experimental values. (4) The steady-state solutions are shown to be independent of initial conditions. This establishes convergence of the numerical integration of a differential equation of motion of the bearing elements. (5) Churning and drag losses in the circulating LOX constitutes the majority of bearing heat generation under the experimental test conditions. (6) The heat dissipated in frictional interactions in hybrid bearings is significantly lower in hybrid bearings in comparison to that simulated in all-steel bearings, while the churning and drag losses are relatively unchanged.

Gupta, P. K.↗
NASA NTRS1976

Tabulation of hybrid theory calculated e-N2 vibrational and rotational cross sections

Vibrational excitation cross sections of N2 by electron impact are tabulated. Integrated cross sections are given for transitions v yields v prime where o=or v=or 8 in the energy range 0.1 eV=or E=or 10 eV. The energy grid is chosen to be most dense in the resonance region (2 to 4 eV) so that the substructure is present in the numerical results. Coefficients in the angular distribution formula (differential scattering cross section) for transitions v=0 yields v prime = or 8 are also numerically given over the same grid of energies. Simultaneous rotation-vibration coefficients are also given for transitions v=o,j=o; 1 yields v prime=o, j=o,2,4; 1,3,5. All results are obtained from the hybrid theory.

Chandra, N.↗
NASA NTRS1999

A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications

In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.

Povitsky, Alex↗
NASA NTRS1977

Prediction of hydrodynamics and chemistry of confined turbulent methane-air flames with attention to formation of oxides of nitrogen

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a tubular combustor is presented. A numerical procedure for the solution of the governing differential equations is described, and models for chemical equilibrium and chemical kinetics calculations are presented. The chemical equilibrium model is used to characterize the hydrocarbon reactions. The chemical kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor consists of a cylindrical duct of varying cross sections with concentric streams of gaseous fuel and air entering the duct at one end. Four sample cases with specified inlet and boundary conditions are considered, and the results are discussed

Elghobashi, S.↗
NASA NTRS1989

Personal computer study of finite-difference methods for the transonic small disturbance equation

Calculation of unsteady flow phenomena requires careful attention to the numerical treatment of the governing partial differential equations. The personal computer provides a convenient and useful tool for the development of meshes, algorithms, and boundary conditions needed to provide time accurate solution of these equations. The one-dimensional equation considered provides a suitable model for the study of wave propagation in the equations of transonic small disturbance potential flow. Numerical results for effects of mesh size, extent, and stretching, time step size, and choice of far-field boundary conditions are presented. Analysis of the discretized model problem supports these numerical results. Guidelines for suitable mesh and time step choices are given.

Bland, Samuel R.↗
NASA NTRS1990

Personal computer study of finite-difference methods for the transonic small disturbance equation

Calculation of unsteady flow phenomena requires careful attention to the numerical treatment to the governing partial differential equations. The personal computer provides a convenient and useful tool for the development of meshes, algorithms, and boundary conditions needed to provide time accurate solution of these equations. The one-dimensional equation considered provides a suitable model for the study of wave propagation in the equations of transonic small disturbance potential flow. Numerical results for effects of mesh size, extent, and stretching, time step size, and choice of far-field boundary conditions are presented. Analysis of the discretized model problem supports these numerical results. Guidelines for suitable mesh and time step choices are given.

Bland, Samuel R.↗
NASA NTRS1991

Distributed Minimal Residual (DMR) method for acceleration of iterative algorithms

A new method for enhancing the convergence rate of iterative algorithms for the numerical integration of systems of partial differential equations was developed. It is termed the Distributed Minimal Residual (DMR) method and it is based on general Krylov subspace methods. The DMR method differs from the Krylov subspace methods by the fact that the iterative acceleration factors are different from equation to equation in the system. At the same time, the DMR method can be viewed as an incomplete Newton iteration method. The DMR method was applied to Euler equations of gas dynamics and incompressible Navier-Stokes equations. All numerical test cases were obtained using either explicit four stage Runge-Kutta or Euler implicit time integration. The formulation for the DMR method is general in nature and can be applied to explicit and implicit iterative algorithms for arbitrary systems of partial differential equations.

Lee, Seungsoo↗
NASA NTRS1987

Legendre-tau approximation for functional differential equations. II - The linear quadratic optimal control problem

The numerical scheme based on the Legendre-tau approximation is proposed to approximate the feedback solution to the linear quadratic optimal control problem for hereditary differential systems. The convergence property is established using Trotter ideas. The method yields very good approximations at low orders and provides an approximation technique for computing closed-loop eigenvalues of the feedback system. A comparison with existing methods (based on averaging and spline approximations) is made.

Ito, Kazufumi↗
NASA NTRS1977

The separated turbulent boundary layer over a wavy wall

A study and application of the fourth order spline collocation procedure, numerical solution of boundary layer like differential equations, is presented. A simple inversion algorithm for the simultaneous solution of the resulting difference equations is given. Particular attention is focused on the boundary condition representation for the spline second derivative approximations. Solutions using the spline procedure, as well as the three point finite difference method, are presented for several model problems in order to assess and improve the spline numerical scheme. Application of the resulting algorithm to the incompressible laminar self similar boundary layer equations is presented.

Polak, A.↗
NASA NTRS1985

Characteristics pertaining to a stiffness cross-coupled Jeffcott model

Rotordynamic studies of complex systems utilizing multiple degree-of-freedom analysis have been performed to understand response, loads, and stability. In order to understand the fundamental nature of rotordynamic response, the Jeffcott rotor model has received wide attention. The purpose of this paper is to provide a generic rotordynamic analysis of a stiffness cross-coupled Jeffcott rotor model to illustrate characteristics of a second order stiffness-coupled linear system. The particular characteristics investigated were forced response, force vector diagrams, response orbits, and stability. Numerical results were achieved through a fourth order Runge-Kutta method for solving differential equations and the Routh Hurwitz stability criterion. The numerical results were verified to an exact mathematical solution for the steady state response.

Spanyer, K. L.↗
NASA NTRS1991

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗