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.
Report discusses methods for numerical simulation of Navier-Stokes equations of viscous, incompressible flow, with emphasis on pseudocompressibility method. Advanced computational techniques reviewed contribute to improved designs by rapidly providing detailed information about conditions at all points in flow fields of proposed designs.
The purpose of this research is to construct accurate finite difference schemes for incompressible unsteady flow simulations such as LES (large-eddy simulation) or DNS (direct numerical simulation). In this report, conservation properties of the continuity, momentum, and kinetic energy equations for incompressible flow are specified as analytical requirements for a proper set of discretized equations. Existing finite difference schemes in staggered grid systems are checked for satisfaction of the requirements. Proper higher order accurate finite difference schemes in a staggered grid system are then proposed. Plane channel flow is simulated using the proposed fourth order accurate finite difference scheme and the results compared with those of the second order accurate Harlow and Welch algorithm.
In this paper we introduce a high order discontinuous Galerkin method for two dimensional incompressible flow in vorticity streamfunction formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method The streamfunction is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability The method is suitable for inviscid or high Reynolds number flows. Optimal error estimates are proven and verified by numerical experiments.
This paper presents viewgraphs of time-accurate simulations of unsteady incompressible flow with moving boundaries.
A new free surface flow model, RIPPLE, is summarized. RIPPLE obtains finite difference solutions for incompressible flow problems having strong surface tension forces at free surfaces of arbitrarily complex topology. The key innovation is the continuum surface force model which represents surface tension as a (strongly) localized volume force. Other features include a higher-order momentum advection model, a volume-of-fluid free surface treatment, and an efficient two-step projection solution method. RIPPLE's unique capabilities are illustrated with two example problems: low-gravity jet-induced tank flow, and the collision and coalescence of two cylindrical rods.
The problem of lifting surfaces and complex aircraft configurations in steady incompressible flow is considered. For lifting surfaces the problem is formulated in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential on the lifting surface. For complex configurations the problem is formulated in terms of an integral equation relating the potential to its normal derivative on the surface of the aircraft. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential (or the potential discontinuity) and its normal derivative to be constant within each element. The wake geometry is obtained by iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines.
An approximate indicial lift function associated with circulation was developed for tapered, swept wings in incompressible flow. The function is derived by representing the wings with a simple vortex system. The results from the derived equations compare well with the limited available results from more rigorous and complex methods. The equations, as derived, are not very convenient for calculating the dynamic response of aircraft, parameter extraction, or for determining frequency-response curves for wings. Therefore, an expression is developed to convert the indicial response function to an exponential form which is more convenient for these purposes.
Theoretical study of effect of pivot location on aerodynamic balance movement of variable sweep wings in incompressible flow
A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow.
Prandtl first order boundary layer equations for two dimensional laminar incompressible flow past circulation controlled circular lifting rotor
A new approach to panel methods is explored for two-dimensional steady incompressible flows. The method uses linear distributions of sources and vortices on straight-line panels, but satisfies boundary conditions on the actual body surface, at nodes that are also end points of the panels. The result is continuity in body-surface velocity distribution, without recourse to numerical quadrature for the velocity influence coefficients. The method is unusually sensitive to the distribution of the nodes. For example, it almost always fails to give acceptable results when the nodes are distributed randomly. However, the continuity of the velocity distribution makes possible a unique node redistribution scheme, which may be iterated to give accurate results reliably.
A new method for the numerical simulation of three-dimensional incompressible flows is described. The vortex-in-cell (VIC) method presented traces the motion of the vortex filaments in the velocity field which these filaments create. The velocity field is not calculated directly by the Biot-Savart law of interaction but by creating a mesh-record of the vorticity field, then integrating a Poisson's equation via the fast Fourier transform to generate a mesh-record of the velocity field. The computed scales of motion are assumed to be essentially inviscid. Viscous of subgrid-scale effects are incorporated into a filtering procedure in wave vector space. Results of tracing a periodic array of single vortex rings are compared with a Green's function calculation. The agreement is very good.
A method is presented for estimating the lift due to flap deflection at low angles of attack in incompressible flow. In this method provision is made for the use of incremental section-lift data for estimating the effectiveness of high-lift flaps. The method is applicable to swept wings of any aspect ratio or taper ratio. The present method differs from other current methods mainly in its ease of application and its more general application. Also included is a simplified method of estimating the lift-curve slope throughout the subsonic speed range.
A new method is presented for the solution of incompressible flow in generalized coordinates. This method is based on the substitution of the pressure weighted form of the momentum equations into the continuity equation. The algorithm is rigorously derived and a Fourier analysis is used to assess its suitability to act as an error smoother. Linear stability analysis results indicate that the performance of the new pressure substitution method (PSM) and the pressure correction method (PCM) is about the same at low Reynolds numbers, with no significant pressure gradient. At high Reynolds numbers the PSM shows much faster convergence. Likewise prediction of various flows indicate that the PSM has better accuracy for high Reynolds number flows with significant pressure gradients. Since most practical aerodynamic flows have significant pressure gradients, the PSM seems to be attractive for such flows. Solutions for both laminar and turbulent flow are compared with the experimental data. A two-equation low Reynolds number turbulence model is used to resolve the turbulent flowfield.
The calculation of pressures when the penalty-function approximation is used in finite-element solutions of laminar incompressible flows is addressed. A Poisson equation for the pressure is formulated that involves third derivatives of the velocity field. The second derivatives appearing in the weak formulation of the Poisson equation are calculated from the C0 velocity approximation using a least-squares method. The present scheme is shown to be efficient, free of spurious oscillations, and accurate. Examples of applications are given and compared with results obtained using mixed formulations.
A mathematical model is developed to analyze the suppression of rotating stall in an incompressible flow centrifugal compressor with a vaned diffuser, thereby addressing the important need for centrifugal compressor rotating stall and surge control. In this model, the precursor to to instability is a weak rotating potential velocity perturbation in the inlet flow field that eventually develops into a finite disturbance. To suppress the growth of this potential disturbance, a rotating control vortical velocity disturbance is introduced into the impeller inlet flow. The effectiveness of this control is analyzed by matching the perturbation pressure in the compressor inlet and exit flow fields with a model for the unsteady behavior of the compressor. To demonstrate instability control, this model is then used to predict the control effectiveness for centrifugal compressor geometries based on a low speed research centrifugal compressor. These results indicate that reductions of 10 to 15 percent in the mean inlet flow coefficient at instability are possible with control waveforms of half the magnitude of the total disturbance at the inlet.
A method was devised for estimating the incompressible-flow pressure distribution over compressor blade sections at design angle of attack. The theoretical incremental velocities due to camber and thickness of the section as an isolated airfoil are assumed proportional to the average passage velocity and are modified by empirically determined interference factors. Comparisons were made between estimated and test pressure distributions of NACA 65-series sections for typical conditions. Good agreement was obtained.
A method of computing span loads and the resulting rolling moments for sideslipping wings of arbitrary plan form in incompressible flow is derived. The method requires that the span load at zero sideslip be known for the wing under consideration. Because this information is available for a variety of wings, this requirement should not seriously restrict the application of the present method. The basic method derived herein requires a mechanical differentiation and integration to obtain the rolling moment for the general wing in sideslip. For wings having straight leading and trailing edges over each semispan, the rolling moment due to sideslip is given by a simple equation in terms of plan-form parameters and the lateral center of pressure of the lift due to angle of attack. Calculated span loads and rolling-moment parameters are compared with experimental values. The comparison indicates good agreement between calculations and available experimental data.