Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

The stabilizing effect of compressibility in turbulent shear flow

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers relevant to homogeneous shear flow are the turbulent Mach number M(t) and the gradient Mach number M(g). Two series of simulations are performed where the initial values of M(g) and M(t) are increased separately. The growth rate of turbulent kinetic energy is observed to decrease in both series of simulations. This 'stabilizing' effect of compressibility on the turbulent energy growth rate is observed to be substantially larger in the DNS series where the initial value of M(g) is changed. A systematic companion of the different DNS cues shows that the compressibility effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence production and not due to explicit dilatational effects. The reduced turbulence production is not a mean density effect since the mean density remains constant in compressible homogeneous shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to increase with the gradient Mach number M(g) in the homogeneous shear flow DNS. Estimates of M(g) for the mixing and the boundary layer are obtained. These estimates show that the parameter M(g) becomes much larger in the high-speed mixing layer relative to the high-speed boundary layer even though the mean flow Mach numbers are the same in the two flows. Therefore, the inhibition of turbulent energy production and consequent 'stabilizing' effect of compressibility on the turbulence (over and above that due to the mean density variation) is expected to be larger in the mixing layer relative to the boundary layer in agreement with experimental observations.

Sarkar, S.↗
NASA NTRS1994

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗
NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗
NASA NTRS1994

Nonlinear Instability of a Uni-directional Transversely Sheared Mean Flow

It is well known that the presence of a weak cross flow in an otherwise two-dimensional shear flow results in a spanwise variation in the mean streamwise velocity profile that can lead to an amplification of certain three-dimensional disturbances through a kind of resonant-interaction mechanism (Goldstein and Wundrow 1994). The spatial evolution of an initially linear, finite-growth-rate, instability wave in such a spanwise-varying shear flow is considered, The base flow, which is governed by the three-dimensional parabolized Navier-Stokes equations, is initiated by imposing a spanwise- periodic cross-flow velocity on an otherwise two-dimensional shear flow at some fixed streamwise location. The resulting mean-flow distortion initially grows with increasing streamwise distance, reaches a maximum and eventually decays through the action of viscosity. This decay, which coincides with the viscous spread of of the shear layer, means that the local growth rate of the instability wave will eventually decrease as the wave propagates downstream. Nonlinear effects can then become important within a thin spanwise-modulated critical layer once the local instability-wave amplitude and growth rate become sufficiently large and small, respectively. The amplitude equation that describes this stage of evolution is shown to be a generalization of the one obtained by Goldstein and Choi (1989) who considered the related problem of the interaction of two oblique modes in a two-dimensional shear layer.

Wundrow, David W.↗
NASA NTRS1994

Three Dimensional Compressible Turbulent Flow Computations for a Diffusing S-Duct With/Without Vortex Generators

Numerical investigations on a diffusing S-duct with/without vortex generators and a straight duct with vortex generators are presented. The investigation consists of solving the full three-dimensional unsteady compressible mass averaged Navier-Stokes equations. An implicit finite volume lower-upper time marching code (RPLUS3D) has been employed and modified. A three-dimensional Baldwin-Lomax turbulence model has been modified in conjunction with the flow physics. A model for the analysis of vortex generators in a fully viscous subsonic internal flow is evaluated. A vortical structure for modeling the shed vortex is used as a source term in the computation domain. The injected vortex paths in the straight duct are compared with the analysis by two kinds of prediction models. The flow structure by the vortex generators are investigated along the duct. Computed results of the flow in a circular diffusing S-duct provide an understanding of the flow structure within a typical engine inlet system. These are compared with the experimental wall static-pressure, static- and total-pressure field, and secondary velocity profiles. Additionally, boundary layer thickness, skin friction values, and velocity profiles in wall coordinates are presented. In order to investigate the effect of vortex generators, various vortex strengths are examined in this study. The total-pressure recovery and distortion coefficients are obtained at the exit of the S-duct. The numerical results clearly depict the interaction between the low velocity flow by the flow separation and the injected vortices.

Cho, Soo-Yong↗
NASA NTRS1994

An oil flow study of the reference Shuttle-C configuration

An oil flow study of the reference Shuttle-C configuration is presented. The Shuttle-C vehicle was a proposed shuttle derived vehicle where the orbiter was to be replaced by an unmanned cargo carrier element. Oil flows are shown for the range of Mach numbers from Mach 1.10 to 3.48 at various angles-of-attack and roll angles. The major flow field phenomena over the Shuttle-C reference configuration are shown in these oil flows. Using the oil flows, a qualitative understanding of the flow around the vehicle can be determined, aiding the quantitative definition of aerodynamic data from theoretical analyses and test results. The oil flows presented in this study were obtained from configurations tested in the NASA Marshall Space Flight Center's 14-Inch Trisonic Wind Tunnel from October 1988 through February 1989.

Pokora, Darlene C.↗
NASA NTRS1994

Local dynamic subgrid-scale models in channel flow

The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Cabot, William H.↗
NASA NTRS1994

Large-eddy simulation of flow through a plane, asymmetric diffuser

A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.

Kaltenbach, Hans-Jakob↗
NASA NTRS1995

Experimental results for a hypersonic nozzle/afterbody flow field

This study was conducted to experimentally characterize the flow field created by the interaction of a single-expansion ramp-nozzle (SERN) flow with a hypersonic external stream. Data were obtained from a generic nozzle/afterbody model in the 3.5 Foot Hypersonic Wind Tunnel at the NASA Ames Research Center, in a cooperative experimental program involving Ames and McDonnell Douglas Aerospace. The model design and test planning were performed in close cooperation with members of the Ames computational fluid dynamics (CFD) team for the National Aerospace Plane (NASP) program. This paper presents experimental results consisting of oil-flow and shadow graph flow-visualization photographs, afterbody surface-pressure distributions, rake boundary-layer measurements, Preston-tube skin-friction measurements, and flow field surveys with five-hole and thermocouple probes. The probe data consist of impact pressure, flow direction, and total temperature profiles in the interaction flow field.

Spaid, Frank W.↗
NASA NTRS1994

Holographic interferometric tomography for reconstructing flow fields

Holographic interferometric tomography is a technique for instantaneously capturing and quantitatively reconstructing three-dimensional flow fields. It has a very useful application potential for high-speed aerodynamics. However, three major challenging tasks need to be accomplished before its practical applications. First, fluid flows are mostly unsteady or at least non repeatable. Consequently, a means for Instantaneously recording three-dimensional flow fields, that is, a simple holographic technique for simultaneously recording multi-directional projections, needs to be developed. Second, while holographic interferometry provides enormous data storage capabilities, expeditious data extraction from complicated interferograms is very important for timely near real-time applications. Third, unlike medical applications, flow tomography does not provide complete data sets but instead involves ill-posed reconstruction problems of incomplete projection and limited angular scanning. During this summer research period, new experimental techniques and corresponding hardware were developed and tested to address the above mentioned tasks. The first task was achieved by diffuser illumination. This concept allows instantaneous capture of many projections with a conventional setup for single-projection recording. For the second task, a phase-shifting technique was incorporated. This technique allows one to acquire multiple phase-stepped interferograms for a single projection and thus to extract phase information from intensity data almost at real-time. For the third task, the research that has been extensively conducted previously was utilized. In this research period, a complete experimental setup that provides the above three major capabilities was designed, built, and tested by integrating all the techniques. A simple laboratory experiment for simulating wind-tunnel testing was then conducted. A test flow was produced by employing a relatively simple device that generated a gravity-driven flow. The flow was then experimentally investigated to check the viability of the holographic interferometric tomographic technique before wind-tunnel application.

Cha, Soyoung S.↗
NASA NTRS1994

High-lift flow-physics flight experiments on a subsonic civil transport aircraft (B737-100)

As part of the subsonic transport high-lift program, flight experiments are being conducted using NASA Langley's B737-100 to measure the flow characteristics of the multi-element high-lift system at full-scale high-Reynolds-number conditions. The instrumentation consists of hot-film anemometers to measure boundary-layer states, an infra-red camera to detect transition from laminar to turbulent flow, Preston tubes to measure wall shear stress, boundary-layer rakes to measure off-surface velocity profiles, and pressure orifices to measure surface pressure distributions. The initial phase of this research project was recently concluded with two flights on July 14. This phase consisted of a total of twenty flights over a period of about ten weeks. In the coming months the data obtained in this initial set of flight experiments will be analyzed and the results will be used to finalize the instrumentation layout for the next set of flight experiments scheduled for Winter and Spring of 1995. The main goal of these upcoming flights will be: (1) to measure more detailed surface pressure distributions across the wing for a range of flight conditions and flap settings; (2) to visualize the surface flows across the multi-element wing at high-lift conditions using fluorescent mini tufts; and (3) to measure in more detail the changes in boundary-layer state on the various flap elements as a result of changes in flight condition and flap deflection. These flight measured results are being correlated with experimental data measured in ground-based facilities as well as with computational data calculated with methods based on the Navier-Stokes equations or a reduced set of these equations. Also these results provide insight into the extent of laminar flow that exists on actual multi-element lifting surfaces at full-scale high-life conditions. Preliminary results indicate that depending on the deflection angle, the slat and flap elements have significant regions of laminar flow over a wide range of angles of attack. Boundary-layer transition mechanisms that were observed include attachment-line contamination on the slat and inflectional instability on the slat and fore flap. Also, the results agree fairly well with the predictions reported in a paper presented at last year's AIAA Fluid Dynamics Conference. The fact that extended regions of laminar flow are shown to exist on the various elements of the high-lift system raises the question what the effect is of loss of laminar flow as a result of insect contamiantion, rain or ice accumulation on high-life performance.

Vandam, Cornelis P.↗
NASA NTRS1993

Unsteady Flow Simulations in Support of the SSME HEX Turning Vane Cracking Investigation with the ATD HPOTP

Unsteady flow computations are being performed with the P&W (ATD) and the Rocketdyne baseline configurations of the SSME LO2 turbine turnaround duct (TAD) and heat exchanger (HEX). The work is in support of the HEX inner turning vane cracking investigation. Fatigue cracking has occurred during hot firings with the P&W configuration on the HEX inner vane, and it appears the fix will involve changes to the TAD splitter vane position and to the TAD inner wall curvature to reduce the dynamic loading on the inner vane. Unsteady flow computations on the P&W baseline and fix and on the Rocketdyne baseline reference follow steady-flow screening computations done by MSFC/ED32 on several trial configurations arriving at the fix. The P&W TAD inlet velocity profile has a strong radial velocity component that directs the flow toward the inner wall and raises the local velocity a factor of two and the dynamic pressure a factor, of four. The fix is intended to redistribute the flow more evenly across the HEX inner and outer vanes like the Rocketdyne baseline reference. Vane buffeting at frequencies around 4,000 Hz is the leading suspected cause of the problem. Our simulations (work in progress) are being done with the USA 2D axisymmetric code approximating the flow as axisymmetric u+v 2D (axial, u, and radial, v, components only). The HEX coils are included in the model to make sure the fix does not adversely affect the HEX environment. Turbulent kinetic energy, k, levels where k = 1/2 v' rms2 are locally as high as 10,000 ft2/sec2 for the P&W baseline at the engine interface (between the TAD and HEX) at the HEX inner vane location. However, k is less than 8,000 on the HEX outer vane and only about 4,500 on the HEX inner vane for the Rocketdyne baseline. Unsteady turbulence intensity, v'rms/v, and pressure, p', are being computed in the present computations to compare with steady-flow Reynolds-averaged computations where p'rms = const (pk) for overall rms random turbulence from 0.1 to 12,000 Hz frequency. Random overall static, p'rms fluctuations as large as 1.7 psi are estimated from k on the HEX inner vane for the P&W baseline configuration but only about 0.7 psi for the Rocketdyne configuration.

Dougherty, N. S.↗
NASA NTRS1995

Calculation of three-dimensional (3-D) internal flow by means of the velocity-vorticity formulation on a staggered grid

A method has been developed to accurately compute the viscous flow in three-dimensional (3-D) enclosures. This method is the 3-D extension of a two-dimensional (2-D) method developed for the calculation of flow over airfoils. The 2-D method has been tested extensively and has been shown to accurately reproduce experimental results. As in the 2-D method, the 3-D method provides for the non-iterative solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body fitted computational mesh incorporating a staggered grid methodology. In the staggered grid method, the three components of vorticity are defined at the centers of the computational cell sides, while the velocity components are defined as normal vectors at the centers of the computational cell faces. The staggered grid orientation provides for the accurate definition of the vorticity components at the vorticity locations, the divergence of vorticity at the mesh cell nodes and the conservation of mass at the mesh cell centers. The solution is obtained by utilizing a fractional step solution technique in the three coordinate directions. The boundary conditions for the vorticity and velocity are calculated implicitly as part of the solution. The method provides for the non-iterative solution of the flow field and satisfies the conservation of mass and divergence of vorticity to machine zero at each time step. To test the method, the calculation of simple driven cavity flows have been computed. The driven cavity flow is defined as the flow in an enclosure driven by a moving upper plate at the top of the enclosure. To demonstrate the ability of the method to predict the flow in arbitrary cavities, results will he shown for both cubic and curved cavities.

Stremel, Paul M.↗
NASA NTRS1995

Numerical Simulation of Flow in a Whirling Annular Seal and Comparison with Experiments

The turbulent flow field in a simulated annular seal with a large clearance/radius ratio (0.015) and a whirling rotor was simulated using an advanced 3D CFD code SCISEAL. A circular whirl orbit with synchronous whirl was imposed on the rotor center. The flow field was rendered quasi-steady by making a transformation to a totaling frame. Standard k-epsilon model with wall functions was used to treat the turbulence. Experimentally measured values of flow parameters were used to specify the seal inlet and exit boundary conditions. The computed flow-field in terms of the velocity and pressure is compared with the experimental measurements inside the seal. The agreement between the numerical results and experimental data with correction is fair to good. The capability of current advanced CFD methodology to analyze this complex flow field is demonstrated. The methodology can also be extended to other whirl frequencies. Half- (or sub-) synchronous (fluid film unstable motion) and synchronous (rotor centrifugal force unbalance) whirls are the most unstable whirl modes in turbomachinery seals, and the flow code capability of simulating the flows in steady as well as whirling seals will prove to be extremely useful in the design, analyses, and performance predictions of annular as well as other types of seals.

Athavale, M. M.↗
NASA NTRS1995

The unique radar scattering properties of silicic lava flows and domes

Silicic (silica-rich) lava flows, such as rhyolite, rhyodacite, and dacite, possess unique physical properties primarily because of the relatively high viscosity of the molten lava. Silicic flows tend to be thicker than basaltic flows, and the resulting large-scale morphology is typically a steep-sided dome or flow lobe, with aspect ratios (height/length) sometimes approaching unity. The upper surfaces of silicic domes and flows are normally emplaced as relatively cool, brittle slabs that fracture as they are extruded from the central vent areas, and are then rafted away toward the flow margin as a brittle carapace above a more ductile interior layer. This mode of emplacement results in a surface with unique roughness characteristics, which can be well-characterized by multiparameter synthetic aperture radar (SAR) observations. In this paper, we examine the scattering properties of several silicic domes in the Inyo volcanic chain in the Eastern Sierra of California, using AIRSAR and TOPSAR data. Field measurements of intermediate-scale (cm to tens of m) surface topography and block size are used to assess the mechanisms of the scattering process, and to quantify the unique roughness characteristics of the flow surfaces.

Plaut, Jeffrey J.↗
NASA NTRS1994

Mapping the dark matter in the NGC 5044 group with ROSAT: Evidence for a nearly homogeneous cooling flow with a cooling wake

The NGC 5044 group of galaxies was observed by the ROSAT Position Sensitive Proportional Counter (PSPC) for 30 ks during its reduced pointed phase (1991 July). Due to the relatively cool gas temperature in the group (kT = 0.98 +/- 0.02 keV) and the excellent photon statistics (65,000 net counts), we are able to determine precisely a number of fundamental properties of the group within 250 kpc of the central galaxy. In particular, we present model-independent measurements of the total gravitating mass, the temperature and abundance profiles of the gas, and the mass accretion rate. Between 60 and 250 kpc, the gas is nearly isothermal with T varies as r(exp (-0.13 +/- 0.03)). The total gravitating mass of the group can be unambiguously determined from the observed density and temperature profiles of the gas using the equation of hydrostatic equilibrium. Within 250 kpc, the gravitating mass is 1.6 x 10(exp 13) solar mass, yielding a mass-to-light ratio of 130 solar mass/solar luminosity. The baryons (gas and stars) comprise 12% of the total mass within this radius. At small radii, the temperature clearly increases outward and attains a maximum value at 60 kpc. The positive temperature gradient in the center of the group confirms the existence of a cooling flow. The cooling flow region extends well beyond the temperature maximum with a cooling radius between 100 and 150 kpc. There are two distinct regions in the cooling flow separated by the temperature maximum. In the outer region, the gas is nearly isothermal with a unifor m Fe abundance of approximately 80% solar, the flow is nearly homogeneous with dot-M= 20 to 25 solar mass/year, the X-ray contours are spherically symmetric, and rho(sub gas) varies as r(exp -1.6). In the inner region, the temperature profile has a positive gradient, the mass accretion rate decreases rapidly inward, the gas density profile is steeper, and the X-ray image shows some substrucutre. NGC 5044 is offset from the centroid of the outer X-ray contours indicating that the central galaxy may have a residual velocity with respect to the center of the group potential. There is also a linear X-ray feature with an extent of approximately 30 kpc with one end coincident with NGC 5044. The X-ray emission from this feature is softer than the ambient gas. We interpret this feature as a 'cooling wake' formed by the accreting gas as it is gravitationally focused into the wake of NGC 5044. One of the most surprising results of our PSPC observation is the discovery of a nearly homogeneous cooling flow. Prior results concerning the mass accretion profile in cooling flows indicate that dot-M varies as r. This relation implies that significant mass deposition occurs at large radii which generates an inhomogeneous flow. The mass accretion rate in the NGC 5044 group is essentially a constant beyond 40 kpc (well within the cooling radius). Significant mass deposition (a declining dot-M) does not commence until the gas accretes to within 40 kpc of the group center where the radiative cooling time is approximately equals 10(exp 9) year. Th is radius also corresponds to the temperature maximum, the break in gas density profile, and the onset of structure in the X-ray image. A Hubble constant of H(sub 0) = 50 km/sec/Mpc is used throughout the paper.

David, Laurence P.↗
NASA NTRS1994

Coupling of volatile transport and internal heat flow on Triton

Recently Brown et al. (1991) showed that Triton's internal heat source could amount to 5-20% of the absorbed insolation on Triton, thus significantly affecting volatile transport and atmospheric pressure. Subsequently, Kirk and Brown (1991a) used simple analytical models of the effect of internal heat on the distribution of volatiles on Triton's surface, confirming the speculation of Brown et al. that Triton's internal heat flow could strongly couple to the surface volatile distribution. To further explore this idea, we present numerical models of the permanent distribution of nitrogen ice on Triton that include the effects of sunlight, the two-dimensional distribution of internal heat flow, the coupling of internal heat flow to the surface distribution of nitrogen ice, and the finite viscosity of nitrogen ice. From these models we conclude that: (1) The strong vertical thermal gradient induced in Triton's polar caps by internal heat-flow facilitates viscous spreading to lower latitudes, thus opposing the poleward transport of volatiles by sunlight, and, for plausible viscosities and nitrogen inventories, producing permanent caps of considerable latitudinal extent; (2) It is probable that there is a strong coupling between the surface distribution of nitrogen ice on Triton and internal heat flow; (3) Asymmetries in the spatial distribution of Triton's heat flow, possibly driven by large-scale, volcanic activity or convection in Triton's interior, can result in permanent polar caps of unequal latitudinal extent, including the case of only one permanent polar cap; (4) Melting at the base of a permanent polar cap on Triton caused by internal heat flow can significantly enhance viscous spreading, and, as an alternative to the solid-state greenhouse mechanism proposed by Brown et al. (1990), could provide the necessary energy, fluids, and/or gases to drive Triton's geyser-like plumes; (5) The atmospheric collapse predicted to occur on Triton in the next 20 years (Spencer, 1990) may be plausibly avoided because of the large latitudinal extent expected for permanent polar caps on Triton.

Brown, Robert H.↗
NASA NTRS1994

The flow in the cometosheath of comet Halley

A quasi-one-dimensional, single-fluid model is used to examine the cometosheath flow along the Sun-comet axis at comet, Halley. The effects of mass loading, dissociative recombination, magnetic field line curvature, ion-neutral drag, and flow divergence on the submagnetosonic flow in this region are discussed. In particular, the inclusion of the flow divergence, which is achieved by allowing the central flow tube area to increase toward the comet, is necessary to bring the cometosheath flow to stagnation at the subsolar point of the ionopause. The cometosheath flow and field profiles are determined for two assumed values of the electron-ion dissociative recombination rate. The model calculations performed using the higher recombination rate are shown to be in better agreement with the Giotto data at comet Halley.

Damas, M. C.↗