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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Theoretical kinetic computations in complex reacting systems

Nasa Lewis' studies of complex reacting systems at high temperature are discussed. The changes which occur are the result of many different chemical reactions occurring at the same time. Both an experimental and a theoretical approach are needed to fully understand what happens in these systems. The latter approach is discussed. The differential equations which describe the chemical and thermodynamic changes are given. Their solution by numerical techniques using a detailed chemical mechanism is described. Several different comparisons of computed results with experimental measurements are also given. These include the computation of (1) species concentration profiles in batch and flow reactions, (2) rocket performance in nozzle expansions, and (3) pressure versus time profiles in hydrocarbon ignition processes. The examples illustrate the use of detailed kinetic computations to elucidate a chemical mechanism and to compute practical quantities such as rocket performance, ignition delay times, and ignition lengths in flow processes.

Bittker, David A.↗

NASA NTRS1987

Damage growth in composite laminates with interleaves

The influence of placing interleaves between fiber reinforced plies in multilayered composite laminates is investigated. The geometry of the composite is idealized as two dimensional, isotropic, linearly elastic media made of a damaged layer bonded between two half planes and separated by thin interleaves of low extensional and shear moduli. The damage in the layer is taken in the form of a symmetric crack perpendicular to the interface and may extend up to the interface. The case of an H-shaped crack in the form of a broken layer with delamination along the interface is also analyzed. The interleaves are modeled as distributed shear and tension springs. Fourier integral transform techniques are used to develop solutions in terms of singular integral equations. An asymptotic analysis of the integral equations based on Muskhelishvili's techniques reveals logarithmically singular axial stresses in the half plane at the crack tips for the broken layer. For the H shaped crack, similar singularities are found to exist in the axial stresses at the interface crack tips in the layer and the half plane. The solution of the equations is found numerically for the stresses and displacements by using the Hadamard's concept of direct differentiation of Cauchy integrals as well as Gaussian integration techniques.

Goree, James G.↗

NASA NTRS1988

Nonlinear Analysis Of Rotor Dynamics

Study explores analytical consequences of nonlinear Jeffcott equations of rotor dynamics. Section 1: Summary of previous studies. Section 2: Jeffcott Equations. Section 3: Proves two theorems that provide inequalities on coefficients of differential equations and magnitude of forcing function in absence of side force. Section 4: Numerical investigation of multiple-forcing-function problem by introducing both side force and mass imbalance. Section 5: Examples of numberical solutions of complex generalized Jeffcott equation with two forcing functions of different frequencies f1 and f2. Section 6: Boundedness and stability of solutions.Section 7: Concludes report reviewing analytical results and significance.

Day, William B.↗

NASA NTRS1989

Numerical studies of convective heat transfer in an inclined semiannular enclosure

Natural convection heat transfer in a two-dimensional differentially heated semiannular enclosure is studied. The enclosure is isothermally heated and cooled at the inner and outer walls, respectively. A commercial software based on the SIMPLER algorithm was used to simulate the velocity and temperature profiles. Various parameters that affect the momentum and heat transfer processes were examined. These parameters include the Rayleigh number, Prandtl number, radius ratio, and the angle of inclination. A flow regime extending from conduction-dominated to convection-dominated flow was examined. The computed results of heat transfer are presented as a function of flow parameter and geometric factors. It is found that the heat transfer rate attains a minimum when the enclosure is tilted about +50 deg with respect to the gravitational direction.

Wang, Lin-Wen↗

NASA NTRS1989

Nonlinear mechanical behavior of thermoplastic matrix materials for advanced composites

Two recent theories of nonlinear mechanical response are quantitatively compared and related to experimental data. Computer techniques are formulated to handle the numerical integration and iterative procedures needed to solve the associated sets of coupled nonlinear differential equations. Problems encountered during these formulations are discussed and some open questions described. Bearing in mind these cautions, the consequences of changing parameters that appear in the formulations on the resulting engineering properties are discussed. Hence, engineering approaches to the analysis of thermoplastic matrix material can be suggested.

Arenz, R. J.↗

NASA NTRS1996

Fast Multigrid Techniques in Total Variation-Based Image Reconstruction

Existing multigrid techniques are used to effect an efficient method for reconstructing an image from noisy, blurred data. Total Variation minimization yields a nonlinear integro-differential equation which, when discretized using cell-centered finite differences, yields a full matrix equation. A fixed point iteration is applied with the intermediate matrix equations solved via a preconditioned conjugate gradient method which utilizes multi-level quadrature (due to Brandt and Lubrecht) to apply the integral operator and a multigrid scheme (due to Ewing and Shen) to invert the differential operator. With effective preconditioning, the method presented seems to require Omicron(n) operations. Numerical results are given for a two-dimensional example.

Oman, Mary Ellen↗

NASA NTRS2004

Continuation Methods for Qualitative Analysis of Aircraft Dynamics

A class of numerical methods for constructing bifurcation curves for systems of coupled, non-linear ordinary differential equations is presented. Foundations are discussed, and several variations are outlined along with their respective capabilities. Appropriate background material from dynamical systems theory is presented.

Cummings, Peter A.↗

NASA NTRS2007

Enhanced Elliptic Grid Generation

An enhanced method of elliptic grid generation has been invented. Whereas prior methods require user input of certain grid parameters, this method provides for these parameters to be determined automatically. "Elliptic grid generation" signifies generation of generalized curvilinear coordinate grids through solution of elliptic partial differential equations (PDEs). Usually, such grids are fitted to bounding bodies and used in numerical solution of other PDEs like those of fluid flow, heat flow, and electromagnetics. Such a grid is smooth and has continuous first and second derivatives (and possibly also continuous higher-order derivatives), grid lines are appropriately stretched or clustered, and grid lines are orthogonal or nearly so over most of the grid domain. The source terms in the grid-generating PDEs (hereafter called "defining" PDEs) make it possible for the grid to satisfy requirements for clustering and orthogonality properties in the vicinity of specific surfaces in three dimensions or in the vicinity of specific lines in two dimensions. The grid parameters in question are decay parameters that appear in the source terms of the inhomogeneous defining PDEs. The decay parameters are characteristic lengths in exponential- decay factors that express how the influences of the boundaries decrease with distance from the boundaries. These terms govern the rates at which distance between adjacent grid lines change with distance from nearby boundaries. Heretofore, users have arbitrarily specified decay parameters. However, the characteristic lengths are coupled with the strengths of the source terms, such that arbitrary specification could lead to conflicts among parameter values. Moreover, the manual insertion of decay parameters is cumbersome for static grids and infeasible for dynamically changing grids. In the present method, manual insertion and user specification of decay parameters are neither required nor allowed. Instead, the decay parameters are determined automatically as part of the solution of the defining PDEs. Depending on the shape of the boundary segments and the physical nature of the problem to be solved on the grid, the solution of the defining PDEs may provide for rates of decay to vary along and among the boundary segments and may lend itself to interpretation in terms of one or more physical quantities associated with the problem.

Kaul, Upender K.↗

NASA NTRS2016

Catching Constrains on the Parent Body Genesis of Mesosiderites and a Possible Link to HED (Howardite-Eucrite-Diogenite) Meteorites - A New Hope?

Mesosiderites (MES) are a group of enigmatic stony-iron meteorites exhibiting fragmental matrix breccias and irregular textures; e.g. [1-3]. Mesosiderites contain roughly equal volumes metal (Fe-Ni) and silicates often intimately mixed together (Fig.1). The silicates mostly consist of basaltic, gabbroic, and pyroxenitic components, and appear similar to eucrites and howardites; [4-8]. But unlike HEDs - and other differentiated parent body meteorite groups e.g. ureilites - mesosiderites contain high metal abundances. Several studies have been published to reveal the processes leading to the formation of mesosiderites and attempt to classifiy them [1], [2], [10-15]. Because the silicate inclusions in mesosiderites are often strongly metamorphosed after formation, it is difficult to assess the origin of the silicates and implications for the differentiation process of their parent body [15-17]. Several workers have advanced a formation hypothesis for the mesosiderites where an impact between differentiated bodies occurred prior to 4.47 Ga ago (e.g. [13,18], which could explain the possible incomplete dispersal of the colliding bodies due to their low cosmic ray exposure ages and their special thermal history. However, [13] discuss and favor the model for formation of mesosiderites with the collision of two differentiated bodies, along with disruption events and gravitational re-assembly. The mesosiderites have numerous gabbroid melt clasts with anomalous rare-earth- element (REE) - especially positive Eu - values [19, 20]. HEDs do not show the same. However, the heating mechanisms of both mesosiderites and HED's are puzzling. Mesosiderites are remarkable, they consist of a mix of basalts, which are only found on or near planetary surfaces and undifferentiated metal [1,2]. The probable model is that an asteroid containing a metallic magma impacted onto a second asteroid covered with basalt [18,21]. The mix was then buried under an insulating regolith, and cooled slowly. During cooling and at low temperatures the redox reactions continued to occur and proceed (J.T. Wasson; in pers. comm. 2015).

Eucrite↗

NASA NTRS2018

A Singular Value Decomposition Framework for Retrievals with Vertical Distribution Information from Greenhouse Gas Column Absorption Spectroscopy Measurements

We describe a variation of the Optimal Estimation (OE) method for greenhouse gas remote sensing retrievals using a singular value decomposition (SVD) and an uninformative prior. The SVD method is capable of discerning vertical information in column integrated absorption measurements. While traditional Bayesian optimal estimation (OE) assumes a prior distribution in order to regularize the inversion problem, the SVD approach identifies principal components that can be retrieved from the measurement without explicitly specifying a prior mean and prior covariance matrix. We discuss the method, illustrate its use on an integrated path differential absorption CO2 lidar measurement model, and compare it to traditional optimal estimation using numerical simulations. In the absence of forward model error, the mathematics behind the SVD method guarantee it to be bias-free, which is confirmed by the numerical simulations. In contrast, traditional OE retrievals exhibit bias when the prior mean used in the retrieval differs from the true mean. While the SVD approach can be used for most trace gas retrievals, it is particularly useful for situations where prior knowledge of the trace gas profile is poor. The SVD analysis is also useful in illustrating how vertical information is treated by the traditional OE approach.

Ramanathan, Anand K.↗

NASA NTRS2019

Evaluation of a Multizone Impedance Eduction Method

A computational study is used to evaluate the PyCHE impedance eduction method developed at the NASA Langley Research Center. This method combines an aeroacoustic duct propagation code based on numerical solution to the convected Helmholtz equation with a global optimizer that uses the Differential Evolution algorithm. The efficacy of this method is evaluated with acoustic pressure data simulated to represent that measured with one-zone, two-zone, and three-zone liners mounted in the NASA Langley Grazing Flow Impedance Tube. The PyCHE method has a normalized impedance error of approximately 0.2 for (uniform) one-zone liners with a length of at least 5”, and produces quite reasonable results for liners as short as 2”. Whereas the impedance of the liner has an effect on eduction accuracy, the amount of attenuation is shown to be the dominant parameter. Similar results are observed for two-zone liners, for which the impedance of each zone is unique. The two-zone results also indicate it is more difficult to accurately educe resistance than reactance, and a zone length of at least 6” (slightly longer than for uniform liners) is needed to limit the normalized error to 0.2. The PyCHE method is also demonstrated to successfully educe the impedances for each zone of a three-zone liner. These results are sufficiently encouraging to warrant the continued usage of the PyCHE impedance eduction method for single and multizone liners.

Jones, M. G.↗

NASA NTRSDate not supplied

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗

NASA NTRS1980

An analytical technique for approximating unsteady aerodynamics in the time domain

An analytical technique is presented for approximating unsteady aerodynamic forces in the time domain. The order of elements of a matrix Pade approximation was postulated, and the resulting polynomial coefficients were determined through a combination of least squares estimates for the numerator coefficients and a constrained gradient search for the denominator coefficients which insures stable approximating functions. The number of differential equations required to represent the aerodynamic forces to a given accuracy tends to be smaller than that employed in certain existing techniques where the denominator coefficients are chosen a priori. Results are shown for an aeroelastic, cantilevered, semispan wing which indicate a good fit to the aerodynamic forces for oscillatory motion can be achieved with a matrix Pade approximation having fourth order numerator and second order denominator polynomials.

Dunn, H. J.↗

NASA NTRS1991

Asynchronous and corrected-asynchronous numerical solutions of parabolic PDES on MIMD multiprocessors

A major problem in achieving significant speed-up on parallel machines is the overhead involved with synchronizing the concurrent process. Removing the synchronization constraint has the potential of speeding up the computation. The authors present asynchronous (AS) and corrected-asynchronous (CA) finite difference schemes for the multi-dimensional heat equation. Although the discussion concentrates on the Euler scheme for the solution of the heat equation, it has the potential for being extended to other schemes and other parabolic partial differential equations (PDEs). These schemes are analyzed and implemented on the shared memory multi-user Sequent Balance machine. Numerical results for one and two dimensional problems are presented. It is shown experimentally that the synchronization penalty can be about 50 percent of run time: in most cases, the asynchronous scheme runs twice as fast as the parallel synchronous scheme. In general, the efficiency of the parallel schemes increases with processor load, with the time level, and with the problem dimension. The efficiency of the AS may reach 90 percent and over, but it provides accurate results only for steady-state values. The CA, on the other hand, is less efficient, but provides more accurate results for intermediate (non steady-state) values.

Amitai, Dganit↗

NASA NTRS1984

Differential detection of Gaussian MSK in a mobile radio environment

Minimum shift keying with Gaussian shaped transmit pulses is a strong candidate for a modulation technique that satisfies the stringent out-of-band radiated power requirements of the mobil radio application. Numerous studies and field experiments have been conducted by the Japanese on urban and suburban mobile radio channels with systems employing Gaussian minimum-shift keying (GMSK) transmission and differentially coherent reception. A comprehensive analytical treatment is presented of the performance of such systems emphasizing the important trade-offs among the various system design parameters such as transmit and receiver filter bandwidths and detection threshold level. It is shown that two-bit differential detection of GMSK is capable of offering far superior performance to the more conventional one-bit detection method both in the presence of an additive Gaussian noise background and Rician fading.

Simon, M. K.↗

NASA NTRS2001

End-to-End QoS for Differentiated Services and ATM Internetworking

The Internet was initially design for non real-time data communications and hence does not provide any Quality of Service (QoS). The next generation Internet will be characterized by high speed and QoS guarantee. The aim of this paper is to develop a prioritized early packet discard (PEPD) scheme for ATM switches to provide service differentiation and QoS guarantee to end applications running over next generation Internet. The proposed PEPD scheme differs from previous schemes by taking into account the priority of packets generated from different application. We develop a Markov chain model for the proposed scheme and verify the model with simulation. Numerical results show that the results from the model and computer simulation are in close agreement. Our PEPD scheme provides service differentiation to the end-to-end applications.

Su, Hongjun↗

NASA NTRS1994

Modeling jets in cross flow

Various approaches to the modeling of jets in cross flow are reviewed. These are grouped into four classes, namely: empirical models, integral models, perturbation models, and numerical models. Empirical models depend largely on the correlation of experimental data and are mostly useful for first-order estimates of global properties such as jet trajectory and velocity and temperature decay rates. Integral models are based on some ordinary-differential form of the conservation laws, but require substantial empirical calibration. They allow more details of the flow field to be obtained; simpler versions have to assume similarity of velocity and temperature profiles, but more sophisticated ones can actually calculate these profiles. Perturbation models require little empirical input, but the need for small parameters to ensure convergent expansions limits their application to either the near-field or the far-field. Therefore, they are mostly useful for the study of flow physics. Numerical models are based on conservation laws in partial-differential form. They require little empirical input and have the widest range of applicability. They also require the most computational resources. Although many qualitative and quantitative features of jets in cross flow have been predicted with numerical models, many issues affecting accuracy such as grid resolution and turbulence model are not completely resolved.

Demuren, A. O.↗

NASA NTRS1990

NASA-UCLA Workshop on Computational Techniques in Identification and Control of Flexible Flight Structures, Lake Arrowhead, CA, Nov. 2-4, 1989, Proceedings

The present conference on the identification and control of flexible flight structures with computational techniques encompasses existing and planned testbeds, modeling techniques, adaptive control, and numerical computation and simulation. Specific issues addressed include the NASA/MFSC ground test facility, an active control test on the NASA Minimast, the Astrex testbed for large precision space structures, the dynamic analysis of truss-beam modeling, combined structure-controls-integrated optimization, distributed parameter modeling applicable to the control of flexible flight structures, and approaches to nonclassically damped dynamic systems. Also addressed are adaptive control concepts, techniques for vibration and decentralized slew maneuver control, the discretization of second-order-in-time differential equations which are applicable to nonlinear wave problems, the finite-dimensional approximations of Riccati equations, and numerical simulations of actively controlled space structures.

Balakrishnan, A. V.↗