SEARCH · Engineering Papers
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.
Mean velocity and mean turbulent energy field methods of calculating boundary layer behavior based on numerical solutions of equations of motions
Computer program for numerical solution of unsteady Navier-Stokes equations and application to flow in rectangular cavity with moving wall
The cubically nonlinear Schroedinger equation is an important model of nonlinear phenomena in fluids and plasmas. Numerical solutions in a spatially periodic system commonly involve truncation to a finite number of Fourier modes. These solutions are found to be stochastic in the sense that the largest Liapunov exponent is positive. As the number of modes is increased, the size of this exponent appears to converge to zero, in agreement with the recent demonstration of the integrability of the spatially periodic case.
The structure and decay of a trailing vortex were analyzed during the numerical solutions of the full Navier-Stokes equations. Unsteady forms of the governing equations were recast in terms of circulation, vorticity, and stream function as dependent variables, and a second upwind finite difference scheme was used to integrate them with prescribed initial and boundary conditions. The boundary conditions at the outer edge and at the outflow section of the trailing vortex were considered. Different models of the flow were postulated, and solutions were obtained describing the development of the flow as integration proceeds in time. A parametric study was undertaken with a view to understanding the various phenomena that may possibly occur in the trailing vortex. Using the Hoffman and Joubert law of circulation at the inflow section, the results of this investigation were compared with experimental data for a Convair 990 wind model and a rectangular wing. With an exponentially decaying law of circulation at the inflow section and an adverse pressure gradient at the outer edge of the trailing vortex, solutions depict vortex bursting through the sudden expansion of the core and/or through the stagnation and consequent reversal of the flow on the axis. It was found that this bursting takes place at lower values of the swirl ratio as the Reynolds number increases.
Problems encoutered in investigations of high Reynolds number, incompressible flow are reviewed. A numerical solution computer program was modified to solve the stream function-vorticity form of the Navier-Stokes equations. Using a body fitted coordinate system with a U-shaped outer boundary, a simulation of incompressible flow at a Reynolds number of one million and a body angle of attack of zero was achieved.
The steady separated flow due to an obstruction in a two-dimensional channel is studied through the numerical solution of the Navier-Stokes equations using finite difference techniques. The characteristics of the separated flow behind the obstruction are obtained up to a Reynolds number of 1500. The flow patterns, vorticity contours and total pressure contours are presented. The effects of Reynolds number, height of obstruction, thickness of obstruction and the inlet velocity profile are investigated.
The equations of motion are developed for a perfectly flexible, inelastic tether with a satellite at its extremity. The tether is attached to a space vehicle in orbit. The tether is allowed to possess electrical conductivity. A numerical solution algorithm to provide the motion of the tether and satellite system is presented. The resulting differential equations can be solved by various existing standard numerical integration computer programs. The resulting differential equations allow the introduction of approximations that can lead to analytical, approximate general solutions. The differential equations allow more dynamical insight of the motion.
A six-stage low-storage Runge-Kutta time-marching method is presented and shown to be an efficient method for use with high-accuracy spatial difference operators for wave propagation problems. The accuracy of the method for inhomogeneous ordinary differential equations is demonstrated through numerical solutions of the linear convection equation with forced boundary conditions. Numerical experiments are presented simulating a sine wave and a Gaussian pulse propagating into and through the domain. For practical levels of mesh refinement corresponding to roughly ten points per wavelength, the six-stage Runge-Kutta method is more accurate than the popular fourth-order Runge-Kutta method. Further numerical experiments are presented which show that the numerical boundary scheme at an inflow boundary can be a significant source of error when high-accuracy spatial discretizations are used.
This paper presents a method of computing the laminar flow field of a body with stern-mounted jet. Problems involved are the numerical solutions of boundary layer and Navier-Stokes equations as well as the viscous-inviscid interaction. The adequacy of a simple viscous-inviscid interaction approach is investigated and a technique for connecting hybrid mesh systems without overlapping of meshes is presented. The numerical results indicate that while the external wake flow affects the development of the jet flow, the presence of a jet also influences the flow in the stern region of the body.
Analytical solutions for inviscid supersonic corner flows are virtually nonexistent due to the complexity of the interference geometry. In view of this, numerical solutions for compressive-compressive and expansive-compressive corner flows are obtained. The governing equations are written in strong conservation-law form and are solved iteratively in nonorthogonal conical coordinates by use of a second-order, shock-capturing, finite-difference technique. The computed wave structure and surface pressure distributions are compared with high Reynolds number (Re greater than 2,500,000 ft) experimental data. The results clearly show that the wave structure in the corner is dominated by the inviscid field.
Analytical solutions for inviscid supersonic corner flows are virtually nonexistent due to the complexity of the interference geometry. In view of this, numerical solutions for swept-compressive and swept-expansive corner flows are obtained. The governing equations are written in strong conservation-law form and are solved iteratively in nonorthogonal conical coordinates by use of a second-order, shock-capturing, finite-difference technique. The computed wave structure and surface pressure distributions are compared with high Reynolds number experimental data and show very good agreement. The results clearly show that supersonic corner flow at reasonably high Reynolds numbers including the effect of sweep is dominated by the inviscid field.
The technique of body-fitted coordinate systems is applied in numerical solutions of the complete time-dependent compressible and incompressible Navier-Stokes equations for laminar flow and to the time-dependent mean turbulent equations closed by modified Kolmogorov hypotheses for turbulent flow. Coordinate lines are automatically concentrated near to the bodies at higher Reynolds number so that accurate resolution of the large gradients near the solid boundaries is achieved. Two-dimensional bodies of arbitrary shapes are treated, the body contour(s) being simply input to the program. The complication of the body shape is thus removed from the problem.
Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.
Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem. Previously announced in STAR as N83-30109
The time dependent Navier-Stokes equations in mass averaged variables are solved for transonic flow over axisymmetric boattail plume simulator configurations. Numerical solution of these equations is accomplished with the unsplit explict finite difference algorithm of MacCormack. A grid subcycling procedure and computer code vectorization are used to improve computational efficiency. The two layer algebraic turbulence models of Cebeci-Smith and Baldwin-Lomax are employed for investigating turbulence closure. Two relaxation models based on these baseline models are also considered. Results in the form of surface pressure distribution for three different circular arc boattails at two free stream Mach numbers are compared with experimental data. The pressures in the recirculating flow region for all separated cases are poorly predicted with the baseline turbulence models. Significant improvements in the predictions are usually obtained by using the relaxation models.
A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.
A supersonic three-dimensional viscous forward-marching computer design code called PEPSIS is used to obtain a numerical solution of the three-dimensional problem of the interaction of a glancing sidewall oblique shock wave and a turbulent boundary layer. Very good results are obtained for a test case that was run to investigate the use of the wall-function boundary-condition approximation for a highly complex three-dimensional shock-boundary layer interaction. Two additional test cases (coarse mesh and medium mesh) are run to examine the question of near-wall resolution when no-slip boundary conditions are applied. A comparison with experimental data shows that the PEPSIS code gives excellent results in general and is practical for three-dimensional supersonic inlet calculations.
Laminar channel flow over a backward-facing step is investigated. The finite analytic (FA) method is used to obtain the numerical solution. The FA solutions predict the recirculation zone lengths and the recirculated mass flow rates for Reynolds numbers, Re, of 25, 50, 73, 125, 191 and 229 which correlate well with experimental measurements. The general flow patterns of the recirculation region flows for the Reynolds numbers considered in this study are similar to each other.