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

Simulation of three-dimensional liquid sloshing flows using a strongly implicit calculation procedure

A coupled strongly implicit solution strategy for unsteady three-dimensional free surface flows has been developed based on an artificial compressibility formulation for the incompressible Navier-Stokes equations. A pseudotime term has been used in the continuity equation to permit time accurate calculations to be achieved. The scheme appears capable of tracking the free surface reasonably accurately inside a partially-filled spherical container undergoing a general rotating motion characteristic of that experienced by a spin-stabilized satellite. Five different free surface calculations have been presented. Some of the results exhibit an interesting Reynolds number dependent oscillatory behavior which is believed to be physical although no experimental results appears to be available for verification to date.

Chen, K.-H.↗
NASA NTRS2000

Statistically Reliable 'Atomistic' Simulation of Sub 100 nm MOSFETs

A 3D 'atomistic' simulation technique to study random impurity induced threshold voltage lowering and fluctuations in sub 0. 1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level-Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen↗
NASA NTRS1998

Efficient 3D 'Atomistic' Simulation Technique for Studying of Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Decanano MOSFETs

A 3D 'atomistic' simulation technique to study random dopant induced threshold voltage lowering and fluctuations in sub 0.1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level. Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen↗
NASA NTRS1998

Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Sub-0.1 (micron)meter MOSFET's: A 3-D 'Atomistic' Simulation Study

A three-dimensional (3-D) "atomistic" simulation study of random dopant induced threshold voltage lowering and fluctuations in sub-0.1 microns MOSFET's is presented. For the first time a systematic analysis of random dopant effects down to an individual dopant level was carried out in 3-D on a scale sufficient to provide quantitative statistical predictions. Efficient algorithms based on a single multigrid solution of the Poisson equation followed by the solution of a simplified current continuity equation are used in the simulations. The effects of various MOSFET design parameters, including the channel length and width, oxide thickness and channel doping, on the threshold voltage lowering and fluctuations are studied using typical samples of 200 atomistically different MOSFET's. The atomistic results for the threshold voltage fluctuations were compared with two analytical models based on dopant number fluctuations. Although the analytical models predict the general trends in the threshold voltage fluctuations, they fail to describe quantitatively the magnitude of the fluctuations. The distribution of the atomistically calculated threshold voltage and its correlation with the number of dopants in the channel of the MOSFET's was analyzed based on a sample of 2500 microscopically different devices. The detailed analysis shows that the threshold voltage fluctuations are determined not only by the fluctuation in the dopant number, but also in the dopant position.

Asenov, Asen↗
NASA NTRS1981

The effect of microscopic turbulence on magnetosphere-ionosphere coupling

The effect of turbulence on the coupling of the magnetosphere and ionosphere has been investigated by including an effective collision frequency in the electron equation of motion. When this term is combined with the continuity equation, the ion equation of motion and Maxwell's equations, a dispersion relation for the kinetic Alfven wave including effective collisions is found. The wave-particle interaction leads to a strong damping of the wave. Inclusion of the effects of plasma sheet kinetics yields a scale size transverse to the magnetic field which corresponds to the size of visual auroral arcs.

Lysak, R. L.↗
NASA NTRS1986

A numerical model of magnetosphere-ionosphere coupling Preliminary results

A three-dimensional simulation model was developed to study magnetosphere-ionosphere coupling in the auroral region. One-fluid MHD equations are adopted to model the magnetosphere, and current density continuity equations were solved consistently to model the ionosphere. In the preliminary simulation runs described here, the electrodynamics of region 1 field-aligned currents were modeled. Initially, the electric field is taken to be in the magnetospheric equatorial plane. Alfven waves then propagate down to the ionosphere, accompanying a field-aligned current and exciting the electrostatic potential in the ionosphere by electrodynamic coupling. The results of the preliminary runs directly correspond with the fundamental characteristics of global magnetosphere-ionosphere coupling. These characteristics include an ionospheric electrostatic potential which varies in its development in time, depending upon the ratio of ionospheric resistance to magnetospheric impedance, the flowing of field-aligned currents into the ionosphere on the dawnside and out on the duskside, and the distribution along the geomagnetic field of the amplitude of the field-aligned current density in proportion to the field intensity.

Watanabe, K.↗
NASA NTRS1996

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗
NASA NTRS2000

Descending Layer Variability Over Arecibo

Descending layers of ionization over Arecibo exhibit very diverse behavior from night to night that does not appear to be strongly correlated to geomagnetic activity, solar forcing, or average semidiurnal tidal winds. On some nights, three or more distinct layers are observed to form near 170 km over timescales of approx. 2 hours. Rather than descending smoothly over periods of several hours, these layers stall, abruptly disappear, or even reverse direction in the midst of their descent. The time scales for their disappearance are examined and compared to loss rates arising from diffusion and recombination. Diffusion alone is found to be too slow to account for the observations, but recombination is fast enough provided that the convergent wind shear that forms the layer is relatively weak coincident with their disappearance. The continuity equation is solved in conjunction with a time sequence of radar profiles to estimate the vertical drift and horizontal neutral wind consistent with the observed behavior. The resultant wind field is northward, has an average speed of approx. 80 m/ s, and varies significantly near the altitude where the layers are observed. These inferred winds are consistent with the presence of the observed layers, and their magnitudes as obtained from the classical continuity and momentum equations are reasonable for this altitude range.

Earle, G. D.↗
NASA NTRS1982

A composite velocity procedure for the compressible Navier-Stokes equations

A new boundary-layer relaxation procedure is presented. In the spirit of the theory of matched asymptotic expansions, a multiplicative composite of the appropriate velocity representations for the inviscid and viscous regions is prescribed. The resulting equations are structured so that far from the surface of the body the momentum equations lead to the Bernoulli relation for the pressure, while the continuity equation reduces to the familiar compressible potential equation. Close to the body surface, the governing equations and solution techniques are characteristic of those describing interacting boundary-layers; although, the full Navier-Stokes equations are considered here. Laminar flow calculations for the subsonic flow over an axisymmetric boattail simulator geometry are presented for a variety of Reynolds and Mach numbers. A strongly implicit solution method is applied for the coupled velocity components.

Khosla, P. K.↗
NASA NTRS1987

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗
NASA NTRS1990

Prediction of vortical flows on wings using incompressible Navier-Stokes equations

Numerical solutions to the incompressible Navier-Stokes equations are in greater demand than ever before as the field of computational fluid dynamics (CFD) increases its impact as an engineering tool. An implicit upwind-relation finite-difference scheme is developed for solving unsteady incompressible Navier-Stokes equations incorporating an artificial time derivative of the pressure to the continuity equation. The essential features observed in experiments for complicated vortical flows past three round-edged low-aspect-ratio wings at high angles of attack are reproduced both qualitatively and quantitatively.

Hsu, C.-H.↗
NASA NTRS1989

Numerical simulation of thermocapillary bubble migration under microgravity for large Reynolds and Marangoni numbers

A numerical procedure in which the Navier-Stokes equations are discretized using tightly coupled discretizations of pressure derivatives and continuity equations is used here to extend the range of known terminal velocities of gaseous bubbles in liquids well beyond that in previous investigations. Computations performed for Reynolds numbers up to 2000 and Marangoni numbers up to 1000 show only a modest variation of the scaled bubble velocity between 0.16 and 0.5. The bubble velocity is influenced more by the Marangoni number than by the Reynolds number.

Balasubramaniam, R.↗
NASA NTRS1994

Numerical modeling of a cryogenic fluid within a fuel tank

The computational method developed to study the cryogenic fluid characteristics inside a fuel tank in a hypersonic aircraft is presented. The model simulates a rapid draining of the tank by modeling the ullage vapor and the cryogenic liquid with a moving interface. A mathematical transformation was developed and applied to the Navier-Stokes equations to account for the moving interface. The formulation of the numerical method is a transient hybrid explicit-implicit technique where the pressure term in the momentum equations is approximated to first order in time by combining the continuity equation with an ideal equation of state.

Greer, Donald S.↗
NASA NTRS1987

Force and moment rotordynamic coefficients for pump-impeller shroud surfaces

Governing equations of motion are derived for a bulk-flow model of the leakage path between an impeller shroud and a pump housing. The governing equations consist of a path-momentum, a circumferential - momentum, and a continuity equation. The fluid annulus between the impeller shroud and pump housing is assumed to be circumferentially symmetric when the impeller is centered; i.e., the clearance can vary along the pump axis but does not vary in the circumferential direction. A perturbation expansion of the governing equations in the eccentricity ratio yields a set of zeroth and first-order governing equations. The zeroth-order equations define the leaking rate and the circumferential and path velocity distributions and pressure distributions for a centered impeller position. The first-order equations define the perturbations in the velocity and pressure distributions due to either a radial-displacement perturbation or a tilt perturbation of the impeller. Integration of the perturbed pressure and shear-stress distribution acting on the rotor yields the reaction forces and moments acting on the impeller face.

Childs, Dara W.↗
NASA NTRS1993

A multigrid solver for the semiconductor equations

We present a multigrid solver for the exponential fitting method. The solver is applied to the current continuity equations of semiconductor device simulation in two dimensions. The exponential fitting method is based on a mixed finite element discretization using the lowest-order Raviart-Thomas triangular element. This discretization method yields a good approximation of front layers and guarantees current conservation. The corresponding stiffness matrix is an M-matrix. 'Standard' multigrid solvers, however, cannot be applied to the resulting system, as this is dominated by an unsymmetric part, which is due to the presence of strong convection in part of the domain. To overcome this difficulty, we explore the connection between Raviart-Thomas mixed methods and the nonconforming Crouzeix-Raviart finite element discretization. In this way we can construct nonstandard prolongation and restriction operators using easily computable weighted L(exp 2)-projections based on suitable quadrature rules and the upwind effects of the discretization. The resulting multigrid algorithm shows very good results, even for real-world problems and for locally refined grids.

Bachmann, Bernhard↗
NASA NTRS1982

'Scaling' analysis of the ice accretion process on aircraft surfaces

A comprehensive set of scaling parameters is developed for the ice accretion process by analyzing the energy equations of the dynamic freezing zone and the already frozen ice layer, the continuity equation associated with supercooled liquid droplets entering into and impacting within the dynamic freezing zone, and energy equation of the ice layer. No initial arbitrary judgments are made regarding the relative magnitudes of each of the terms. The method of intrinsic reference variables in employed in order to develop the appropriate scaling parameters and their relative significance in rime icing conditions in an orderly process, rather than utilizing empiricism. The significance of these parameters is examined and the parameters are combined with scaling criteria related to droplet trajectory similitude.

Keshock, E. G.↗