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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1949

Compression Shocks in Two-Dimensional Gas Flows

The following are arguments on the compression shocks in gas flow start with a simplified representation of the results of the study made by Th. Meyer as published in the Forschungsheft 62 of the VDI, supplemented by several amplifications for the application.In the treatment of compression shocks, the equation of energy, the equation of continuity, the momentum equation, the equation of state of the particular gas, as well as the condition Of the second law of thermodynamics that no decrease of entropy is possible in an isolated system, must be taken into consideration. The result is that, in those cases where the sudden change of state according to the second law of thermodynamics is possible, there always occurs a compression of the gas which is uniquely determined by the other conditions.

Busemann, A.↗
DOE OSTI2022

Scaling patch analysis of planar turbulent mixing layers

Proper scales for the mean flow and Reynolds shear stress in planar turbulent mixing layers are determined from a scaling patch analysis of the mean continuity and momentum equations. By seeking an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is determined as V ref =(dδ/dx)U ref , where dδ/dx is the growth rate of the mixing layer width and U ref =U h –U l is the difference between the velocity of the high speed stream U h and the velocity of the low speed stream U l . By seeking an admissible scaling for the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is determined as R uv,ref =U avg V ref =[$\frac{1}{2A_u}$$\frac{dδ}{dx}$]$U^{2}_{ref}$ where A u $_{=}^{def}$(U h –U l )/(U h +U l ) is the normalized velocity difference that emerges naturally in the admissible scaling of the mean momentum equation. Self-similar equations for the scaled mean transverse flow V* and Reynolds shear stress $R^{*}_{uv}$=R uv /R uv,ref are derived from the mean continuity and mean momentum equations. Finally, approximate equations for V* and $R^{*}_{uv}$ are developed and found to agree well with experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Ion loss processes in the stratosphere

Small ions consisting of aggregates of a few molecules determine th stratospheric electrical parameters such as conductivity, mobility, etc. The small ion density is controlled by the ionizing mechanisms for the production of ions and electrons and the loss processes for these charged particles. Ion production in stratosphere is chiefly due to galactic cosmic rays, and the loss processes due to recombination and attachment. Free electrons do not exist at stratospheric heights. The primary positive ion O2+ and the electrons are converted into complex clusters of positive and negative ions. The equilibrium ion density is governed by the equation of continuity for the production and loss of these ions. In the generalized equation of continuity, the gain or loss of ions due to transport is neglected. When aerosols or particulates are not present in the atmosphere, ions are loss due to ion-ion recombination. In the presence of aerosols, small ions can also be lost by attachment to the aerosols, and thus aerosols are likely to cause perturbations in the stratosphere electrification. A simplified model approach is adopted to study the effect of aerosols on the equilibrium ion density. The results of the analysis for the equatorial station Thumba (8.5 deg N) are presented.

Prasad, B. S. N.↗
NASA NTRS1980

On the numerical solution of time-dependent viscous incompressible fluid flows involving solid boundaries

An inherent numerical problem associated with the fully explicit pseudospectral numerical simulation of the incompressible Navier-Stokes equation for viscous flows with no-slip walls is described. A semi-implicit scheme which circumvents this numerical difficulty is presented. In this algorithm the equation of continuity rather than the Poisson equation for pressure is solved directly. Pseudospectral formulation of the channel flow problem using Fourier series and Chebyshev polynomials expansions is given for this scheme. An example demonstrating the applicability of the method is given.

Moin, P.↗
NASA NTRS1993

Three-dimensional numerical simulation of near-surface flows over the Martian north polar cap

Measurements made by Viking Lander VL-2 (48 N) have shown that the near-surface wind and temperature regime on Mars displays striking similarities to terrestrial counterparts. The diurnal radiative cycle is responsible for establishment of a well-defined thermal circulation in which downslope (Katabatic) flows prevail during the nighttime hours and weak upslope (anabatic) conditions prevail during the daytime. Previous work has indicated that the slope flows are much like those found on Earth, particularly the Katabatic winds, which show striking similarities to drainage flows observed over Antarctica. The low-level wind regime appears to be an important factor in the scouring of the martian landscape. The north polar cap shows evidence of eolian features such as dunes, frost streaks, and grooves from Viking imagery. The direction of the prevailing wind can in cases be inferred from the eolian features. We examine the thermally induced flows that result from the radiative heating and cooling of the martian north polar region using a comprehensive three-dimensional atmospheric mesoscale numerical model. The same model has been used previously for simulation of Antarctic Katabatic winds. The model equations are written in terrain-following coordinates to allow for irregular terrain; prognostic equations include the flux forms of the horizontal momentum equations, temperature, continuity. A surface energy budget equation is also incorporated in which the surface temperature is determined. Explicit parameterization of both terrestrial (longwave) and solar (shortwave) radiation is included. Turbulent transfer of heat and momentum in the martian atmosphere is assumed to follow the similarity expressions in the surface boundary layer on Earth. The terrain heights for the martian north polar region have been obtained from the U.S. Geological Survey map and digitized onto a 57x57 grid with a spacing of 75 km. The resulting terrain map is shown in Fig. 1. The vertical grid consists of 15 levels (delta equals 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere is necessary to capture details of the boundary layer flows. The lowest level corresponds to a height of approximately 20 m above the ground, the seconds level 100 m.

Parish, Thomas R.↗
NASA NTRS2022

Stability Analysis of Streaks Induced by Optimized Vortex Generators

Numerical computations are performed to investigate the potential for transition control in an axisymmetric boundary layer via fully realizable, streamwise stationary streaks induced by an azimuthally periodic array of surface mounted vortex generators (VGs). Previous work has shown that suitable streaks of this type can significantly reduce the growth of Mack’s second mode instabilities, but large streak amplitudes can make the flow susceptible to previously absent streak instabilities that can become the leading cause of transition. Here, we use the adjoint capabilities of the SU2 flow solver to optimize the VG shape to maximize the reduction in the growth of second-mode disturbances while also preventing the streak amplitudes from reaching large enough values to precipitate an earlier onset of transition via streak instabilities. The geometry and the freestream flow conditions are selected to match a relevant trajectory lo-cation from the HIFiRE-1 flight experiment. Results show that the optimized VGs can increase the mean streak amplitude by 117% with respect to a manually developed baseline design. The stability of this optimized basic state is analyzed via the plane-marching parabolized stability equations, predicting a fully laminar flow over the entire cone, or equivalently, yielding transition delay of 130% versus the 17% for the baseline VGs.

Boundary layer transition↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗
NASA NTRS1975

A prefiltering version of the Kalman filter with new numerical integration formulas for Riccati equations

A prefiltering version of the Kalman filter is derived for both discrete and continuous measurements. The derivation consists of determining a single discrete measurement that is equivalent to either a time segment of continuous measurements or a set of discrete measurements. This prefiltering version of the Kalman filter easily handles numerical problems associated with rapid transients and ill-conditioned Riccati matrices. Therefore, the derived technique for extrapolating the Riccati matrix from one time to the next constitutes a new set of integration formulas which alleviate ill-conditioning problems associated with continuous Riccati equations. Furthermore, since a time segment of continuous measurements is converted into a single discrete measurement, Potter's square root formulas can be used to update the state estimate and its error covariance matrix. Therefore, if having the state estimate and its error covariance matrix at discrete times is acceptable, the prefilter extends square root filtering with all its advantages, to continuous measurement problems.

Womble, M. E.↗
DOE OSTI2022

The radial action from probe amplitudes to all orders

We extract the relativistic classical radial action from scattering amplitudes, to all orders in perturbation theory, in the probe limit. Our sources include point charges and monopoles, as well as the Schwarzschild and pure-NUT gravitational backgrounds. A characteristic relativistic effect, that scattering trajectories may wind around these sources any number of times, can be recovered when all-order amplitudes are available. We show that the amplitude for scattering a probe off a pure NUT is given by the solution of a transcendental equation involving continued fractions, and explain how to solve this equation to any desired loop order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1965

Atmospheric oscillations

Motion, continuity, and adiabatic equations for upper atmospheric oscillation

ADIABATIC EQUATION↗
DOE OSTI2023

MFiX: Fractional-Step Method Implementation

A comprehensive, multiphase computational fluid dynamics (CFD) simulation solves several coupled transport equations including continuity, momentum, species, and energy. Chemical reactions further couple these equations through heats of reaction and rates of formation of products and rates of destruction of reactants. A fractional-step method separates changes attributed to chemical reactions from transport phenomena like convection and diffusion. When the governing equations are split into the transport and reacting components, efficient and independent methodologies can be exploited to solve the different systems. Specifically, discretization of field variable transport equations results in large, sparse matrices which are loosely coupled. These systems are solved in succession using iterative techniques that take advantage of the matrix structure. In contrast, chemical reactions tightly couple field variables locally within the domain (e.g., within a single computational cell) resulting in low dimensional but dense, nonlinear systems that are better solved using direct integration techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

A note on the effect of bending stiffness of stringers attached to a plate with a crack

An infinite stringer partially bonded to a plate through a layer of adhesive is considered. The stringer had bending as well as longitudinal stiffness. The effect of the stringer's bending rigidity on the stress intensity factor at the tip of the crack is illustrated. Shear stress distribution between the plate and the stringer and the stress intensity factors are obtained from the solution of a system of Fredholm integral equations. The continuity of displacements along the line of bond is represented by the equations.

Arin, K.↗