Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary value problem”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Nonlinear singularly perturbed optimal control problems with singular arcs

A third order, nonlinear, singularly perturbed optimal control problem is considered under assumptions which assure that the full problem is singular and the reduced problem is nonsingular. The separation between the singular arc of the full problem and the optimal control law of the reduced one, both of which are hypersurfaces in state space, is of the same order as the small parameter of the problem. Boundary layer solutions are constructed which are stable and reach the outer solution in a finite time. A uniformly valid composite solution is then formed from the reduced and boundary layer solutions. The value of the approximate solution is that it is relatively easy to obtain and does not involve singular arcs. To illustrate the utility of the results, the technique is used to obtain an approximate solution of a simplified version of the aircraft minimum time-to-climb problem. A numerical example is included.

Ardema, M. D.↗
NASA NTRS1979

Nonlinear singularly perturbed optimal control problems with singular arcs

Singular perturbation techniques are studied for dealing with singular arc problems by analyzing a relatively low-order but otherwise general system. This system encompasses many flight mechanic problems including Goddard's problem and a version of the minimum time-to-climb problem. Boundary layer solutions are constructed which are stable and reach the outer solution in a finite time. A uniformly valid composite solution is then formed from the reduced and boundary layer solutions. The value of the approximate solution is that it is relatively easy to obtain and does not involve singular arcs. To illustrate the utility of the results, the technique is used to obtain an approximate solution of a simplified version of the aircraft minimum time-to-climb problem.

Ardema, M. D.↗
NASA NTRS1977

Stress-intensity factors for cracks emanating from the loaded fastener hole

Using a series approach and the Muskhelishvili formulation in the two-dimensional theory of elasticity, stress-intensity factors K are derived for problems in which cracks emanate radially from the boundary of an arbitrarily loaded internal circular hole in an infinite plate. Numerical values are obtained for K(I) and K(II) for radial cracks from a hole containing a loose-fitted pin or rivet that is pulled perpendicular to the crack direction in the plane of the plate. The method is a general one for determining K for a set of symmetrically emanating radial cracks for a variety of concentrated or distributed tractions on the circular hole.

Shivakumar, V.↗
NASA NTRS1947

Boundary-layer and stalling characteristics of two symmetrical NACA low-drag airfoil sections

Two symmetrical airfoils, an NACA 633-018 and an NACA 631-012, were investigated for the purpose of determining their stalling and boundary-layer characteristics with a view toward the eventual application of this information to the problem of boundary-layer control. Force measurements, pressure distributions, tuft studies, and boundary-layer-profile measurements were made at a value of 5,800,000 Reynolds number. It was found that the 18-percent-thick airfoil stalled progressively from the trailing edge because of separation of the turbulent boundary layer. In contrast, the12-percent-thick airfoil stalled abruptly from a separation of flow near the leading edge before the turbulent boundary layer became subject to separation. From this it was concluded that if high values of lift are to be obtained with thin, high-critical-speed sections by means of boundary-layer control, the work must be directed toward delaying the separation of flow near the leading edge. It was found that the presence of a nose flap on the 12-percent-thick section caused the airfoil to stall in a manner similar to that of the 18-percent-thick section.

Mccullough, George B↗
NASA NTRS1989

A finite difference treatment of Stokes-type flows: Preliminary report

The equations Laplacian operator omega = 0, (1.1a) and omega = Laplacian operator Chi, (1.1b) describe, in suitable units, 2-D Stokes flow of an incompressible fluid occupying a domain D in which omega is the vorticity and Chi is the stream function. The flow is uniquely determined by specifying the velocity on the boundary B of D, a condition which leads to specifying the stream function Chi and its normal derivative Chi sub n on B. A mathematically similar problem arises in describing the equilibrium of a flat plate in structural mechanics where a related 1-D problem by finite difference or finite element methods is to introduce effective methods for imposing the boundary conditions through which (1.1a) is coupled to (1.1b). These models thus provide a simple starting point for examining the general treatment of boundary conditions for more general time dependent Navier-Stokes incompressible flows. For the purpose of discussion it is assumed that D is a square domain. A standard finite difference method to solve (1.1) is to introduce a uniform grid and then use standard five point finite difference operators to express each equation in (1.1). At any point on the boundary B a value of Chi is specified by the boundary conditions but a value of omega at the same boundary mesh point will also be required to complete the computation. Methods are discussed which overcome the difficulty in solving these problems.

Rose, M. E.↗
NASA NTRS1981

Analyses of shuttle orbiter approach and landing conditions

A study of one shuttle orbiter approach and landing conditions are summarized. Causes of observed PIO like flight deficiencies are identified and potential cures are examined. Closed loop pilot/vehicle analyses are described and path/attitude stability boundaries defined. The latter novel technique proved of great value in delineating and illustrating the basic causes of this multiloop pilot control problem. The analytical results are shown to be consistent with flight test and fixed base simulation. Conclusions are drawn relating to possible improvements of the shuttle orbiter/digital flight control system.

Teper, G. L.↗
NASA NTRS1993

Flutter analysis using transversality theory

A new method of calculating flutter boundaries of undamped aeronautical structures is presented. The method is an application of the weak transversality theorem used in catastrophe theory. In the first instance, the flutter problem is cast in matrix form using a frequency domain method, leading to an eigenvalue matrix. The characteristic polynomial resulting from this matrix usually has a smooth dependence on the system's parameters. As these parameters change with operating conditions, certain critical values are reached at which flutter sets in. Our approach is to use the transversality theorem in locating such flutter boundaries using this criterion: at a flutter boundary, the characteristic polynomial does not intersect the axis of the abscissa transversally. Formulas for computing the flutter boundaries and flutter frequencies of structures with two degrees of freedom are presented, and extension to multi-degree of freedom systems is indicated. The formulas have obvious applications in, for instance, problems of panel flutter at supersonic Mach numbers.

Afolabi, D.↗
NASA NTRS1994

Evaluation of water vapor distribution in general circulation models using satellite observations

This paper presents a comparison of the water vapor distribution obtained from two general circulation models, the European Centre for Medium-Range Weather Forecasts (ECMWF) model and the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM), with satellite observations of total precipitable water (TPW) from Special Sensor Microwave/Imager (SSM/I) and upper tropospheric relative humidity (UTH) from GOES. Overall, both models are successful in capturing the primary features of the observed water vapor distribution and its seasonal variation. For the ECMWF model, however, a systematic moist bias in TPW is noted over well-known stratocumulus regions in the eastern subtropical oceans. Comparison with radiosonde profiles suggests that this problem is attributable to difficulties in modeling the shallowness of the boundary layer and large vertical water vapor gradients which characterize these regions. In comparison, the CCM is more successful in capturing the low values of TPW in the stratocumulus regions, although it tends to exhibit a dry bias over the eastern half of the subtropical oceans and a corresponding moist bias in the western half. The CCM also significantly overestimates the daily variability of the moisture fields in convective regions, suggesting a problem in simulating the temporal nature of moisture transport by deep convection. Comparison of the monthly mean UTH distribution indicates generally larger discrepancies than were noted for TPW owing to the greater influence of large-scale dynamical processes in determining the distribution of UTH. In particular, the ECMWF model exhibits a distinct dry bias along the Intertropical Convergence Zone (ITCZ) and a moist bias over the subtropical descending branches of the Hadley cell, suggesting an underprediction in the strength of the Hadley circulation. The CCM, on the other hand, demonstrates greater discrepancies in UTH than are observed for the ECMWF model, but none that are as clearly correlated with well-known features of the large-scale circulation.

Soden, Brian J.↗
NASA NTRS1986

Kelvin-Helmholtz instability at the magnetopause boundary

An MHD model is defined for the convective Kelvin-Helmholtz instability (CKHI) which satellite data indicates produces vortices in the near-earth region of the plasma sheet through turbulence caused by shear flow. The MHD equations are solved as an initial value problem, which reveals that the CKHI growth rate is the same as that of the periodic KHI (PKHI). However, the CKHI is not as limited in the range of amplitude as is the PKHI, which experiences stability in nonlinear growth because of the periodic characteristics of the growth. A perpendicularly moving CKHI can produce a large vortex flow and shocks off the vortex boundary. The satellite data indicate that the vortices are in the size range of 10-20 earth radii. Additional effects can include the expulsion of magnetic flux and the formation of the boundary layer. The success of the model up to the point of formation of the boundary layer is taken as an incentive for further model development using more realistic magnetospheric data.

Wu, C. C.↗
NASA NTRS1973

Laminar film boiling on inclined isothermal flat plates.

Laminar film boiling from an inclined flat plate has been investigated analytically. Using the singular perturbation scheme, the complete set of Navier-Stokes equations is solved. The zeroth-order perturbation coinciding with the boundary-layer equations for vertical flat plates governs the problem. The higher-order perturbations become important near the leading edge and for large values of the inclination angle from the vertical. The assumption of zero interfacial velocity shows that, except for fluids having large (rho x mu) ratios, the results can be predicted using the vertical flat plate results by defining a modified Grashof parameter containing a cos phi term. When the interfacial shear is considered, the solutions indicate that for fluids having large (rho x mu) ratios, the heat transfer rates will be larger (approximately 15% maximum) than those predicted by the simplified model using zero interfacial velocity. In general, the inclination decreases the rate of heat transfer as well as the rate of evaporation.

Nagendra, H. R.↗
NASA NTRS1982

Analyses of Shuttle Orbiter approach and landing

A study of the Shuttle Orbiter approach and landing conditions is summarized. The causes of observed PIO-like flight deficiencies are listed, and possible corrective measures are examined. Closed-loop pilot/vehicle analyses are described, and a description is given of path-attitude stability boundaries. The latter novel approach is found to be of great value in delineating and illustrating the basic causes of this multiloop pilot control problem. It is shown that the analytical results are consistent with flight test and fixed-base simulation. Conclusions are drawn concerning possible improvements in the Shuttle Orbiter/Digital Flight Control System.

Ashkenas, I. L.↗
NASA NTRS1982

Chemorheology of reactive systems: Finite element analysis

The equations which govern the nonisothermal flow of reactive fluids are outlined, and the means by which finite element analysis is used to solve these equations for the sort of arbitrary boundary conditions encountered in industrial practice are described. The performance of the computer code is illustrated by several trial problems, selected more for their value in providing insight to polymer processing flows than as practical production problems. Although a good deal remains to be learned as to the performance and proper use of this numerical technique, it is undeniably useful in providing better understanding of today's complicated polymer processing problems.

Douglas, C.↗
NASA NTRS1983

Boundary conditions at grain boundaries

Grain boundaries in polycrystalline semiconductor material are frequently idealized as surfaces possessing a characteristic surface recombination velocity. Minority carrier densities generated by external sources (electron beams, light, etc.) in polycrystalline p-n junctions must satisfy certain boundary conditions at these grain boundaries which will be derived. It will also be shown that a 'folding technique' introduced recently to deal with this problem is of limited value.

Von Roos, O.↗
NASA NTRS1986

Turbulent relaxation processes in magnetohydrodynamics

Competing processes previously called 'selective decay' and 'dynamic alignment' are studied numerically for two-dimensional magnetohydrodynamic turbulence. In selective decay, the energy decays relatively to mean-square vector potential, and in dynamic alignment, the energy decays relatively to cross helicity. In the former case, the kinetic (fluid) energy decays to zero and the magnetic energy occupies the largest scales allowed by the boundary conditions. In the latter case, the velocity field and magnetic field become aligned and energetically equipartitioned. An extensive study of the initial value problem, with viscous and resistive dissipation, indicates that four distinct regimes of behavior are possible: a magnetically dominated regime, a velocity dominated regime, a dynamic alignment dominated regime, and a transition regime separating the others in parameter space. An analytical variational problem predicts several features seen in the computations, including geometrical alignment of velocity and magnetic fields, constraints on global quantities in the long-time limit, and in some cases, condensation to long wavelength modes.

Ting, Antonio C.↗
NASA NTRS1986

Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.

Gottlieb, D.↗
NASA NTRS1976

A new method for solving problems of sound radiation from a duct for the case of low Mach number

A method for rapidly solving duct end-plane impedance problems is presented which is based on deriving a form of the Helmholtz integral formula expressing the normal derivative of the acoustic pressure at a boundary point. A set of integral expressions is obtained by satisfying the existing boundary conditions and then solved using a general collocation method. The results obtained for the Levine-Schwinger (1948) problem agree well with the exact values. In the case of a duct flow having a temperature mismatch with the surroundings, the results show an increase in the magnitude of the reflection coefficient for the case of a cold core, and a decrease in its value for the case of a hot core, as compared with the case of exhaust into a uniform medium. A phase change of between pi/2 and pi is nearly always achieved, which indicates the tendency towards maintaining a constant pressure at the exit plane. The method can easily be extended to arbitrary duct shapes as long as they are axisymmetric.

Wahbah, M. M.↗