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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

A Systematic Solution Approach for Neutron Transport Problems in Diffuse Regimes

A systematic solution approach for the neutron transport equation, based on a least-squares finite-element discretization, is presented. This approach includes the theory for the existence and uniqueness of the analytical as well as of the discrete solution, bounds for the discretization error, and guidance for the development of an efficient multigrid solver for the resulting discrete problem. To guarantee the accuracy of the discrete solution for diffusive regimes, a scaling transformation is applied to the transport operator prior to the discretization. The key result is the proof of the V-ellipticity and continuity of the scaled least-squares bilinear form with constants that are independent of the total cross section and the absorption cross section. For a variety of least-squares finite-element discretizations this leads to error bounds that remain valid in diffusive regimes. Moreover, for problems in slab geometry a full multigrid solver is presented with V(1, 1)-cycle convergence rates approximately equal to 0.1, independent of the size of the total cross section and the absorption cross section.

Manteuffel, T. A.↗
NASA NTRS1981

Mesh generation by conformal and quasiconformal mappings

It is pointed out that many recent advances in the finite-difference solution of elliptic equations have been limited to regions whose boundary contours coincide with coordinate lines of the Cartesian coordinate system. The reason for this is related to the fact that in the case of an arbitrary curvilinear coordinate system the original equation becomes much more complex. However, there is no added complexity if an orthogonal coordinate system is generated from a conformal mapping. In the present investigation, a finite difference method developed for the construction of conformal mappings has been generalized to construct quasi-conformal mappings. It is expected that the use of more sophisticated numerical algorithms could lead to improvements in both speed and accuracy. Quasi-conformal mappings have applications not only in the solution of elliptic equations but also in other areas such as orthogonal mesh generation on surfaces and the solution of certain fluid flow problems.

Mastin, C. W.↗
NASA NTRS1985

Kalman-like estimation for static distributed systems Antenna shape from radiation measurements

This paper advances an approach to the determination of shape of static distributed systems. It also illustrates the application of the approach to the problems of surface diagnosis of large parabolic reflectors. The estimation methods developed combine in an optimal sense the information from an elliptic model of the structure and from measurements of the structural deflection and of the far-field pattern changes due to the structural deformation. The estimators have a predictor-corrector structure, quite similar to that of a Kalman filter. The system model is first used to obtain a predicted estimate. A correction term is then added to the prediction to obtain the final state estimate. The relative weighting between prediction and correction terms is determined by an estimator gain. As in a Kalman filter, the estimator gain can be expressed in terms of the state estimation error covariance.

Rodriguez, G.↗
NASA NTRS1988

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

Three-dimensional internal viscous flow problems are presently addressed by a consistent and efficient set of boundary conditions for the multisweep space-marching, pressure-elliptic Reduced Navier-Stokes scheme. Since continuity is by this means directly satisfied at all points in the flow domain, the first-order momentum equations are directly solvable for pressure without the requirement for a Poisson pressure-correction equation. Incompressible flow solutions are obtained for straight and curved ducts of square cross section, in order to validate the procedure. Usefulness is demonstrated for internal flows with strong interactions, as would be found in turbomachine geometries.

Reddy, D. R.↗
NASA NTRS2005

Perturbed Equations of Motion for Formation Flight Near the Sun-Earth L2 Point

This Memorandum Report consists of a compilation of three individual reports, of increasing complexity, describing investigations of formation flight of spacecraft in the vicinity of the L2 Sun-Earth 1ibration point. The individual reports form the following parts of this compilation: - Introduction to the relative motion of spacecraft about the Sun-Earth L2 Point - Linear and quadratic modelling and solution of the relative motion - Modelling the Perturbations - Elliptical Earth Orbit, Lunar Gravity, Solar Radiation Pressure, Thrusters. The three parts are self-contained, with somewhat, varying notation and terminology. After fair1y significant literature searches: this new work (of Parts 2 and 3) is deemed to be unique because it describes the primary perturbations to the description of relative motion between nearby spacecraft. The effect of the elliptical motion of the Earth about the Sun was verified to be the dominant perturbation to the circular restricted three body problem. Contributions due to lunar gravity and solar radiation pressure are seen to have much smaller effect.

Segerman, Alan M.↗
NASA NTRS2008

Limits to Open Class Performance?

This presentation discusses open or unlimited class aircraft performance limitations and design solutions. Limitations in this class of aircraft include slow climbing flight which requires low wing loading, high cruise speed which requires high wing loading, gains in induced or viscous drag alone which result in only half the gain overall and other structural problems (yaw inertia and spins, flutter and static loads integrity). Design solutions include introducing minimum induced drag for a given span (elliptical span load or winglets) and introducing minimum induced drag for a bell shaped span load. It is concluded that open class performance limits (under current rules and technologies) is very close to absolute limits, though some gains remain to be made from unexplored areas and new technologies.

Bowers, Albion H.↗
NASA NTRS1980

Predicted vortex shedding from noncircular bodies in supersonic flow

A method to predict nose vortex shedding from circular and noncircular bodies in supersonic flow at angles of attack and roll is presented. The body is represented by a supersonic panel method, and the lee side vortex wake is modeled by discrete vortices in crossflow planes. The three-dimensional steady flow problem is reduced to a two-dimensional, unsteady, separated flow problem for solution. Comparison of measured and predicted surface pressure distributions, flow field surveys, and aerodynamic characteristics are presented for circular and elliptic bodies.

Mendenhall, M. R.↗
NASA NTRS1995

On the consistency of Reynolds stress turbulence closures with hydrodynamic stability theory

The consistency of second-order closure models with results from hydrodynamic stability theory is analyzed for the simplified case of homogeneous turbulence. In a recent study, Speziale, Gatski, and MacGiolla Mhuiris showed that second-order closures are capable of yielding results that are consistent with hydrodynamic stability theory for the case of homogeneous shear flow in a rotating frame. It is demonstrated in this paper that this success is due to the fact that the stability boundaries for rotating homogeneous shear flow are not dependent on the details of the spatial structure of the disturbances. For those instances where they are -- such as in the case of elliptical flows where the instability mechanism is more subtle -- the results are not so favorable. The origins and extent of this modeling problem are examined in detail along with a possible resolution based on rapid distortion theory (RDT) and its implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1996

The Effects of Dissipation and Coarse Grid Resolution for Multigrid in Flow Problems

The objective of this paper is to investigate the effects of the numerical dissipation and the resolution of the solution on coarser grids for multigrid with the Euler equation approximations. The convergence is accomplished by multi-stage explicit time-stepping to steady state accelerated by FAS multigrid. A theoretical investigation is carried out for linear hyperbolic equations in one and two dimensions. The spectra reveals that for stability and hence robustness of spatial discretizations with a small amount of numerical dissipation the grid transfer operators have to be accurate enough and the smoother of low temporal accuracy. Numerical results give grid independent convergence in one dimension. For two-dimensional problems with a small amount of numerical dissipation, however, only a few grid levels contribute to an increased speed of convergence. This is explained by the small numerical dissipation leading to dispersion. Increasing the mesh density and hence making the problem over resolved increases the number of mesh levels contributing to an increased speed of convergence. If the steady state equations are elliptic, all grid levels contribute to the convergence regardless of the mesh density.

Eliasson, Peter↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗
NASA NTRS1974

Transonic conical flow

The problem of inviscid, steady transonic conical flow, formulated in terms of the small disturbance theory, is studied. The small disturbance equation and similarity rules are presented, and a boundary value problem is formulated for the case of a supersonic freestream Mach number. The equation for the perturbation potential is solved numerically using an elliptic finite difference system. The difference equations are solved with a point relaxation algorithm that is also capable of capturing the shock wave during the iteration procedure by using the boundary conditions at the shock. Numerical calculations, for shock location, pressure distribution and drag coefficient, are presented for a family of nonlifting conical wings. The theory of slender wings is also presented and analytical results for pressure and drag coefficients are obtained.

Agopian, K. G.↗
NASA NTRS1974

Finite element solution theory for three-dimensional boundary flows

A finite element algorithm is derived for the numerical solution of a three-dimensional flow field described by a system of initial-valued, elliptic boundary value partial differential equations. The familiar three-dimensional boundary layer equations belong to this description when diffusional processes in only one coordinate direction are important. The finite element algorithm transforms the original description into large order systems of ordinary differential equations written for the dependent variables discretized at node points of an arbitrarily irregular computational lattice. The generalized elliptic boundary conditions is piecewise valid for each dependent variable on boundaries that need not explicitly coincide with coordinate surfaces. Solutions for sample problems in laminar and turbulent boundary flows illustrate favorable solution accuracy, convergence, and versatility.

Baker, A. J.↗
NASA NTRS1986

Bodies with noncircular cross sections and bank-to-turn missiles

An evaluation is made of prospective missile applications for noncircular cross section bodies, and of recent developments in bank-to-turn missile configuration aerodynamics. The discussion encompasses cross-flow analysis techniques, as well as study results obtained for bodies with elliptical and square cross sections and with variable cross sections. Attention is given to both the performance advantages and the stability and control problems of bank-to-turn missile configurations; the aerodynamic data presented for monoplanar configurations extend to those incorporating airbreathing propulsion systems.

Jackson, C. M., Jr.↗
NASA NTRS2002

Analytical Modal Analysis for Thin-Film Flat Lenses

Due to strong potential applications and more demanding requirements imposed upon thin-film structures for space deployable, there has been increasing research and development activities during recent years in the field of vibration analysis of these types of structures. Moreover, interests in employing these structural components have received renewed emphasis in recent years within NASA and the Air Force. This is due to their inherent lightweight, low packaging and launch volume, and relative simplicity of deployment. Among the potential mission concepts for which these structural elements are included, one can mention solar sails, space solar power generation systems, solar thermal propulsion vehicles, large space telescopes, and inflatable communication antennas. This paper presents analytical procedures to determine vibration and physical characteristics of thin film lenses with circular and elliptical shapes membranes considered in design of a solar concentrator. In general, three methods are used to obtain approximate solutions of Helmholtz boundary value problems. One method requires that solution satisfy the differential equation exactly and the boundary condition approximately. Another method demands a solution that satisfies the boundary conditions exactly and the governing equations approximately. The third method sees a solution that satisfies both the governing equation and boundary conditions approximately. Extensive reviews of vibrations of membrane and plates are provided by Leissa and Mazumdar.

Hamid R. Hamidzadeh↗
NASA NTRS1980

Doubly-periodic orbits in the Sun-Earth-Moon system

A series of periodic orbits in the Earth-Moon circular restricted problem of three bodies was found which is ideally suited for exploring the Earth's geomagnetic tail. The mean apsidal motion of the basic highly elliptical Earth orbit was maintained at about one degree per day by a sequence of lunar swingbys, keeping the apogees in the anti-Sun direction. The orbits were periodic in reference frames rotating at both lunar and solar rates. Apogee distances were alternately raised and lowered by the lunar swingby maneuvers. Several categories of these Sun-synchronous double lunar swingby orbits were identified. The strength and flexibility of this trajectory concept was demonstrated with real world simulations.

Farohar, R.↗
NASA NTRS1979

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Periaux, J.↗