Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Incompressible Navier-Stokes”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

A finite element solution algorithm for the Navier-Stokes equations

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the steady-state kinematics and thermodynamics of a variable viscosity, compressible multiple-species fluid. For an incompressible fluid, the motion may be transient as well. The primitive dependent variables are replaced by a vorticity-streamfunction description valid in domains spanned by rectangular, cylindrical and spherical coordinate systems. Use of derived variables provides a uniformly elliptic partial differential equation description for the Navier-Stokes system, and for which the finite element algorithm is established. Explicit non-linearity is accepted by the theory, since no psuedo-variational principles are employed, and there is no requirement for either computational mesh or solution domain closure regularity. Boundary condition constraints on the normal flux and tangential distribution of all computational variables, as well as velocity, are routinely piecewise enforceable on domain closure segments arbitrarily oriented with respect to a global reference frame.

Baker, A. J.↗
NASA NTRS1976

A numerical study of flow phenomena within a vortex sink rate sensor

A numerical investigation of a viscous incompressible fluid with constant properties within a vortex sink rate sensor is undertaken. A computer program was developed, for obtaining the numerical solution of the Navier-Stokes equations, assuming laminar flow, having general prescribed inlet conditions and axisymmetric boundary conditions. Computational results for various viscous flows and assorted boundary conditions have been obtained. In the vortex chamber, the radial momentum is conserved, whereas the axial momentum is negligible. In the sink region, the radial and tangential velocities, as well as the tangential vorticities, have significantly high values. In the sink tube, the radial momentum is negligible. The radial momentum within the vortex chamber is converted to the axial momentum in the sink tube. A change in the flow rate had an appreciable effect within the sensor and particularly near the sink tube entrance. As the flow rate was increased, both the tangential velocity and tangential vorticity increased rapidly. At the higher flow rates, vortices were produced at the corner of the entrance section of the sink tube.

Patel, D. K.↗
NASA NTRS1975

Remarks on the theory of aerodynamic noise

The accuracy of approximations employed in the Lighthill theory in current aerodynamic noise research is critically evaluated. Based on the method of matched asymptotic expansions, the full Navier-Stokes equations are expanded for small Mach numbers. To the first order, the near-field is generally nonisentropic. The pressure field (pseudo-sound) is generated by the incompressible Reynolds stresses in the turbulent flow and the velocity, pressure perturbation, and their derivatives on the boundaries. In the far-field, the first-order pressure (acoustic) field satisfying a linear wave equation is obtained by matching with the pseudo-sound field. A uniformly valid solution for the pressure field in a stationary or uniformly moving medium is obtained. The solution shows that the generation of the first-order aerodynamic noise does not depend on the viscous, thermal, or entropy effects in the adiabatic flow nor on the shear stress on a smooth rigid boundary.

Pan, Y. S.↗
NASA NTRS1974

Steady simple shear flow past a circular cylinder at moderate Reynolds numbers - A numerical solution

The two-dimensional steady flow of an incompressible viscous fluid past a circular cylinder, placed symmetrically in a simple shear field, has been studied for both the stationary and the freely rotating case by solving numerically the Navier-Stokes equations for values of the Reynolds number R in the range from 0.047 to 70. At R = 0.047, the results obtained are in substantial agreement with the analytic small-R perturbation solution given by Robertson and Acrivos (1970). Inertia effects were found, however, to play a significant role even at R = 1, and hence the calculated flow pattern for R greater than or equal to 1 differs significantly from that of the creeping-flow solution. Specifically, for the freely rotating case, the region of closed streamlines decreases rapidly in extent with increasing R, two symmetrically placed wakes are formed on either side of the cylinder, and the dimensionless rotational speed of the freely suspended cylinder decreases as the reciprocal of the square root of R.

Kossack, C. A.↗
DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
NASA NTRS1977

A second-order closure analysis of turbulent diffusion flames

A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.

Varma, A. K.↗
NASA NTRS1974

Numerical solution of the time-dependent compressible Navier-Stokes equations in inlet regions

The results of a study to determine the effects of compressibility on the viscous flow through channels that have straight, parallel walls are presented. Two channel configurations are considered, the flow between two semi-infinite flat plates with uniform flow prescribed at the inlet plane and a cascade of semi-infinite flat plates with uniform flow introduced upstream. The flow field is modeled by using the time dependent, compressible Navier-Stokes equations. Time dependent solutions are obtained by using an explicit finite difference technique which advances the pressure on near field subsonic boundaries such that accurate steady state solutions are obtained. Steady state results at Reynolds number 20 and 150 are presented for Mach numbers between 0.09 and 0.36 and compared with the incompressible solutions of previous studies.

Olson, L. E.↗
NASA NTRS1975

Investigation A: Numerical solutions of driven vortices of binary fluid in a cylindrical geometry

The steady state laminar motion of a viscous, incompressible and binary fluid is studied for a rotating flow in a cylindrical geometry. The mathematical model employed is a cold flow simulation of the fluid mechanics of the light-bulb concept of the gaseous core nuclear engine. A numerical treatment is developed for the rotating flow which includes a description of the nuclear fuel addition. The problem is formulated with the complete Navier-Stokes equations in order to show the interaction between the fuel addition, the main flow, and the boundary layer flow in an accurate manner. The results presented show holdup of the nuclear fuel for the case of steady fuel addition.

Shin, Y. W.↗
NASA NTRS1975

Perturbation solution of the Navier-Stokes equations and its relation to the Lighthill-Curle solution of aerodynamic sound

The aerodynamic sound described by the Lighthill-Curle solution is reexamined using a method of matched asymptotic expansions. The governing Navier-Stokes equations written in nondimensional form are expanded for a small Mach number. First- and second-order solutions for the pressure field are obtained, and the singular nature of the expansion at large distances is indicated. The nearfield pressure is governed by the Poisson equation, whereas the farfield equations describe a linear wave system in a dissipative medium. The pseudosound is related to the incompressible Reynolds stresses associated with a solenoidal velocity field, the velocity, the pressure perturbation, and their derivatives on the boundaries. A uniformly valid first-order solution for the pressure is obtained. It is shown that viscosity, thermal conductivity, and entropy in the flow do not contribute to the first-order noise generation, while the viscous stress contributes to noise only from some boundaries. The application of the proposed perturbation method to a subsonically moving surface and a hot jet is discussed.

Pan, Y. S.↗
NASA NTRS1978

The computation of incompressible, three-dimensional flows in containers with extensible walls

The theoretical study of flow in biological vascular systems is made very difficult in connection with local and temporal changes of the cross section. Experimental investigations with the aid of model tests are, however, not enough for a solution of the problems, and numerical solutions are more and more employed for a description of such flows. A description is presented of a difference procedure for the solution of the complete Navier-Stokes equations with curvilinear coordinates for three-dimensional flows in containers. The integration of the Navier-Stokes equations for flows in containers with rigid and moving walls is considered.

Liu, N. S.↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗
NASA NTRS1974

Approach to the origin of turbulence on the basis of two-point kinetic theory

Equations for the fluctuation correlation in an incompressible shear flow are derived on the basis of kinetic theory, utilizing the two-point distribution function which obeys the BBGKY hierarchy equation truncated with the hypothesis of 'ternary' molecular chaos. The step from the molecular to the hydrodynamic description is accomplished by a moment expansion which is a two-point version of the thirteen-moment method, and which leads to a series of correlation equations, viz., the two-point counterparts of the continuity equation, the Navier-Stokes equation, etc. For almost parallel shearing flows the two-point equation is separable and reduces to two Orr-Sommerfeld equations with different physical implications.

Tsuge, S.↗
NASA NTRS1974

Solutions for laminar boundary layers with separation and reattachment

Numerical solutions of the laminar, incompressible boundary-layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified, and the pressure is deduced from the solution. A stream function/vorticity formulation of the boundary-layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite-difference equations, in which the flow direction is taken into account, are obtained by a successive column iteration scheme. Results are compared with Klineberg and Steger's (1974) separated boundary-layer calculations, and with Briley's (1971) solution of the Navier-Stokes equations for a separated region.

Carter, J. E.↗
NASA NTRS1975

Inverse solutions for laminar boundary-layer flows with separation and reattachment

Numerical solutions of the laminar, incompressible boundary layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified; the pressure is deduced from the solution. A vorticity-stream-function formulation of the boundary layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite difference equations, in which the flow direction is taken into account, are obtained by several global iteration schemes which are stable and have unconditional diagonal dominance. Results are compared with Klineberg and Steger's separated boundary layer calculations, and with Briley's solution of Navier-Stokes equations for a separated region. In addition, an approximate technique is presented in which the streamwise convection of vorticity is set equal to zero in the reversed flow region; such a technique results in a quick forward marching procedure for separated flows.

Carter, J. E.↗
NASA NTRS1952

Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System

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.

Mager, Arthur↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗