Search NASA⌕ Search

SEARCH · Search NASA

Results for “Elliptic problems”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Compressible flows with periodic vortical disturbances around lifting airfoils

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method that is presented fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. The numerical scheme is validated through extensive comparisons with known solutions to unsteady vortical flow problems. In general, it is seen that the agreement between the numerical and analytical results is very good for reduced frequencies ranging from 0 to 4, and for Mach numbers ranging from .1 to .8. Numerical results are also presented for a wide variety of flow configurations for the purpose of determining the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. It is seen that each of these parameters can have a significant effect on the unsteady airfoil response to the incident disturbances, and that the effect depends strongly upon the reduced frequency and the dimensionality of the gust. For a one-dimensional (transverse) or two-dimensional (transverse and longitudinal) gust, the results indicate that airfoil thickness increases the unsteady lift and moment at the low reduced frequencies but decreases it at the high reduced frequencies. The results show that an increase in airfoil Mach number leads to a significant increase in the unsteady lift and moment for the low reduced frequencies, but a significant decrease for the high reduced frequencies.

Scott, James R.↗
NASA NTRS1991

A finite-difference, frequency-domain numerical scheme for the solution of the linearized unsteady Euler equations

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method, that is presented, fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. Extensive comparisons are made with known solutions to unsteady vortical flow problems, and it is seen that the agreement is generally very good for reduced frequencies ranging from 0 up to 4.

Scott, James R.↗
NASA NTRS1990

On the three-dimensional instability of strained vortices

The three-dimensional (3-D) instability of a two-dimensional (2-D) flow with elliptical streamlines has been proposed as a generic mechanism for the breakdown of many 2-D flows. A physical interpretation for the mechanism is presented together with an analytical treatment of the problem. It is shown that the stability of an elliptical flow is governed by an Ince equation. An analytical representation for a localized solution is given and establishes a direct link with previous computations and experiments.

Waleffe, Fabian↗
NASA NTRS1987

Development and application of unified algorithms for problems in computational science

A framework is presented for developing computationally unified numerical algorithms for solving nonlinear equations that arise in modeling various problems in mathematical physics. The concept of computational unification is an attempt to encompass efficient solution procedures for computing various nonlinear phenomena that may occur in a given problem. For example, in Computational Fluid Dynamics (CFD), a unified algorithm will be one that allows for solutions to subsonic (elliptic), transonic (mixed elliptic-hyperbolic), and supersonic (hyperbolic) flows for both steady and unsteady problems. The objectives are: development of superior unified algorithms emphasizing accuracy and efficiency aspects; development of codes based on selected algorithms leading to validation; application of mature codes to realistic problems; and extension/application of CFD-based algorithms to problems in other areas of mathematical physics. The ultimate objective is to achieve integration of multidisciplinary technologies to enhance synergism in the design process through computational simulation. Specific unified algorithms for a hierarchy of gas dynamics equations and their applications to two other areas: electromagnetic scattering, and laser-materials interaction accounting for melting.

Shankar, Vijaya↗
NASA NTRS2005

How Does Abundance Affect the Strength of UV Emission in Elliptical Galaxies?

This program used the Far Ultraviolet Spectroscopic Explorer (FUSE) to observe elliptical galaxies with the intention of measuring the chemical abundances in their hot stellar populations. It was designed to complement an earlier FUSE program that observed elliptical galaxies with strong UV emission. The current program originally planned observations of two ellipticals with weak UV emission (M32 and M49). Once FUSE encountered pointing control problems in certain regions of the sky (particularly Virgo, which is very unfortunate for the study of ellipticals in general), M49 was replaced with the bulge of M31, which has a similar UV-to-optical flux ratio as the center of M49. As the closest elliptical galaxy and the one with the weakest UV-to-optical flux ratio, M32 was an obvious choice of target, but M49 was the ideal complementary target, because it has a very low reddening (unlike M32). With the inability of FUSE to point at Virgo, nearly all of the best elliptical galaxies (bright galaxies with low foreground extinction) were also lost, and this severely hampered three FUSE programs of the PI, all focused on the hot stellar populations of ellipticals. M31 was the best replacement for M49, but like M32, it suffers from significant foreground reddening. Strong Galactic ISM lines heavily contaminate the FUSE spectra of M31 and M32. These ISM lines are coincident with the photospheric lines from the stellar populations (whereas M49, with little foreground ISM and significant redshift, would not have suffered from this problem). We have reduced the faint (and thus difficult) data for M31 and M32, producing final co-added spectra representing all of the exposures, but we have not yet finished our analysis, due to the complication of the contaminating ISM. The silver lining here is the set of CHI lines at 1175 Angstroms, which are not significantly contaminated by the ISM. A comparison of the M31 spectrum with other galaxies observed by FEE showed a surprising result: the hot stars in M31 seem to have a similar carbon abundance to those stars in galaxies with much brighter UV emission. The fraction of these hot stars in a population should be a strong function of chemical abundances, so this finding warrants further exploration, and we are proceeding with our analysis. Because the UV emission in these galaxies comes from a population of extreme horizontal branch stars, the PI (Brown) presented this result at a June 2003 conference on such stars.

Sonneborn, George↗
NASA NTRS1986

Survey of the status of finite element methods for partial differential equations

The finite element methods (FEM) have proved to be a powerful technique for the solution of boundary value problems associated with partial differential equations of either elliptic, parabolic, or hyperbolic type. They also have a good potential for utilization on parallel computers particularly in relation to the concept of domain decomposition. This report is intended as an introduction to the FEM for the nonspecialist. It contains a survey which is totally nonexhaustive, and it also contains as an illustration, a report on some new results concerning two specific applications, namely a free boundary fluid-structure interaction problem and the Euler equations for inviscid flows.

Temam, Roger↗
NASA NTRS1991

Calculating satellite umbra/penumbra entry and exit positions and times

The problem of calculating earth satellite entry/exit positions and times through the earth's umbra and penumbra for elliptical satellite orbits is solved in closed form; i.e., without iteration, by reducing the problem to that of finding the roots of a quartic polynomial in a cartesian coordinate. Typical results from an algorithm constructed from the method are shown.

Mullins, Larry D.↗
NASA NTRS1987

Problem size, parallel architecture and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗
NASA NTRS1988

Problem size, parallel architecture, and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗
NASA NTRS1983

Preconditioned Minimal Residual Methods for Chebyshev Spectral Caluclations

The problem of preconditioning the pseudospectral Chebyshev approximation of an elliptic operator is considered. The numerical sensitiveness to variations of the coefficients of the operator are investigated for two classes of preconditioning matrices: one arising from finite differences, the other from finite elements. The preconditioned system is solved by a conjugate gradient type method, and by a DuFort-Frankel method with dynamical parameters. The methods are compared on some test problems with the Richardson method and with the minimal residual Richardson method.

Canuto, C.↗
NASA NTRS1985

Preconditioned minimal residual methods for Chebyshev spectral calculations

The problem of preconditioning the pseudospectral Chebyshev approximation of an elliptic operator is considered. The numerical sensitiveness to variations of the coefficients of the operator are investigated for two classes of preconditioning matrices: one arising from finite differences, the other from finite elements. The preconditioned system is solved by a conjugate gradient type method, and by a Dufort-Frankel method with dynamical parameters. The methods are compared on some test problems with the Richardson method and with the minimal residual Richardson method.

Canuto, C.↗
NASA NTRS1982

A minimal residual method for transonic potential flows

For transonic flow calculations, a combination of the successive line over-relaxation (SLOR) and the preconditioned conjugate gradient (CG) method has been suggested by Wong and Hafez (1981). This paper studies the method of minimal residual (MR) which avoids a combined iteration. This method is closely related to the CG method, may be regarded as a first-order gradient method, and is applicable to symmetric and nonsymmetric matrices. The problem is formulated as a nonlinear mixed elliptic-hyperbolic partial differential equation which includes an artificial viscosity and a switching function which is zero in subsonic regions and nonzero in supersonic regions. Alternatives to the SLOR method which provide faster convergence rates are introduced. The preconditioned MR algorithm is developed, and transonic potential flows around NACA 0012 airfoil are calculated for different Mach numbers and angles of attack. Preliminary results are presented, demonstrating that the MR algorithm requires no parameter estimation and rapidly converges for subsonic flows.

Wong, Y. S.↗
NASA NTRS1987

A new method for the aerodynamic analysis of lifting surfaces

The use of the original Ffowcs Williams and Hawkings equation to improve the solution to simple problems of two-dimensional and incompressible flow is studied. It is shown that a linear formulation, with the quadrupole term neglected, leads to better correlation with the results of potential theory. A solution to steady and two-dimensional problems is developed and applied to families of elliptic cylinders and symmetric airfoils.

Brandao, Mauricio Pazini↗
NASA NTRS1992

Towards a theory of automated elliptic mesh generation

The theory of elliptic mesh generation is reviewed and the fundamental problem of constructing computational space is discussed. It is argued that the construction of computational space is an NP-Complete problem and therefore requires a nonstandard approach for its solution. This leads to the development of graph-theoretic, combinatorial optimization and integer programming algorithms. Methods for the construction of two dimensional computational space are presented.

Cordova, J. Q.↗
NASA NTRS1994

An experiment in hurricane track prediction using parallel computing methods

The barotropic model is used to explore the advantages of parallel processing in deterministic forecasting. We apply this model to the track forecasting of hurricane Elena (1985). In this particular application, solutions to systems of elliptic equations are the essence of the computational mechanics. One set of equations is associated with the decomposition of the wind into irrotational and nondivergent components - this determines the initial nondivergent state. Another set is associated with recovery of the streamfunction from the forecasted vorticity. We demonstrate that direct parallel methods based on accelerated block cyclic reduction (BCR) significantly reduce the computational time required to solve the elliptic equations germane to this decomposition and forecast problem. A 72-h track prediction was made using incremental time steps of 16 min on a network of 3000 grid points nominally separated by 100 km. The prediction took 30 sec on the 8-processor Alliant FX/8 computer. This was a speed-up of 3.7 when compared to the one-processor version. The 72-h prediction of Elena's track was made as the storm moved toward Florida's west coast. Approximately 200 km west of Tampa Bay, Elena executed a dramatic recurvature that ultimately changed its course toward the northwest. Although the barotropic track forecast was unable to capture the hurricane's tight cycloidal looping maneuver, the subsequent northwesterly movement was accurately forecasted as was the location and timing of landfall near Mobile Bay.

Song, Chang G.↗
NASA NTRS1996

Preconditioned Mixed Spectral Element Methods for Elasticity and Stokes Problems

Preconditioned iterative methods for the indefinite systems obtained by discretizing the linear elasticity and Stokes problems with mixed spectral elements in three dimensions are introduced and analyzed. The resulting stiffness matrices have the structure of saddle point problems with a penalty term, which is associated with the Poisson ratio for elasticity problems or with stabilization techniques for Stokes problems. The main results of this paper show that the convergence rate of the resulting algorithms is independent of the penalty parameter, the number of spectral elements Nu and mildly dependent on the spectral degree eta via the inf-sup constant. The preconditioners proposed for the whole indefinite system are block-diagonal and block-triangular. Numerical experiments presented in the final section show that these algorithms are a practical and efficient strategy for the iterative solution of the indefinite problems arising from mixed spectral element discretizations of elliptic systems.

Pavarino, Luca F.↗