SEARCH · Engineering Papers
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 method was developed for determining the fundamental tidal frequencies for closed basins of water, by means of an eigenvalue analysis. The mathematical model employed, was the Laplace tidal equations.
A series representation is presented of the dynamic behavior of incompressible non-viscous liquids contained in partially filled elastic shells of revolution. Using a series in which each term satisfies exactly Laplace's equation, the liquid portion of the problem becomes an integral equation on the wetted wall of the shell. Physical insight required in selecting the terms of the series is such that selection may be readily programed on a computer. Liquid-shell vibration interaction is investigated by modifying an existing shell of revolution vibration program with the series expansion of the liquid pressure loading on the wetted shell wall. Application of the modified program is demonstrated on a hemispherical tank and a realistic scaled model of a launch vehicle tank having a more complex geometry. Comparison of the present analysis with a NASTRAN finite element analysis for a hemispherical tank, and comparison with test results for a launch vehicle tank shows excellent agreement and rapid convergence of the series representation of the liquid to the true solution.
Program performs general analysis of linear and continuous, discrete and sampled-data systems using state-variable techniques. Program is especially suited analysis of linearized control system problems. It also can be used to model system described by combination of differential equations and Laplace transform blocks, such as aircraft control system.
A general technique for expanding an unknown potential field in terms of a linear summation of weighted dipole or quadrupole fields is described. Computational methods were developed for the iterative addition of dipole fields. Various solution potentials were compared inside the boundary with a more precise calculation of the potential to derive optimal schemes for locating the singularities of the dipole fields. Then, the problem of determining solutions to Laplace's equation on an unbounded domain as constrained by pertinent electron trajectory data was considered.
The geomagnetic field in the near-earth magnetotail is made up of three major contributions related to the dipole field, the magnetotail current system, and the magnetopause current system exclusive of the return currents in the magnetotail current system. The present investigation is concerned with the third contribution. The calculation of the magnetospheric field is discussed, taking into account the calculation of the field inside the magnetopause on the basis of the numerical integration of the Biot Savart integral. The representation of the field for distances greater than 10 earth radii on the dark side of the magnetosphere is considered. The field internal to the magnetopause is expressed as the negative gradient of a scalar potential which is expressed as a sum of solutions to Laplace's equation in cylindrical coordinates. The origin of coordinates is the earth.
Consideration is being given to representing the magnetic field in the area 40 deg N to 83 deg N by means of functions in spherical coordinates. A solution to Laplace's equation for the magnetic potential over a restricted area was found, and programming and testing are currently being carried out. Magnetic anomaly modelling is proceeding. The program SPHERE, which was adapted to function correctly on the Cyber computer, is now operational, for deriving gravity and magnetic models in a spherical coordinate system.
An algebraic decomposition of the Jacobian matrix which relates physical and computational variables is presented. This invertible decomposition parameterizes the mesh by the physically intuitive qualities of cell orientation, cell orthogonality, cell volume, and cell aspect ratio. The decomposition can be used to analyze numerically generated curvilinear coordinate meshes and to assess the contribution of the mesh to the truncation error for any specific differential operator and algorithm. This is worked out in detail for Laplace's equation in nonconservative and conservative forms. The analysis is applied to the solution of the full potential code TAIR, showing grid plots, carpet plots, and truncation error for a NACA 0012 airfoil.
The classical methods of initial orbit determination are brought together within a larger viewpoint. This new synthesis stresses that all such techniques follow one of three approaches. Either they seek to compute the orbital element set, or its equivalent, by attacking the differential equations of motion (Laplace), the first integrals of the equations of motion (Taff), or the solution itself (Gauss). The particular technique pursued within a given type of approach should depend upon the nature of the observational data, the amount of a priori information one is willing to presume, and the object of the exercise. This might be a binary star system, a moon, a minor planet, or an artificial satellite. The efficacy of some algorithms for each approach is discussed briefly. Unfortunately, none of them work very well. Extensions of these techniques to radars or laser radars are trivial and have provided no new insights into the overall problem.
A solution for the tunnel wall boundary layer effects for three-dimensional subsonic tunnels is presented. The model potentials are represented with simple singularities placed on the centerline of the tunnel and Laplace's equation in cylindrical coordinates is solved for either the conventional homogeneous slotted-wall boundary condition, the solid-wall viscous boundary condition, or a combination of them. The most pronounced wall boundary layer effect is on solid blockage for completely closed wind tunnels. Boundary layers on the wall reduce the blockage from the solid-wall, no-boundary-layer case in a manner similar to opening slots in a solid wall. Additionally, for solid-wall tunnel configurations, the streamline curvature interference factor is reduced by a significant amount, whereas the lift interference factor at the model station does not depend on the boundary layer parameter. For combination wall configurations, the slot effect of the horizontal walls dominates the viscous effect of the solid sidewalls.
A sequential modular computational scheme has been proposed for the calculation of temperature profiles in the melt and crystal including the melt-crystal interface shape and the crystal pulling rate for a fixed crystal radius. The heat transfer in the melt is assumed to occur by conduction only. The shape of melt-gas meniscus, described by Laplace-Young equation, is incorporated in the model and its effects on interface shape and pulling rate are examined. Further, the model incorporates the detailed radiation interaction among the various surfaces in the puller using the Gebhart enclosure theory which accounts for both direct and reflected radiation. The effects of various process parameters on the interface shape and the pulling rate have been investigated. The results of this study indicate that the shape of the melt-gas meniscus has a significant effect on the pulling rate and the interface shape in conformity with earlier studies. Further, it is shown that a simple radiation model (Stefan's model) is inadequate to model the radiation heat exchange in the Czochralski puller apparatus and the detailed Gebhart analysis is necessary for accurate calculation of both the temperature profile in the crystal as well as the interface shape.
The effect of surface tension on steady-state rotating fluids in a low gravity environment is studied. All the values of the physical parameters used in these calculations, except in the low gravity environments, are based on the measurements carried out by Leslie (1985) in the low gravity environment of a free-falling aircraft. The profile of the interface of two fluids is derived from Laplace's equation relating the pressure drop across an interface to the radii of curvature which has been applied to a low gravity rotating bubble that contacts the container boundary. The interface shape depends on the ratio of gravity to surface tension forces, the ratio of centrifugal to surface tension forces, the contact radius of the interface to the boundary, and the contact angle. The shape of the bubble is symmetric about its equator in a zero-gravity environment. This symmetry disappears and gradually shifts to parabolic profiles as the gravity environment becomes non-zero. The location of the maximum radius of the bubble moves upward from the center of the depth toward the top boundary of the cylinder as gravity increases. The contact radius of interface to the boundary r0 at the top side of cylinder increases and r0 at the bottom side of the cylinder decreases as the gravity environment increases from zero to 1 g.
Theoretical models of planetary-atmosphere tidal fields are examined analytically, comparing models based on geometric-height coordinates with those employing log-pressure coordinates. The relationship between the linearized meteorological variables in the two systems is explored for classical tidal theory and its reduction to the Laplace tidal equation, and it is shown that identical horizontal and vertical equations are obtained. Also considered are the tidal zonal-mean bilinear flux convergences and their Eliassen-Palm formulations. Numerical results for problems involving the earth and Mars atmospheres are presented in graphs and discussed.
An analytic solution illustrating the behavior of incompressible irrotational plasma outflow in two dimensions is given for application to the Jovian magnetosphere. Plasma is assumed to emanate from the Alfven surface and is channeled downtail by the confining magnetopause surface. Conformal transformation methods are used to map this complex magnetotail geometry into a rectangular box for which the solution to Laplace's equation is expressed in simple form. The model suggests that a dipole layer is formed on the dawn side in the vicinity of the Alfven surface where the flow bifurcates into clockwise and counterclockwise motion. An electric field discontinuity of this type fixed in local time may be able to account for several anomalous spin-periodic phenomena, namely, the modulation of Io-independent radio emissions and the injection of energetic electrons into interplanetary space.
The dynamic pole tide is determined by solving Laplace tide equations which take into account the presence of continents in oceans, oceanic self-gravitation and loading, and mantle elasticity. Dynamical effects are found to be only mild. It is shown that the dynamical pole tide contributes about one day more to the Chandler period than a static pole tide would, and dissipates wobble energy at a very weak rate. It is noted that, depending on the wobble period predicted for an oceanless elastic earth, mantle anelasticity at low frequencies may nevertheless contribute negligibly to the Chandler period.
The results of the analysis of nonzonal perturbations (m = 1, 2, 3) of the geopotential field at a 30 mb level are presented. A long period modulation of the harmonics' amplitude is discovered. Calculations of eigenfunctions and eigennumbers of the Laplace tidal equation are carried out for a real latitudinal wind profile. The observed first zonal harmonic in different years is caused by the same mode. Thus, the difference in the wave amplitudes could not be accounted for by the difference in stratospheric zonal circulation in different years and should be related to tropospheric processes.
An important goal of geodesy is to determine the anomalous potential and its derivatives outside of the earth. Representing the surface anomalies by a series of spherical harmonics is useful since it is then possible to do a term by term solution of Laplace's equation and upward continuation. The problem of finding such a spherical harmonic series for anomaly values given on an equiangular surface grid is addressed. (This is a first step toward the more complicated problem of finding a function such that locally averaged values fit a grid of mean anomalies.) Three approaches to this fitting problem are discussed and compared: the discrete Fourier technique, the discrete integral technique, and a new approach. The peculiar nature of the equiangular grid, with its increasing density of (noisy) data toward the poles, causes each method to exhibit a different type of difficulty. The new method is shown to be practical as well as precise since the numerical conditioning problems which appear can be successfully handled by such well-known techniques as a (simple) Kalman filter.
This User's Manual describes the Hypercluster Parallel Processing Library, composed of FORTRAN-callable subroutines which enable a FORTRAN programmer to manipulate and transfer information throughout the Hypercluster at NASA Lewis Research Center. Each subroutine and its parameters are described in detail. A simple heat flow application using Laplace's equation is included to demonstrate the use of some of the library's subroutines. The manual can be used initially as an introduction to the parallel features provided by the library. Thereafter it can be used as a reference when programming an application.
The exact potential corresponding to confining fields inside a linear rf quadrupole particle trap of finite length is presented. The analytic expressions for the trapping potential is derived by introducing a linear trap employing a relatively simple cylindrical geometry and solving Laplace's equation for the trap electrodes. The finite length of linear traps results in field distortion near the trap ends. An exact analytic determination of the fields is useful because the profile of the trapped ion cloud is highly dependent on the fields confining it. It is shown that near the ends of the trap, the effective potential arising from the rf fields acts to propel particles out of the trap, and further, that the addition of a dc bias generates an inhomogeneous in the trap that influences the particles both perpendicularly to and along the trap's long axis.