Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELLIPTIC EQUATION”

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 235 records · Page 13
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 NTRS1989

Some fast elliptic solvers on parallel architectures and their complexities

The discretization of separable elliptic partial differential equations leads to linear systems with special block triangular matrices. Several methods are known to solve these systems, the most general of which is the Block Cyclic Reduction (BCR) algorithm which handles equations with nonconsistant coefficients. A method was recently proposed to parallelize and vectorize BCR. Here, the mapping of BCR on distributed memory architectures is discussed, and its complexity is compared with that of other approaches, including the Alternating-Direction method. A fast parallel solver is also described, based on an explicit formula for the solution, which has parallel computational complexity lower than that of parallel BCR.

Gallopoulos, E.↗
NASA NTRS1991

Effect of nose shape on three-dimensional stagnation region streamlines and heating rates

A new method for calculating the three-dimensional inviscid surface streamlines and streamline metrics using Cartesian coordinates and time as the independent variable of integration has been developed. The technique calculates the streamline from a specified point on the body to a point near the stagnation point by using a prescribed pressure distribution in the Euler equations. The differential equations, which are singular at the stagnation point, are of the two point boundary value problem type. Laminar heating rates are calculated using the axisymmetric analog concept for three-dimensional boundary layers and approximate solutions to the axisymmetric boundary layer equations. Results for elliptic conic forebody geometries show that location of the point of maximum heating depends on the type of conic in the plane of symmetry and the angle of attack, and that this location is in general different from the stagnation point. The new method was found to give smooth predictions of heat transfer in the nose region where previous methods gave oscillatory results.

Hassan, Basil↗
NASA NTRS1992

A regularization method for extrapolation of solar potential magnetic fields

The mathematical basis of a Tikhonov regularization method for extrapolating the chromospheric-coronal magnetic field using photospheric vector magnetograms is discussed. The basic techniques show that the Cauchy initial value problem can be formulated for potential magnetic fields. The potential field analysis considers a set of linear, elliptic partial differential equations. It is found that, by introducing an appropriate smoothing of the initial data of the Cauchy potential problem, an approximate Fourier integral solution is found, and an upper bound to the error in the solution is derived. This specific regularization technique, which is a function of magnetograph measurement sensitivities, provides a method to extrapolate the potential magnetic field above an active region into the chromosphere and low corona.

Gary, G. A.↗
NASA NTRS1992

Applications of multigrid software in the atmospheric sciences

Elliptic partial differential equations from different areas in the atmospheric sciences are efficiently and easily solved utilizing the multigrid software package named MUDPACK. It is demonstrated that the multigrid method is more efficient than other commonly employed techniques, such as Gaussian elimination and fixed-grid relaxation. The efficiency relative to other techniques, both in terms of storage requirement and computational time, increases quickly with grid size.

Adams, J.↗
NASA NTRS1951

Flow through cascades in tandem

An exact treatment of the problem of finding the incompressible, inviscid two-dimensional flow around two cascades in tandem is presented. The analysis includes solutions of both the direct and the inverse problems. These problems are solved by conformally mapping the tandem cascade onto the region between two concentric circles in which region there are suitably placed flow singularities. Formulas for the velocity and the potential in the annular region are presented in a closed form by means of elliptic functions. The equations are presented in a form suitable for computation.

Spraglin, William E↗
NASA NTRS1993

MGGHAT: Elliptic PDE software with adaptive refinement, multigrid and high order finite elements

MGGHAT (MultiGrid Galerkin Hierarchical Adaptive Triangles) is a program for the solution of linear second order elliptic partial differential equations in two dimensional polygonal domains. This program is now available for public use. It is a finite element method with linear, quadratic or cubic elements over triangles. The adaptive refinement via newest vertex bisection and the multigrid iteration are both based on a hierarchical basis formulation. Visualization is available at run time through an X Window display, and a posteriori through output files that can be used as GNUPLOT input. In this paper, we describe the methods used by MGGHAT, define the problem domain for which it is appropriate, illustrate use of the program, show numerical and graphical examples, and explain how to obtain the software.

Mitchell, William F.↗
NASA NTRS1996

Preconditioning Operators on Unstructured Grids

We consider systems of mesh equations that approximate elliptic boundary value problems on arbitrary (unstructured) quasi-uniform triangulations and propose a method for constructing optimal preconditioning operators. The method is based upon two approaches: (1) the fictitious space method, i.e., the reduction of the original problem to a problem in an auxiliary (fictitious) space, and (2) the multilevel decomposition method, i.e., the construction of preconditioners by decomposing functions on hierarchical meshes. The convergence rate of the corresponding iterative process with the preconditioner obtained is independent of the mesh step. The preconditioner has an optimal computational cost: the number of arithmetic operations required for its implementation is proportional to the number of unknowns in the problem. The construction of the preconditioning operators for three dimensional problems can be done in the same way.

Nepomnyaschikh, S. V.↗
NASA NTRS1996

MGLab: An Interactive Multigrid Environment

MGLab is a set of Matlab functions that defines an interactive environment for experimenting with multigrid algorithms. The package solves two-dimensional elliptic partial differential equations discretized using either finite differences or finite volumes, depending on the problem. Built-in problems include the Poisson equation, the Helmholtz equation, a convection-diffusion problem, and a discontinuous coefficient problem. A number of parameters controlling the multigrid V-cycle can be set using a point-and-click mechanism. The menu-based user interface also allows a choice of several Krylov subspace methods, including CG, GMRES(k), and Bi-CGSTAB, which can be used either as stand-alone solvers or as multigrid acceleration schemes. The package exploits Matlab's visualization and sparse matrix features and has been structured to be easily extensible.

Bordner, James↗
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 NTRS2011

Extracting Zero-Gravity Surface Figure of a Mirror

The technical innovation involves refinement of the classic optical technique of averaging surface measurements made in different orientations with respect to gravity, so the effects of gravity cancel in the averaged image. Particularly for large, thin mirrors subject to substantial deformation, the further requirement is that mount forces must also cancel when averaged over measurement orientations. The zerogravity surface figure of a mirror in a hexapod mount is obtained by analyzing the summation of mount forces in the frame of the optic as surface metrology is averaged over multiple clockings. This is illustrated with measurements taken from the Space Interferometry Mission (SIM) PT-Ml mirror for both twofold and threefold clocking. The positive results of these measurements and analyses indicate that, from this perspective, a lighter mirror could be used; that is, one might place less reliance on the damping effects of the elliptic partial differential equations that describe the propagation of forces through glass. The advantage over prior art is relaxing the need for an otherwise substantial thickness of glass that might be needed to ensure accurate metrology in the absence of a detailed understanding and analysis of the mount forces. The general insights developed here are new, and provide the basic design principles on which mirror mount geometry may be chosen.

Bloemhof, Eric E.↗
NASA NTRS2016

Cluster Observations of Non-Time Continuous Magnetosonic Waves

Equatorial magnetosonic waves are normally observed as temporally continuous sets of emissions lasting from minutes to hours. Recent observations, however, have shown that this is not always the case. Using Cluster data, this study identifies two distinct forms of these non temporally continuous use missions. The first, referred to as rising tone emissions, are characterized by the systematic onset of wave activity at increasing proton gyroharmonic frequencies. Sets of harmonic emissions (emission elements)are observed to occur periodically in the region +/- 10 off the geomagnetic equator. The sweep rate of these emissions maximizes at the geomagnetic equator. In addition, the ellipticity and propagation direction also change systematically as Cluster crosses the geomagnetic equator. It is shown that the observed frequency sweep rate is unlikely to result from the sideband instability related to nonlinear trapping of suprathermal protons in the wave field. The second form of emissions is characterized by the simultaneous onset of activity across a range of harmonic frequencies. These waves are observed at irregular intervals. Their occurrence correlates with changes in the spacecraft potential, a measurement that is used as a proxy for electron density. Thus, these waves appear to be trapped within regions of localized enhancement of the electron density.

Walker, Simon N.↗
NASA NTRS2009

Formation Flying Control Implementation in Highly Elliptical Orbits

The Tschauner-Hempel equations are widely used to correct the separation distance drifts between a pair of satellites within a constellation in highly elliptical orbits [1]. This set of equations was discretized in the true anomaly angle [1] to be used in a digital steady-state hierarchical controller [2]. This controller [2] performed the drift correction between a pair of satellites within the constellation. The objective of a discretized system is to develop a simple algorithm to be implemented in the computer onboard the satellite. The main advantage of the discrete systems is that the computational time can be reduced by selecting a suitable sampling interval. For this digital system, the amount of data will depend on the sampling interval in the true anomaly angle [3]. The purpose of this paper is to implement the discrete Tschauner-Hempel equations and the steady-state hierarchical controller in the computer onboard the satellite. This set of equations is expressed in the true anomaly angle in which a relation will be formulated between the time and the true anomaly angle domains.

Capo-Lugo, Pedro A.↗
NASA NTRS1975

Advances in Application of Fast Semidirect Computational Methods in Transonic Flow

This paper is intended as a review and summary of the advances made in a recently developed approach for rapid numerical solution of the equations of inviscid transonic aerodynamics. The investigation has been limited to two-dimensional, steady, inviscid flow over airfoils in a subsonic free stream, with emphasis on development of a rapid computational technique, rather than on generality of application. The approach uses finite-difference algorithms called "fast direct elliptic solvers" within an iteration scheme. "Direct" means that the entire computation field is solved at once, rather than in successive traverses over the field as in a point- or line-relaxation method. Such an iterative method is referred to as "semidirect." The iterative convergence can be faster than in other relaxation methods because changes are felt simultaneously at all points in each succeeding iteration. Direct elliptic solvers and semidirect methods have restrictions, but these are gradually being removed. Direct solvers were first developed for solving Poisson's equation on a rectangle without interior boundaries. A method to treat first-order systems, a direct Cauchy-Riemann solver has also been developed. Numerical treatment of part of a system of nonlinear equations by a Poisson solver has been reported. Also Poisson solvers in semidirect methods were used for nonseparable elliptic equations. The semidirect method was extended to the solution of a problem of mixed type, where the improved Murman-Cole transonic small-disturbance difference equations were solved. A slightly supercritical flow over a biconvex airfoil was treated successfully, but the iterations did not converge for more strongly supercritical conditions In another work the addition of terms ot both sides of the difference equations stabilized the iteration for supercritical conditions with large supersonic zones. For this, the Cauchy-Riemann solver was revised to incl,ude the needed terms. Most recently, the evaluation of parameters for rapid convergence and comparisons, with Murman's line-relaxation method was described. The method was extended to full second order accuracy in a fully conservative formulation in another work.

Martin, E. Dale↗
NASA NTRS1985

A space-marching method for incompressible Navier-Stokes equations

This paper deals with the development of a space-marching method for incompressible flows. The method solves the continuity and momentum equations as a coupled system at each streamwise station. The character of the system of equations has been changed from elliptic to hyperbolic/parabolic in order to enable the equations to be marched in space. The present method has many advantages compared to the existing parabolic or space-marching methods for incompressible flow: (1) it avoids the solution of Poisson equations, (2) it conserves the mass flow with no additional computation, (3) it does not require the specification of an assumed pressure field when used in the prediction of duct flows. The present method can capture strong secondary velocities and strong transverse pressure gradients. Predictions of the flow through straight and curved ducts are in good agreement with analytical and experimental results.

Pouagare, M.↗
NASA NTRS1977

Supersonic flow about circular cones at large angles of attack - A floating discontinuity approach

The technique of floating shock fitting is adapted to the computation of the inviscid flowfield about circular cones in a supersonic free stream at angles of attack that exceed the cone half-angle. In those regions in which the governing conical equations are mixed elliptic-hyperbolic, the fully hyperbolic form is obtained by the addition of the temporal derivative. The resulting equations are applicable over the complete range of free-stream Mach numbers, angles of attack and cone half-angles for which the bow shock is attached. An explicit finite-difference algorithm is used to obtain the solution by an unsteady relaxation approach. The bow shock, embedded crossflow shock, and vortical singularity in the leeward symmetry plane are all treated as floating discontinuities in a fixed computational mesh. The method yields excellent results for the bow and embedded shocks, however, the solution in the leeward symmetry plane exhibits viscous-like effects and does not appear to adequately predict the behavior of the vortical singularity.

Daywitt, J.↗