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 NTRS1996

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗
NASA NTRS1947

Cones in Supersonic Flow

In the case of cones in axially symmetric flow of supersonic velocity, adiabatic compression takes place between shock wave and surface of the cone. Interpolation curves betwen shock polars and the surface are therefore necessary for the complete understanding of this type of flow. They are given in the present report by graphical-numerical integration of the differential equation for all cone angles and airspeeds.

Hantzsche, W.↗
NASA NTRS2005

A Study of Multigrid Preconditioners Using Eigensystem Analysis

The convergence properties of numerical schemes for partial differential equations are studied by examining the eigensystem of the discrete operator. This method of analysis is very general, and allows the effects of boundary conditions and grid nonuniformities to be examined directly. Algorithms for the Laplace equation and a two equation model hyperbolic system are examined.

Roberts, Thomas W.↗
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
NASA NTRS2011

Collocation and Galerkin Time-Stepping Methods

We study the numerical solutions of ordinary differential equations by one-step methods where the solution at tn is known and that at t(sub n+1) is to be calculated. The approaches employed are collocation, continuous Galerkin (CG) and discontinuous Galerkin (DG). Relations among these three approaches are established. A quadrature formula using s evaluation points is employed for the Galerkin formulations. We show that with such a quadrature, the CG method is identical to the collocation method using quadrature points as collocation points. Furthermore, if the quadrature formula is the right Radau one (including t(sub n+1)), then the DG and CG methods also become identical, and they reduce to the Radau IIA collocation method. In addition, we present a generalization of DG that yields a method identical to CG and collocation with arbitrary collocation points. Thus, the collocation, CG, and generalized DG methods are equivalent, and the latter two methods can be formulated using the differential instead of integral equation. Finally, all schemes discussed can be cast as s-stage implicit Runge-Kutta methods.

Huynh, H. T.↗
NASA NTRS2010

High-Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

It has been known for some time that Taylor series (TS) integration is among the most efficient and accurate numerical methods in solving differential equations. However, the full benefit of the method has yet to be realized in calculating spacecraft trajectories, for two main reasons. First, most applications of Taylor series to trajectory propagation have focused on relatively simple problems of orbital motion or on specific problems and have not provided general applicability. Second, applications that have been more general have required use of a preprocessor, which inevitably imposes constraints on computational efficiency. The latter approach includes the work of Berryman et al., who solved the planetary n-body problem with relativistic effects. Their work specifically noted the computational inefficiencies arising from use of a preprocessor and pointed out the potential benefit of manually coding derivative routines. In this Engineering Note, we report on a systematic effort to directly implement Taylor series integration in an operational trajectory propagation code: the Spacecraft N-Body Analysis Program (SNAP). The present Taylor series implementation is unique in that it applies to spacecraft virtually anywhere in the solar system and can be used interchangeably with another integration method. SNAP is a high-fidelity trajectory propagator that includes force models for central body gravitation with N X N harmonics, other body gravitation with N X N harmonics, solar radiation pressure, atmospheric drag (for Earth orbits), and spacecraft thrusting (including shadowing). The governing equations are solved using an eighth-order Runge-Kutta Fehlberg (RKF) single-step method with variable step size control. In the present effort, TS is implemented by way of highly integrated subroutines that can be used interchangeably with RKF. This makes it possible to turn TS on or off during various phases of a mission. Current TS force models include central body gravitation with the J2 spherical harmonic, other body gravitation, thrust, constant atmospheric drag from Earth's atmosphere, and solar radiation pressure for a sphere under constant illumination. The purpose of this Engineering Note is to demonstrate the performance of TS integration in an operational trajectory analysis code and to compare it with a standard method, eighth-order RKF. Results show that TS is 16.6 times faster on average and is more accurate in 87.5% of the cases presented.

Scott, James R.↗
NASA NTRSDate not supplied

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗
NASA NTRSDate not supplied

Responses of Microbes to Modeled Space Radiation

The built environment of spaceships is host to a microbial community that affects crew and craft alike. While the static composition of this community has been characterized and its temporal dynamics examined, the mechanisms controlling its make-up and evolutionary trajectory are not understood. Systematic analyses of microbial diversity show consistent patterns in community composition and function. Understanding these patterns' ecological origins remains a significant challenge, as it requires connecting processes at varying temporal and spatial scales. However, it is clear that the state and trajectories of microbial communities are in-part determined by their physical environment. In this regard, the spaceflight environment includes numerous interacting factors that differentiate it from Earth environments, including an altered atmospheric composition, reduced gravity (and thus altered fluid dynamics), and increased ionizing radiation. These factors impart selective pressures on microbial communities that affect their evolutionary trajectories and thus the risks and benefits these communities represent to crew and craft. The radiation environment of space leads to chronic exposure to low doses and is difficult to mimic on Earth. Thus, little is known about how microbial communities in spacecraft will respond and evolve. Therefore, given the limitations of existing studies, we aim to empirically determine how exposure to low doses of ionizing radiation for thousands of cell divisions affects rates of mutation accumulation in bacteria and the trajectory of their evolution. In this way, we will provide a critical set of data for designing safe and robust space missions. Here we discuss our progress towards this aim, including the construction of exposure facilities, our culturing and analysis approach, and preliminary data.

radiation↗
NASA NTRSDate not supplied

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗
NASA NTRSDate not supplied

Derivatives Across Adaptive-Step Integration Accurately and Efficiently

The numerical integration of ordinary differential equations (ODEs) is common in trajectory optimization. Previous studies have shown that the use of adaptive-step Runge-Kutta integration schemes can introduce errors in the derivatives due to non-differentiable behavior within the integration scheme itself. The authors previously demonstrated the use of an implicit adaptation to the Runge-Kutta methods to mitigate this issue, but at a significant computational cost. In this work, the authors demonstrate an explicit adaptation to the Runge-Kutta methods to address the issue that does not incur a significant computational cost.

Kaushik Ponnapalli↗
NASA NTRS1980

Newton's method: A link between continuous and discrete solutions of nonlinear problems

Newton's method for nonlinear mechanics problems replaces the governing nonlinear equations by an iterative sequence of linear equations. When the linear equations are linear differential equations, the equations are usually solved by numerical methods. The iterative sequence in Newton's method can exhibit poor convergence properties when the nonlinear problem has multiple solutions for a fixed set of parameters, unless the iterative sequences are aimed at solving for each solution separately. The theory of the linear differential operators is often a better guide for solution strategies in applying Newton's method than the theory of linear algebra associated with the numerical analogs of the differential operators. In fact, the theory for the differential operators can suggest the choice of numerical linear operators. In this paper the method of variation of parameters from the theory of linear ordinary differential equations is examined in detail in the context of Newton's method to demonstrate how it might be used as a guide for numerical solutions.

Thurston, G. A.↗
NASA NTRS1988

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: Stationary modes

Spatially discrete difference approximations for hyperbolic initial-boundary-value problems (IBVPs) require numerical boundary conditions in addition to the analytical boundary conditions specified for the differential equations. Improper treatment of a numerical boundary condition can cause instability of the discrete IBVP even though the approximation is stable for the pure initial-value or Cauchy problem. In the discrete IBVP stability literature there exists a small class of discrete approximations called borderline cases. For nondissipative approximations, borderline cases are unstable according to the theory of the Gustafsson, Kreiss, and Sundstrom (GKS) but they may be Lax-Richtmyer stable or unstable in the L sub 2 norm on a finite domain. It is shown that borderline approximation can be characterized by the presence of a stationary mode for the finite-domain problem. A stationary mode has the property that it does not decay with time and a nontrivial stationary mode leads to algebraic growth of the solution norm with mesh refinement. An analytical condition is given which makes it easy to detect a stationary mode; several examples of numerical boundary conditions are investigated corresponding to borderline cases.

Warming, Robert F.↗
NASA NTRS2008

Analytic Formulation and Numerical Implementation of an Acoustic Pressure Gradient Prediction

Two new analytical formulations of the acoustic pressure gradient have been developed and implemented in the PSU-WOPWOP rotor noise prediction code. The pressure gradient can be used to solve the boundary condition for scattering problems and it is a key aspect to solve acoustic scattering problems. The first formulation is derived from the gradient of the Ffowcs Williams-Hawkings (FW-H) equation. This formulation has a form involving the observer time differentiation outside the integrals. In the second formulation, the time differentiation is taken inside the integrals analytically. This formulation avoids the numerical time differentiation with respect to the observer time, which is computationally more efficient. The acoustic pressure gradient predicted by these new formulations is validated through comparison with available exact solutions for a stationary and moving monopole sources. The agreement between the predictions and exact solutions is excellent. The formulations are applied to the rotor noise problems for two model rotors. A purely numerical approach is compared with the analytical formulations. The agreement between the analytical formulations and the numerical method is excellent for both stationary and moving observer cases.

Lee, Seongkyu↗
NASA NTRS1976

Studies of aircraft differential maneuvering. Report 75-27: Calculating of differential-turning barrier surfaces. Report 75-26: A user's guide to the aircraft energy-turn and tandem-motion computer programs. Report 75-7: A user's guide to the aircraft energy-turn hodograph program

The numerical analysis of composite differential-turn trajectory pairs was studied for 'fast-evader' and 'neutral-evader' attitude dynamics idealization for attack aircraft. Transversality and generalized corner conditions are examined and the joining of trajectory segments discussed. A criterion is given for the screening of 'tandem-motion' trajectory segments. Main focus is upon the computation of barrier surfaces. Fortunately, from a computational viewpoint, the trajectory pairs defining these surfaces need not be calculated completely, the final subarc of multiple-subarc pairs not being required. Some calculations for pairs of example aircraft are presented. A computer program used to perform the calculations is included.

Kelley, H. J.↗
NASA NTRS1999

A Textbook for a First Course in Computational Fluid Dynamics

This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Zingg, D. W.↗
DOE OSTI2022

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗
NASA NTRS1986

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme.

Solomonoff, A.↗