SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1996

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗

DOE OSTI2015

Experimental study of a linear/non-linear flux rope

Not Available

DeHaas, Timothy↗

NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗

NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗

NASA NTRS1986

Supersonic full-potential method applied to missile bodies

The NCOREL full-potential method with an entropy correction is presently applied to supersonic missile flowfield problems. After defining the salient characteristics of the method, a combination of linear theory with NCOREL and experimental data is used to isolate the nonlinear features of the supersonic flow so that the influence of geometry and flow conditions on the development of such flow nonlinearities can be appreciated. Comparisons of experimental longitudinal force and moment data with NCOREL and various linear theory predictions are presented for several generic missile airframe configurations of circular and elliptic cross section. The NCOREL code solves the nonconservative full potential equation in a spherical coordinate system; exact boundary conditions are defined on the missile surface.

Pittman, J. L.↗

NASA NTRS1990

Paired and interacting galaxies: Conference summary

The author gives a summary of the conference proceedings. The conference began with the presentation of the basic data sets on pairs, groups, and interacting galaxies with the latter being further discussed with respect to both global properties and properties of the galactic nuclei. Then followed the theory, modelling and interpretation using analytic techniques, simulations and general modelling for spirals and ellipticals, starbursts and active galactic nuclei. Before the conference the author wrote down the three questions concerning pairs, groups and interacting galaxies that he hoped would be answered at the meeting: (1) How do they form, including the role of initial conditions, the importance of subclustering, the evolution of groups to compact groups, and the fate of compact groups; (2) How do they evolve, including issues such as relevant timescales, the role of halos and the problem of overmerging, the triggering and enhancement of star formation and activity in the galactic nuclei, and the relative importance of dwarf versus giant encounters; and (3) Are they important, including the frequency of pairs and interactions, whether merging and interactions are very important aspects of the life of a normal galaxy at formation, during its evolution, in forming bars, shells, rings, bulges, etc., and in the formation and evolution of active galaxies? Where possible he focuses on these three central issues in the summary.

Norman, Colin A.↗

DOE OSTI2017

Quantifying Environmental and Line-of-sight Effects in Models of Strong Gravitational Lens Systems

Matter near a gravitational lens galaxy or projected along the line of sight (LOS) can affect strong lensing observables by more than contemporary measurement errors. We simulate lens fields with realistic three-dimensional mass configurations (self-consistently including voids), and then fit mock lensing observables with increasingly complex lens models to quantify biases and uncertainties associated with different ways of treating the lens environment (ENV) and LOS. We identify the combination of mass, projected offset, and redshift that determines the importance of a perturbing galaxy for lensing. Foreground structures have a stronger effect on the lens potential than background structures, due to nonlinear effects in the foreground and downweighting in the background. There is dramatic variation in the net strength of ENV/LOS effects across different lens fields; modeling fields individually yields stronger priors for H {sub 0} than ray tracing through N-body simulations. Models that ignore mass outside the lens yield poor fits and biased results. Adding external shear can account for tidal stretching from galaxies at redshifts z⩾z{sub lens}, but it requires corrections for external convergence and cannot reproduce nonlinear effects from foreground galaxies. Using the tidal approximation is reasonable for most perturbers as long as nonlinear redshift effects are included. Even then, the scatter in H {sub 0} is limited by the lens profile degeneracy. Asymmetric image configurations produced by highly elliptical lens galaxies are less sensitive to the lens profile degeneracy, so they offer appealing targets for precision lensing analyses in future surveys like LSST and Euclid.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1989

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Stremel, Paul M.↗

DOE OSTI2014

Co-evolution of galactic nuclei and globular cluster systems

We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in the understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M{sub NSC}/M{sub ∗}≈0.0025 M{sub ∗,11}{sup −0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10{sup 11} M {sub ☉}. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach ∼10{sup 5} M {sub ☉} within ≲ 1 Gyr of the beginning of globular cluster formation.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1974

Solution of an eigenvalue problem for the Laplace operator on a spherical surface

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (also referred to as the membrane eigenvalue problem for the vibration equation) on the unit spherical surface are developed. Two specific types of spherical surface domains are considered: (1) the interior of a spherical triangle, i.e., the region bounded by arcs of three great circles, and (2) the exterior of a great circle arc extending for less than pi radians on the sphere (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is derived. The symmetric (generally five-point) finite difference equations that develop are written in matrix form and then solved by the iterative method of point successive overrelaxation. Upon convergence of this iterative method, the fundamental eigenvalue is approximated by iteration utilizing the power method as applied to the finite Rayleigh quotient.

Walden, H.↗

NASA NTRS1988

Viscous induced drag

A fundamentally new approach to the aircraft minimum induced drag problem is presented. The method, a 'viscous lifting line', is based on the minimum entropy production principle and does not require the planar wake assumption. An approximate, closed form solution is obtained and compared with several classical results. In addition, the problem of optimizing in-plane wing sweep with constant wing root bending moment is considered. Like the classical lifting line theory, this theory predicts that induced drag is proportional to the square of the lift coefficient and inversely proportional to the wing aspect ratio. Unlike the classical theory, it predicts that in-plane wing sweep may significantly reduce induced drag, that induced drag is Reynolds number dependent, and that the optimum spanwise circulation distribution is non-elliptic.

Greene, George C.↗

NASA NTRS1976

The log N-log S curve for 3CR radio galaxies and the problem of identifying faint radio galaxies

The log N-log S slope of optically identified galaxies in the 3CR catalog is analyzed for galaxies with flux levels greater than 10 Jy and located at galactic latitudes greater than 7 deg. For the 86 galaxies with known redshifts, the slope is found to be about 1.16; for all 119 galaxies, the slope is found to be about 1.50. If the 113 unidentified sources are included, the slope steepens to about 1.81. It is shown that this steepness is caused by the unidentified sources. Assuming that the slope is due to evolution at redshifts of 1 to 3 and that the unidentified sources are bright ellipticals, the mean redshifts and apparent brightnesses of these galaxies are calculated. The results indicate that it is impossible to observe the unidentified sources with existing ground-based telescopes, implying that it is presently impossible to establish directly that evolution is responsible for the steep log N-log S curve of the 3CR galaxies.

Burbidge, G. R.↗

DOE OSTI2013

Optimal Control of a Parabolic Equation with Dynamic Boundary Condition

We investigate a control problem for the heat equation. The goal is to find an optimal heat transfer coefficient in the dynamic boundary condition such that a desired temperature distribution at the boundary is adhered. To this end we consider a function space setting in which the heat flux across the boundary is forced to be an L{sup p} function with respect to the surface measure, which in turn implies higher regularity for the time derivative of temperature. We show that the corresponding elliptic operator generates a strongly continuous semigroup of contractions and apply the concept of maximal parabolic regularity. This allows to show the existence of an optimal control and the derivation of necessary and sufficient optimality conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Parallel time integration software

This package implements an optimal-scaling multigrid solver for the (non) linear systems that arise from the discretization of problems with evolutionary behavior. Typically, solution algorithms for evolution equations are based on a time-marching approach, solving sequentially for one time step after the other. Parallelism in these traditional time-integrarion techniques is limited to spatial parallelism. However, current trends in computer architectures are leading twards system with more, but not faster. processors. Therefore, faster compute speeds must come from greater parallelism. One approach to achieve parallelism in time is with multigrid, but extending classical multigrid methods for elliptic poerators to this setting is a significant achievement. In this software, we implement a non-intrusive, optimal-scaling time-parallel method based on multigrid reduction techniques. The examples in the package demonstrate optimality of our multigrid-reduction-in-time algorithm (MGRIT) for solving a variety of parabolic equations in two and three sparial dimensions. These examples can also be used to show that MGRIT can achieve significant speedup in comparison to sequential time marching on modern architectures.

Source record↗

NASA NTRS1992

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗

DOE OSTI2014

Parallel time integration software

Source record↗

NASA NTRS1979

The Inertia Coefficients of an Airship in a Frictionless Fluid

The apparent inertia of an airship hull is examined. The exact solution of the aerodynamical problem is studied for hulls of various shapes with special attention given to the case of an ellipsoidal hull. So that the results for the ellipsoidal hull may be readily adapted to other cases, they are expressed in terms of the area and perimeter of the largest cross section perpendicular to the direction of motion by means of a formula involving a coefficient kappa which varies only slowly when the shape of the hull is changed, being 0.637 for a circular or elliptic disk, 0.5 for a sphere, and about 0.25 for a spheroid of fineness ratio. The case of rotation of an airship hull is investigated and a coefficient is defined with the same advantages as the corresponding coefficient for rectilinear motion.

Bateman, H.↗

NASA NTRS1980

Baffle design for earth radiation rejection in the Cryogenic Limb-Scanning Interferometer/Radiometer

The Cryogenic Limb-Scanning Interferometer/Radiometer (CLIR) is being developed to observe infrared emissions of the earth's upper atmosphere from space. The earth's surface is an extended source of intense background radiation with a small angular separation from the desired scene. The CLIR employs an off-axis Gregorian Telescope whose primary mirror and baffles are cooled by an open-cycle cryogen system. A system of specular annular baffles has been developed to minimize both stray light problems and cryogen consumption by retro-mapping the aperture into itself. Each off-axis ray which enters the aperture and strikes the specular baffle surface is reflected so that it passes out of the aperture and is not absorbed on a cryogenic surface. The specular baffle which lies closest to the aperture is an ellipsoid whose foci trace out the circular aperture on revolution about the axis. Its theoretical 'ray trace' efficiency is 100 percent. A subsequent baffle has an elliptical cross section whose near focus traces out the central hole in the ellipsoidal baffle and whose far focus traces out the aperture. Its theoretical efficiency is about 90 percent. These baffles reduce the earth radiation heat load on the cryogenic cooler by an order of magnitude, changing it from the dominant cause of cryogen consumption to a relatively small effect. An aperture shield is also desirable to reduce cryogen consumption, stray light, and contamination.

Bremer, J. C.↗