Conductive heating of the solar wind. III - Breakdown of the Navier-Stokes equations, 4 June - 4 December 1965
Navier-Stokes equation in predicting solar wind flow pattern
SEARCH · Search NASA
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.
Navier-Stokes equation in predicting solar wind flow pattern
Explore the source record for details and available documents.
The profiles and thicknesses of normal shock waves in argon at Mach numbers of 1.335, 1.454, 1.576, and 1-713 were determined experimentally by means of a free-molecule probe whose equilibrium temperature is related by kinetic theory to the local flow properties and their gradients. Comparisons were made between the experimental shock profiles and the theoretical profiles calculated from the Navier-Stokes equations, the Grad 13-moment equations, and the Burnett equations. New, very accurate numerical integrations of the Burnett equations were obtained for this purpose with results quite different from those found by Zoller, to whom the solution of this problem is frequently attributed. The experimental shock profiles were predicted with approximately equal success by the Navier-Stokes and Burnett theories, while the 13-moment method was definitely less satisfactory. A surprising feature of the theoretical results is the relatively small difference in predictions between the Navier-Stokes and Burnett theories in the present range of shock strengths and the contrastingly large difference between predictions of Burnett and the 13-moment theories. It is concluded that the Navier-Stokes equations are correct for weak shocks and that within the present shock strength range the Burnett equations make no improvement which merits the trouble of solving them. For shocks of noticeably greater strength, say with a shock Mach number of more than 2.5, it remains fundamentally doubtful that any of these theories can be correct.
Numerical solutions to Navier-Stokes equations which describe streaming solar atmosphere
Pressure distribution in hydrostatic bearing of multi-wells obtained by use of Navier-Stokes equations
Two-point, two-time correlation equations are obtained by considering the Navier-Stokes equations for two points in a fluid at two time. By neglecting the triple correlations in the equations, a solution is obtained for the final period of decay. The analysis is extended to earlier times by considering three points at three different times. The set of equations is made determinate by neglecting the quadruple correlations in comparison with the triple correlations. The diffusion of particles from a source in a decaying turbulent field is calculated approximately by assuming that the velocity fluctuations are small.
Fluid amplifier dynamic characteristics and Fortran program for numerical solution of time- dependent two-dimensional Navier-Stokes equation for viscous jet in arbitrary flow field
The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.
Following the law of stress adopted in the Navier-Stokes equations, the configuration of the viscous flow in planes at right angles to the axis of an infinite cylinder is found to be independent of the axial motion of the cylinder. In the limiting case of a yawed or swept wing of very high aspect ratio, certain boundary-layer and separation phenomena are thus determined independently by the crosswise component of velocity.
Following a law of stress adopted in the Navier-Stokes equations, the configuration of the viscous flow in planes at right angles to the axis of an infinite cylinder is found to be independent of the axial motion of the cylinder. In the limiting case of a yawed or swept wing of very high aspect ratio, certain boundary-layer and separation phenomena are thus determined independently by the crosswise component of velocity. It follows that the effect of sweepback is to increase the area of stable laminar flow and to decrease the lift coefficient at which flow separation occurs.
By an analysis of the Navier-Stokes equations it is shown that the aerodynamic coefficients of an infinite rectangular swept wing in an isothermal or adiabatic flow of a compressible gas can be determined from the aerodynamic coefficients of the unswept wing. When the flow is neither isothermal nor adiabatic, a three-dimensional boundary layer theory is developed and applied to the special case of a swept flat plate.
A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.
The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.
Equations for inhomogeneous turbulence constructed from Navier-Stokes and energy equations
Iterative solution to Krook-Boltzmann kinetic equation noting Navier-Stokes numerical solution, computation and analysis of distribution function within shock wave
Creeping flow solution of Leidenfrost phenomenon by use of Navier-Stokes, continuity, and energy equations
Shock wave profile computation using krook collision model equation in iteration scheme starting from navier-stokes solution
Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.