SEARCH · Search NASA
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.
A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.
The basic crack problem which is essential for the study of subcritical crack propagation and fracture of layered structural materials is considered. Because of the apparent analytical difficulties, the problem is idealized as one of plane strain or plane stress. An additional simplifying assumption is made by restricting the formulation of the problem to crack geometries and loading conditions which have a plane of symmetry perpendicular to the interface. The general problem is formulated in terms of a coupled system of four integral equations. For each relevant crack configuration of practical interest, the singular behavior of the solution near and at the ends and points of intersection of the cracks is investigated and the related characteristic equations are obtained. The edge crack terminating at and crossing the interface, the T-shaped crack consisting of a broken layer and a delamination crack, the cross-shaped crack which consists of a delamination crack intersecting a crack which is perpendicular to the interface, and a delamination crack initiating from a stress-free boundary of the bonded layers are some of the practical crack geometries considered.
An algorithm was developed and incorporated into a computer program for solving the damped, Timoshenko beam equations with free-free boundary conditions and prescribed initial data. A number of special cases are considered and, where appropriate, comparisons are made with known results.
Two theoretical problems are presented to illustrate the equilibrium structure of a gravitating gas threaded by magnetic fields. The gas is heated everywhere at a rate proportional to the local density and cooled by thermal conduction channeled along the magnetic field. Emphasis is placed on the influence of anisotropic thermal conduction on the distribution of the gas. In the first problem, the magnetic field has straight field lines that are nearly parallel. The simple field line geometry allows us to calculate, in closed form, the equilibria configurations of a gas slab which is cooled nonuniformly by thermal conduction along field lines which leave the gas slab at varying angles. The problem involves the gas having a free boundary, which is one of the unknowns to be determined. A discussion is given of the set of all possible equilibria so constructed. In the second problem there is no electric current flowing in the gas but there is a potential magnetic field having curved field lines. Exact solutions of equilibriums are presented to illustrate the heating of stellar and galactic coronas.
Axisymmetric equilibrium shapes of conducting drops and bubbles, (1) pendant or sessile on one face of a circular parallel-plate capacitor or (2) free and surface-charged, are found by solving simultaneously the free boundary problem consisting of the augmented Young-Laplace equation for surface shape and the Laplace equation for electrostatic field, given the surface potential. The problem is nonlinear and the method is a finite element algorithm employing Newton iteration, a modified frontal solver, and triangular as well as quadrilateral tessellations of the domain exterior to the drop in order to facilitate refined analysis of sharply curved drop tips seen in experiments. The stability limit predicted by this computer-aided theoretical analysis agrees well with experiments.
The quasi-steady evolution of solar magnetic fields in response to gradual photospheric changes is considered, with particular attention given to the threshold of a sudden eruption in the solar atmosphere. The formal model of an evolving, force-free field dependent on two Cartesian coordinates is extended to a field which is not force free but in static equilibrium with plasma pressure and gravity. The basic physics is illustrated through the evolution of a loop-shaped electric current sheet enclosing a potential bipolar field with footpoints rooted in the photosphere. A free-boundary problem is posed and then solved for the equilibrium configuration of the current sheet in a hydrostatically supported isothermal atmosphere. As the footpoints move apart to spread a constant photospheric magnetic flux over a larger region, the equilibria available extend the field to increasing heights.
A description is presented of the results of a research program directed at an improved understanding of condensate deposition rate phenomena in combustion systems. The conducted experiments make use of real-time optical laser reflectance-interference-polarization techniques in flame environments. The obtained new data and the results of previous gravimetric experiments are employed as a basis for the development of a comprehensive convective diffusion deposition theory, taking into account the assumption of a multicomponent vapor or multisize class particles 'source-free' boundary layer. The theory makes it possible to provide self-consistent salt/ash/soot deposition rate predictions over a wide variety of environmental conditions.
The computer code for calculating web temperature distribution was expanded to provide a graphics output in addition to numerical and punch card output. The new code was used to examine various modifications of the J419 configuration and, on the basis of the results, a new growth geometry was designed. Additionally, several mathematically defined temperature profiles were evaluated for the effects of the free boundary (growth front) on the thermal stress generation. Experimental growth runs were made with modified J419 configurations to complement the modeling work. A modified J435 configuration was evaluated.
Internal and edge cracks were studied by using Reissner's transverse shear theory. The effect of stress-free boundaries on the stress intensity factors in plates under bending were investigated. Among the results found, particularly interesting are those relating to the limiting cases of the crack geometries. The numerical results are given for a single internal crack, two collinear cracks, and two edge cracks. The effect of Poisson's ratio on the stress intensity factors was studied.
Rayleigh's integral formula is evaluated numerically for planar radiators of any shape, with any specified velocity in the source plane using the fast Fourier transfrom algorithm. The major advantage of this technique is its speed of computation - over 400 times faster than a straightforward two-dimensional numerical integration. The technique is developed for computation of the radiated pressure in the nearfield of the source and can be easily extended to provide, with little computation time, the vector intensity in the nearfield. Computations with the FFT of the nearfield pressure of baffled rectangular plates with clamped and free boundaries are compared with the 'exact' solution to illuminate any errors. The bias errors, introduced by the FFT, are investigated and a technique is developed to significantly reduce them.
The basic crack problem which is essential for the study of subcritical crack propagation and fracture of layered structural materials is considered. Because of the apparent analytical difficulties, the problem is idealized as one of plane strain or plane stress. An additional simplifying assumption is made by restricting the formulation of the problem to crack geometries and loading conditions which have a plane of symmetry perpendicular to the interface. The general problem is formulated in terms of a coupled systems of four integral equations. For each relevant crack configuration of practical interest, the singular behavior of the solution near and at the ends and points of intersection of the cracks is investigated and the related characteristic equations are obtained. The edge crack terminating at and crossing the interface, the T-shaped crack consisting of a broken layer and a delamination crack, the cross-shaped crack which consists of a delamination crack intersecting a crack which is perpendicular to the interface, and a delamination crack initiating from a stress-free boundary of the bonded layers are some of the practical crack geometries considered. Previously announced in STAR as N80-18428 and N80-18429
In this paper the title problem is studied by using Reissner's transverse shear theory. The main purpose of the paper is to investigate the effect of stress-free boundaries on the stress intensity factors in plates under bending. Among the results found particularly interesting are those relating to the limiting cases of the crack geometries. The numerical results are given for a single internal crack, two collinear cracks, and two edge cracks. Also studied is the effect of Poisson's ratio on the stress intensity factors.
Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.
It is observed that a complete and explicit example of Hopf bifurcation appears not to be known in fluid mechanics. Such an example is presented for the rotating Benard problem with free boundary conditions on the upper and lower faces, and horizontally periodic solutions. Normal modes are found for the linearization, and the Veronis computation of the wave numbers is modified to take into account the imposed horizontal periodicity. An invariant subspace of the phase space is found in which the hypotheses of the Joseph-Sattinger theorem are verified, thus demonstrating the Hopf bifurcation. The criticality calculations are carried through to demonstrate rigorously, that the bifurcation is subcritical for certain cases, and to demonstrate numerically that it is subcritical for all the cases in the paper.
Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.
Design principles suggested here aim toward high conversion efficiency (greater than 15 percent) in polysilicon cells. The principles seek to decrease the liabilities of both intragranular and grain-boundary-surface defects. The advantages of a phosphorus atom concentration gradient in a thin (less than 50 microns) base of a p(+)/n(x)/n(+) drift-field solar cell, which produces favorable gradients in chemical potential, minority-carrier mobility and diffusivity, and recombination lifetime (via phosphorus gettering) are suggested. The degrading effects of grain boundaries are reduced by these three gradients and by substituting atoms (P, H, F or Li) for vacancies on the grain-boundary surface. From recent experiments comes support for the benefits of P diffusion down grain boundaries and, for quasi-grain-boundary-free and related structures. New analytic solutions for the n(x)-base include the effect of a power-law dependence between P concentration and lifetime. These provide an upper-bound estimate on the open circuit voltage. Finite-difference numerical solutions of the six Shockley equations furnish complete information about all solar-cell parameters and add insight concerning design.
A general adaptive grid algorithm is developed on triangular grids. The adaptivity is provided by a combination of node addition, dynamic node connectivity and a simple node movement strategy. While the local restructuring process and the node addition mechanism take place in the physical plane, the nodes are displaced on a monitor surface, constructed from the salient features of the physical problem. An approximation to mean curvature detects changes in the direction of the monitor surface, and provides the pulling force on the nodes. Solutions to the axisymmetric Grad-Shafranov equation demonstrate the capturing, by triangles, of the plasma-vacuum interface in a free-boundary equilibrium configuration.
An inverse swept wing code is described that is based on the widely used transonic flow program FLO22. The new code incorporates a free boundary algorithm permitting the pressure distribution to be prescribed over a portion of the wing surface. A special routine is included to calculate the wave drag, which can be minimized in its dependence on the pressure distribution. An alternate formulation of the boundary condition at infinity was introduced to enhance the speed and accuracy of the code. A FORTRAN listing of the code and a listing of a sample run are presented. There is also a user's manual as well as glossaries of input and output parameters.