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.
The breakdown of a vortex flow in a tube is studied for a slightly diverging tube by means of a numerical solution of the Navier-Stokes equations for a quasi-stationary, incompressible, laminar flow. Numerical results are compared to experiments of Faler and Leibovich. The numerical data display features similar to the experimental results concerning breakdown conditions, location, and structure.
The measurements of part V (reference 1) of this series of reports, which concerned comparatively long ring profiles, are supplemented by measurements on shorter rings as they are used for shrouded propellers and cowlings of ring-shaped radiators. Mass-flow coefficients and profile drags are given. Furthermore, it has to be determined how far the potential theory describes the flow phenomenon with sufficient accuracy and whether the present theory for the calculation of thin annular profile yields useful profile forms and is suitable for determination of the mass flow for thick profiles.
A computational study of the separated flow through a 2-D asymmetric subsonic diffuser has been performed. The Wind Computational Fluid Dynamics code is used to predict the separation and reattachment behavior for an incompressible diffuser flow. The diffuser inlet flow is a two-dimensional, turbulent, and fully-developed channel flow with a Reynolds number of 20,000 based on the centerline velocity and the channel height. Wind solutions computed with the Menter SST, Chien k-epsilon, Spalart-Allmaras and Explicit Algebraic Reynolds Stress turbulence models are compared with experimentally measured velocity profiles and skin friction along the upper and lower walls. In addition to the turbulence model study, the effects of grid resolution and use of wall functions were investigated. The grid studies varied the number of grid points across the diffuser and varied the initial wall spacing from y(sup +) = 0.2 to 60. The wall function study assessed the applicability of wall functions for analysis of separated flow. The SST and Explicit Algebraic Stress models provide the best agreement with experimental data, and it is recommended wall functions should only be used with a high level of caution.
The theoretical and experimental internal flow characteristics of a 13.97-cm-diam inlet with centerbody retracted and extended are presented at STOL takeoff and approach operating conditions. The theoretical results were obtained from incompressible potential flow corrected for compressibility and boundary layer. Comparisons between theoretical internal surface static-pressure distributions and experimental data are presented for free-stream velocities of 0, 24, 32, and 45 m/sec for a range of inlet incidence angles from 0 to 50 deg. Surface static-pressure distributions are illustrated at circumferential locations of 0, 60, 120, and 180 deg. Surface Mach number distributions from the stagnation point to the diffuser exit are presented along with turbulent boundary-layer shape factors. In general, good agreement was found between the theoretical and experimental surface static pressure distributions.
The objective is to determine by analytical means how drag on a plane wall may be modified favorably using a minimal amount of flow information - preferably only information at the wall. What quantities should be measured? How should that information be assimilated in order to arrive at effective control? As a prototypical problem, incompressible, viscous flow, governed by the Navier-Stokes equations, past a plane wall at which the no-slip condition was modified was considered. The streamwise and spanwise velocity components are required to be zero, but the normal component is to be specified according to some control law. The challenge is to choose the wall-normal velocity component based on flow conditions at the wall so that the mean drag is as small as possible. There can be no net mass flux through the wall, and the total available control energy is constrained. A turbulent flow is highly unsteady and has detailed spatial structure. The mean drag on the wall is the integral over the wall of the local shear forces exerted by the fluid, which is then averaged in time; it is a 'macroscopic' property of the flow. It is not obvious how unsteady boundary control is to be applied in order to modify the mean flow most effectively, especially in view of the non- self-adjoint nature of the governing equations. An approximate analytical solution to the suboptimal scheme is pursued.
In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.
A numerical approach for computing unsteady compressible viscous flows was developed. This approach offers the capability of confining the region of computation to the viscous region of the flow. The viscous region is defined as the region where the vorticity is nonnegligible and the difference in dilatation between the potential flow and the real flow around the same geometry is also nonnegligible. The method was developed and tested. Also, an application of the procedure to the solution of the steady Navier-Stokes equations for incompressible internal flows is presented.
The computation of time-dependent flows has inspired a new, higher-dimensional formulation of the conservation equations of fluid dynamics in which time is treated as a fourth coordinate. The formulation is derived for a constant-density flow, and then extended to a variable-density flow by introducing a fifth, fictitious coordinate. This new coordinate can also act as a source coordinate, so that external source terms can be included. The analysis is carried out for both incompressible, stratified flow, and compressible equilibrium flow. The results are then extended to non-equilibrium and magnetohydrodynamic flows. Several applications of the new formulation to the computation of time-dependent flows are discussed.
A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.
In order to reduce the computational difficulty associated with a single grid (SG) solution procedure, the multigrid (MG) technique was identified as a useful means for improving the convergence rate of iterative methods. A full MG full approximation storage (FMG/FAS) algorithm is used to solve the incompressible recirculating flow problems in complex geometries. The algorithm is implemented in conjunction with a pressure correction staggered grid type of technique using the curvilinear coordinates. In order to show the performance of the method, two flow configurations, one a square cavity and the other a channel, are used as test problems. Comparisons are made between the iterations, equivalent work units, and CPU time. Besides showing that the MG method can yield substantial speed-up with wide variations in Reynolds number, grid distributions, and geometry, issues such as the convergence characteristics of different grid levels, the choice of convection schemes, and the effectiveness of the basic iteration smoothers are studied. An adaptive grid scheme is also combined with the MG procedure to explore the effects of grid resolution on the MG convergence rate as well as the numerical accuracy.
Results of four thin, unswept, flutter airfoils attached to two freely falling bodies are reported. Two airfoils fluttered at a Mach number of 0.85, a third airfoil fluttered at a Mach number of 1.03, and a fourth fluttered at a Mach number of 1.07. Results of calculations of flutter speed using incompressible and compressible air-force coefficients, including a Mach number of 1.0, are presented.
An attempt was made to model the so called leading edge vortex which forms over the leading edge of delta wings at high angles of attack. A simplified model was considered, namely that of a two-dimensional, inviscid, incompressible steady flow around a flat plate at an angle of attack with a stationary vortex detached on top, as well as a sink to simulate the strong spanwise flow. The results appear to agree qualitatively with experiments. A comparison was also made between the lift and the drag of this model and the corresponding results for two classical solutions: (1) that of totally attached flow over the plate with the Kutta condition satisfied at the trailing edge only; and (2) the Helmholtz solution of totally separated flow over the plate.
An attempt was made to model the so called leading edge vortex which forms over the leading edge of delta wings at high angles of attack. A simplified model was considered, namely that of a two-dimensional, inviscid, incompressible steady flow around a flat plate at an angle of attack with a stationary vortex detached on top, as well as a sink to simulate the strong spanwise flow. The results appear to agree qualitatively with experiments. A comparison was also made between the lift and the drag of this model and the corresponding results for two classical solutions: (1) that of totally attached flow over the plate with the Kutta condition satisfied at the trailing edge only: and (2) the Helmholtz solution of totally separated flow over the plate.
In order to study the effects of ion wind on the plane Poiseuille flow in a channel and Blasius flow on a flat plate, the equations of electrogasdynamics are solved numerically under the assumptions that the flow is incompressible, that the electric field is steady, and that the fluid velocity is negligible compared to ion drift velocity. The results show that ion wind strongly affects the skin-friction distribution for these flows.
Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.
The potential of realizable, stationary streaks undergoing non-modal growth to stabilize a hypersonic boundary-layer flow and, subsequently, delay the laminar-turbulent transition onset, is studied via numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach number of 5.30, freestream unit Reynolds number equal to 13.42 x 10(exp 6)/m, and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the test article. This paper investigates flow modifications induced by wall-mounted vortex generators (VGs), followed by an analysis of the modal instability of the perturbed, streaky boundary-layer flow. Results are presented both for a single array of VGs that was designed on the basis of optimal growth theory and for a VG configuration involving two separate arrays with opposite orientations that ware designed to provide staged control of flow instabilities while simultaneously reducing the amplification of streak instabilities resulting from the control devices. Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment, which did not include any control devices, correlated with an amplification factor of N = 14.7 for the planar Mack modes. If one assumes that the transition N -factor is not affected by the introduction of the VGs, then the control configurations based on a single array of VGs and two separate arrays would result in a transition delay of 17% and 40%, respectively. These findings suggest a passive flow control s to induce streaks that would delay transition in hypersonic boundary dominated by Mack-mode instabilities.
This paper describes the simulation of turbulent boundary layers by direct numerical solution of the three-dimensional, time-dependent Navier-Stokes equations, using a spectral method. The flow is incompressible, with Re sub delta = 1000. The equations are written in the self-similar coordinate system and periodic streamwise and spanwise boundary conditions are imposed. A family of nine 'equilibrium' boundary layers, from the strongly accelerated 'sink' flow to Stratford's separating boundary layer is treated. Good general agreement with experiments is observed. The effects of pressure gradients on the structures and statistics, both in the wall and wake regions are discussed.
An implicit flux-difference splitting scheme is used to compute incompressible laminar vortical flows about a thin round-edged double-delta wing with 80 deg and 60 deg leading-edged sweep for the strake and the main wing respectively. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. It is second-order accurate spatially by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. The interaction between the two primary vortices emanating from the apex and kink leading-edges is successfully simulated. The computed trajectories of vortical cores compare well with experimental data.