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 325 records · Page 18

Multigrid techniques for the solution of the passive scalar advection-diffusion equation

The solution of elliptic passive scalar advection-diffusion equations is required in the analysis of many turbulent flow and convective heat transfer problems. The accuracy of the solution may be affected by the presence of regions containing large gradients of the dependent variables. The multigrid concept of local grid refinement is a method for improving the accuracy of the calculations in these problems. In combination with the multilevel acceleration techniques, an accurate and efficient computational procedure is developed. In addition, a robust implementation of the QUICK finite-difference scheme is described. Calculations of a test problem are presented to quantitatively demonstrate the advantages of the multilevel-multigrid method.

Phillips, R. E.↗

Spline-based Rayleigh-Ritz methods for the approximation of the natural modes of vibration for flexible beams with tip bodies

Rayleigh-Ritz methods for the approximation of the natural modes for a class of vibration problems involving flexible beams with tip bodies using subspaces of piecewise polynomial spline functions are developed. An abstract operator-theoretic formulation of the eigenvalue problem is derived and spectral properties investigated. The existing theory for spline-based Rayleigh-Ritz methods applied to elliptic differential operators and the approximation properties of interpolatory splines are used to argue convergence and establish rates of convergence. An example and numerical results are discussed.

Rosen, I. G.↗

Reynolds stress closure in jet flows using wave models

A wave model was developed for the two dimensional shear layer. This configuration is being used as a test case for the closure schemes. Numerical methods are under development to solve the nonseparable Rayleigh equation. A model problem is being used to assist in the algorithm development. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique are being examined.

Morris, P. J.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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.↗

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↗

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.↗

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.↗

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.↗

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.↗

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↗