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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A complete spherical harmonic approach to luni-solar tides

In this work a spherical harmonic theory of ocean tides is presented. The theory is based on Laplace tide equations modified to include turbulence with constant eddy viscosity, linearized bottom friction, and oceanic loading and self-gravitation. Variable bathymetry is also treated in harmonic terms, and no-flow boundary conditions are applied at continental coastlines. The tide and boundary constraint equations are reduced to matrix form and solved by a weighted least-squares procedure. Five zonal luni-solar tides, ranging in period from 14 days to 18.6 yr, are investigated using the theory; such tides have typically been difficult to compute using traditional numerical approaches. The polar motion and changes in the length of day induced by these long-period tides are calculated. Tidal solutions are compared extensively with results from other tidal theories and from recent satellite and sea-level observations. The greatest limitation to accurate prediction of zonal tides - for any theory - appears to be the marginal failure of all tide theories to conserve mass globally; the use of additional mass constraints may be warranted.

Dickman, S. R.↗
NASA NTRS1990

Eigenfrequencies and horizontal structure of divergent barotropic instability originating in tropical latitudes

A systematic investigation of eigenfunctions of the generalized Laplace tidal equation, for the case of tropical-latitude monotonic mean zonal wind profiles with a single narrow reversed vorticity gradient region, has revealed that, in the limit of low planetary zonal wavenumber, the barotropic instability modes bifurcate into weakly divergent modes of hemispheric scale, on the one hand, and on the other strongly divergent 'internal' modes trapped about the source region. These latter disturbances penetrate into the deep tropics as a Kelvin wave-type phenomenon not previously seen in this context. These results suggest that hemispheric barotropic instability need not be purely nondivergent, and that the growth of weakly divergent modes is preferred.

Dunkerton, Timothy J.↗
NASA NTRS1989

Construction of surface pressure field from scatterometer wind field

An account of the construction of surface pressure fields from Seasat-A satellite scatterometer (SASS) winds as carried out by different methods, and the comparison of these pressure fields with those derived from in situ ship observations is presented. On the assumption that the pressure adjusts itself instantaneously to the motion field, it may be computed by various methods. One of these makes use of planetary boundary theory, and of the possible techniques in this category a two-layer iterative scheme admitting of the parametrization of diabatic and baroclinic effects and of secondary flow was chosen. A second method involves the assumption of zero two-dimensional divergence, leading to a Laplace's equation (the balance equation) in pressure, with the wind field serving as a forcing function. This method does not accommodate adiabatic or baroclinic effects, and requires a knowledge of the pressure at all boundary points. Two comparison fields are used for validation: the conventional operational analyses of the US National Meteorological Center (NMC), and the special analyses of the Gulf of Alaska Experiment (GOASEX), which were done by hand. The results of the computations were as follows: (1) The pressure fields, as computed from the SASS winds alone, closely approximated the NMC fields in regions where reasonable in situ coverage was available (typically, one or two mb differences over most of the chart, three to four mb in extreme cases); (2) In some cases the SASS-derived pressure fields displayed high-resolution phenomena not detected by the NMC fields, but evident in the GOASEX data; and, (3) As expected, the pressure fields derived from the balance equation were much smoother and less well resolved than the SASS-derived or NMC fields. The divergence as measured from the SASS winds is smaller than, but of the same order of magnitude as, the vorticity.

Wurtele, Morton G.↗
NASA NTRS1991

Analytic potential in a linear radio-frequency quadrupole trap with cylindrical electrodes

An analytical expression is derived for a radio-frequency ion trap of novel configuration consisting of a four-sectored hollow cylinder enclosed between two end caps. The cylindrical geometry of the sectored trap provides shielding against the buildup of charge and also makes it possible to calculate the potential within the trap by solving Laplace's equation for given boundary conditions. Equations are presented for calculating the time-averaged potential generated by the RF fields, the end-cap potential, and the potential arising from the application of a dc bias on two of the four electrode sectors. It is shown that, near the ends of this trap, the effective potential arising from the RF fields acts to propel ions out of the trap and that the addition of a dc bias on two neighboring sectors generates an inhomogeneous field in the trap which produces a force on the ions along the trap's long axis in a direction dependent on the sign of the bias.

Melbourne, R. K.↗
NASA NTRS1992

Ion collection from a plasma by a pinhole

Ion focusing by a biased pinhole is studied numerically. Laplace's equation is solved in 3-D for cylindrical symmetry on a constant grid to determine the potential field produced by a biased pinhole in a dielectric material. Focusing factors are studied for ions of uniform incident velocity with a 3-D Maxwellian distribution superimposed. Ion currents to the pinhole are found by particle tracking. The focusing factor of positive ions as a function of initial velocity, temperature, injection radius, and hole size is reported. For a typical Space Station Freedom environment (oxygen ions having a 4.5 eV ram energy, 0.1 eV temperature, and a -140 V biased pinhole), a focusing factor of 13.35 is found for a 1.5 mm radius pinhole.

Snyder, David B.↗
NASA NTRS1990

The pole tide in deep oceans

The fluid-dynamical theory of the pole tide is examined by describing the oceanic response to the Chandler wobble and assessing its implications for mantle anelasticity and low-frequency ocean dynamics. The Laplace tide equations accounting for bottom friction are given, and a spherical harmonic approach is delineated in which the time-independent portion of the tide height is expanded. Pole-tide height and related inertia products are linearly proportional to wobble amplitude, and the final equations are modified to account for mantle elasticity and oceanic loading. Results for pole tide effects are given for various earth models with attention to the role of boundary constraints. A dynamic effect is identified which lengthens the Chandler period by about 1 day more than static lengthening, a contribution that suggests a vigorous low-frequency response. The values derived are shown to agree with previous models that do not incorporate the effects of the pole tide.

Dickman, S. R.↗
NASA NTRS1992

Antarctic marine gravity field from high-density satellite altimetry

High-density (about 2-km profile spacing) Geosat/GM altimetry profiles were obtained for Antarctic waters (6-deg S to 72 deg S) and converted to vertical gravity gradient, using Laplace's equation to directly calculate gravity gradient from vertical deflection grids and Fourier analysis to construct gravity anomalies from two vertical deflection grids. The resultant gravity grids have resolution and accuracy comparable to shipboard gravity profiles. The obtained gravity maps display many interesting and previously uncharted features, such as a propagating rift wake and a large 'leaky transform' along the Pacific-Antarctic Rise.

Sandwell, David T.↗
NASA NTRS1991

Global ocean tide mapping using TOPEX/Poseidon altimetry

The investigation's main goals are to produce accurate tidal maps of the main diurnal, semidiurnal, and long-period tidal components in the world's deep oceans. This will be done by the application of statistical estimation techniques to long time series of altimeter data provided by the TOPEX/POSEIDON mission, with additional information provided by satellite tracking data. In the prelaunch phase, we will use in our simulations and preliminary work data supplied by previous oceanographic missions, such as Seasat and Geosat. These results will be of scientific interest in themselves. The investigation will also be concerned with the estimation of new values, and their uncertainties, for tidal currents and for the physical parameters appearing in the Laplace tidal equations, such as bottom friction coefficients and eddy viscosity coefficients. This will be done by incorporating the altimetry-derived charts of vertical tides as boundary conditions in the integration of those equations. The methodology of the tidal representation will include the use of appropriate series expansions such as ocean-basin normal modes and spherical harmonics. The results of the investigation will be space-determined tidal models of coverage and accuracy superior to that of the present numerical models of the ocean tides, with the concomitant benefits to oceanography and associated disciplinary fields.

Sanchez, Braulio V.↗
NASA NTRS1993

Vortex/surface interaction

This paper considers the interaction of a vortex generated upstream in a flow field with a downstream aerodynamic surface that possesses a large chord. The flow is assumed to be steady, incompressible, inviscid and irrotational, and the surface to be semiinfinite. The vortex is considered to be a straight vortex filament. To lowest order the problem is modeled using potential theory, where the 3D Laplace's equation for the velocity potential on the surface is solved exactly. The closed-form equation for pressure distribution obtained from this theory is found to have a square root singularity at the leading-edge. It also converges, as x goes to infinity, to the solution of the 2D point-vortex/infinite plane problem. The pressure coefficient presents an anti-symmetric behavior, near the leading-edge and a symmetric behavior as x goes to infinity.

Bodstein, G. C. R.↗
NASA NTRS1936

Potential flow about elongated bodies of revolution

This report presents a method of solving the problem of axial and transverse potential flows around arbitrary elongated bodies of revolution. The solutions of Laplace's equation for the velocity potentials of the axial and transverse flows, the system of coordinates being an elliptic one in a meridian plane, are given. The theory is applied to a body of revolution obtained from a symmetrical Joukowsky profile, a shape resembling an airship hull. The pressure distribution and the transverse-force distribution are calculated and serve as examples of the procedure to be followed in the case of an actual airship. A section on the determination of inertia coefficients is also included in which the validity of some earlier work is questioned.

Kaplan, Carl↗
NASA NTRS1993

A study of metal-ceramic wettability in SiC-Al using dynamic melt infiltration of SiC

Pressure-assisted infiltration with a 2014 Al alloy of plain and Cu-coated single crystal platelets of alpha silicon carbide was used to study particulate wettability under dynamic conditions relevant to pressure casting of metal-matrix composites. The total penetration length of infiltrant metal in porous compacts was measured at the conclusion of solidification as a function of pressure, infiltration time, and SiC size for both plain and Cu-coated SiC. The experimental data were analyzed to obtain a threshold pressure for the effect of melt intrusion through SiC compacts. The threshold pressure was taken either directly as a measure of wettability or converted to an effective wetting angle using the Young-Laplace capillary equation. Cu coating resulted in partial but beneficial improvements in wettability as a result of its dissolution in the melt, compared to uncoated SiC.

Asthana, R.↗
NASA NTRS1994

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗
NASA NTRS1995

Coarsening strategies for unstructured multigrid techniques with application to anisotropic problems

Over the years, multigrid has been demonstrated as an efficient technique for solving inviscid flow problems. However, for viscous flows, convergence rates often degrade. This is generally due to the required use of stretched meshes (i.e., the aspect-ratio AR = delta y/delta x is much less than 1) in order to capture the boundary layer near the body. Usual techniques for generating a sequence of grids that produce proper convergence rates on isotopic meshes are not adequate for stretched meshes. This work focuses on the solution of Laplace's equation, discretized through a Galerkin finite-element formulation on unstructured stretched triangular meshes. A coarsening strategy is proposed and results are discussed.

Morano, E.↗
NASA NTRS1996

Dynamical Modeling of Surface Tension

In a recent review it is said that free-surface flows 'represent some of the difficult remaining challenges in computational fluid dynamics'. There has been progress with the development of new approaches to treating interfaces, such as the level-set method and the improvement of older methods such as the VOF method. A common theme of many of the new developments has been the regularization of discontinuities at the interface. One example of this approach is the continuum surface force (CSF) formulation for surface tension, which replaces the surface stress given by Laplace's equation by an equivalent volume force. Here, we describe how CSF formulation might be made more useful. Specifically, we consider a derivation of the CSF equations from a minimization of surface energy as outlined by Jacqmin (1996). This reformulation suggests that if one eliminates the computation of curvature in terms of a unit normal vector, parasitic currents may be eliminated. For this reformulation to work, it is necessary that transition region thickness be controlled. Various means for this, in addition to the one discussed by Jacqmin (1996), are discussed.

Brackbill, Jeremiah U.↗
NASA NTRS1996

Coarsening Strategies for Unstructured Multigrid Techniques with Application to Anisotropic Problems

Over the years, multigrid has been demonstrated as an efficient technique for solving inviscid flow problems. However, for viscous flows, convergence rates often degrade. This is generally due to the required use of stretched meshes (i.e. the aspect-ratio AR = (delta)y/(delta)x much less than 1) in order to capture the boundary layer near the body. Usual techniques for generating a sequence of grids that produce proper convergence rates on isotropic meshes are not adequate for stretched meshes. This work focuses on the solution of Laplace's equation, discretized through a Galerkin finite-element formulation on unstructured stretched triangular meshes. A coarsening strategy is proposed and results are discussed.

Morano, E.↗
NASA NTRS1995

NAS Experiences of Porting CM Fortran Codes to HPF on IBM SP2 and SGI Power Challenge

Current Connection Machine (CM) Fortran codes developed for the CM-2 and the CM-5 represent an important class of parallel applications. Several users have employed CM Fortran codes in production mode on the CM-2 and the CM-5 for the last five to six years, constituting a heavy investment in terms of cost and time. With Thinking Machines Corporation's decision to withdraw from the hardware business and with the decommissioning of many CM-2 and CM-5 machines, the best way to protect the substantial investment in CM Fortran codes is to port the codes to High Performance Fortran (HPF) on highly parallel systems. HPF is very similar to CM Fortran and thus represents a natural transition. Conversion issues involved in porting CM Fortran codes on the CM-5 to HPF are presented. In particular, the differences between data distribution directives and the CM Fortran Utility Routines Library, as well as the equivalent functionality in the HPF Library are discussed. Several CM Fortran codes (Cannon algorithm for matrix-matrix multiplication, Linear solver Ax=b, 1-D convolution for 2-D datasets, Laplace's Equation solver, and Direct Simulation Monte Carlo (DSMC) codes have been ported to Subset HPF on the IBM SP2 and the SGI Power Challenge. Speedup ratios versus number of processors for the Linear solver and DSMC code are presented.

Saini, Subhash↗
NASA NTRS1998

Parallelization of Rocket Engine Simulator Software (PRESS)

We have outlined our work in the last half of the funding period. We have shown how a demo package for RESSAP using MPI can be done. However, we also mentioned the difficulties with the UNIX platform. We have reiterated some of the suggestions made during the presentation of the progress of the at Fourth Annual HBCU Conference. Although we have discussed, in some detail, how TURBDES/PUMPDES software can be run in parallel using MPI, at present, we are unable to experiment any further with either MPI or PVM. Due to X windows not being implemented, we are also not able to experiment further with XPVM, which it will be recalled, has a nice GUI interface. There are also some concerns, on our part, about MPI being an appropriate tool. The best thing about MPr is that it is public domain. Although and plenty of documentation exists for the intricacies of using MPI, little information is available on its actual implementations. Other than very typical, somewhat contrived examples, such as Jacobi algorithm for solving Laplace's equation, there are few examples which can readily be applied to real situations, such as in our case. In effect, the review of literature on both MPI and PVM, and there is a lot, indicate something similar to the enormous effort which was spent on LISP and LISP-like languages as tools for artificial intelligence research. During the development of a book on programming languages [12], when we searched the literature for very simple examples like taking averages, reading and writing records, multiplying matrices, etc., we could hardly find a any! Yet, so much was said and done on that topic in academic circles. It appears that we faced the same problem with MPI, where despite significant documentation, we could not find even a simple example which supports course-grain parallelism involving only a few processes. From the foregoing, it appears that a new direction may be required for more productive research during the extension period (10/19/98 - 10/18/99). At the least, the research would need to be done on Windows 95/Windows NT based platforms. Moreover, with the acquisition of Lahey Fortran package for PC platform, and the existing Borland C + + 5. 0, we can do work on C + + wrapper issues. We have carefully studied the blueprint for Space Transportation Propulsion Integrated Design Environment for the next 25 years [13] and found the inclusion of HBCUs in that effort encouraging. Especially in the long period for which a map is provided, there is no doubt that HBCUs will grow and become better equipped to do meaningful research. In the shorter period, as was suggested in our presentation at the HBCU conference, some key decisions regarding the aging Fortran based software for rocket propellants will need to be made. One important issue is whether or not object oriented languages such as C + + or Java should be used for distributed computing. Whether or not "distributed computing" is necessary for the existing software is yet another, larger, question to be tackled with.

Cezzar, Ruknet↗
NASA NTRS2011

Existence and Stability of Menisci in Detached Bridgman Growth

Detached growth, also referred to as dewetted growth, is a Bridgman crystal growth process in which the melt is in contact with the crucible wall but the crystal is not. A meniscus bridges the gap between the top of the crystal and the crucible wall. The Young-Laplace capillary equation was used to calculate the crystal radii of detached states as a function of the pressure differential across the meniscus. The detached states depend on the contact angle of the melt with the crucible wall, the growth angle of the melt with respect to the solidifying crystal, and the Bond number. A static stability analysis was performed on the calculated detached states. The stability criterion was the sign of the second variation of the potential energy upon admissible meniscus shape perturbations. The conditions considered corresponded to the growth of Ge and InSb, in both terrestrial and microgravity conditions. Stability was found to depend significantly on whether the interior surface was considered to be microscopically rough or smooth, corresponding to pinned or unpinned states. It was also found that all meniscus shapes which are single-valued functions of the radius are statically stable in a microgravity environment.

Volz, M. P.↗