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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Solution of partial differential equations on vector and parallel computers

The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed.

Ortega, J. M.↗
NASA NTRS1985

Progress and supercomputing in computational fluid dynamics; Proceedings of U.S.-Israel Workshop, Jerusalem, Israel, December 1984

Current developments and future trends in the application of supercomputers to computational fluid dynamics are discussed in reviews and reports. Topics examined include algorithm development for personal-size supercomputers, a multiblock three-dimensional Euler code for out-of-core and multiprocessor calculations, simulation of compressible inviscid and viscous flow, high-resolution solutions of the Euler equations for vortex flows, algorithms for the Navier-Stokes equations, and viscous-flow simulation by FEM and related techniques. Consideration is given to marching iterative methods for the parabolized and thin-layer Navier-Stokes equations, multigrid solutions to quasi-elliptic schemes, secondary instability of free shear flows, simulation of turbulent flow, and problems connected with weather prediction.

Murman, E. M.↗
NASA NTRS1994

Theoretical study of the incompressible Navier-Stokes equations by the least-squares method

Usually the theoretical analysis of the Navier-Stokes equations is conducted via the Galerkin method which leads to difficult saddle-point problems. This paper demonstrates that the least-squares method is a useful alternative tool for the theoretical study of partial differential equations since it leads to minimization problems which can often be treated by an elementary technique. The principal part of the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation consists of two div-curl systems, so the three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy variable and by using the least-squares method, this paper shows that the div-curl system is properly determined and elliptic, and has a unique solution. The same technique then is employed to prove that the Stokes equations are properly determined and elliptic, and that four boundary conditions on a fixed boundary are required for three-dimensional problems. This paper also shows that under four combinations of non-standard boundary conditions the solution of the Stokes equations is unique. This paper emphasizes the application of the least-squares method and the div-curl method to derive a high-order version of differential equations and additional boundary conditions. In this paper, an elementary method (integration by parts) is used to prove Friedrichs' inequalities related to the div and curl operators which play an essential role in the analysis.

Jiang, Bo-Nan↗
NASA NTRS1976

Calculations of polarization and radiance in the atmosphere

Two different methods for the solution of radiative transfer problems, matrix operator and Monte Carlo, are discussed. As an example of the use of the method, results are given for the radiance and polarization of the radiation scattered from haze layers as well as from models of the real atmosphere. The variation of the radiance, polarization, and ellipticity with the aerosol amount in the real atmosphere model is presented. The Monte Carlo method is applied to the problem of calculation of the radiance and polarization of the photons in the atmosphere when there is an ocean as the lower surface. It is found that the turbidity of the ocean can be deduced from the upwelling radiance measured at some height in the atmosphere. The downwelling radiation just beneath the ocean surface is elliptically polarized at those angles where it is derived from the total internal reflection of the upwelling radiation at the ocean surface.

Plass, G. N.↗
NASA NTRS1977

Downstream boundary conditions for viscous flow problems

The problem of the specification of artificial outflow conditions in flow problems is studied. It is shown that for transport type equations incorrect outflow conditions will adversely affect the solution only in a small region near the outflow boundary, while for elliptic equations, e.g. those governing the streamfunction or pressure, a correct boundary specification is essential. In addition, integral outflow boundary conditions for fluid dynamical problems are considered. It is shown that such conditions are well posed, and their effect on the solutions of the Navier-Stokes equations is also considered.

Fix, G.↗
NASA NTRS1991

Simplified computational methods for elastic and elastic-plastic fracture problems

This paper presents an overview of some of the recent (1984-1991) developments in computational/analytical methods in the mechanics of fracture: (i) analytical solutions for elliptical or circular cracks embedded in isotropic or transversely isotropic solids (the crack-plane being at an arbitrary angle to the axis of transverse isotropy), with crack-faces being subjected to arbitrary tractions; (ii) finite-element or boundary-element alternating methods for two- and three-dimensional crack problems; (iii) a "direct-stiffness" method for stiffened panels with flexible fasteners and with multiple cracks, using the alternating method; (iv) multiplesite-damage near a row of fastener holes; (v) analysis of cracks with bonded repair patches; (vi) methods for generation of weight-functions for 2- and 3-D crack problems, and (vii) domain-integral methods for elastic-plastic or inelastic crack mechanics.

Aging Aircraft↗
NASA NTRS1978

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗
NASA NTRS1978

Computerized symbolic manipulation in structural mechanics Progress and potential

Status and recent applications of computerized symbolic manipulation to structural mechanics problems are summarized. The applications discussed include; (1) generation of characteristic arrays of finite elements; (2) evaluation of effective stiffness and mass coefficients of continuum models for repetitive lattice structures; and (3) application of Rayleigh-Ritz technique to free vibration analysis of laminated composite elliptic plates. The major advantages of using computerized symbolic manipulation in each of these applications are outlined. A number of problem areas which limit the realization of the full potential of computerized symbolic manipulation in structural mechanics are examined and some of the means of alleviating them are discussed.

Noor, A. K.↗
NASA NTRS1977

Supercritical cascade design

The use of the method of complex extension to achieve better aerodynamic designs for supercritical cascades applicable to transonic turbomachinery is discussed. The method permits the computation of analytical solutions to elliptic, hyperbolic or mixed second-order partial equations in two dimensions. Boundary value problems formulated to develop an airfoil shape having a prescribed speed distribution for subsonic flow and a nearby speed distribution in the transonic case are also considered. Computing times necessary to run the blade design program are described as acceptably short.

Korn, D. G.↗
NASA NTRS1983

State and model error estimation for elliptic systems: Applications to large antenna static shape determination

This paper outlines the application of various estimation approaches to the problem of static shape determination for large antenna systems. The problem consists of estimating the shape of an antenna surface from measurements of its static deflection. The estimation schemes are based on any one of the modeling options of a single PDE for early insight and understanding, coarse-resolution multiple-PDE models for parametric studies and fine-resolution piecewise-continuum models for detailed design. For any one of these three models, estimator design can be developed using an infinite-dimensional approach, where the necessary finite-element truncation and approximation is conducted after the analytical design has taken place, or it can be based on a finite-dimensional approach, where the model is truncated before the estimation problem is formulated. One of the main objectives of the paper is to develop both approaches while simultaneously investigating their differences and similarities. Simulation results of an application of the finite-dimensional approach to a large parabolic reflector are presented.

Rodriguez, G.↗
NASA NTRS1993

An adaptive grid with directional control

An adaptive grid generator for adaptive node movement is here derived by combining a variational formulation of Winslow's (1981) variable-diffusion method with a directional control functional. By applying harmonic-function theory, it becomes possible to define conditions under which there exist unique solutions of the resulting elliptic equations. The results obtained for the grid generator's application to the complex problem posed by the fluid instability-driven magnetic field reconnection demonstrate one-tenth the computational cost of either a Eulerian grid or an adaptive grid without directional control.

Brackbill, J. U.↗
NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗
NASA NTRS2004

On the Transition from Two-Dimensional to Three-Dimensional MHD Turbulence

We report a theoretical investigation of the robustness of two-dimensional inviscid MHD flows at low magnetic Reynolds numbers with respect to three-dimensional perturbations. We analyze three model problems, namely flow in the interior of a triaxial ellipsoid, an unbounded vortex with elliptical streamlines, and a vortex sheet parallel to the magnetic field. We demonstrate that motion perpendicular to the magnetic field with elliptical streamlines becomes unstable with respect to the elliptical instability once the velocity has reached a critical magnitude whose value tends to zero as the eccentricity of the streamlines becomes large. Furthermore, vortex sheets parallel to the magnetic field, which are unstable for any velocity and any magnetic field, are found to emit eddies with vorticity perpendicular to the magnetic field and with an aspect ratio proportional to N(sup 1/2). The results suggest that purely two-dimensional motion without Joule energy dissipation is a singular type of flow which does not represent the asymptotic behaviour of three-dimensional MHD turbulence in the limit of infinitely strong magnetic fields.

Thess, A.↗
NASA NTRS1976

On the computational cost of approximating and recognizing noise-perturbed straight lines and quadratic arcs in the plane

Recognition of underlying straight lines and quadratic arcs in line drawings, and approximation of very noisy data by such line/curve segments, is addressed as a subproblem of the more general problem of optimum recognition of complicated line/curve drawings. Some specific algorithms are presented with extensions and interpretations for more complicated applications, and a data generation model is developed for the problem. Data are generated as a perturbation of a single underlying straight line or an elliptic or hyperbolic arc. Recursive estimation techniques, minimization of central processing unit time, decision making with controlled error probabilities, and modeling and recognition of pictures consisting of noisy curves are dealt with. Applications envisaged include: picture data compression, contour line representation in maps, intelligent data searches, and ballistic missile decoy tracking.

Cooper, D. B.↗
NASA NTRS1985

The role of service areas in the optimization of FSS orbital and frequency assignments

A relationship is derived, on a single-entry interference basis, for the minimum allowable spacing between two satellites as a function of electrical parameters and service-area geometries. For circular beams, universal curves relate the topocentric satellite spacing angle to the service-area separation angle measured at the satellite. The corresponding geocentric spacing depends only weakly on the mean longitude of the two satellites, and this is true also for alliptical antenna beams. As a consequence, if frequency channels are preassigned, the orbital assignment synthesis of a satellite system can be formulated as a mixed-integer programming (MIP) problem or approximated by a linear programming (LP) problem, with the interference protection requirements enforced by constraints while some linear function is optimized. Possible objective-function choices are discussed and explicit formulations are presented for the choice of the sum of the absolute deviations of the orbital locations from some prescribed ideal location set. A test problem is posed consisting of six service areas, each served by one satellite, all using elliptical antenna beams and the same frequency channels. Numerical results are given for the three ideal location prescriptions for both the MIP and LP formulations. The resulting scenarios also satisfy reasonable aggregate interference protection requirements.

Levis, C. A.↗
NASA NTRS1988

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗
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↗