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 37 records · Page 2

Determination of the dispersion constant in a constrained vapor bubble thermosyphon

The isothermal profiles of the extended meniscus in a quartz cuvette were measured in a gravitational field using an image analyzing interferometer which is based on computer enhanced video microscopy of the naturally occurring interference fringes. The experimental results for heptane and pentane menisci were analyzed using the extended Young-Laplace equation. These isothermal results characterized the interfacial force field in-situ at the start of the heat transfer experiments by quantifying the dispersion constant, which is a function of the liquid-solid system and cleaning procedures. The experimentally obtained values of the disjoining pressure and the dispersion constants were compared to that predicted from the DLP theory and good agreements were obtained. The measurements are critical to the subsequent non-isothermal experiments because one of the major variables in the heat sink capability of the CVBT is the dispersion constant. In all previous studies of micro heat pipes the value of the dispersion constant has been 'guesstimated'. One of the major advantages of the current glass cell is the ability to view the extended meniscus at all times. Experimentally. we find that the extended Young-Laplace equation is an excellent model for the force field at the solid-liquid-vapor interfaces.

Dasgupta, SUNANDO.↗

Determination of the dispersion constant in a constrained vapor bubble thermosyphon

The isothermal profiles of the extended meniscus in a quartz cuvette were measured in a gravitational field using an image analyzing interferometer which is based on computer enhanced video microscopy of the naturally occurring interference fringes. The experimental results for heptane and pentane menisci were analyzed using the extended Young Laplace Equation. These isothermal results characterized the interfacial force field in-siru at the start of the heat transfer experiments by quantifying the dispersion constant, which is a function of the liquid-solid system and cleaning procedures. The experimentally obtained values of the disjoining pressure and the dispersion constants were compared to that predicted from the DLP theory and good agreements were obtained. The measurements are critical to the subsequent non-isothermal experiments because one of the major variables in the heat sink capability of the Constrained Vapor Bubble Thermosyphon, CVBT, is the dispersion constant. In all previous studies of micro heat pipes the value of the dispersion constant has been 'estimated'. One of the major advantages of the current glass cell is the ability to view the extended meniscus at all times. Experimentally, we find that the extended Young-Laplace Equation is an excellent model for the force field at the solid-liquid-vapor interfaces.

Dasgupta, Sunando↗

Evaluation of 3 numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finitie volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, S. F.↗

Evaluation of three numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finite volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, Steven F.↗

Premixed-Gas Flame Propagation in Hele-Shaw Cells

It is well known that buoyancy and thermal expansion affect the propagation ra and shapes of premixed gas flames. The understanding of such effects is complicated by the large density ratio between the reactants and products, which induces a baroclinic production of vorticity due to misalignment of density and pressure gradients at the front, which in turn leads to a complicated multi-dimensional flame/flow interaction. The Hele-Shaw cell, i.e., the region between closely-spaced flat parallel plates, is probably the simplest system in which multi-dimensional convection is presents consequently, the behavior of fluids in this system has been studied extensively (Homsy, 1987). Probably the most important characteristic of Hele-Shaw flows is that when the Reynolds number based on gap width is sufficiently small, the Navier-Stokes equations averaged over the gap reduce to a linear relation, namely a Laplace equation for pressure (Darcy's law). In this work, flame propagation in Hele-Shaw cells is studied to obtain a better understanding of buoyancy and thermal expansion effects on premixed flames. This work is also relevant to the study of unburned hydrocarbon emissions produced by internal combustion engines since these emissions are largely a result of the partial burning or complete flame quenching in the narrow, annular gap called the "crevice volume" between the piston and cylinder walls (Heywood, 1988). A better understanding of how flames propagate in these volumes through experiments using Hele-Shaw cells could lead to identification of means to reduce these emissions.

Sharif, J.↗

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗

The unified acoustic and aerodynamic prediction theory of advanced propellers in the time domain

This paper presents some numerical results for the noise of an advanced supersonic propeller based on a formulation published last year. This formulation was derived to overcome some of the practical numerical difficulties associated with other acoustic formulations. The approach is based on the Ffowcs Williams-Hawkings equation and time domain analysis is used. To illustrate the method of solution, a model problem in three dimensions and based on the Laplace equation is solved. A brief sketch of derivation of the acoustic formula is then given. Another model problem is used to verify validity of the acoustic formulation. A recent singular integral equation for aerodynamic applications derived from the acoustic formula is also presented here.

Farassat, F.↗

Subsonic and Supersonic Effects in Bose-Einstein Condensate

A paper presents a theoretical investigation of subsonic and supersonic effects in a Bose-Einstein condensate (BEC). The BEC is represented by a time-dependent, nonlinear Schroedinger equation that includes terms for an external confining potential term and a weak interatomic repulsive potential proportional to the number density of atoms. From this model are derived Madelung equations, which relate the quantum phase with the number density, and which are used to represent excitations propagating through the BEC. These equations are shown to be analogous to the classical equations of flow of an inviscid, compressible fluid characterized by a speed of sound (g/Po)1/2, where g is the coefficient of the repulsive potential and Po is the unperturbed mass density of the BEC. The equations are used to study the effects of a region of perturbation moving through the BEC. The excitations created by a perturbation moving at subsonic speed are found to be described by a Laplace equation and to propagate at infinite speed. For a supersonically moving perturbation, the excitations are found to be described by a wave equation and to propagate at finite speed inside a Mach cone.

Zak, Michail↗

A Study of Multigrid Preconditioners Using Eigensystem Analysis

The convergence properties of numerical schemes for partial differential equations are studied by examining the eigensystem of the discrete operator. This method of analysis is very general, and allows the effects of boundary conditions and grid nonuniformities to be examined directly. Algorithms for the Laplace equation and a two equation model hyperbolic system are examined.

Roberts, Thomas W.↗

Stability analysis of intermediate boundary conditions in approximate factorization schemes

In many cases, approximate factorization schemes have provided a significant increase in efficiency over previously used solution methods in certain problems. The present investigation is concerned with the importance of intermediate boundary conditions in approximate factorization schemes, taking into account a specific example regarding a boundary-induced stability restriction in a scheme for the transonic full-potential equation. The considered scheme has been discussed by Holst (1979). Holst's scheme is a variation of the AF2 schemes described by Ballhaus and Steger (1975). The application of the AF2 scheme to the two-dimensional Laplace equation in a rectangle is studied giving attention to the stability of the scheme in connection with various boundary conditions for the intermediate variable.

South, J. C., Jr.↗

Validating corrosion models: Influence of physical properties

In chloride containing environments, two metals in physical contact can undergo galvanic corrosion limiting the lifetime of components. Being able to accurately predict galvanic corrosion damage distributions over time has not been widely presented in literature. Therefore, Finite Element Method (FEM) corrosion models were experimentally validated for two galvanic couples as a function of governing equations, environment, anode material, anode:cathode ratio, and time. Carbon steel/stainless steel (CS/SS) and zinc/stainless steel (Zn/SS) galvanic couples were exposed to NaCl solutions at room temperatures for up to 14 days. For the galvanic couples and environments, the Laplace equation with variable conductivity and reactions was sufficient to model the experimental corrosion damage. Increasing the chloride concentration for the CS/SS galvanic couple, regardless of the anode:cathode ratio, decreased the observed and modeled corrosion damage. Decreasing the anode:cathode ratio (i.e., increasing the cathode length), increased the experimentally observed and modeled corrosion damage. For small anode:cathode ratios, the governing equations deviate over long time periods with the Nernst-Plank equation being conservative. For the Zn/SS galvanic couple, the cathode length controls the dominant cathodic reduction reaction. For small cathode lengths, the hydrogen evolution reaction is dominant. For larger cathode lengths, the oxygen reduction reaction is dominant. The results are discussed with regard to the influence of solution chemistry, ohmic drop, and governing reactions. Overall, validated FEM models were presented, and the resultant models, physics, and mechanisms can be applied to other corrosion scenarios with confidence.

Carbon steel↗

Application of a viscous-inviscid interaction procedure to predict separated flows with heat transfer

A viscous-inviscid interaction procedure is described for predicting heat transfer in separated flows. The separating flow in a rearward-facing step/asymmetric channel expansion is considered. For viscous regions, the boundary-layer momentum and continuity equations are solved inversely in a coupled manner by a finite-difference numerical scheme. The streamwise convective term is altered to permit marching the solution through regions of reversed flow. The inviscid flow is computed by numerically solving the Laplace equation for stream-function in the region bounded by the displacement surfaces used in the inverse boundary-layer solution. The viscous and inviscid solutions are repeated iteratively until the edge velocities obtained from both solutions are in agreement. Predictions using this method compare favorably with experimental data and other predictions.

Hall, E. J.↗

On the normal modes of Laplace's tidal equations for zonal wavenumber zero

The characteristic differences between two different rotational modes of Laplace's tidal equations for wavenumber m = 0, called the K- and the S-modes, are compared in their energy ratio and structures. It is shown that the K-mode representation captures most of the observed zonal energy with a few terms, whereas the S-mode representation requires many terms. For small vertical scale components, the K-mode series converges faster than the S-mode series. Attention is also given to the differences between the energy spectra projected upon the K- and S-modes and the merits of each set as expansion functions for the zonal atmospheric motions.

Tanaka, H. L.↗

Measurement of the Surface Dilatational Viscosity of an Insoluble Surfactant Monolayer at the Air/Water Interface Using a Pendant Drop Apparatus

When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Lorenzo, Jose↗

An extended numerical manifold method for unsaturated soil-water interaction analysis at micro-scale

To investigate unsaturated soil-water interaction at micro-scale, this work extends the numerical manifold method (NMM) by incorporating a soil-water coupling model considering specific capillary water distribution and capillary force calculation. The soil skeleton is constructed by a soil skeleton generation algorithm with random polygons. To more realistically capture the interaction between soil grains and capillary water, a capillary mechanics-based geometric algorithm is proposed to iteratively calculate the capillary water distribution. The capillary forces corresponding to the capillary water distribution are calculated based on the Young-Laplace equation. The proposed capillary water solving framework is first verified by reproducing the soil-water characteristic curve and the capillary water distribution of an ideal contact-disk model against analytical solutions. To further validate the ability of the capillary water solving framework to predict hydraulic behavior of the real soil, a laboratory test on the Toyoura sand is reproduced numerically. Then an ideal direct shear test is performed to further validate the two-way soil-water coupling procedure, in which a comparison between the numerical and analytical results regarding the shear strength and matric suction is presented. Finally, microscopic hydraulic and compression tests are conducted on two soil specimens with the same porosity and mean grain diameter but different uniformity coefficients. The results elucidate that the extended method is a potential tool to explore unsaturated soil behaviors at micro-scale.

58 GEOSCIENCES↗