Search NASA⌕ Search

SEARCH · Search NASA

Results for “DISCRETE FUNCTION”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method↗
NASA NTRS1990

Quadrature imposition of compatibility conditions in Chebyshev methods

Often, in solving an elliptic equation with Neumann boundary conditions, a compatibility condition has to be imposed for well-posedness. This condition involves integrals of the forcing function. When pseudospectral Chebyshev methods are used to discretize the partial differential equation, these integrals have to be approximated by an appropriate quadrature formula. The Gauss-Chebyshev (or any variant of it, like the Gauss-Lobatto) formula can not be used here since the integrals under consideration do not include the weight function. A natural candidate to be used in approximating the integrals is the Clenshaw-Curtis formula, however it is shown that this is the wrong choice and it may lead to divergence if time dependent methods are used to march the solution to steady state. The correct quadrature formula is developed for these problems. This formula takes into account the degree of the polynomials involved. It is shown that this formula leads to a well conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.

Gottlieb, D.↗
NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗
NASA NTRS1971

Computational requirements for a discrete Kalman filter.

Computational requirements - i.e., computing time per cycle (iteration) and required storage - which determine minimum sampling rates and computer memory size, were obtained as functions of the dimensions of the important system matrices for a discrete Kalman filter. Two types of measurement processing are discussed: simultaneous and sequential. It is shown that it is often better to process statistically independent measurements in more than one batch and then use sequential processing than to process them together via simultaneous processing.

Mendel, J. M.↗
NASA NTRS1985

Improving the Accuracy of the Boundary Integral Method Based on the Helmholtz Integral

Several recent papers in the literature have been based on various forms of the Helmholtz integral to compute the radiation fields of vibrating bodies. The surface integral form is given. The symbols of P,R micron, rho,G,R,V, and S micron are acoustic pressure, source coordinate, angular frequency, fluid density, Green function, field coordinate, surface velocity and body surface respectively. A discretized form of the surface integral is also given. Solutions to the surface integral are complicated with the singularity of the Green function at R=R micron and with the uniqueness problem at interior eigen frequencies of the enclosed space. The use of the interior integral circumvents the singularity problem since the field points are chosen in the interior space of the vibrating body where a zero pressure condition exists. The interior integral form is given. The method to improve the accuracy is detailed. Examples of the method is presented for a variety of radiators.

Koopmann, G. H.↗
NASA NTRS1989

Extension-torsion coupling behavior of advanced composite tilt-rotor blades

An analytic model was developed to study the extension-bend-twist coupling behavior of an advanced composite helicopter or tilt-rotor blade. The outer surface of the blade is defined by rotating an arbitrary cross section about an initial twist axis. The cross section can be nonhomogeneous and composed of generally anisotropic materials. The model is developed based upon a three dimensional elasticity approach that is recast as a coupled two-dimensional boundary value problem defined in a curvilinear coordinate system. Displacement solutions are written in terms of known functions that represent extension, bending, and twisting and unknown functions for local cross section deformations. The unknown local deformation functions are determined by applying the principle of minimum potential energy to the discretized two-dimensional cross section. This is an application of the Ritz method, where the trial function family is the displacement field associated with a finite element (8-node isoparametric quadrilaterals) representation of the section. A computer program was written where the cross section is discretized into 8-node quadrilateral subregions. Initially the program was verified using previously published results (both three-dimensional elasticity and technical beam theory) for pretwisted isotropic bars with an elliptical cross section. In addition, solid and thin-wall multi-cell NACA-0012 airfoil sections were analyzed to illustrate the pronounced effects that pretwist, initial twist axis location, and spar location has on coupled behavior. Currently, a series of advanced composite airfoils are being modeled in order to assess how the use of laminated composite materials interacts with pretwist to alter the coupling behavior of the blade. These studies will investigate the use of different ply angle orientations and the use of symmetric versus unsymmetric laminates.

Kosmatka, J. B.↗
NASA NTRS1972

An innovative approach to compensator design

The primary goal is to present for a control system a computer-aided-compensator design technique from a frequency domain point of view. The thesis for developing this technique is to describe the open loop frequency response by n discrete frequency points which result in n functions of the compensator coefficients. Several of these functions are chosen so that the system specifications are properly portrayed; then mathematical programming is used to improve all of these functions which have values below minimum standards. In order to do this several definitions in regard to measuring the performance of a system in the frequency domain are given. Next, theorems which govern the number of compensator coefficients necessary to make improvements in a certain number of functions are proved. After this a mathematical programming tool for aiding in the solution of the problem is developed. Then for applying the constraint improvement algorithm generalized gradients for the constraints are derived. Finally, the necessary theory is incorporated in a computer program called CIP (compensator improvement program).

Mitchell, J. R.↗
NASA NTRS1973

An innovative approach to compensator design

The design is considered of a computer-aided-compensator for a control system from a frequency domain point of view. The design technique developed is based on describing the open loop frequency response by n discrete frequency points which result in n functions of the compensator coefficients. Several of these functions are chosen so that the system specifications are properly portrayed; then mathematical programming is used to improve all of these functions which have values below minimum standards. To do this, several definitions in regard to measuring the performance of a system in the frequency domain are given, e.g., relative stability, relative attenuation, proper phasing, etc. Next, theorems which govern the number of compensator coefficients necessary to make improvements in a certain number of functions are proved. After this a mathematical programming tool for aiding in the solution of the problem is developed. This tool is called the constraint improvement algorithm. Then for applying the constraint improvement algorithm generalized, gradients for the constraints are derived. Finally, the necessary theory is incorporated in a Computer program called CIP (compensator Improvement Program). The practical usefulness of CIP is demonstrated by two large system examples.

Mitchell, J. R.↗
NASA NTRS1974

Computer programs for pressurization (RAMP) and pressurized expulsion from a cryogenic liquid propellant tank

An analysis to predict the pressurant gas requirements for the discharge of cryogenic liquid propellants from storage tanks is presented, along with an algorithm and two computer programs. One program deals with the pressurization (ramp) phase of bringing the propellant tank up to its operating pressure. The method of analysis involves a numerical solution of the temperature and velocity functions for the tank ullage at a discrete set of points in time and space. The input requirements of the program are the initial ullage conditions, the initial temperature and pressure of the pressurant gas, and the time for the expulsion or the ramp. Computations are performed which determine the heat transfer between the ullage gas and the tank wall. Heat transfer to the liquid interface and to the hardware components may be included in the analysis. The program output includes predictions of mass of pressurant required, total energy transfer, and wall and ullage temperatures. The analysis, the algorithm, a complete description of input and output, and the FORTRAN 4 program listings are presented. Sample cases are included to illustrate use of the programs.

Masters, P. A.↗
NASA NTRS1975

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗
NASA NTRS1976

Computer measurement and representation of the heart in two and three dimensions

Methods for the measurement and display by minicomputer of cardiac images obtained from fluoroscopy to permit an accurate assessment of functional changes are discussed. Heart contours and discrete points can be digitized automatically or manually, with the recorded image in a video, cine, or print format. As each frame is digitized it is assigned a code name identifying the data source, experiment, run, view, and frame, and the images are filed for future reference in any sequence. Two views taken at the same point in the heart cycle are used to compute the spatial position of the ventricle apex and the midpoint of the aortic valve. The remainder of the points on the chamber border are corrected for the linear distortion of the X-rays by projection to a plane containing the chord between the apex and the aortic valve center and oriented so that lines perpendicular to the chord are parallel to the image intensifier face. The image of the chamber surface is obtained by generating circular cross sections with diameters perpendicular to the major chord. The transformed two- and three-dimensional imagery can be displayed in either static or animated form using a graphics terminal.

Rasmussen, D.↗
NASA NTRS1978

Cockpit displays and the growing role of the pilot in the ATC system

It is pointed out that advances in avionics and computer technology have reached the point where serious attention may be given to structuring a new controlled air traffic environment. The principal feature of such an environment is related to the capability of pilots to perform certain active ATC functions from the cockpit. The main function is concerned with the ability to maintain safe separation from other aircraft under all meteorological conditions. Technical features available to permit such a functional change include such concepts as the Discrete Address Beacon System, a full capability Beacon Collision Avoidance System, data links, computer generated displays, and flight management computers.

Lyman, E. G.↗
NASA NTRS1980

Aliased noise in radiometric measurements

The magnitude of aliased noise that degrades the accuracy of continuous reconstructions of discrete radiometric measurements was evaluated as a function of the spatial response and sampling intervals of the radiometer, and of the resolution of the reconstructed measurements. A Wiener spectrum, representative of a wide range of scenes, was used to characterize the radiance fluctuations.

Huck, F. O.↗
NASA NTRS1988

Spectral synthesis in the ultraviolet. II - Stellar populations and star formation in blue compact galaxies

An initial attempt to apply optimizing spectral synthesis techniques to the far-UV spectra of blue compact galaxies (BCGs) is presented. The far-UV absorption-line spectra of the galaxies are clearly composite, with the signatures of the main-sequence types between O3 and mid-A. Most of the low-ionization absorption lines have a stellar origin. The Si IV and C IV features in several objects have P Cygni profiles. In Haro I the strength of Si IV indicates a significant blue supergiant population. The metal-poor blue compact dwarf Mrk 209 displays weak absorption lines, evidence that the stellar component has the same low metallicity as observed in the ionized gas. Good fits to the data are obtained the technique of optimizing population synthesis. The solutions yield stellar luminosity functions which display large discontinuities, indicative of discrete star formation episodes or bursts. The amount of UV extinction is low.

Fanelli, Michael N.↗
NASA NTRS1989

Nonlinear finite element formulation for the large displacement analysis in multibody system dynamics

A total Lagrangian finite element formulation for the deformable bodies in multibody mechanical systems that undergo finite relative rotations is developed. The deformable bodies are discretized using finite element methods. The shape functions that are used to describe the displacement field are required to include the rigid body modes that describe only large translational displacements. This does not impose any limitations on the technique because most commonly used shape functions satisfy this requirement. The configuration of an element is defined using four sets of coordinate systems: Body, Element, Intermediate element, Global. The body coordinate system serves as a unique standard for the assembly of the elements forming the deformable body. The element coordinate system is rigidly attached to the element and therefore it translates and rotates with the element. The intermediate element coordinate system, whose axes are initially parallel to the element axes, has an origin which is rigidly attached to the origin of the body coordinate system and is used to conveniently describe the configuration of the element in undeformed state with respect to the body coordinate system.

Rismantab-Sany, J.↗
NASA NTRS1992

Discrete shaped strain sensors for intelligent structures

Design of discrete, highly distributed sensor systems for intelligent structures has been studied. Data obtained indicate that discrete strain-averaging sensors satisfy the functional requirements for distributed sensing of intelligent structures. Bartlett and Gauss-Hanning sensors, in particular, provide good wavenumber characteristics while meeting the functional requirements. They are characterized by good rolloff rates and positive Fourier transforms for all wavenumbers. For the numerical integration schemes, Simpson's rule is considered to be very simple to implement and consistently provides accurate results for five sensors or more. It is shown that a sensor system that satisfies the functional requirements can be applied to a structure that supports mode shapes with purely sinusoidal curvature.

Andersson, Mark S.↗
NASA NTRS1996

Heat and Momentum Transfer Studies in High Reynolds Number Wavy Films at Normal and Reduced Gravity Conditions

We examined the effect of the gas flow on the liquid film when the gas flows in the countercurrent direction in a vertical pipe at normal gravity conditions. The most dramatic effect of the simultaneous flow of gas and liquid in pipes is the greatly increased transport rates of heat, mass, and momentum. In practical situations this enhancement can be a benefit or it can result in serious operational problems. For example, gas-liquid flow always results in substantially higher pressure drop and this is usually undesirable. However, much higher heat transfer coefficients can be expected and this can obviously be of benefit for purposes of design. Unfortunately, designers know so little of the behavior of such two phase systems and as a result these advantages are not utilized. Due to the complexity of the second order boundary model as well as the fact that the pressure variation across the film is small compared to the imposed gas phase pressure, the countercurrent gas flow affect was studied for the standard boundary layer model. A different stream function that can compensate the shear stress affect was developed and this stream function also can predict periodic solutions. The discretized model equations were transformed to a traveling wave coordinate system. A stability analysis of these sets of equations showed the presence of a Hopf bifurcation for certain values of the traveling wave velocity and the shear stress. The Hopf celerity was increased due to the countercurrent shear. For low flow rate the increases of celerity are more than for the high flow rate, which was also observed in experiments. Numerical integration of a traveling wave simplification of the model also predicts the existence of chaotic large amplitude, nonperiodic waves as observed in the experiments. The film thickness was increased by the shear.

Balakotaiah, V.↗
NASA NTRS1997

LQR Control of Thin Shell Dynamics: Formulation and Numerical Implementation

A PDE-based feedback control method for thin cylindrical shells with surface-mounted piezoceramic actuators is presented. Donnell-Mushtari equations modified to incorporate both passive and active piezoceramic patch contributions are used to model the system dynamics. The well-posedness of this model and the associated LQR problem with an unbounded input operator are established through analytic semigroup theory. The model is discretized using a Galerkin expansion with basis functions constructed from Fourier polynomials tensored with cubic splines, and convergence criteria for the associated approximate LQR problem are established. The effectiveness of the method for attenuating the coupled longitudinal, circumferential and transverse shell displacements is illustrated through a set of numerical examples.

delRosario, R. C. H.↗