Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAPLACE EQUATION”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

PC-CUBE: A Personal Computer Based Hypercube

PC-CUBE is an ensemble of IBM PCs or close compatibles connected in the hypercube topology with ordinary computer cables. Communication occurs at the rate of 115.2 K-band via the RS-232 serial links. Available for PC-CUBE is the Crystalline Operating System III (CrOS III), Mercury Operating System, CUBIX and PLOTIX which are parallel I/O and graphics libraries. A CrOS performance monitor was developed to facilitate the measurement of communication and computation time of a program and their effects on performance. Also available are CXLISP, a parallel version of the XLISP interpreter; GRAFIX, some graphics routines for the EGA and CGA; and a general execution profiler for determining execution time spent by program subroutines. PC-CUBE provides a programming environment similar to all hypercube systems running CrOS III, Mercury and CUBIX. In addition, every node (personal computer) has its own graphics display monitor and storage devices. These allow data to be displayed or stored at every processor, which has much instructional value and enables easier debugging of applications. Some application programs which are taken from the book Solving Problems on Concurrent Processors (Fox 88) were implemented with graphics enhancement on PC-CUBE. The applications range from solving the Mandelbrot set, Laplace equation, wave equation, long range force interaction, to WaTor, an ecological simulation.

Ho, Alex↗
NASA NTRS1988

A numerical method for computing unsteady 2-D boundary layer flows

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results.

Krainer, Andreas↗
NASA NTRS1990

Cascade aerodynamic gust response including steady loading effects

To predict the unsteady convected gust aerodynamic response of a cascade comprised of arbitrary thick and cambered aerofoils in an incompressible, inviscid, flow field, a complete first-order model is formulated. The flow is analyzed by considering a periodic flow channel. The velocity potential is separated into steady and unsteady harmonic components, each described by a Laplace equation. The strong dependence of the unsteady aerodynamics on the steady effects of aerofoil and cascade geometry and incidence angle is manifested in the coupling of the unsteady and steady flow fields through the unsteady boundary conditions. Analytical solutions in individual grid elements of a body-fitted computational grid are then determined, with the complete solution obtained by assembly of these local solutions. The validity and capabilities of this model and solution technique are then demonstrated by analyzing the steady and unsteady aerodynamics of both theoretical and experimental cascade configurations.

Chiang, Hsiao-Wei D.↗
NASA NTRS1991

Electron distribution functions in electric field environments

The amount of current carried by an electric discharge in its early stages of growth is strongly dependent on its geometrical shape. Discharges with a large number of branches, each funnelling current to a common stem, tend to carry more current than those with fewer branches. The fractal character of typical discharges was simulated using stochastic models based on solutions of the Laplace equation. Extension of these models requires the use of electron distribution functions to describe the behavior of electrons in the undisturbed medium ahead of the discharge. These electrons, interacting with the electric field, determine the propagation of branches in the discharge and the way in which further branching occurs. The first phase in the extension of the referenced models , the calculation of simple electron distribution functions in an air/electric field medium, is discussed. Two techniques are investigated: (1) the solution of the Boltzmann equation in homogeneous, steady state environments, and (2) the use of Monte Carlo simulations. Distribution functions calculated from both techniques are illustrated. Advantages and disadvantages of each technique are discussed.

Rudolph, Terence H.↗
NASA NTRS1990

Keyhole and weld shapes for plasma arc welding under normal and zero gravity

A first order study of the interfacial (keyhole) shape between a penetrating argon plasma arc jet and a stationary liquid metal weld pool is presented. The interface is determined using the Young-Laplace equation by assuming that the plasma jet behaves as a one-dimensional ideal gas flow and by neglecting flow within the weld pool. The solution for the keyhole shape allows an approximate determination of the liquid-solid metal phase boundary location based on the assumption that the liquid melt is a stagnant thermal boundary layer. Parametric studies examine the effect of plasma mass flow rate, initial plasma enthalpy, liquid metal surface tension, and jet shear on weldment shape under both normal and zero gravity. Among the more important findings of this study is that keyhole and weld geometries are minimally affected by gravity, suggesting that data gathered under gravity can be used in planning in-space welding.

Keanini, R. G.↗
NASA NTRS1993

An augmented Young-Laplace model of an evaporating meniscus in a micro-channel with high heat flux

High flux evaporations from a steady meniscus formed in a 2 micron channel is modeled using the augmented Young-Laplace equation. The heat flux is found to be a function of the long range van der Waals dispersion force which represents interfacial conditions between heptane and various substrates. Heat fluxes of (1.3-1.6) x 10(exp 6) W/m(exp 2) based on the width of the channel are obtained for heptane completely wetting the substrate at 100 C. Small channels are used to obtain these large fluxes. Even though the real contact angle is 0 deg, the apparent contact angle is found to vary between 24.8 deg and 25.6 deg. The apparent contact angle, which represents viscous losses near the contact line, has a large effect on the heat flow rate because of its effect on capillary suction and the area of the meniscus. The interfacial heat flux is modeled using kinetic theory for the evaporation rate. The superheated state depends on the temperature and the pressure of the liquid phase. The liquid pressure differs from the pressure of the vapor phase due to capillarity and long range van der Waals dispersion forces which are relevant in the ultra think film formed at the leading edge of the meniscus. Important pressure gradients in the thin film cause a substantial apparent contact angle for a complete wetting system. The temperature of the liquid is related to the evaporation rate and to the substrate temperature through the steady heat conduction equation. Conduction in the liquid phase is calculated using finite element analysis except in the vicinity of the thin film. A lubrication theory solution for the thin film is combined with the finite element analysis by the method of matched asymptotic expansions.

Wayner, P. C., Jr.↗
NASA NTRS1994

Solution of the Chapman-Ferraro problem with an arbitrary magnetopause

We present a global model of the magnetic field of the magnetosphere that includes the effects of the Chapman-Ferraro currents at the magnetopause. In contrast to ealier models, the magnetopause shape is arbitrary, thus allowing the use of more realistic geometries. The internal magnetospheric field model of Hilmer and Voigt (1993), is completely shielded within the magnetopause by solving the Laplace equation with Neumann boundary conditions using a finite difference method on a non-orthogonal, curvilinear grid. The resulting model magnetosphere is perfectly closed although the method can also be applied with more general boundary conditions, to generate a set of open models based on the approach of Toffoletto and Hill (1989, 1993). The purpose of this paper is to demonstrate the feasibility of a purely numerical approach to solving the Chapman-Ferraro problem with arbitrary magnetopause shape and boundary conditions.

Toffoletto, F. R.↗
NASA NTRS1995

Analysis of discretization errors in LES

All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Ghosal, Sandip↗
NASA NTRS1996

Electrohydrodynamic Interaction of a Pair of Spherical Drops

The axisymmetric electrohydrodynamic interaction between two spherical emulsion drops has been examined, using the leaky dielectric model to represent the constitutive behavior of the liquid phases. The results follow from the general solutions in bispherical coordinates to the Laplace equation for the electric potential and the Stokes equations for the velocity field. For drops of similar composition, the electrical interactions induced between the drops by the imposition of the electric field are always attractive, meaning they favor coalescence of the drop pair. The hydrodynamic interactions, however, are not always favorable and, indeed, are shown in certain circumstances to drive the drops apart.

Erker, J. A.↗
NASA NTRS1995

Interfacial Force Field Characterization in a Constrained Vapor Bubble Thermosyphon

Isothermal profiles of the extended meniscus in a quartz cuvette were measured in the earth's gravitational field using an image-analyzing interferometer that is based on computer-enhanced video microscopy of the naturally occurring interference fringes. These profiles are a function of the stress field. Experimentally, the augmented Young-Laplace equation is an excellent model for the force field at the solid-liquid-vapor interfaces for heptane and pentane menisci on quartz and tetradecane on SFL6. The effects of refractive indices of the solid and liquid on the measurement techniques were demonstrated. Experimentally obtained values of the disjoining pressure and dispersion constants were compared to those predicted from the Dzyaloshinskii - Lifshitz - Pilaevskii theory for an ideal surface and reasonable agreements were obtained. A parameter introduced gives a quantitative measurement of the closeness of the system to equilibrium. The nonequilibrium behavior of this parameter is also presented

DasGupta, Sunando↗
NASA NTRS1996

Development of an Automatic Grid Generator for Multi-Element High-Lift Wings

The procedure to generate the grid around a complex wing configuration is presented in this report. The automatic grid generation utilizes the Modified Advancing Front Method as a predictor and an elliptic scheme as a corrector. The scheme will advance the surface grid one cell outward and the newly obtained grid is corrected using the Laplace equation. The predictor-corrector step ensures that the grid produced will be smooth for every configuration. The predictor-corrector scheme is extended for a complex wing configuration. A new technique is developed to deal with the grid generation in the wing-gaps and on the flaps. It will create the grids that fill the gap on the wing surface and the gap created by the flaps. The scheme recognizes these configurations automatically so that minimal user input is required. By utilizing an appropriate sequence in advancing the grid points on a wing surface, the automatic grid generation for complex wing configurations is achieved.

Eberhardt, Scott↗
NASA NTRS2004

On the Inversion for Mass (Re)Distribution from Global (Time-Variable) Gravity Field

The well-known non-uniqueness of the gravitational inverse problem states the following: The external gravity field, even if completely and exactly known, cannot Uniquely determine the density distribution of the body that produces the gravity field. This is an intrinsic property of a field that obeys the Laplace equation, as already treated in mathematical as well as geophysical literature. In this paper we provide conceptual insight by examining the problem in terms of spherical harmonic expansion of the global gravity field. By comparing the multipoles and the moments of the density function, we show that in 3-S the degree of knowledge deficiency in trying to inversely recover the density distribution from external gravity field is (n+l)(n+2)/2 - (2n+l) = n(n-1)/2 for each harmonic degree n. On the other hand, on a 2-D spherical shell we show via a simple relationship that the inverse solution of the surface density distribution is unique. The latter applies quite readily in the inversion of time-variable gravity signals (such as those observed by the GRACE space mission) where the sources over a wide range of the scales largely come from the Earth's Surface.

Chao, Benjamin F.↗
NASA NTRS2013

Zero-gravity Mean Free Surface Curvature of a Confined Liquid in a Radially-Vaned Container

A variety of increasingly intricate container geometries are under consideration for the passive manipulation of liquids aboard spacecraft where the impact of gravity may be neglected. In this study we examine the mean curvature of a liquid volume confined within a radial array of disconnected vanes of infinite extent. This particular geometry possesses a number of desirable characteristics relevant to waste water treatment aboard spacecraft for life support. It is observed that under certain conditions the slender shape of the free surface approaches an asymptote, which can be predicted analytically using new hybrid boundary conditions proposed herein. This contribution represents possibly the final extension of what has been referred to as the method of de Lazzer et al. (1996). The method enables the integration of the Young-Laplace equation over a domain with its boundaries, including the wetted portion of the solid boundaries, symmetry planes, and circular arcs representing free surfaces at the center plane of the liquid body. Asymptotic solutions at several limits are obtained and the analysis is confirmed with numerical computations.

Chen, Yongkang↗
NASA NTRS2017

Acoustic Effects in Classical Nucleation Theory

The effect of sound wave oscillations on the rate of nucleation in a parent phase can be calculated by expanding the free energy of formation of a nucleus of the second phase in powers of the acoustic pressure. Since the period of sound wave oscillation is much shorter than the time scale for nucleation, the acoustic effect can be calculated as a time average of the free energy of formation of the nucleus. The leading non-zero term in the time average of the free energy is proportional to the square of the acoustic pressure. The Young-Laplace equation for the surface tension of the nucleus can be used to link the time average of the square of the pressure in the parent phase to its time average in the nucleus of the second phase. Due to the surface tension, the pressure in the nuclear phase is higher than the pressure in the parent phase. The effect is to lower the free energy of formation of the nucleus and increase the rate of nucleation.

Baird, J. K.↗
NASA NTRS1976

The theory of the gravitational potential applied to orbit prediction

A complete derivation of the geopotential function and its gradient is presented. Also included is the transformation of Laplace's equation from Cartesian to spherical coordinates. The analytic solution to Laplace's equation is obtained from the transformed version, in the classical manner of separating the variables. A cursory introduction to the method devised by Pines, using direction cosines to express the orientation of a point in space, is presented together with sample computer program listings for computing the geopotential function and the components of its gradient. The use of the geopotential function is illustrated.

Kirkpatrick, J. C.↗
NASA NTRS1991

Nonlinear oscillations of inviscid free drops

The present analysis of free liquid drops' inviscid oscillations proceeds through solution of Bernoulli's equation to obtain the free surface shape and of Laplace's equation for the velocity potential field. Results thus obtained encompass drop-shape sequences, pressure distributions, particle paths, and the temporal evolution of kinetic and surface energies; accuracy is verified by the near-constant drop volume and total energy, as well as the diminutiveness of mass and momentum fluxes across drop surfaces. Further insight into the nature of oscillations is provided by Fourier power spectrum analyses of mode interactions and frequency shifts.

Patzek, T. W.↗