SEARCH · Search NASA

Results for “collocation”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Using spatial coherence to retrieve cloud properties - ERBE experience and CERES applications

Cloud properties for individual ERBE scanner fields of view are obtained from simultaneous, collocated imagery data. The cloud properties are used to assess the performance of the ERBE maximum likelihood estimate of cloud conditions which ultimately influences estimates of the radiative fluxes at the top of the atmosphere. During EOS the analysis of MODIS observations within the fields of view of the CERES instruments promises to yield, in addition, estimates of surface radiative fluxes as well as the vertical distribution of radiative heating within the atmosphere.

Coakley, James A., Jr.↗

NASA NTRS1990

Quadrature imposition of compatibility conditions in Chebyshev methods

This paper handles a problem encountered in the solution by way of Chebyshev spectral collocation discretization of the Poisson equation with homogeneous Neumann boundary conditions. When pseudospectral Chebyshev methodology is utilized to discretize the partial differential equation these integrals have to be approximated by an appropriate quadrature formula. The correct quadrature formula is developed for these problems, taking 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 NTRS1990

Improved operation of magnetic bearings for flywheel energy storage system

Analysis and operation of prototype 500-Wh flywheel at low speeds have shown that many factors affect the correct functioning of the magnetic bearings. An examination is made of a number of these, including magnetic bearing control system nonlinearities and displacement transducer positioning, and their effects upon the successful operation of the suspension system. It is observed that the bearing control system is extremely sensitive to actuator parameters such as coil inductance. As a consequence of the analysis of bearing relaxation oscillations, the bearing actuator design methodology which has previously been used, where coil parameter selection is based upon static considerations, has been revised. Displacement transducer sensors which overcome the collocation problem are discussed.

Zmood, R. B.↗

NASA NTRS1991

On the viscous modes of instability of a trailing line vortex

A viscous linear stability analysis of a trailing line (Batchelor) vortex is presented. Employing a staggered Chebyshev spectral collocation technique, very accurate results were obtained. The destabilizing role of viscous forces has been shown to produce two types of viscous instability modes. These viscous disturbances consist of an axisymmetric mode and an asymmetric mode. Both disturbances are long-wave instabilities with maximum growth rates which are orders of magnitude smaller than the inviscid modes which have been found by others. Comparison with experimental results and condensation trail observations are found to be in good qualitative agreement with the present study.

Khorrami, Mehdi R.↗

NASA NTRS1990

Second mode interactions in supersonic boundary layers

The nonlinear evolution of a two-dimensional second mode unstable wave in a Mach 4.5 wall-bounded flow is computed by solving the full time-dependent compressible Navier-Stokes equations. A highly accurate solution is obtained using spectral collocation methods. It is shown that departure from linearity first occurs in the critical layer due to the cubic nonlinearities in the momentum equation. This is a direct result of the large density perturbations in this regime. Time evolution studies of the growth rate as a function of normal distance from the plate suggests that the mode is evolving toward a nonlinear saturated state, and that this problem is possibly amenable to standard weakly nonlinear perturbation methods.

Erlebacher, Gordon↗

NASA NTRS1990

Enstrophy budget in decaying compressible turbulence

The interaction of vorticity and shock waves in decaying compressible turbulence is investigated by direct numerical simulation of the Navier-Stokes equation using a spectral method with 64-cubed collocation points. The Taylor microscale Reynolds number decreases from 36 to 17 and the rms Mach number from 0.8 to 0.4 during the decay. It is observed that vorticity is created inside curved shocks through baroclinic interaction and intensified by the compression of fluid elements. Vortex stretching is, on the other hand, dominant outside shocks and is responsible for enhancing the strength of vorticity. A dipole of vorticity is typically created at the vertex of curved shocks, the structure of which is explained in terms of jump relations of velocity derivatives across a curved shock.

Kida, Shigeo↗

NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗

NASA NTRS1991

Spectral methods for the Euler equations - The blunt body problem revisited

The present use of the Chebyshev spectral collocation method, in conjunction with shock-fitting, to solve the blunt-body problem gives attention to the boundary and the shock-acceleration equations. The crux of these procedures is the use of the characteristic compatibility relations to compute the body pressure and shock velocity. It is shown that converged solutions are obtainable without artificial smoothing, and that spectral accuracy is achieved.

Kopriva, David A.↗

NASA NTRS1991

Application of the Global Positioning System to crustal deformation measurement. II - The influence of errors in orbit determination networks

The errors introduced into GPS measurements of a geodetic network by means of adopting different sets of stations in central and southern California as fixed points (termed fiducial stations) were investigated. Solutions for a crustal deformation network spanning distances up to 350 km were computed with a variety of fiducial networks, using fiducial coordinates determined from VLBI. Results showed that three receivers collocated with VLBI sites provide a stable reference system at a few millimeter level in the horizontal components and at 10 mm in the vertical components, indicating that local survey errors between VLBI and GPS markers at fiducial sites do not severely affect the estimates of interstation vectors.

Larson, Kristine M.↗

NASA NTRS1991

Multidomain spectral solution of the Euler gas-dynamics equations

The present interfacial treatments for Euler gasdynamic equation computations via the multidomain Chebyshev spectral collocation method are applicable both at subdomain corners and in overlapping or patched subdomains, for interfaces located in the sub-, super-, or transonic regions of a flow. The results thus obtained are found to be (1) spectrally accurate, (2) both more accurate and more efficient than a single-domain calculation, and (3) potentially more accurate and efficient than a single-domain calculation. Interfacial wave reflection is insignificant.

Kopriva, David A.↗

NASA NTRS1991

Estimation of geopotential differences over intercontinental locations using satellite and terrestrial measurements

An error analysis study was conducted in order to assess the current accuracies and the future anticipated improvements in the estimation of geopotential differences over intercontinental locations. An observation/estimation scheme was proposed and studied, whereby gravity disturbance measurements on the Earth's surface, in caps surrounding the estimation points, are combined with corresponding data in caps directly over these points at the altitude of a low orbiting satellite, for the estimation of the geopotential difference between the terrestrial stations. The mathematical modeling required to relate the primary observables to the parameters to be estimated, was studied for the terrestrial data and the data at altitude. Emphasis was placed on the examination of systematic effects and on the corresponding reductions that need to be applied to the measurements to avoid systematic errors. The error estimation for the geopotential differences was performed using both truncation theory and least squares collocation with ring averages, in case observations on the Earth's surface only are used. The error analysis indicated that with the currently available global geopotential model OSU89B and with gravity disturbance data in 2 deg caps surrounding the estimation points, the error of the geopotential difference arising from errors in the reference model and the cap data is about 23 kgal cm, for 30 deg station separation.

Pavlis, Nikolaos K.↗

NASA NTRS1991

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗

NASA NTRS1991

Development of a CFD Code for Analysis of Fluid Dynamic Forces in Seals

The aim is to develop a 3-D computational fluid dynamics (CFD) code for the analysis of fluid flow in cylindrical seals and evaluation of the dynamic forces on the seals. This code is expected to serve as a scientific tool for detailed flow analysis as well as a check for the accuracy of the 2D industrial codes. The features necessary in the CFD code are outlined. The initial focus was to develop or modify and implement new techniques and physical models. These include collocated grid formulation, rotating coordinate frames and moving grid formulation. Other advanced numerical techniques include higher order spatial and temporal differencing and an efficient linear equation solver. These techniques were implemented in a 2D flow solver for initial testing. Several benchmark test cases were computed using the 2D code, and the results of these were compared to analytical solutions or experimental data to check the accuracy. Tests presented here include planar wedge flow, flow due to an enclosed rotor, and flow in a 2D seal with a whirling rotor. Comparisons between numerical and experimental results for an annular seal and a 7-cavity labyrinth seal are also included.

Athavale, Mahesh M.↗

NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗

NASA NTRS1992

Vibration attenuation of the NASA Langley evolutionary structure experiment using H(sub infinity) and structured singular value (micron) robust multivariable control techniques

The use is studied of active control to attenuate structural vibrations of the NASA Langley Phase Zero Evolutionary Structure due to external disturbance excitations. H sub infinity and structured singular value (mu) based control techniques are used to analyze and synthesize control laws for the NASA Langley Controls Structures Interaction (CSI) Evolutionary Model (CEM). The CEM structure experiment provides an excellent test bed to address control design issues for large space structures. Specifically, control design for structures with numerous lightly damped, coupled flexible modes, collocated and noncollocated sensors and actuators and stringent performance specifications. The performance objectives are to attenuate the vibration of the structure due to external disturbances, and minimize the actuator control force. The control design problem formulation for the CEM Structure uses a mathematical model developed with finite element techniques. A reduced order state space model for the control design is formulated from the finite element model. It is noted that there are significant variations between the design model and the experimentally derived transfer function data.

Balas, Gary J.↗

NASA NTRS1992

Stability and performance of notch filter control for unbalance response

Many current applications of magnetic bearings for rotating machinery employ notch filters in the feedback control loop to reduce the synchronous forces transmitted through the bearings. The capabilities and limitations of notch filter control are investigated. First, a rigid rotor is examined with some classical root locus techniques. Notch filter control is shown to result in conditional stability whenever complete synchronous attenuation is required. Next, a nondimensional parametric symmetric flexible three mass rotor model is constructed. An examination of this model for several test cases illustrates the limited attenuation possible with notch filters at and near the system critical speeds when the bearing damping is low. The notch filter's alteration of the feedback loop is shown to cause stability problems which limits performance. Poor transient response may also result. A high speed compressor is then examined as a candidate for notch filter control. A collocated 22 mass station model with lead-lag control is used. The analysis confirms the reduction in stability robustness that can occur with notch filter control. It is concluded that other methods of synchronous vibration control yield greater performance without compromising stability.

Knospe, C. R.↗

NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

NASA NTRS1990

Shear buckling of specially orthotropic plates with centrally located cutouts

There is significant industry demand for a method of analyzing the shear buckling of composite plates with cutouts. This method should be able to easily accommodate frequent changes in model design; inflexibility is the major drawback to current finite element methods. The approach taken here is broken into two problems, prebuckling and buckling. To solve for the prebuckling stresses, complex variable equations are used in conjunction with boundary collocation. The least squares approach is utilized to improve the accuracy of the results. The buckling problem is solved using the Ritz method. A product of the Ritz method is a complicated integral equation which is solved using numerical integration. To date, the aforementioned method of determining the prebuckling stresses was verified against infinite plate theory and finite elements.

Owen, Vicki L.↗