Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity 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 271 records · Page 15

A case study using kinematic quantities derived from a triangle of VHF Doppler wind profilers

Horizontal divergence, relative vorticity, kinematic vertical velocity, and geostrophic and ageostrophic winds are computed from Colorado profiler network data to investigate an upslope snowstorm in northeastern Colorado. Horizontal divergence and relative vorticity are computed using the Gauss and Stokes theorems, respectively. Kinematic vertical velocities are obtained from the surface to 9 km by vertically integrating the continuity equation. The geostrophic and ageostrophic winds are computed by applying a finite differencing technique to evaluate the derivatives in the horizontal equations of motion. Comparison of the synoptic-scale data with the profiler network data reveals that the two datasets are generally consistent. Also, the profiler-derived quantities exhibit coherent vertical and temporal patterns consistent with conceptual and theoretical flow fields of various meteorological phenomena. It is suggested that the profiler-derived quantities are of potential use to weather forecasters in that they enable the dynamic and kinematic interpretation of weather system structure to be made and thus have nowcasting and short-term forecasting value.

Carlson, Catherine A.↗

Large velocity shears and associated electrostatic waves and turbulence in the auroral F region

Broadband electrostatic waves at 10-1000 Hz have been observed with very large shears in the plasma flow velocity transverse to the ambient magnetic field in the auroral F region. The shears were detected through their perpendicular electric field signatures, which changed by as much as 200 mV/m over distances of only a few hundred meters. Transverse shears can be uniquely related to field-aligned currents through the current continuity equation, and the resulting field-aligned drift exceeds the threshold for excitation of current-driven electrostatic ion-cyclotron waves. A numerical simulation of this instability has been used to generate electric-field spectra in the rocket frame of reference, and these spectra are similar to the spectra generated form the actual rocket data.

Earle, G. D.↗

Multigrid Calculations of a Jet in Crossflow

A series of numerical calculations are made of a jet in crossflow using progressive mesh refinement up to 2.4 million grid points. A solution algorithm that combined the features of a point-coupled solution of the momentum and continuity equations with multigrid acceleration demonstrated efficient convergence of the governing equations. In comparisons with experimental data, some features of the jet flow field were found to be predicted reasonably well, but the results still remain affected by numerical diffusion.

Claus, R. W.↗

Stratospheric constituent trends from ER-2 profile data

Trend analysis of the Airborne Arctic Stratospheric Expedition's ER-2 profile data reveals an average decrease in N2O on potential temperature isentropes which can be attributed to diabatic cooling of inner vortex air. This conclusion is independently supported by radiative transfer computations. Trends in ozone and water vapor over the same period are not consistent with the magnitude of the diabatic descent. After accounting for the diabatic motion (estimated from N2O), an additional 0.44 + or - 0.3 percent/day average anomalous O3 decrease above 440 K (about 20 km) is needed to balance the continuity equation. This ozone decrease suggests additional photochemical destruction of ozone in the presence of the high amounts of ClO observed during the mission. A 0.4 + or - 0.3 percent/day average anomalous increase in H2O is also observed near 420 K (about 18 km) which may be due to the evaporation of ice crystals falling from higher, colder stratospheric layers.

Schoeberl, Mark R.↗

Aerodynamic stability and heating analyses for the aeroassist flight experiment vehicle

Since ground based flow simulations are presently unable to model flight conditions expected for AOTVs (Aeroassist Orbital Transfer Vehicle) and other hypersonic space vehicles, computer codes are being developed to provide design parameters necessary for structure, guidance, and control aspects. Over the past four years, VRFLO (Viscous Reactive Flow) has been written to model finite-rate chemistry and viscous effects for a variety of aerobrake bodies. VRFLO includes a number of unique features that are summarized as follows: (1) Grid generation is an integral part of the code for several aerobrake configurations which includes the wake flow region; (2) The formulation is valid for three air chemical models; (3) An ADI central difference technique is used to solve the Navier-Stokes and species continuity equations in split groups; and (4) Grid density and numerical damping are minimized by shock-fitting and conformal mapping of body points.

Mcgary, J.↗

Surface and implantation effects on p-n junctions

The contribution of the graded region of implanted p-n junctions is analyzed using an exponential profile. Though previously neglected, it was recently shown that this contribution to the saturation current of HgCdTe diodes is significant. Assuming a dominant Auger recombination, an analytical solution to the continuity equation is obtained. An expression for the current generation by the graded region is presented for both ohmic and reflecting boundary conditions. A revised condition for a wide region is derived. When the region is narrow, the current differs drastically from that of the zero-gradient case. The effects of the junction depth and the substrate and surface concentrations on the current are investigated. It is shown that the reverse current does not saturate.

Schacham, Samuel E.↗

Modeling of flow systems for implementation under KATE

The modeling of flow systems is a task currently being investigated at Kennedy Space Center in parallel with the development of the KATE artificial intelligence system used for monitoring diagnosis and control. Various aspects of the modeling issues are focussed on with particular emphasis on a water system scheduled for demonstration within the KATE environment in September of this year. LISP procedures were written to solve the continuity equations for three internal pressure nodes using Newton's method for simultaneous nonlinear equations.

Whitlow, Jonathan E.↗

A simple, analytical, axisymmetric microburst model for downdraft estimation

A simple analytical microburst model was developed for use in estimating vertical winds from horizontal wind measurements. It is an axisymmetric, steady state model that uses shaping functions to satisfy the mass continuity equation and simulate boundary layer effects. The model is defined through four model variables: the radius and altitude of the maximum horizontal wind, a shaping function variable, and a scale factor. The model closely agrees with a high fidelity analytical model and measured data, particularily in the radial direction and at lower altitudes. At higher altitudes, the model tends to overestimate the wind magnitude relative to the measured data.

Vicroy, Dan D.↗

Plasma density enhancement at the Comet Halley diamagnetic cavity boundary

A one-dimensional single-fluid MHD model was used by Cravens (1989) to predict the existence of a narrow (50 km) layer of enhanced plasma density at the boundary of the Comet Halley diamagnetic cavity. The existence of such a layer was confirmed by measurements made by the Giotto ion mass spectrometer (Goldstein et al., 1989). Here, the time-dependent coupled continuity equations are solved for several species including H3O(+), H2O(+), OH(+), O(+), NH4(+), NH3(+), CH4(+), and CH3(+). For several ion species, factors affecting the magnitude of this enhancement and its relationship to the thickness of the transition layer are investigated. For example, ion species with short chemical lifetimes are shown to have smaller density enhancements than species with longer lifetimes. The cometocentric distance of the cavity boundary also strongly affects the magnitude of the density enhancement; the enhancement increases with increasing distance.

Keller, C. N.↗

Vertical velocities and momentum fluxes derived from wind measurements in the dusk auroral oval

Results are presented on a chemical release experiment carried out on March 20, 1985 in connection with NASA's Cooperative Observations of Polar Electrodynamics I campaign. Simultaneous neutral wind measurements in E region were carried out at three separate locations over a triangular region with sides of about 150 km in the dusk auroral oval, and the three wind profiles were used to calculate values of divergence and vorticity over the area. The vertical velocity over the hight range was calculated using the mass continuity equation, and the instantaneous vertical momentum fluxes in the E region were derived using the combination of horizontal-wind measurements and calculated vertical velocities. Results show that there is strong coupling between layers in the E region and that the momentum-flux-induced accelerations are at least comparable in magnitude to the pressure gradient and Coriolis force accelerations.

Larsen, M. F.↗

Analytical methods for the development of Reynolds-stress closures in turbulence

The derivation of Reynolds-stress models for viscous incompressible turbulent flow on the basis of the Navier-Stokes and continuity equations is explored in an analytical review. The formulation of the basic equations is outlined, and particular attention is given to zero-equation and one-equation models based on eddy viscosity, two-equation (k-l, k-epsilon, and k-omega) models, and second-order closure models. The superior performance of the second-order models is demonstrated by numerical results for (1) the return to isotropy from anisotropic homogeneous turbulence, (2) homogeneous turbulent shear flow in a rotating frame, (3) and fully developed inhomogeneous turbulent channel flow in a rotating frame. The inherent limitations of the Reynolds-stress approach and areas for further improvement are discussed.

Speziale, Charles G.↗

Stochastic estimation of coherent structures in turbulent boundary layers

The stochastic estimation technique is described as well as its extension to two-point and space-time conditional estimates. Based on the experimental data of Guezennec (1985), a comparison is made between the stochastic estimates and the conventional ensemble averaging results. As only the normal velocity and streamwise components were measured at each point on the sampling grid, the assumption of statistical symmetry in the spanwise direction about the centerline was used and the spanwise velocity component was reconstructed by integrating the continuity equation from the centerline outward.

Guezennec, Y. G.↗

The turbomachine blading design using S2-S1 approach

The boundary conditions corresponding to the design problem when the blades being simulated by the bound vorticity distribution are presented. The 3D flow is analyzed by the two steps S2 - S1 approach. In the first step, the number of blades is supposed to be infinite, the vortex distribution is transformed into an axisymmetric one, so that the flow field can be analyzed in a meridional plane. The thickness distribution of the blade producing the flow channel striction is taken into account by the modification of metric tensor in the continuity equation. Using the meridional stream function to define the flow field, the mass conservation is satisfied automatically. The governing equation is deduced from the relation between the azimuthal component of the vorticity and the meridional velocity. The value of the azimuthal component of the vorticity is provided by the hub to shroud equilibrium condition. This step leads to the determination of the axisymmetric stream sheets as well as the approximate camber surface of the blade. In the second step, the finite number of blades is taken into account, the inverse problem corresponding to the blade to blade flow confined in each stream sheet is analyzed. The momentum equation implies that the free vortex of the absolute velocity must be tangential to the stream sheet. The governing equation for the blade to blade flow stream function is deduced from this condition. At the beginning, the upper and the lower surfaces of the blades are created from the camber surface obtained from the first step with the assigned thickness distribution. The bound vorticity distribution and the penetrating flux conservation applied on the presumed blade surface constitute the boundary conditions of the inverse problem. The detection of this flux leads to the rectification of the geometry of the blades.

Luu, T. S.↗

Viscous-shock-layer analysis of hypersonic flows over long slender vehicles

An efficient and accurate method for solving the viscous shock layer equations for hypersonic flows over long slender bodies is presented. The two first order equations, continuity and normal momentum, are solved simultaneously as a coupled set. The flow conditions included are from high Reynolds numbers at low altitudes to low Reynolds numbers at high altitudes. For high Reynolds number flows, both chemical nonequilibrium and perfect gas cases are analyzed with surface catalytic effects and different turbulence models, respectively. At low Reynolds number flow conditions, corrected slip models are implemented with perfect gas case. Detailed comparisons are included with other predictions and experimental data.

Lee, Kam-Pui↗

Users manual for updated computer code for axial-flow compressor conceptual design

An existing computer code that determines the flow path for an axial-flow compressor either for a given number of stages or for a given overall pressure ratio was modified for use in air-breathing engine conceptual design studies. This code uses a rapid approximate design methodology that is based on isentropic simple radial equilibrium. Calculations are performed at constant-span-fraction locations from tip to hub. Energy addition per stage is controlled by specifying the maximum allowable values for several aerodynamic design parameters. New modeling was introduced to the code to overcome perceived limitations. Specific changes included variable rather than constant tip radius, flow path inclination added to the continuity equation, input of mass flow rate directly rather than indirectly as inlet axial velocity, solution for the exact value of overall pressure ratio rather than for any value that met or exceeded it, and internal computation of efficiency rather than the use of input values. The modified code was shown to be capable of computing efficiencies that are compatible with those of five multistage compressors and one fan that were tested experimentally. This report serves as a users manual for the revised code, Compressor Spanline Analysis (CSPAN). The modeling modifications, including two internal loss correlations, are presented. Program input and output are described. A sample case for a multistage compressor is included.

Glassman, Arthur J.↗

Pressure oscillation in the leakage annulus between a shrouded impeller and its housing due to impeller-discharge-pressure disturbances

The perturbed flow in the leakage path between a shrouded-pump impeller and its housing is analyzed using experiences with the Space Shuttle Main Engine (SSME), high pressure fuel turbopump (HPFTP) wearing-ring seals. Analysis is based on a bulk-flow model which consists of the path-momentum, circumferential momentum, and continuity equations. The pressure oscillations in the leakage annulus are driven by a circumferential variation of the impeller discharge pressure. It is shown that the occurrence and nature of the pressure oscillations depend on the tangential-velocity ratio of the fluid entering the seal, the order of the Fourier coefficient, the closeness of the precessional frequency of the rotating pressure field to the first natural frequency of the fluid annulus, and the clearance of the wearing-ring seal. The results obtained may explain the internal melting observed on SSME HPFTP seal parts.

Childs, D. W.↗

A modeling study of the time-averaged electric currents in the vicinity of isolated thunderstorms

A thorough examination of the results of a time-dependent computer model of a dipole thunderstorm revealed that there are numerous similarities between the time-averaged electrical properties and the steady-state properties of an active thunderstorm. Thus, the electrical behavior of the atmosphere in the vicinity of a thunderstorm can be determined with a formulation similar to what was first described by Holzer and Saxon (1952). From the Maxwell continuity equation of electric current, a simple analytical equation was derived that expresses a thunderstorm's average current contribution to the global electric circuit in terms of the generator current within the thundercloud, the intracloud lightning current, the cloud-to-ground lightning current, the altitudes of the charge centers, and the conductivity profile of the atmosphere. This equation was found to be nearly as accurate as the more computationally expensive numerical model, even when it is applied to a thunderstorm with a reduced conductivity thundercloud, a time-varying generator current, a varying flash rate, and a changing lightning mix.

Driscoll, Kevin T.↗

Modeling ice streams: Derived quantities

The model addressed is a finite-element, map-plane, time-dependent, column-averaged continuity equation solver. The key to the fitting process involves the balance between ice motion dominated by flow in internal layers, and ice motion dominated by sliding at the bed. The fitting process involves an iterative process carried out in the time domain. Beginning with the portion of the ice sheet being modeled identical to the present ice sheet with uniform flow, sliding, and fraction specified at nominal values, the model monitors each nodal point surface elevation. As the calculated surface elevation deviates from the present surface, a correction proportional to the difference is applied to selected parameter sets. This correction is in a sense that would tend to improve the fit at the particular nodal point. A calculated surface elevation that was higher than the present surface would result in an increased fraction, which would tend to lower the calculated surface (if the flow or sliding constant were being used as the fitting parameter, they would be lowered to improve the fit). This process is allowed to proceed as long as is necessary for the situation to stabilize. Typically, this takes tens of thousands of model years, but the rate is dependent on other external forcings such as the accumulation rate. The primary result is that while a typical sample of ice streams from around Antarctica can be fitted quite reasonably using only the fraction of the velocity due to sliding, a different mechanism seems to be in play along the Siple Coast, where reduced sliding constants are required to attain a reasonable fit. Flow is more strongly channelized in this region, and velocities are, in general, higher than are observed in other regions. It is unlikely that the mechanism that controls the ice movement along the Siple Coast is exactly similar to the mechanisms in the other ice streams. The concept of deformable sediments and their contribution to the fast flow along the Siple Coast may have limited applicability to other Antarctic ice streams.

Fastook, James↗