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 523 records · Page 29

Ionospheric models for Venus and Mars

The ways in which different models employ the multifluid continuity, momentum, and energy equations to characterize the ionospheres of Mars and Venus on the basis of in situ neutral and ion compositions are discussed. The Venus ionosphere is permeated by large-scale magnetic fields when the solar wind dynamic pressure is high, but is free from such fields under conditions of low solar wind dynamic pressure. Large-scale magnetic fields are nearly always present in the Mars ionosphere. Plasma moves downward on both planets for magnetized conditions, and magnetic flux is carried from the magnetosheath down into the ionosphere. Ohmic dissipation of the currents responsible for the magnetic field takes place deep in the ionospheres of both planets.

Cravens, T. E.↗

Equations, tables, and charts for compressible flow

This report, which is a revision and extension of NACA-TN-1428, presents a compilation of equations, tables, and charts useful in the analysis of high-speed flow of a compressible fluid. The equations provide relations for continuous one-dimensional flow, normal and oblique shock waves, and Prandtl-Meyer expansions for both perfect and imperfect gases. The tables present useful dimensionless ratios for continuous one-dimensional flow and for normal shock waves as functions of Mach number for air considered as a perfect gas. One series of charts presents the characteristics of the flow of air (considered a perfect gas) for oblique shock waves and for cones in a supersonic air stream. A second series shows the effects of caloric imperfections on continuous one-dimensional flow and on the flow through normal and oblique shock waves. (author)

Source record↗

Transport phenomena in the micropores of plug-type phase separators

This study numerically investigates the transport phenomena within and across a porous-plug phase separator. The effect of temperature differential across a single pore and of the sidewall boundary conditions, i.e., isothermal or linear thermal gradient, are presented and discussed. The effects are quantified in terms of the evaporation mass flux across the boundary and the mean surface temperature. A two-dimensional finite element model is used to solve the continuity, momentum, and energy equations for the liquid. Temperature differentials across the pore interface of 1.0, and 1.5 K are examined and their effect on evaporation flux and mean surface temperature is shown. For isothermal side boundary conditions, the evaporation flux across the pore is directly proportional and linear with Delta T. For the case of an imposed linear thermal gradient on the side boundaries, Biot numbers of 0.0, 0.15, and 0.5 are examined. The most significant effect of Biot number is to lower the overall surface temperature and evaporation flux.

Fazah, M. M.↗

Two-stream modeling of plasmaspheric refilling

Plasmaspheric refilling on an L = 4 flux tube was studied by using a time-dependent, hydrodynamic plasmaspheric flow model in which the ion streams from the two hemispheres are treated as distinct fluids. In the model the continuity, momentum, and energy equations of a two-ion (O(+) and H(+)), quasi-neutral, currentless plasma are solved along a closed geomagnetic field line; diffusive equilibrium is not assumed. collisions between all stream pairs and with neutral species are included. The model includes a corotating, tilted dipole magnetic field and neutral winds. Ionospheric sources and sinks are accounted for in a self-consistent manner. Electrons are assumed to be heated by photoelectrons. The model flux tube extends from a 200-km altitude in one hemisphere to a 200-km altitude in the other hemisphere. Initially, the upwelling streams pass through each other practically unimpeded. When the streams approach the boundary in the conjugate ionosphere, a shock develops there, which moves upward and dissipates slowly; at about the same time a reverse shock develops in the hemisphere of origin, which moves upward. After about 1 hour, large shocks develop in each stream near the equator; these shocks move toward the equator and downward after crossing the equator. However, these shocks are probably artificial, because counterstreaming flows occur in each H(+) fluid, which the model can only handle by creating shocks.

Guiter, S. M.↗

Plasmasphere Modeling with Ring Current Heating

Coulomb collisions between ring current ions and the thermal plasma in the plasmasphere will heat the plasmaspheric electrons and ions. During a storm such heating would lead to significant changes in the temperature and density of the thermal plasma. This was modeled using a time- dependent, one-stream hydrodynamic model for plasmaspheric flows, in which the model flux tube is connected to the ionosphere. The model simultaneously solves the coupled continuity, momentum, and energy equations of a two-ion (H(+) and O(+) quasineutral, currentless plasma. Heating rates due to collisions with ring current ions were calculated along the field line using a kinetic ring current model. First, diurnally reproducible results were found assuming only photoelectron heating of the thermal electrons. Then results were found with heating of the H(+) ions by the ring current during the recovery phase of a magnetic storm.

Guiter, S. M.↗

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗

Simulation of Gas Mixtures in RF Discharges for Nitride Deposition

Tungsten and titanium nitride films have long been grown by chemical vapor deposition (CVD) methods. However, there has been recent interest in low temperature growth using plasma enhanced CVD. For the present work, we focus on the radio frequency (BE) discharge characteristics of gas mixtures used in nitride deposition (for example, WF6 and ammonia). Because the radial variations for a standard 200 mm, parallel plate reactor are limited to a small zone near the edges of the electrodes, a 1-D (one-dimensional) analysis is considered. This model consists of a self-consistent, 3-D (three-dimensional) moment fluid simulation that solves the continuity, momentum, and energy equations for neutral and charged species. The results in terms of plasma structure, radical concentrations, and local deposition rate will be presented. We will also compare the 1-D results with those obtained from a 2-D hybrid plasma equipment model (HPEM) developed at the University of Illinois.

Hwang, Helen H.↗

Noise and Spreading of a Subsonic Coannular Jet - Comparison with Single Equivalent Jet

The issue of scaling of noise as well as spreading of subsonic coannular jets is revisited. Far-field noise and centerline Pitot-static pressure surveys are conducted with concentric, circular nozzles having an outer-to-inner diameter ratio of 1.42. Both the inner nozzle and the outer annular passage are convergent. Outer-to-inner Mach number ratio (R) is varied over a large range from 0 to approximately 10. Results are examined on the basis of single equivalent jet parameters calculated by satisfying continuity, momentum and energy equations. The results confirm that coannular jets with normal velocity profiles are noisier than the single equivalent jet. Jets with "inverted" velocity profiles are also found to be noisier except in a narrow R-range of 1-1.5. In the latter range, contrasting the inference in previous studies of IVP jets, the present data do not exhibit a clear noise reduction. When normalized with equivalent jet parameters the asymptotic Mach number decay rate, as well as potential core length, are found to be comparable to those of a single jet. However, an abrupt shift in the virtual origin is noted across R=1.

Zaman, K. B. M. Q.↗

Small Scale Mass Flow Plug Calibration

A simple control volume model has been developed to calculate the discharge coefficient through a mass flow plug (MFP) and validated with a calibration experiment. The maximum error of the model in the operating region of the MFP is 0.54%. The model uses the MFP geometry and operating pressure and temperature to couple continuity, momentum, energy, an equation of state, and wall shear. Effects of boundary layer growth and the reduction in cross-sectional flow area are calculated using an in- integral method. A CFD calibration is shown to be of lower accuracy with a maximum error of 1.35%, and slower by a factor of 100. Effects of total pressure distortion are taken into account in the experiment. Distortion creates a loss in flow rate and can be characterized by two different distortion descriptors.

mass-flow measurements↗

Numerical Calibration of Mass Flow Plug for Inlet Testing

A simple control volume model has been developed to calculate the discharge coefficient through a mass flow plug (MFP) and validated with a calibration experiment. The maximum error of the model within the operating region of the MFP is 0.54%. The control volume analysis developed work is comprised of a sequence of flow calculations through the MFP. The model uses the MFP geometry and operating pressure and temperature to couple continuity, momentum, energy, an equation of state, and wall shear. The discharge coefficient calculation also includes the effects of boundary layer growth, including the reduction in cross-sectional flow area as characterized by the boundary layer displacement thickness. The last calculation in the sequence uses an integral method to calculate the growth of the boundary layer, from which the displacement thickness is then determined. The result of these successive calculations is an accurate one-dimension model of the velocity, pressure, and temperature through the MFP. For comparison, a computational fluid dynamic (CFD) calibration is shown, which when compared to the presented numerical model, had a lower accuracy with a maximum error of 1.35% in addition to being slower by a factor of 100."

Flow Measurement↗

A Higher-Order Space-Time Finite-Element Method for Moving-Body and Fluid-Structure Interaction Problems

We present a high-order finite-element method for moving body and fluid/structure interaction problems. Our solution strategy is based on a space-time discontinuous Galerkin (DG) spectral-element discretization which extends to arbitrary order of accuracy. The space-time DG discretization is a natural choice for moving body and fluid-structure interaction problems as moving surfaces are incorporated simply by considering curved space-time elements whose space-time faces align with the moving body. We present a discontinuous-Galerkin in time discretization for six-degree of motion modeling of rigid bodies, and a continuous-Galerkin discretization for equations of linear elasticity to generate curved space-time meshes. Numerical results for several simple 2D test cases are presented in order to verify the implementation of the different models. Finally we present a preliminary dynamic simulation of a parachute.

High-Order↗

Translational control of a graphically simulated robot arm by kinematic rate equations that overcome elbow joint singularity

An operator commands a robot hand to move in a certain direction relative to its own axis system by specifying a velocity in that direction. This velocity command is then resolved into individual joint rotational velocities in the robot arm to effect the motion. However, the usual resolved-rate equations become singular when the robot arm is straightened. To overcome this elbow joint singularity, equations were developed which allow continued translational control of the robot hand even though the robot arm is (or is nearly) fully extended. A feature of the equations near full arm extension is that an operator simply extends and retracts the robot arm to reverse the direction of the elbow bend (difficult maneuver for the usual resolved-rate equations). Results show successful movement of a graphically simulated robot arm.

Barker, L. K.↗

Implementation and Testing of Turbulence Models for the F18-HARV Simulation

This report presents three methods of implementing the Dryden power spectral density model for atmospheric turbulence. Included are the equations which define the three methods and computer source code written in Advanced Continuous Simulation Language to implement the equations. Time-history plots and sample statistics of simulated turbulence results from executing the code in a test program are also presented. Power spectral densities were computed for sample sequences of turbulence and are plotted for comparison with the Dryden spectra. The three model implementations were installed in a nonlinear six-degree-of-freedom simulation of the High Alpha Research Vehicle airplane. Aircraft simulation responses to turbulence generated with the three implementations are presented as plots.

Yeager, Jessie C.↗

A new flux-conserving numerical scheme for the steady, incompressible Navier-Stokes equations

This paper is concerned with the continued development of a new numerical method, the space-time solution element (STS) method, for solving conservation laws. The present work focuses on the two-dimensional, steady, incompressible Navier-Stokes equations. Using first an integral approach, and then a differential approach, the discrete flux conservation equations presented in a recent paper are rederived. Here a simpler method for determining the flux expressions at cell interfaces is given; a systematic and rigorous derivation of the conditions used to simulate the differential form of the governing conservation law(s) is provided; necessary and sufficient conditions for a discrete approximation to satisfy a conservation law in E2 are derived; and an estimate of the local truncation error is given. A specific scheme is then constructed for the solution of the thin airfoil boundary layer problem. Numerical results are presented which demonstrate the ability of the scheme to accurately resolve the developing boundary layer and wake regions using grids which are much coarser than those employed by other numerical methods. It is shown that ten cells in the cross-stream direction are sufficient to accurately resolve the developing airfoil boundary layer.

Scott, James R.↗

Simple spline-function equations for fracture mechanics calculations

The paper presents simple spline-function equations for fracture mechanics calculations. A spline function is a sequence of piecewise polynomials of degree n greater than 1 whose coefficients are such that the function and its first n-1 derivatives are continuous. Second-degree spline equations are presented for the compact, three point bend, and crack-line wedge-loaded specimens. Some expressions can be used directly, so that for a cyclic crack propagation test using a compact specimen, the equation given allows the cracklength to be calculated from the slope of the load-displacement curve. For an R-curve test, equations allow the crack length and stress intensity factor to be calculated from the displacement and the displacement ratio.

Orange, T. W.↗

Electron temperature anisotropy in the polar wind

The steady state flow of a fully ionized H(+)-O(+)-electron plasma along geomagnetic field lines in the high-latitude topside ionosphere is investigated theoretically, with emphasis on the electron temperature anisotropy and heat flow in the polar wind. The 13-moment system of transport equations developed by Schunk (1975, 1977), which contains a continuity, momentum, internal energy, stress tensor and heat flow equation for each species, is employed to represent the electrons, with a simplified set of transport equations used for the ions, and the transport equations are solved at altitudes from 1500 to 12,000 km for both subsonic and supersonic H(+) outflows. For subsonic H(+) outflows, the electron gas is found to remain collision dominated to high altitudes. For supersonic H(+) outflows, the electron gas is also found to be collision dominated at altitudes below 2500 km, with an anisotropy in electron temperature distribution developing such that the temperature of the perpendicular electrons is greater than that of the parallel electrons at higher altitudes. In addition, the magnitude of downward electron heat flux at 1500 km is shown to have a dramatic effect on the individual parallel and perpendicular electron temperatures.

Schunk, R. W.↗

Second-order discrete Kalman filtering equations for control-structure interaction simulations

A general form for the first-order representation of the continuous, second-order linear structural dynamics equations is introduced in order to derive a corresponding form of first-order Kalman filtering equations (KFE). Time integration of the resulting first-order KFE is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete KFE involving only symmetric, N x N solution matrix.

Park, K. C.↗

A Localized Tau Method PDE Solver

In this paper we present a new form of the collocation method that allows one to find very accurate solutions to time marching problems without the unwelcome appearance of Gibb's phenomenon oscillations. The basic method is applicable to any partial differential equation whose solution is a continuous, albeit possibly rapidly varying function. Discontinuous functions are dealt with by replacing the function in a small neighborhood of the discontinuity with a spline that smoothly connects the function segments on either side of the discontinuity. This will be demonstrated when the solution to the inviscid Burgers equation is discussed.

Cottam, Russell↗