The development of a solution-adaptive 3D Navier-Stokes solver for turbomachinery
NASA Astrophysics Data System (ADS)
Dawes, W. N.
1991-06-01
This paper describes the early stages in the development of a solution-adaptive fully three-dimensional Navier-Stokes solver. The compressible, Navier-Stokes equations, closed with k-epsiton turbulence modeling, are discretized on an unstructured mesh formed from tetrahedral computational control volumes. At the mesh generation stage and at stages during the solution process itself, mesh refinement is carried out by flagging cells which satisfy particular critera. These criteria include geometric features such as proximity to wetted surfaces and features associated with the particular flowfield, such as fractional variation of a flow variable over cell faces. Solutions are presented for the highly three-dimensional flows associated with a truncated cylinder in a cross flow, a three-dimensional swept transonic bump, and the corner stall and secondary flow in a transonic compressor cascade.
Calculations of separated 3-D flows with a pressure-staggered Navier-Stokes equations solver
NASA Technical Reports Server (NTRS)
Kim, S.-W.
1991-01-01
A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.
NASA Technical Reports Server (NTRS)
Loh, Ching Y.; Himansu, Ananda; Hultgren, Lennart S.
2003-01-01
A 3-D space-time CE/SE Navier-Stokes solver using an unstructured hexahedral grid is described and applied to a circular jet screech noise computation. The present numerical results for an underexpanded jet, corresponding to a fully expanded Mach number of 1.42, capture the dominant and nonaxisymmetric 'B' screech mode and are generally in good agreement with existing experiments.
NASA Astrophysics Data System (ADS)
Dawes, W. N.
This paper describes some recent improvements made to an unstructed mesh, solution-adaptive three-dimensional Navier-Stokes solver aimed at extending the range of geometric complexity which can be handled in the general context of turbomachinery. The methodology involves generation of a topologically cuboidal mesh, and then the detetion of cells which are not required to allow the formation of relatively complex geometries. This comparatively simple approach permits much of the benefits of an unstructured solution environment to be achieved with minimal complication. Solutions are presented for the highly three-dimensional flows associated with a truncated cylinder in a cross flow, a periodic array of coolant ejection holes, and the overtip leakage flow in an annular cascade of turbine blades.
NASA Technical Reports Server (NTRS)
Lee-Rausch, Elizabeth M.; Hammond, Dana P.; Nielsen, Eric J.; Pirzadeh, S. Z.; Rumsey, Christopher L.
2010-01-01
FUN3D Navier-Stokes solutions were computed for the 4th AIAA Drag Prediction Workshop grid convergence study, downwash study, and Reynolds number study on a set of node-based mixed-element grids. All of the baseline tetrahedral grids were generated with the VGRID (developmental) advancing-layer and advancing-front grid generation software package following the gridding guidelines developed for the workshop. With maximum grid sizes exceeding 100 million nodes, the grid convergence study was particularly challenging for the node-based unstructured grid generators and flow solvers. At the time of the workshop, the super-fine grid with 105 million nodes and 600 million elements was the largest grid known to have been generated using VGRID. FUN3D Version 11.0 has a completely new pre- and post-processing paradigm that has been incorporated directly into the solver and functions entirely in a parallel, distributed memory environment. This feature allowed for practical pre-processing and solution times on the largest unstructured-grid size requested for the workshop. For the constant-lift grid convergence case, the convergence of total drag is approximately second-order on the finest three grids. The variation in total drag between the finest two grids is only 2 counts. At the finest grid levels, only small variations in wing and tail pressure distributions are seen with grid refinement. Similarly, a small wing side-of-body separation also shows little variation at the finest grid levels. Overall, the FUN3D results compare well with the structured-grid code CFL3D. The FUN3D downwash study and Reynolds number study results compare well with the range of results shown in the workshop presentations.
NASA Astrophysics Data System (ADS)
Bruno, Oscar P.; Cubillos, Max
2016-02-01
This paper introduces alternating-direction implicit (ADI) solvers of higher order of time-accuracy (orders two to six) for the compressible Navier-Stokes equations in two- and three-dimensional curvilinear domains. The higher-order accuracy in time results from 1) An application of the backward differentiation formulae time-stepping algorithm (BDF) in conjunction with 2) A BDF-like extrapolation technique for certain components of the nonlinear terms (which makes use of nonlinear solves unnecessary), as well as 3) A novel application of the Douglas-Gunn splitting (which greatly facilitates handling of boundary conditions while preserving higher-order accuracy in time). As suggested by our theoretical analysis of the algorithms for a variety of special cases, an extensive set of numerical experiments clearly indicate that all of the BDF-based ADI algorithms proposed in this paper are "quasi-unconditionally stable" in the following sense: each algorithm is stable for all couples (h , Δt)of spatial and temporal mesh sizes in a problem-dependent rectangular neighborhood of the form (0 ,Mh) × (0 ,Mt). In other words, for each fixed value of Δt below a certain threshold, the Navier-Stokes solvers presented in this paper are stable for arbitrarily small spatial mesh-sizes. The second-order formulation has further been rigorously shown to be unconditionally stable for linear hyperbolic and parabolic equations in two-dimensional space. Although implicit ADI solvers for the Navier-Stokes equations with nominal second-order of temporal accuracy have been proposed in the past, the algorithms presented in this paper are the first ADI-based Navier-Stokes solvers for which second-order or better accuracy has been verified in practice under non-trivial (non-periodic) boundary conditions.
Development of advanced Navier-Stokes solver
NASA Technical Reports Server (NTRS)
Yoon, Seokkwan
1994-01-01
The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.
INS3D: An incompressible Navier-Stokes code in generalized three-dimensional coordinates
NASA Technical Reports Server (NTRS)
Rogers, S. E.; Kwak, D.; Chang, J. L. C.
1987-01-01
The operation of the INS3D code, which computes steady-state solutions to the incompressible Navier-Stokes equations, is described. The flow solver utilizes a pseudocompressibility approach combined with an approximate factorization scheme. This manual describes key operating features to orient new users. This includes the organization of the code, description of the input parameters, description of each subroutine, and sample problems. Details for more extended operations, including possible code modifications, are given in the appendix.
Application of Aeroelastic Solvers Based on Navier Stokes Equations
NASA Technical Reports Server (NTRS)
Keith, Theo G., Jr.; Srivastava, Rakesh
2001-01-01
The propulsion element of the NASA Advanced Subsonic Technology (AST) initiative is directed towards increasing the overall efficiency of current aircraft engines. This effort requires an increase in the efficiency of various components, such as fans, compressors, turbines etc. Improvement in engine efficiency can be accomplished through the use of lighter materials, larger diameter fans and/or higher-pressure ratio compressors. However, each of these has the potential to result in aeroelastic problems such as flutter or forced response. To address the aeroelastic problems, the Structural Dynamics Branch of NASA Glenn has been involved in the development of numerical capabilities for analyzing the aeroelastic stability characteristics and forced response of wide chord fans, multi-stage compressors and turbines. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading is available. To obtain the steady and unsteady aerodynamic forces for the complex flows around the engine components, for the flow regimes encountered by the rotor, an advanced compressible Navier-Stokes solver is required. A finite volume based Navier-Stokes solver has been developed at Mississippi State University (MSU) for solving the flow field around multistage rotors. The focus of the current research effort, under NASA Cooperative Agreement NCC3- 596 was on developing an aeroelastic analysis code (entitled TURBO-AE) based on the Navier-Stokes solver developed by MSU. The TURBO-AE code has been developed for flutter analysis of turbomachine components and delivered to NASA and its industry partners. The code has been verified. validated and is being applied by NASA Glenn and by aircraft engine manufacturers to analyze the aeroelastic stability characteristics of modem fans, compressors
Algorithmic Enhancements to the VULCAN Navier-Stokes Solver
NASA Technical Reports Server (NTRS)
Litton, D. K.; Edwards, J. R.; White, J. A.
2003-01-01
VULCAN (Viscous Upwind aLgorithm for Complex flow ANalysis) is a cell centered, finite volume code used to solve high speed flows related to hypersonic vehicles. Two algorithms are presented for expanding the range of applications of the current Navier-Stokes solver implemented in VULCAN. The first addition is a highly implicit approach that uses subiterations to enhance block to block connectivity between adjacent subdomains. The addition of this scheme allows more efficient solution of viscous flows on highly-stretched meshes. The second algorithm addresses the shortcomings associated with density-based schemes by the addition of a time-derivative preconditioning strategy. High speed, compressible flows are typically solved with density based schemes, which show a high level of degradation in accuracy and convergence at low Mach numbers (M less than or equal to 0.1). With the addition of preconditioning and associated modifications to the numerical discretization scheme, the eigenvalues will scale with the local velocity, and the above problems will be eliminated. With these additions, VULCAN now has improved convergence behavior for multi-block, highly-stretched meshes and also can solve the Navier-Stokes equations for very low Mach numbers.
A Navier-Stokes solver for turbomachinery applications
Arnone, A.; Swanson, R.C. )
1993-04-01
A computer code for solving the Reynolds-averaged full Navier-Stokes equations has been developed and applied using H- and C-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady-state calculations. A grid independence analysis is presented for a transonic rotor blade. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications.
Laser engine simulation using pressure based Navier-Stokes solver
NASA Astrophysics Data System (ADS)
Youssef, Hazim Saad
1994-03-01
Analysis of the flow field in a laser engine represents a difficult computational problem involving combinations of complex physical and gas-dynamical processes. Following a brief discussion of these processes a calculation procedure using primitive variables formulation on a nonstaggered grid system is introduced. Based on this procedure, a pressure based Navier-Stokes solver (PBNS) is developed using a generalized curvilinear coordinate system. The solver is first tested in application to a subsonic compressible flow over an insulated flat plate and to a flow in an axisymmetric converging-diverging nozzle. Next, the PBNS code is used to analyze the flowfield and performance of a laser thruster. The physical/numerical model includes the geometric ray tracing for the laser beam, beam power absorption, plasma radiation losses, and plasma thermophysical and optical properties. Equilibrium hydrogen is used as a flowing gas and its properties are calculated using the Hydrogen Properties Calculation (HPC) based on the methods of statistical thermodynamics. Two thrustor configurations, two laser types (CO2 and iodide), various laser power levels, and various injection conditions are tested. The results of these tests include the temperature, pressure, velocity, and Mach number contours, as well as tables of the laser beam power absorbed, radiation losses to the thrustor walls, thrust level, and specific impulse. The maximum specific impulse obtained in these tests is 1537 sec for a CO2 laser thruster and 827 sec for an iodide laser thruster. Up to 100% power absorption can be achieved; however, radiation losses from the hot plasma are quite high disallowing a full conversion of the absorbed power into the thermal energy of the propellant. The PBNS code can be used to study the effects of various design parameters on the performance of a laser thruster and provide guidelines for the preliminary design of a laser engine.
Exponential Mixing of the 3D Stochastic Navier-Stokes Equations Driven by Mildly Degenerate Noises
Albeverio, Sergio; Debussche, Arnaud; Xu Lihu
2012-10-15
We prove the strong Feller property and exponential mixing for 3D stochastic Navier-Stokes equation driven by mildly degenerate noises (i.e. all but finitely many Fourier modes being forced) via a Kolmogorov equation approach.
Navier-Stokes solutions for rotating 3-D duct flows
NASA Astrophysics Data System (ADS)
Srivastava, B. N.
1988-07-01
This paper deals with the computation of three-dimensional viscous turbulent flow in a rotating rectangular duct of low aspect ratio using thin-layer Navier-Stokes equations. Scalar form of an approximate factorization implicit scheme along with a modified q-omega turbulence model has been utilized for mean flow predictions. The predicted mean flow behavior has been favorably compared with the experimental data for mean axial velocity, channel pressure and cross-flow velocities at a flow Mach number of 0.05 and a rotational speed of 300 rpm.
Numerical solution of 3D Navier-Stokes equations with upwind implicit schemes
NASA Technical Reports Server (NTRS)
Marx, Yves P.
1990-01-01
An upwind MUSCL type implicit scheme for the three-dimensional Navier-Stokes equations is presented. Comparison between different approximate Riemann solvers (Roe and Osher) are performed and the influence of the reconstructions schemes on the accuracy of the solution as well as on the convergence of the method is studied. A new limiter is introduced in order to remove the problems usually associated with non-linear upwind schemes. The implementation of a diagonal upwind implicit operator for the three-dimensional Navier-Stokes equations is also discussed. Finally the turbulence modeling is assessed. Good prediction of separated flows are demonstrated if a non-equilibrium turbulence model is used.
Navier-Stokes Neutral and Plasma Fluid Modelling in 3D
Riemann, J; Borchardt, M; Schneider, R; Mutzke, A; Rognlien, T; Umansky, M
2004-05-17
The 3D finite volume transport code BoRiS is applied to a system of coupled plasma and neutral fluid equations in a slab. Demonstrating easy implementation of new equations, a new parallel BoRiS version is tested on three different models for the neutral fluid - diffusive, parallel Navier-Stokes and full Navier-Stokes - and the results are compared to each other. Typical effects like density enhancement by ionization of recycled neutrals in front of a target plate can be seen and differences are linked to the neutral models in use.
Turbomachinery Heat Transfer and Loss Modeling for 3D Navier-Stokes Codes
NASA Technical Reports Server (NTRS)
DeWitt, Kenneth; Ameri, Ali
2005-01-01
This report's contents focus on making use of NASA Glenn on-site computational facilities,to develop, validate, and apply models for use in advanced 3D Navier-Stokes Computational Fluid Dynamics (CFD) codes to enhance the capability to compute heat transfer and losses in turbomachiney.
On the Dynamic Programming Approach for the 3D Navier-Stokes Equations
Manca, Luigi
2008-06-15
The dynamic programming approach for the control of a 3D flow governed by the stochastic Navier-Stokes equations for incompressible fluid in a bounded domain is studied. By a compactness argument, existence of solutions for the associated Hamilton-Jacobi-Bellman equation is proved. Finally, existence of an optimal control through the feedback formula and of an optimal state is discussed.
NASA Technical Reports Server (NTRS)
Biedron, Robert T.; Vatsa, Veer N.; Atkins, Harold L.
2005-01-01
We apply an unsteady Reynolds-averaged Navier-Stokes (URANS) solver for unstructured grids to unsteady flows on moving and stationary grids. Example problems considered are relevant to active flow control and stability and control. Computational results are presented using the Spalart-Allmaras turbulence model and are compared to experimental data. The effect of grid and time-step refinement are examined.
A Navier-Stokes solver using the LU-SSOR TVD algorithm
NASA Technical Reports Server (NTRS)
Yoon, Seokkwan
1987-01-01
A new Navier-Stokes solver is developed by combining the efficiency of the LU-SSOR scheme and the accuracy of the flux-limited dissipation scheme. Application to laminar and turbulent flows and hypersonic flows proves the reliability of the new algorithm.
NASA Technical Reports Server (NTRS)
Mahajan, Aparajit J.; Dowell, Earl H.; Bliss, Donald B.
1991-01-01
A Lanczos procedure is presently applied to a Navier-Stokes (N-S) solver for eigenvalues and eigenvectors associated with the small-perturbation analysis of the N-S equations' finite-difference representation for airfoil flows; the matrix used is very large, sparse, real, and nonsymmetric. The Lanczos procedure is shown to furnish complete spectral information for the eigenvalues, as required for transient-stability analysis of N-S solvers.
Euler/Navier-Stokes Solvers Applied to Ducted Fan Configurations
NASA Technical Reports Server (NTRS)
Keith, Theo G., Jr.; Srivastava, Rakesh
1997-01-01
Due to noise considerations, ultra high bypass ducted fans have become a more viable design. These ducted fans typically consist of a rotor stage containing a wide chord fan and a stator stage. One of the concerns for this design is the classical flutter that keeps occurring in various unducted fan blade designs. These flutter are catastrophic and are to be avoided in the flight envelope of the engine. Some numerical investigations by Williams, Cho and Dalton, have suggested that a duct around a propeller makes it more unstable. This needs to be further investigated. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading be available. Aerodynamic solvers based on unsteady three-dimensional analysis will provide accurate and fast solutions and are best suited for aeroelastic analysis. The Euler solvers capture significant physics of the flowfield and are reasonably fast. An aerodynamic solver Ref. based on Euler equations had been developed under a separate grant from NASA Lewis in the past. Under the current grant, this solver has been modified to calculate the aeroelastic characteristics of unducted and ducted rotors. Even though, the aeroelastic solver based on three-dimensional Euler equations is computationally efficient, it is still very expensive to investigate the effects of multiple stages on the aeroelastic characteristics. In order to investigate the effects of multiple stages, a two-dimensional multi stage aeroelastic solver was also developed under this task, in collaboration with Dr. T. S. R. Reddy of the University of Toledo. Both of these solvers were applied to several test cases and validated against experimental data, where available.
Navier-Stokes cascade analysis with a stiff Kappa-Epsilon turbulence solver
NASA Technical Reports Server (NTRS)
Liu, Jong-Shang; Sockol, Peter M.; Prahl, Joseph M.
1987-01-01
The two dimensional, compressible, thin layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and the kinetic energy-energy dissipation (k-epsilon) model are solved numerically to simulate the flow through a cascade. The governing equations are solved for the entire flow domain, without the boundary layer assumptions. The stiffness of the k-epsilon equations is discussed. A semi-implicit, Runge-Kutta, time-marching scheme is developed to solve the k-epsilon equations. The impact of the k-epsilon solver on the explicit Runge-Kutta Navier-Stokes solver is discussed. Numerical solutions are presented for two dimensional turbulent flow over a flat plate and a double circular arc cascade and compared with experimental data.
Navier-Stokes Solvers and Generalizations for Reacting Flow Problems
Elman, Howard C
2013-01-27
This is an overview of our accomplishments during the final term of this grant (1 September 2008 -- 30 June 2012). These fall mainly into three categories: fast algorithms for linear eigenvalue problems; solution algorithms and modeling methods for partial differential equations with uncertain coefficients; and preconditioning methods and solvers for models of computational fluid dynamics (CFD).
Bypass Transitional Flow Calculations Using a Navier-Stokes Solver and Two-Equation Models
NASA Technical Reports Server (NTRS)
Liuo, William W.; Shih, Tsan-Hsing; Povinelli, L. A. (Technical Monitor)
2000-01-01
Bypass transitional flows over a flat plate were simulated using a Navier-Stokes solver and two equation models. A new model for the bypass transition, which occurs in cases with high free stream turbulence intensity (TI), is described. The new transition model is developed by including an intermittency correction function to an existing two-equation turbulence model. The advantages of using Navier-Stokes equations, as opposed to boundary-layer equations, in bypass transition simulations are also illustrated. The results for two test flows over a flat plate with different levels of free stream turbulence intensity are reported. Comparisons with the experimental measurements show that the new model can capture very well both the onset and the length of bypass transition.
Navier-Stokes calculations on multi-element airfoils using a chimera-based solver
NASA Technical Reports Server (NTRS)
Jasper, Donald W.; Agrawal, Shreekant; Robinson, Brian A.
1993-01-01
A study of Navier-Stokes calculations of flows about multielement airfoils using a chimera grid approach is presented. The chimera approach utilizes structured, overlapped grids which allow great flexibility of grid arrangement and simplifies grid generation. Calculations are made for two-, three-, and four-element airfoils, and modeling of the effect of gap distance between elements is demonstrated for a two element case. Solutions are obtained using the thin-layer form of the Reynolds averaged Navier-Stokes equations with turbulence closure provided by the Baldwin-Lomax algebraic model or the Baldwin-Barth one equation model. The Baldwin-Barth turbulence model is shown to provide better agreement with experimental data and to dramatically improve convergence rates for some cases. Recently developed, improved farfield boundary conditions are incorporated into the solver for greater efficiency. Computed results show good comparison with experimental data which include aerodynamic forces, surface pressures, and boundary layer velocity profiles.
Nearly Interactive Parabolized Navier-Stokes Solver for High Speed Forebody and Inlet Flows
NASA Technical Reports Server (NTRS)
Benson, Thomas J.; Liou, May-Fun; Jones, William H.; Trefny, Charles J.
2009-01-01
A system of computer programs is being developed for the preliminary design of high speed inlets and forebodies. The system comprises four functions: geometry definition, flow grid generation, flow solver, and graphics post-processor. The system runs on a dedicated personal computer using the Windows operating system and is controlled by graphical user interfaces written in MATLAB (The Mathworks, Inc.). The flow solver uses the Parabolized Navier-Stokes equations to compute millions of mesh points in several minutes. Sample two-dimensional and three-dimensional calculations are demonstrated in the paper.
Flowfield Comparisons from Three Navier-Stokes Solvers for an Axisymmetric Separate Flow Jet
NASA Technical Reports Server (NTRS)
Koch, L. Danielle; Bridges, James; Khavaran, Abbas
2002-01-01
To meet new noise reduction goals, many concepts to enhance mixing in the exhaust jets of turbofan engines are being studied. Accurate steady state flowfield predictions from state-of-the-art computational fluid dynamics (CFD) solvers are needed as input to the latest noise prediction codes. The main intent of this paper was to ascertain that similar Navier-Stokes solvers run at different sites would yield comparable results for an axisymmetric two-stream nozzle case. Predictions from the WIND and the NPARC codes are compared to previously reported experimental data and results from the CRAFT Navier-Stokes solver. Similar k-epsilon turbulence models were employed in each solver, and identical computational grids were used. Agreement between experimental data and predictions from each code was generally good for mean values. All three codes underpredict the maximum value of turbulent kinetic energy. The predicted locations of the maximum turbulent kinetic energy were farther downstream than seen in the data. A grid study was conducted using the WIND code, and comments about convergence criteria and grid requirements for CFD solutions to be used as input for noise prediction computations are given. Additionally, noise predictions from the MGBK code, using the CFD results from the CRAFT code, NPARC, and WIND as input are compared to data.
Application of a Scalable, Parallel, Unstructured-Grid-Based Navier-Stokes Solver
NASA Technical Reports Server (NTRS)
Parikh, Paresh
2001-01-01
A parallel version of an unstructured-grid based Navier-Stokes solver, USM3Dns, previously developed for efficient operation on a variety of parallel computers, has been enhanced to incorporate upgrades made to the serial version. The resultant parallel code has been extensively tested on a variety of problems of aerospace interest and on two sets of parallel computers to understand and document its characteristics. An innovative grid renumbering construct and use of non-blocking communication are shown to produce superlinear computing performance. Preliminary results from parallelization of a recently introduced "porous surface" boundary condition are also presented.
Parallelization of Unsteady Adaptive Mesh Refinement for Unstructured Navier-Stokes Solvers
NASA Technical Reports Server (NTRS)
Schwing, Alan M.; Nompelis, Ioannis; Candler, Graham V.
2014-01-01
This paper explores the implementation of the MPI parallelization in a Navier-Stokes solver using adaptive mesh re nement. Viscous and inviscid test problems are considered for the purpose of benchmarking, as are implicit and explicit time advancement methods. The main test problem for comparison includes e ects from boundary layers and other viscous features and requires a large number of grid points for accurate computation. Ex- perimental validation against double cone experiments in hypersonic ow are shown. The adaptive mesh re nement shows promise for a staple test problem in the hypersonic com- munity. Extension to more advanced techniques for more complicated ows is described.
Parallel performance investigations of an unstructured mesh Navier-Stokes solver
NASA Technical Reports Server (NTRS)
Mavriplis, Dimitri J.
2000-01-01
A Reynolds-averaged Navier-Stokes solver based on unstructured mesh techniques for analysis of high-lift configurations is described. The method makes use of an agglomeration multigrid solver for convergence acceleration. Implicit line-smoothing is employed to relieve the stiffness associated with highly stretched meshes. A GMRES technique is also implemented to speed convergence at the expense of additional memory usage. The solver is cache efficient and fully vectorizable, and is parallelized using a two-level hybrid MPI-OpenMP implementation suitable for shared and/or distributed memory architectures, as well as clusters of shared memory machines. Convergence and scalability results are illustrated for various high-lift cases.
A point implicit unstructured grid solver for the Euler and Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Thareja, Rajiv R.; Stewart, James R.; Hassan, Obey; Morgan, Ken; Peraire, Jaime
1988-01-01
An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.
Nonlinear Aeroelastic Analysis Using a Time-Accurate Navier-Stokes Equations Solver
NASA Technical Reports Server (NTRS)
Kuruvila, Geojoe; Bartels, Robert E.; Hong, Moeljo S.; Bhatia, G.
2007-01-01
A method to simulate limit cycle oscillation (LCO) due to control surface freeplay using a modified CFL3D, a time-accurate Navier-Stokes computational fluid dynamics (CFD) analysis code with structural modeling capability, is presented. This approach can be used to analyze aeroelastic response of aircraft with structural behavior characterized by nonlinearity in the force verses displacement curve. A limited validation of the method, using very low Mach number experimental data for a three-degrees-of-freedom (pitch/plunge/flap deflection) airfoil model with flap freeplay, is also presented.
Oasis: A high-level/high-performance open source Navier-Stokes solver
NASA Astrophysics Data System (ADS)
Mortensen, Mikael; Valen-Sendstad, Kristian
2015-03-01
Oasis is a high-level/high-performance finite element Navier-Stokes solver written from scratch in Python using building blocks from the FEniCS project (fenicsproject.org). The solver is unstructured and targets large-scale applications in complex geometries on massively parallel clusters. Oasis utilizes MPI and interfaces, through FEniCS, to the linear algebra backend PETSc. Oasis advocates a high-level, programmable user interface through the creation of highly flexible Python modules for new problems. Through the high-level Python interface the user is placed in complete control of every aspect of the solver. A version of the solver, that is using piecewise linear elements for both velocity and pressure, is shown to reproduce very well the classical, spectral, turbulent channel simulations of Moser et al. (1999). The computational speed is strongly dominated by the iterative solvers provided by the linear algebra backend, which is arguably the best performance any similar implicit solver using PETSc may hope for. Higher order accuracy is also demonstrated and new solvers may be easily added within the same framework.
A least-squares finite element method for 3D incompressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Jiang, Bo-Nan; Lin, T. L.; Hou, Lin-Jun; Povinelli, Louis A.
1993-01-01
The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.
Parallel computation of 3-D Navier-Stokes flowfields for supersonic vehicles
NASA Technical Reports Server (NTRS)
Ryan, James S.; Weeratunga, Sisira
1993-01-01
Multidisciplinary design optimization of aircraft will require unprecedented capabilities of both analysis software and computer hardware. The speed and accuracy of the analysis will depend heavily on the computational fluid dynamics (CFD) module which is used. A new CFD module has been developed to combine the robust accuracy of conventional codes with the ability to run on parallel architectures. This is achieved by parallelizing the ARC3D algorithm, a central-differenced Navier-Stokes method, on the Intel iPSC/860. The computed solutions are identical to those from conventional machines. Computational speed on 64 processors is comparable to the rate on one Cray Y-MP processor and will increase as new generations of parallel computers become available.
Discrete sensitivity derivatives of the Navier-Stokes equations with a parallel Krylov solver
NASA Technical Reports Server (NTRS)
Ajmani, Kumud; Taylor, Arthur C., III
1994-01-01
This paper solves an 'incremental' form of the sensitivity equations derived by differentiating the discretized thin-layer Navier Stokes equations with respect to certain design variables of interest. The equations are solved with a parallel, preconditioned Generalized Minimal RESidual (GMRES) solver on a distributed-memory architecture. The 'serial' sensitivity analysis code is parallelized by using the Single Program Multiple Data (SPMD) programming model, domain decomposition techniques, and message-passing tools. Sensitivity derivatives are computed for low and high Reynolds number flows over a NACA 1406 airfoil on a 32-processor Intel Hypercube, and found to be identical to those computed on a single-processor Cray Y-MP. It is estimated that the parallel sensitivity analysis code has to be run on 40-50 processors of the Intel Hypercube in order to match the single-processor processing time of a Cray Y-MP.
A three dimensional multigrid Reynolds-averaged Navier-Stokes solver for unstructured meshes
NASA Technical Reports Server (NTRS)
Mavriplis, D. J.
1994-01-01
A three-dimensional unstructured mesh Reynolds averaged Navier-Stokes solver is described. Turbulence is simulated using a single field-equation model. Computational overheads are minimized through the use of a single edge-based data-structure, and efficient multigrid solution technique, and the use of multi-tasking on shared memory multi-processors. The accuracy and efficiency of the code are evaluated by computing two-dimensional flows in three dimensions and comparing with results from a previously validated two-dimensional code which employs the same solution algorithm. The feasibility of computing three-dimensional flows on grids of several million points in less than two hours of wall clock time is demonstrated.
Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D
NASA Technical Reports Server (NTRS)
Pao, S. Paul; Abdol-Hamid, Khaled S.
1996-01-01
This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.
Preconditioned implicit solvers for the Navier-Stokes equations on distributed-memory machines
NASA Technical Reports Server (NTRS)
Ajmani, Kumud; Liou, Meng-Sing; Dyson, Rodger W.
1994-01-01
The GMRES method is parallelized, and combined with local preconditioning to construct an implicit parallel solver to obtain steady-state solutions for the Navier-Stokes equations of fluid flow on distributed-memory machines. The new implicit parallel solver is designed to preserve the convergence rate of the equivalent 'serial' solver. A static domain-decomposition is used to partition the computational domain amongst the available processing nodes of the parallel machine. The SPMD (Single-Program Multiple-Data) programming model is combined with message-passing tools to develop the parallel code on a 32-node Intel Hypercube and a 512-node Intel Delta machine. The implicit parallel solver is validated for internal and external flow problems, and is found to compare identically with flow solutions obtained on a Cray Y-MP/8. A peak computational speed of 2300 MFlops/sec has been achieved on 512 nodes of the Intel Delta machine,k for a problem size of 1024 K equations (256 K grid points).
A high-order gas-kinetic Navier-Stokes flow solver
Li Qibing; Xu Kun; Fu Song
2010-09-20
The foundation for the development of modern compressible flow solver is based on the Riemann solution of the inviscid Euler equations. The high-order schemes are basically related to high-order spatial interpolation or reconstruction. In order to overcome the low-order wave interaction mechanism due to the Riemann solution, the temporal accuracy of the scheme can be improved through the Runge-Kutta method, where the dynamic deficiencies in the first-order Riemann solution is alleviated through the sub-step spatial reconstruction in the Runge-Kutta process. The close coupling between the spatial and temporal evolution in the original nonlinear governing equations seems weakened due to its spatial and temporal decoupling. Many recently developed high-order methods require a Navier-Stokes flux function under piece-wise discontinuous high-order initial reconstruction. However, the piece-wise discontinuous initial data and the hyperbolic-parabolic nature of the Navier-Stokes equations seem inconsistent mathematically, such as the divergence of the viscous and heat conducting terms due to initial discontinuity. In this paper, based on the Boltzmann equation, we are going to present a time-dependent flux function from a high-order discontinuous reconstruction. The theoretical basis for such an approach is due to the fact that the Boltzmann equation has no specific requirement on the smoothness of the initial data and the kinetic equation has the mechanism to construct a dissipative wave structure starting from an initially discontinuous flow condition on a time scale being larger than the particle collision time. The current high-order flux evaluation method is an extension of the second-order gas-kinetic BGK scheme for the Navier-Stokes equations (BGK-NS). The novelty for the easy extension from a second-order to a higher order is due to the simple particle transport and collision mechanism on the microscopic level. This paper will present a hierarchy to construct such
A fast Poisson solver for unsteady incompressible Navier-Stokes equations on the half-staggered grid
NASA Technical Reports Server (NTRS)
Golub, G. H.; Huang, L. C.; Simon, H.; Tang, W. -P.
1995-01-01
In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.
Implementation of algebraic stress models in a general 3-D Navier-Stokes method (PAB3D)
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.
1995-01-01
A three-dimensional multiblock Navier-Stokes code, PAB3D, which was developed for propulsion integration and general aerodynamic analysis, has been used extensively by NASA Langley and other organizations to perform both internal (exhaust) and external flow analysis of complex aircraft configurations. This code was designed to solve the simplified Reynolds Averaged Navier-Stokes equations. A two-equation k-epsilon turbulence model has been used with considerable success, especially for attached flows. Accurate predicting of transonic shock wave location and pressure recovery in separated flow regions has been more difficult. Two algebraic Reynolds stress models (ASM) have been recently implemented in the code that greatly improved the code's ability to predict these difficult flow conditions. Good agreement with Direct Numerical Simulation (DNS) for a subsonic flat plate was achieved with ASM's developed by Shih, Zhu, and Lumley and Gatski and Speziale. Good predictions were also achieved at subsonic and transonic Mach numbers for shock location and trailing edge boattail pressure recovery on a single-engine afterbody/nozzle model.
Quasi-3D Navier-Stokes model for a rotating airfoil
Shen, W.Z.; Soerensen, J.N.
1999-04-10
The design of blade shapes for wind turbines is typically based on employing the blade-element momentum-theory (BEM) with lift and drag forces determined from 2D measurements. The results obtained are reasonable in the vicinity of the design point, but in stalled conditions the BEM is known to underpredict the forces acting on the blades. Here, a quasi-3D model of the unsteady Navier-Stokes equations in a rotating frame of reference has been developed. The equations governing the flow past a rotating blade are approximated using an order of magnitude analysis on the spanwise derivatives. The model takes into account rotational effects and spanwise outflow at computing expenses in the order of what is typical for similar 2D calculations. Results are presented for both laminar and turbulent flows past blades in pure rotation. In the turbulent case the influence of small-scale turbulence is modelled by the one-equation Baldwin-Barth turbulence model. The computations demonstrate that the main influence of rotation is to increase the maximum lift.
Manufactured solutions for steady-flow Reynolds-averaged Navier-Stokes solvers
NASA Astrophysics Data System (ADS)
Eça, L.; Hoekstra, M.; Vaz, G.
2012-06-01
This paper presents manufactured solutions (MS's) for code verification of incompressible flow solvers based on the Reynolds-averaged Navier-Stokes (RANS) equations. The proposed solutions mimic statistically steady, two-dimensional or three-dimensional near-wall turbulent flows in a simple domain (rectangle or rectangular box) at a given Reynolds number. The proposed analytical functions cover the mean flow quantities and the dependent variables of several eddy-viscosity turbulence models. Namely, the undamped eddy-viscosity of the Spalart and Allmaras and Menter one-equations models, from the one (SKL) and two-equation (KSKL) models proposed by Menter, the turbulence kinetic energy and the turbulence frequency included in two-equation k - ω models. A basic flow field resembling a turbulent flat plate flow is constructed with the turbulence quantities defined from 'automatic wall functions' that are supposed to reproduce more or less the normal behaviour of these variables. Alternative flow fields are constructed superposing a perturbation flow field that creates a 'recirculation zone'. However, the near-wall solution of the basic flow is kept to avoid zero friction at the wall. Three-dimensional MS's are obtained from the blending of the basic two-dimensional MS's in the transverse direction. All flow fields satisfy mass conservation, i.e. mean velocity fields are divergence-free. The source functions required for the balancing of momentum and turbulence quantities transport equations and all the dependent variables and their derivatives are available in Fortran 90 modules.
NASA Astrophysics Data System (ADS)
Ge, Liang; Sotiropoulos, Fotis
2007-08-01
A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g. the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [A. Gilmanov, F. Sotiropoulos, A hybrid cartesian/immersed boundary method for simulating flows with 3d, geometrically complex, moving bodies, Journal of Computational Physics 207 (2005) 457-492.]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow
On the use of a three-dimensional Navier-Stokes solver for rocket engine pump impeller design
NASA Astrophysics Data System (ADS)
Chen, Wei-Chung; Prueger, George H.; Chan, Daniel C.; Eastland, Anthony H.
1992-07-01
A 3D Reynolds-averaged Navier-Stokes Solver and a Fast Grid Generator (FGG), developed specially for centrifugal impeller design, were incorporated into the pump impeller design process. The impeller performance from the CFD analysis was compared to one-dimensional prediction. Both analyses showed good agreement of the impeller hydraulic efficiency, 94.5 percent, but with an 8 percent discrepancy of Euler head prediction. The impeller blade angle, discharge hub to shroud width, axial length and blade stacking were systematically changed to achieve an optimum impeller design. Impeller overall efficiency, loss distribution, hub-to-tip flow angle distortion and blade-to-blade flow angle change are among those criteria used to evaluate impeller performance. Two grid sizes, one with 10 K grid points and one with 80 K grid points were used to evaluate grid dependency issues. The effects of grid resolution on the accuracy and turnaround time are discussed. In conclusion, it is demonstrated that CFD can be effectively used for design and optimization of rocket engine pump components.
On the use of a three-dimensional Navier-Stokes solver for rocket engine pump impeller design
NASA Technical Reports Server (NTRS)
Chen, Wei-Chung; Prueger, George H.; Chan, Daniel C.; Eastland, Anthony H.
1992-01-01
A 3D Reynolds-averaged Navier-Stokes Solver and a Fast Grid Generator (FGG), developed specially for centrifugal impeller design, were incorporated into the pump impeller design process. The impeller performance from the CFD analysis was compared to one-dimensional prediction. Both analyses showed good agreement of the impeller hydraulic efficiency, 94.5 percent, but with an 8 percent discrepancy of Euler head prediction. The impeller blade angle, discharge hub to shroud width, axial length and blade stacking were systematically changed to achieve an optimum impeller design. Impeller overall efficiency, loss distribution, hub-to-tip flow angle distortion and blade-to-blade flow angle change are among those criteria used to evaluate impeller performance. Two grid sizes, one with 10 K grid points and one with 80 K grid points were used to evaluate grid dependency issues. The effects of grid resolution on the accuracy and turnaround time are discussed. In conclusion, it is demonstrated that CFD can be effectively used for design and optimization of rocket engine pump components.
Local existence and Gevrey regularity of 3-D Navier-Stokes equations with ℓp initial data
NASA Astrophysics Data System (ADS)
Biswas, Animikh
We obtain local existence and Gevrey regularity of 3-D periodic Navier-Stokes equations in case the sequence of Fourier coefficients of the initial data is in ℓp (p<3/2). The ℓp norm of the sequence of Fourier coefficients of the solution and its analogous Gevrey norm remains bounded on a time interval whose length depends only on the size of the body force and the ℓp norm of the Fourier coefficient sequence of the initial data. The control on the Gevrey norm produces explicit estimates on the analyticity radius of the solution as in Foias and Temam (J. Funct. Anal. 87 (1989) 359-369). The results provide an alternate approach in estimating the space-analyticity radius of solutions to Navier-Stokes equations than the one presented by Grujić and Kukavica (J. Funct. Anal. 152 (1998) 447-466).
NASA Technical Reports Server (NTRS)
Kiris, Cetin C.; Kwak, Dochan; Rogers, Stuart E.
2002-01-01
This paper reviews recent progress made in incompressible Navier-Stokes simulation procedures and their application to problems of engineering interest. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. A summary of efforts in flow solver development is given followed by numerical studies of a few example problems of current interest. Both steady and unsteady solution algorithms and their salient features are discussed. Solvers discussed here are based on a structured-grid approach using either a finite -difference or a finite-volume frame work. As a grand-challenge application of these solvers, an unsteady turbopump flow simulation procedure has been developed which utilizes high performance computing platforms. In the paper, the progress toward the complete simulation capability of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of two parallel computing algorithms that have been implemented in the INS3D code. The relative motion of the grid systems for the rotorstator interaction was obtained using overact grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on SCSI Origin 3000 systems at NASA Ames Research Center. The same procedure has been extended to the development of NASA-DeBakey Ventricular Assist Device (VAD) that is based on an axial blood pump. Computational, and clinical analysis of this device are presented.
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Lakshmanan, B.; Carlson, John R.
1995-01-01
A three-dimensional Navier-Stokes solver was used to determine how accurately computations can predict local and average skin friction coefficients for attached and separated flows for simple experimental geometries. Algebraic and transport equation closures were used to model turbulence. To simulate anisotropic turbulence, the standard two-equation turbulence model was modified by adding nonlinear terms. The effects of both grid density and the turbulence model on the computed flow fields were also investigated and compared with available experimental data for subsonic and supersonic free-stream conditions.
On Energy Cascades in the Forced 3D Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Dascaliuc, R.; Grujić, Z.
2016-02-01
We show—in the framework of physical scales and (K_1,K_2) -averages—that Kolmogorov's dissipation law combined with the smallness condition on a Taylor length scale is sufficient to guarantee energy cascades in the forced Navier-Stokes equations. Moreover, in the periodic case we establish restrictive scaling laws—in terms of Grashof number—for kinetic energy, energy flux, and energy dissipation rate. These are used to improve our sufficient condition for forced cascades in physical scales.
On Energy Cascades in the Forced 3D Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Dascaliuc, R.; Grujić, Z.
2016-06-01
We show—in the framework of physical scales and (K_1,K_2)-averages—that Kolmogorov's dissipation law combined with the smallness condition on a Taylor length scale is sufficient to guarantee energy cascades in the forced Navier-Stokes equations. Moreover, in the periodic case we establish restrictive scaling laws—in terms of Grashof number—for kinetic energy, energy flux, and energy dissipation rate. These are used to improve our sufficient condition for forced cascades in physical scales.
Three-dimensional full Navier-Stokes solvers for incompressible flows past arbitrary geometries
NASA Astrophysics Data System (ADS)
Deng, G. B.; Piquet, J.; Queutey, P.; Visonneau, M.
1991-05-01
The computation of the three-dimensional viscous flow past several geometries is investigated. An iterative technique resting on the fully elliptic mode is applied to the Reynolds-Averaged Navier-Stokes Equations written in a nonorthogonal curvilinear body-fitted coordinate system. Results of the computation are compared with available experiments.
On the Global Regularity of a Helical-Decimated Version of the 3D Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Biferale, Luca; Titi, Edriss S.
2013-06-01
We study the global regularity, for all time and all initial data in H 1/2, of a recently introduced decimated version of the incompressible 3D Navier-Stokes (dNS) equations. The model is based on a projection of the dynamical evolution of Navier-Stokes (NS) equations into the subspace where helicity (the L 2-scalar product of velocity and vorticity) is sign-definite. The presence of a second (beside energy) sign-definite inviscid conserved quadratic quantity, which is equivalent to the H 1/2-Sobolev norm, allows us to demonstrate global existence and uniqueness, of space-periodic solutions, together with continuity with respect to the initial conditions, for this decimated 3D model. This is achieved thanks to the establishment of two new estimates, for this 3D model, which show that the H 1/2 and the time average of the square of the H 3/2 norms of the velocity field remain finite. Such two additional bounds are known, in the spirit of the work of H. Fujita and T. Kato (Arch. Ration. Mech. Anal. 16:269-315, 1964; Rend. Semin. Mat. Univ. Padova 32:243-260, 1962), to be sufficient for showing well-posedness for the 3D NS equations. Furthermore, they are directly linked to the helicity evolution for the dNS model, and therefore with a clear physical meaning and consequences.
3-D Navier-Stokes Analysis of Blade Root Aerodynamics for a Tiltrotor Aircraft In Cruise
NASA Technical Reports Server (NTRS)
Romander, Ethan
2006-01-01
The blade root area of a tiltrotor aircraft's rotor is constrained by a great many factors, not the least of which is aerodynamic performance in cruise. For this study, Navier-Stokes CFD techniques are used to study the aerodynamic performance in cruise of a rotor design as a function of airfoil thickness along the blade and spinner shape. Reducing airfoil thickness along the entire blade will be shown to have the greatest effect followed by smaller but still significant improvements achieved by reducing the thickness of root airfoils only. Furthermore, altering the shape of the spinner will be illustrated as a tool to tune the aerodynamic performance very near the blade root.
Dawes, W.N. )
1993-04-01
This paper describes recent developments to a three-dimensional, unstructured mesh, solution-adaptive Navier-Stokes solver. By adopting a simple, pragmatic but systematic approach to mesh generation, the range of simulations that can be attempted is extended toward arbitrary geometries. The combined benefits of the approach result in a powerful analytical ability. Solutions for a wide range of flows are presented, including a transonic compressor rotor, a centrifugal impeller, a steam turbine nozzle guide vane with casing extraction belt, the internal coolant passage of a radial inflow turbine, and a turbine disk cavity flow.
Calculation of a simulated 3-D high speed inlet using the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Knight, D. D.
1983-01-01
A hybrid numerical algorithm, developed to solve the full three-dimensional Navier-Stokes equations, is applied to the computation of the flowfield in a simulated three-dimensional high speed aircraft inlet at a Mach number of 2.5 and Reynolds number of 1.4 x 10 to the 7th based on inlet length. The numerical algorithm incorporates a coordinate transformation in order to handle general flow geometries, and utilizes the algebraic turbulent eddy viscosity model of Baldwin and Lomax. The hybrid algorithm has been vectorized on the CDC CYBER 203 computer using the SL/1 vector programming language developed at NASA Langley. The computed results are compared with experimental measurements of the ramp and cowl static pressures, and boundary layer pitot profiles. The results are also compared with a previous two-dimensional Navier-Stokes computation of the same configuration. The agreement with the experimental data is generally good; however, additional improvements in turbulence modeling are needed.
Development of an Aeroelastic Code Based on an Euler/Navier-Stokes Aerodynamic Solver
NASA Technical Reports Server (NTRS)
Bakhle, Milind A.; Srivastava, Rakesh; Keith, Theo G., Jr.; Stefko, George L.; Janus, Mark J.
1996-01-01
This paper describes the development of an aeroelastic code (TURBO-AE) based on an Euler/Navier-Stokes unsteady aerodynamic analysis. A brief review of the relevant research in the area of propulsion aeroelasticity is presented. The paper briefly describes the original Euler/Navier-Stokes code (TURBO) and then details the development of the aeroelastic extensions. The aeroelastic formulation is described. The modeling of the dynamics of the blade using a modal approach is detailed, along with the grid deformation approach used to model the elastic deformation of the blade. The work-per-cycle approach used to evaluate aeroelastic stability is described. Representative results used to verify the code are presented. The paper concludes with an evaluation of the development thus far, and some plans for further development and validation of the TURBO-AE code.
Implementation of wall boundary conditions for transpiration in F3D thin-layer Navier-Stokes code
NASA Technical Reports Server (NTRS)
Kandula, M.; Martin, F. W., Jr.
1991-01-01
Numerical boundary conditions for mass injection/suction at the wall are incorporated in the thin-layer Navier-Stokes code, F3D. The accuracy of the boundary conditions and the code is assessed by a detailed comparison of the predictions of velocity distributions and skin-friction coefficients with exact similarity solutions for laminar flow over a flat plate with variable blowing/suction, and measurements for turbulent flow past a flat plate with uniform blowing. In laminar flow, F3D predictions for friction coefficient compare well with exact similarity solution with and without suction, but produces large errors at moderate-to-large values of blowing. A slight Mach number dependence of skin-friction coefficient due to blowing in turbulent flow is computed by F3D code. Predicted surface pressures for turbulent flow past an airfoil with mass injection are in qualitative agreement with measurements for a flat plate.
Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Wedan, Bruce W.; Abid, Ridha
1989-01-01
A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.
On the Helicity in 3D-Periodic Navier-Stokes Equations II: The Statistical Case
NASA Astrophysics Data System (ADS)
Foias, Ciprian; Hoang, Luan; Nicolaenko, Basil
2009-09-01
We study the asymptotic behavior of the statistical solutions to the Navier-Stokes equations using the normalization map [9]. It is then applied to the study of mean energy, mean dissipation rate of energy, and mean helicity of the spatial periodic flows driven by potential body forces. The statistical distribution of the asymptotic Beltrami flows are also investigated. We connect our mathematical analysis with the empirical theory of decaying turbulence. With appropriate mathematically defined ensemble averages, the Kolmogorov universal features are shown to be transient in time. We provide an estimate for the time interval in which those features may still be present. Our collaborator and friend Basil Nicolaenko passed away in September of 2007, after this work was completed. Honoring his contribution and friendship, we dedicate this article to him.
NASA Technical Reports Server (NTRS)
Demuren, A. O.; Ibraheem, S. O.
1993-01-01
The convergence characteristics of various approximate factorizations for the 3D Euler and Navier-Stokes equations are examined using the von-Neumann stability analysis method. Three upwind-difference based factorizations and several central-difference based factorizations are considered for the Euler equations. In the upwind factorizations both the flux-vector splitting methods of Steger and Warming and van Leer are considered. Analysis of the Navier-Stokes equations is performed only on the Beam and Warming central-difference scheme. The range of CFL numbers over which each factorization is stable is presented for one-, two-, and three-dimensional flow. Also presented for each factorization is the CFL number at which the maximum eigenvalue is minimized, for all Fourier components, as well as for the high frequency range only. The latter is useful for predicting the effectiveness of multigrid procedures with these schemes as smoothers. Further, local mode analysis is performed to test the suitability of using a uniform flow field in the stability analysis. Some inconsistencies in the results from previous analyses are resolved.
Iterative solvers for Navier-Stokes equations: Experiments with turbulence model
Page, M.; Garon, A.
1994-12-31
In the framework of developing software for the prediction of flows in hydraulic turbine components, Reynolds averaged Navier-Stokes equations coupled with {kappa}-{omega} two-equation turbulence model are discretized by finite element method. Since the resulting matrices are large, sparse and nonsymmetric, strategies based on CG-type iterative methods must be devised. A segregated solution strategy decouples the momentum equation, the {kappa} transport equation and the {omega} transport equation. These sets of equations must be solved while satisfying constraint equations. Experiments with orthogonal projection method are presented for the imposition of essential boundary conditions in a weak sense.
An implicit three-dimensional Navier-Stokes solver for compressible flow
NASA Technical Reports Server (NTRS)
Yoon, Seokkwan; Kwak, Dochan
1991-01-01
A three-dimensional numerical method based on the lower-upper symmetric-Gauss-Seidel implicit scheme in conjunction with the flux-limited dissipation model is developed for solving the compressible Navier-Stokes equations. A new computer code which is based on this method requires only 9 microsec per grid-point per iteration on a single processor of a Cray YMP computer and executes at the sustained rate of 170 MFLOPS. A reduction of 4 orders of magnitude in the residual for a high Reynolds number flow using 230 K grid points is obtained in 24 minutes. The computational results compare well with available experimental data.
A unified multigrid solver for the Navier-Stokes equations on mixed element meshes
NASA Technical Reports Server (NTRS)
Mavriplis, D. J.; Venkatakrishnan, V.
1995-01-01
A unified multigrid solution technique is presented for solving the Euler and Reynolds-averaged Navier-Stokes equations on unstructured meshes using mixed elements consisting of triangles and quadrilaterals in two dimensions, and of hexahedra, pyramids, prisms, and tetrahedra in three dimensions. While the use of mixed elements is by no means a novel idea, the contribution of the paper lies in the formulation of a complete solution technique which can handle structured grids, block structured grids, and unstructured grids of tetrahedra or mixed elements without any modification. This is achieved by discretizing the full Navier-Stokes equations on tetrahedral elements, and the thin layer version of these equations on other types of elements, while using a single edge-based data-structure to construct the discretization over all element types. An agglomeration multigrid algorithm, which naturally handles meshes of any types of elements, is employed to accelerate convergence. An automatic algorithm which reduces the complexity of a given triangular or tetrahedral mesh by merging candidate triangular or tetrahedral elements into quadrilateral or prismatic elements is also described. The gains in computational efficiency afforded by the use of non-simplicial meshes over fully tetrahedral meshes are demonstrated through several examples.
Progress in development of a Navier-Stokes solver for evaluation of iced airfoil performance
NASA Technical Reports Server (NTRS)
Potapczuk, M. G.; Gerhart, P. M.
1985-01-01
A method is being developed for evaluation of the flow field behavior about an airfoil with significant ice accretion on the leading edge. The computer code, being evaluated for this purpose, solves the Navier-Stokes equations in a body-fitted curvilinear coordinate system. This requires the use of a grid generation code to transform the x-y coordinates of the physical space into xi-eta coordinates of the computational space. Evaluation of the suitability of these two codes for predicting iced airfoil performance is presently being carried out in anticipation of use in an overall icing analysis effort. Results of this evaluation to date indicate good correlation with known information on clean airfoils. Preliminary results for rime and glaze, iced airfoil shapes are also presented.
On the role of artificial viscosity in Navier-Stokes solvers
NASA Technical Reports Server (NTRS)
Mahajan, Aparajit J.; Dowell, Earl H.; Bliss, Donald B.
1989-01-01
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over the airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small perturbation about a steady inviscid flow over a NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 degrees. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time marching scheme are also presented. This procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of terms for artificial viscosity.
Application of an unsteady Navier-Stokes solver to transonic turbine design
NASA Technical Reports Server (NTRS)
Rangwalla, Akil A.; Madavan, Nateri K.; Johnson, Paul D.
1991-01-01
This study presents a numerical evaluation of the performance of the first stage of a new-generation turbine design. The numerical method solves the two-dimensional Navier-Stokes equations using a system of patched grids. Three-dimensional effects of stream-tube contraction are also modeled. The study focuses on the effects of axial gap variation on the unsteady rotor-stator interactions and on stage performance. Results are presented for three different axial gaps. The results indicate that the unsteady interactions can be very large in this design. These interactions affect not only the stage efficiency but also substantially alter the time-averaged features of the flow. In particular, for the case of the smallest axial gap, it was found that there was an unsteady shock on the stator suction surface which spanned the gap region and impinged upon the moving rotor airfoils.
NASA Astrophysics Data System (ADS)
Cervone, A.; Manservisi, S.; Scardovelli, R.
2010-09-01
A multilevel VOF approach has been coupled to an accurate finite element Navier-Stokes solver in axisymmetric geometry for the simulation of incompressible liquid jets with high density ratios. The representation of the color function over a fine grid has been introduced to reduce the discontinuity of the interface at the cell boundary. In the refined grid the automatic breakup and coalescence occur at a spatial scale much smaller than the coarse grid spacing. To reduce memory requirements, we have implemented on the fine grid a compact storage scheme which memorizes the color function data only in the mixed cells. The capillary force is computed by using the Laplace-Beltrami operator and a volumetric approach for the two principal curvatures. Several simulations of axisymmetric jets have been performed to show the accuracy and robustness of the proposed scheme.
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.
1990-01-01
The development and applications of multiblock/multizone and adaptive grid methodologies for solving the three-dimensional simplified Navier-Stokes equations are described. Adaptive grid and multiblock/multizone approaches are introduced and applied to external and internal flow problems. These new implementations increase the capabilities and flexibility of the PAB3D code in solving flow problems associated with complex geometry.
Time-Shifted Boundary Conditions Used for Navier-Stokes Aeroelastic Solver
NASA Technical Reports Server (NTRS)
Srivastava, Rakesh
1999-01-01
Under the Advanced Subsonic Technology (AST) Program, an aeroelastic analysis code (TURBO-AE) based on Navier-Stokes equations is currently under development at NASA Lewis Research Center s Machine Dynamics Branch. For a blade row, aeroelastic instability can occur in any of the possible interblade phase angles (IBPA s). Analyzing small IBPA s is very computationally expensive because a large number of blade passages must be simulated. To reduce the computational cost of these analyses, we used time shifted, or phase-lagged, boundary conditions in the TURBO-AE code. These conditions can be used to reduce the computational domain to a single blade passage by requiring the boundary conditions across the passage to be lagged depending on the IBPA being analyzed. The time-shifted boundary conditions currently implemented are based on the direct-store method. This method requires large amounts of data to be stored over a period of the oscillation cycle. On CRAY computers this is not a major problem because solid-state devices can be used for fast input and output to read and write the data onto a disk instead of storing it in core memory.
Role of artificial viscosity in Euler and Navier-Stokes solvers
NASA Technical Reports Server (NTRS)
Mahajan, Aparajit J.; Dowell, Earl H.; Bliss, Donald B.
1991-01-01
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small time-dependent perturbation about a steady inviscid flow over an NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 deg. The method has been applied to inviscid flows here, but as discussed is also applicable to viscous flows. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundaries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time-marching scheme are also presented. The eigenvalue procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of artificial viscosity.
A lattice-Boltzmann scheme of the Navier-Stokes equations on a 3D cuboid lattice
NASA Astrophysics Data System (ADS)
Min, Haoda; Peng, Cheng; Wang, Lian-Ping
2015-11-01
The standard lattice-Boltzmann method (LBM) for fluid flow simulation is based on a square (in 2D) or cubic (in 3D) lattice grids. Recently, two new lattice Boltzmann schemes have been developed on a 2D rectangular grid using the MRT (multiple-relaxation-time) collision model, by adding a free parameter in the definition of moments or by extending the equilibrium moments. Here we developed a lattice Boltzmann model on 3D cuboid lattice, namely, a lattice grid with different grid lengths in different spatial directions. We designed our MRT-LBM model by matching the moment equations from the Chapman-Enskog expansion with the Navier-Stokes equations. The model guarantees correct hydrodynamics. A second-order term is added to the equilibrium moments in order to restore the isotropy of viscosity on a cuboid lattice. The form and the coefficients of the extended equilibrium moments are determined through an inverse design process. An additional benefit of the model is that the viscosity can be adjusted independent of the stress-moment relaxation parameter, thus improving the numerical stability of the model. The resulting cuboid MRT-LBM model is then validated through benchmark simulations using laminar channel flow, turbulent channel flow, and the 3D Taylor-Green vortex flow.
NASA Astrophysics Data System (ADS)
Deng, Liang; Bai, Hanli; Wang, Fang; Xu, Qingxin
2016-06-01
CPU/GPU computing allows scientists to tremendously accelerate their numerical codes. In this paper, we port and optimize a double precision alternating direction implicit (ADI) solver for three-dimensional compressible Navier-Stokes equations from our in-house Computational Fluid Dynamics (CFD) software on heterogeneous platform. First, we implement a full GPU version of the ADI solver to remove a lot of redundant data transfers between CPU and GPU, and then design two fine-grain schemes, namely “one-thread-one-point” and “one-thread-one-line”, to maximize the performance. Second, we present a dual-level parallelization scheme using the CPU/GPU collaborative model to exploit the computational resources of both multi-core CPUs and many-core GPUs within the heterogeneous platform. Finally, considering the fact that memory on a single node becomes inadequate when the simulation size grows, we present a tri-level hybrid programming pattern MPI-OpenMP-CUDA that merges fine-grain parallelism using OpenMP and CUDA threads with coarse-grain parallelism using MPI for inter-node communication. We also propose a strategy to overlap the computation with communication using the advanced features of CUDA and MPI programming. We obtain speedups of 6.0 for the ADI solver on one Tesla M2050 GPU in contrast to two Xeon X5670 CPUs. Scalability tests show that our implementation can offer significant performance improvement on heterogeneous platform.
A compressible Navier-Stokes flow solver using the Newton-Krylov method on unstructured grids
NASA Astrophysics Data System (ADS)
Wong, Peterson
A Newton-Krylov algorithm is presented for the compressible Navier-Stokes equations on hybrid unstructured grids. The Spalart-Allmaras turbulence model is used for turbulent flows. The spatial discretization is based on a finite-volume matrix dissipation scheme. A preconditioned matrix-free generalized minimal residual method is used to solve the linear system that arises in the Newton iterations. The incomplete lower-upper factorization based on an approximate Jacobian is used as the preconditioner after applying the reverse Cuthill-McKee reordering. Various aspects of the Newton-Krylov algorithm are studied to improve efficiency and reliability. The inexact Newton method is studied to avoid over-solving of the linear system to reduce computational cost. The ILU(1) approach is selected in three dimensions, based on a comparison among various preconditioners. Approximate viscous formulations involving only the nearest neighboring terms are studied to reduce the cost of preconditioning. The resulting preconditioners are found to be effective and provide Newton-type convergence. Scaling of the linear system is studied to improve convergence of the inexact matrix-free approach. Numerical studies are performed for two-dimensional cases as well as flows over the ONERA M6 wing and the DLR-F6 wing-body configuration. A ten-order-of-magnitude residual reduction can be obtained with a computing cost equivalent to 4,000 residual function evaluations for two-dimensional cases, while the same convergence can be obtained in 5,500 and 8,000 function evaluations for the wing and wing-body configuration, respectively, on grids with a half million nodes.
Fast solvers for finite difference approximations for the Stokes and Navier-Stokes equations
Shin, D.
1992-01-01
The authors consider several methods for solving the linear equations arising from finite difference discretizations of the Stokes equations. The pressure equation method presented here for the first time, apparently, and the method, presented by Bramble and Pasciak, are shown to have computational effort that grows slowly with the number of grid points. The methods work with second-order accurate discretizations. Computational results are shown for both the Stokes and incompressible Navier-Stokes at low Reynolds number. The inf-sup conditions resulting from three finite difference approximations of the Stokes equations are proven. These conditions are used to prove that the Schur complement Q[sub h] of the linear system generated by each of these approximations is bounded uniformly away from zero. For the pressure equation method, this guarantees that the conjugate gradient method applied to Q[sub h] converges in a finite number of iterations which is independent of mesh size. The fact that Q[sub h] is bounded below is used to prove convergence estimates for the solutions generated by these finite difference approximations. One of the estimates is for a staggered grid and the estimate of the scheme shows that both the pressure and the velocity parts of the solution are second-order accurate. Iterative methods are compared by the use of the regularized central differencing introduced by Strikwerda. Several finite difference approximations of the Stokes equations by the SOR method are compared and the excellence of the approximations by the regularized central differencing over the other finite difference approximation is mentioned. This difference gives rise to a linear equation with a matrix which is slightly non-symmetric. The convergence of the typical steepest descent method and conjugate gradient method, which is almost as same as the typical conjugate gradient method, applied to slightly non-symmetric positive definite matrices are proven.
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.; Van Leer, Bram; Roe, Philip L.
1991-01-01
A new two-dimensional approximate Riemann solver has been developed that obtains fluxes on grid faces via wave decomposition. By utilizing information propagation in the velocity-difference directions rather than in the grid-normal directions, this flux function more appropriately interprets and hence more sharply resolves shock and shear waves when they lie oblique to the grid. The model uses five waves to describe the difference in states at a grid face. Two acoustic waves, one shear wave, and one entropy wave propagate in the direction defined by the local velocity difference vector, while the fifth wave is a shear wave that propagates at a right angle to the other four. Test cases presented include a shock reflecting off a wall, a pure shear wave, supersonic flow over an airfoil, and viscous separated airfoil flow. Results using the new model give significantly sharper shock and shear contours than a grid-aligned solver. Navier-Stokes computations over an aifoil show reduced pressure distortions in the separated region as a result of the grid-independent upwinding.
NASA Astrophysics Data System (ADS)
Zhai, Cuili; Zhang, Ting
2015-09-01
In this article, we consider the global well-posedness to the 3-D incompressible inhomogeneous Navier-Stokes equations with a class of large velocity. More precisely, assuming a 0 ∈ B˙ q , 1 /3 q ( R 3 ) and u 0 = ( u0 h , u0 3 ) ∈ B˙ p , 1 - 1 + /3 p ( R 3 ) for p, q ∈ (1, 6) with sup ( /1 p , /1 q ) ≤ /1 3 + inf ( /1 p , /1 q ) , we prove that if C a↑0↑ B˙q1/3 q α (↑u0 3↑ B˙ p , 1 - 1 + /3 p/μ + 1 ) ≤ 1 , /C μ (↑u0 h↑ B˙ p , 1 - 1 + /3 p + ↑u03↑ B˙ p , 1 - 1 + /3 p 1 - α ↑u0h↑ B˙ p , 1 - 1 + /3 p α) ≤ 1 , then the system has a unique global solution a ∈ C ˜ ( [ 0 , ∞ ) ; B˙ q , 1 /3 q ( R 3 ) ) , u ∈ C ˜ ( [ 0 , ∞ ) ; B˙ p , 1 - 1 + /3 p ( R 3 ) ) ∩ L 1 ( R + ; B˙ p , 1 1 + /3 p ( R 3 ) ) . It improves the recent result of M. Paicu and P. Zhang [J. Funct. Anal. 262, 3556-3584 (2012)], where the exponent form of the initial smallness condition is replaced by a polynomial form.
NASA Technical Reports Server (NTRS)
Yokota, Jeffrey W.
1988-01-01
An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.
NASA Technical Reports Server (NTRS)
Lessard, Victor R.
1993-01-01
Computations of three dimensional vortical flows over a generic High Speed Civil Transport (HSCT) configuration with an aspect ratio of 3.04 are performed using a thin-layer Navier-Stokes solver. The HSCT cruise configuration is modeled without leading or trailing edge flap deflections and without engine nacelles. The flow conditions, which correspond to tests done in the NASA Langley 8-Foot Transonic Pressure Tunnel (TPT), are a subsonic Mach number of 0.3 and Reynolds number of 4.4 million for a range-of-attack (-.23 deg to 17.78 deg). The effects of the farfield boundary location with respect to the body are investigated. The boundary layer is assumed turbulent and simulated using an algebraic turbulence model. The key features of the vortices and their interactions are captured. Grid distribution in the vortex regions is critical for predicting the correct induced lift. Computed forces and surface pressures compare reasonably well with the experimental TPT data.
Navier-Stokes and viscous-inviscid interaction
NASA Technical Reports Server (NTRS)
Steger, Joseph L.; Vandalsem, William R.
1989-01-01
Some considerations toward developing numerical procedures for simulating viscous compressible flows are discussed. Both Navier-Stokes and boundary layer field methods are considered. Because efficient viscous-inviscid interaction methods have been difficult to extend to complex 3-D flow simulations, Navier-Stokes procedures are more frequently being utilized even though they require considerably more work per grid point. It would seem a mistake, however, not to make use of the more efficient approximate methods in those regions in which they are clearly valid. Ideally, a general purpose compressible flow solver that can optionally take advantage of approximate solution methods would suffice, both to improve accuracy and efficiency. Some potentially useful steps toward this goal are described: a generalized 3-D boundary layer formulation and the fortified Navier-Stokes procedure.
NASA Astrophysics Data System (ADS)
Blanchard, M.; Schuller, T.; Sipp, D.; Schmid, P. J.
2015-04-01
The response of a laminar premixed methane-air flame subjected to flow perturbations around a steady state is examined experimentally and using a linearized compressible Navier-Stokes solver with a one-step chemistry mechanism to describe combustion. The unperturbed flame takes an M-shape stabilized both by a central bluff body and by the external rim of a cylindrical nozzle. This base flow is computed by a nonlinear direct simulation of the steady reacting flow, and the flame topology is shown to qualitatively correspond to experiments conducted under comparable conditions. The flame is then subjected to acoustic disturbances produced at different locations in the numerical domain, and its response is examined using the linearized solver. This linear numerical model then allows the componentwise investigation of the effects of flow disturbances on unsteady combustion and the feedback from the flame on the unsteady flow field. It is shown that a wrinkled reaction layer produces hydrodynamic disturbances in the fresh reactant flow field that superimpose on the acoustic field. This phenomenon, observed in several experiments, is fully interpreted here. The additional perturbations convected by the mean flow stem from the feedback of the perturbed flame sheet dynamics onto the flow field by a mechanism similar to that of a perturbed vortex sheet. The different regimes where this mechanism prevails are investigated by examining the phase and group velocities of flow disturbances along an axis oriented along the main direction of the flow in the fresh reactant flow field. It is shown that this mechanism dominates the low-frequency response of the wrinkled shape taken by the flame and, in particular, that it fully determines the dynamics of the flame tip from where the bulk of noise is radiated.
Blanchard, M.; Schuller, T.; Sipp, D.; Schmid, P. J.
2015-04-15
The response of a laminar premixed methane-air flame subjected to flow perturbations around a steady state is examined experimentally and using a linearized compressible Navier-Stokes solver with a one-step chemistry mechanism to describe combustion. The unperturbed flame takes an M-shape stabilized both by a central bluff body and by the external rim of a cylindrical nozzle. This base flow is computed by a nonlinear direct simulation of the steady reacting flow, and the flame topology is shown to qualitatively correspond to experiments conducted under comparable conditions. The flame is then subjected to acoustic disturbances produced at different locations in the numerical domain, and its response is examined using the linearized solver. This linear numerical model then allows the componentwise investigation of the effects of flow disturbances on unsteady combustion and the feedback from the flame on the unsteady flow field. It is shown that a wrinkled reaction layer produces hydrodynamic disturbances in the fresh reactant flow field that superimpose on the acoustic field. This phenomenon, observed in several experiments, is fully interpreted here. The additional perturbations convected by the mean flow stem from the feedback of the perturbed flame sheet dynamics onto the flow field by a mechanism similar to that of a perturbed vortex sheet. The different regimes where this mechanism prevails are investigated by examining the phase and group velocities of flow disturbances along an axis oriented along the main direction of the flow in the fresh reactant flow field. It is shown that this mechanism dominates the low-frequency response of the wrinkled shape taken by the flame and, in particular, that it fully determines the dynamics of the flame tip from where the bulk of noise is radiated.
NASA Technical Reports Server (NTRS)
Jones, Kenneth M.; Biedron, Robert T.; Whitlock, Mark
1995-01-01
A computational study was performed to determine the predictive capability of a Reynolds averaged Navier-Stokes code (CFL3D) for two-dimensional and three-dimensional multielement high-lift systems. Three configurations were analyzed: a three-element airfoil, a wing with a full span flap and a wing with a partial span flap. In order to accurately model these complex geometries, two different multizonal structured grid techniques were employed. For the airfoil and full span wing configurations, a chimera or overset grid technique was used. The results of the airfoil analysis illustrated that although the absolute values of lift were somewhat in error, the code was able to predict reasonably well the variation with Reynolds number and flap position. The full span flap analysis demonstrated good agreement with experimental surface pressure data over the wing and flap. Multiblock patched grids were used to model the partial span flap wing. A modification to an existing patched- grid algorithm was required to analyze the configuration as modeled. Comparisons with experimental data were very good, indicating the applicability of the patched-grid technique to analyses of these complex geometries.
A compressible Navier-Stokes solver with two-equation and Reynolds stress turbulence closure models
NASA Technical Reports Server (NTRS)
Morrison, Joseph H.
1992-01-01
This report outlines the development of a general purpose aerodynamic solver for compressible turbulent flows. Turbulent closure is achieved using either two equation or Reynolds stress transportation equations. The applicable equation set consists of Favre-averaged conservation equations for the mass, momentum and total energy, and transport equations for the turbulent stresses and turbulent dissipation rate. In order to develop a scheme with good shock capturing capabilities, good accuracy and general geometric capabilities, a multi-block cell centered finite volume approach is used. Viscous fluxes are discretized using a finite volume representation of a central difference operator and the source terms are treated as an integral over the control volume. The methodology is validated by testing the algorithm on both two and three dimensional flows. Both the two equation and Reynolds stress models are used on a two dimensional 10 degree compression ramp at Mach 3, and the two equation model is used on the three dimensional flow over a cone at angle of attack at Mach 3.5. With the development of this algorithm, it is now possible to compute complex, compressible high speed flow fields using both two equation and Reynolds stress turbulent closure models, with the capability of eventually evaluating their predictive performance.
On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Ibraheem, S. O.; Demuren, A. O.
1994-01-01
A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.
Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Elmiligui, Alaa; Giriamaji, Sharath S.
2008-01-01
Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.
NASA Technical Reports Server (NTRS)
Rosenfeld, Moshe
1990-01-01
The development, validation and application of a fractional step solution method of the time-dependent incompressible Navier-Stokes equations in generalized coordinate systems are discussed. A solution method that combines a finite-volume discretization with a novel choice of the dependent variables and a fractional step splitting to obtain accurate solutions in arbitrary geometries was previously developed for fixed-grids. In the present research effort, this solution method is extended to include more general situations, including cases with moving grids. The numerical techniques are enhanced to gain efficiency and generality.
NASA Astrophysics Data System (ADS)
Lampson, Alan I.; Plummer, David N.; Erkkila, John H.; Crowell, Peter G.; Helms, Charles A.
1998-05-01
This paper describes a series of analyses using the 3-d MINT Navier-Stokes and OCELOT wave optics codes to calculate beam quality in a COIL laser cavity. To make this analysis tractable, the problem was broken into two contributions to the medium quality; that associated with microscale disturbances primarily from the transverse iodine injectors, and that associated with the macroscale including boundary layers and shock-like effects. Results for both microscale and macroscale medium quality are presented for the baseline layer operating point in terms of single pass wavefront error. These results show that the microscale optical path difference effects are 1D in nature and of low spatial order. The COIL medium quality is shown to be dominated by macroscale effects; primarily pressure waves generated from flow/boundary layer interactions on the cavity shrouds.
NASA Astrophysics Data System (ADS)
Nelson, Daniel A.; Jacobs, Gustaaf B.; Kopriva, David A.
2016-03-01
The effect of curved-boundary representation on the physics of the separated flow over a NACA 65(1)-412 airfoil is thoroughly investigated. A method is presented to approximate curved boundaries with a high-order discontinuous-Galerkin spectral element method for the solution of the Navier-Stokes equations. Multiblock quadrilateral element meshes are constructed with the grid generation software GridPro. The boundary of a NACA 65(1)-412 airfoil, defined by a cubic natural spline, is piecewise-approximated by isoparametric polynomial interpolants that represent the edges of boundary-fitted elements. Direct numerical simulation of the airfoil is performed on a coarse mesh and fine mesh with polynomial orders ranging from four to twelve. The accuracy of the curve fitting is investigated by comparing the flows computed on curved-sided meshes with those given by straight-sided meshes. Straight-sided meshes yield irregular wakes, whereas curved-sided meshes produce a regular Karman street wake. Straight-sided meshes also produce lower lift and higher viscous drag as compared with curved-sided meshes. When the mesh is refined by reducing the sizes of the elements, the lift decrease and viscous drag increase are less pronounced. The differences in the aerodynamic performance between the straight-sided meshes and the curved-sided meshes are concluded to be the result of artificial surface roughness introduced by the piecewise-linear boundary approximation provided by the straight-sided meshes.
NASA Astrophysics Data System (ADS)
Wu, X.; Vahdati, M.; Sayma, A.; Imregun, M.
2005-03-01
This paper describes a large-scale aeroelasticity computation for an aero-engine core compressor. The computational domain includes all 17 bladerows, resulting in a mesh with over 68 million points. The Favre-averaged Navier Stokes equations are used to represent the flow in a non-linear time-accurate fashion on unstructured meshes of mixed elements. The structural model of the first two rotor bladerows is based on a standard finite element representation. The fluid mesh is moved at each time step according to the structural motion so that changes in blade aerodynamic damping and flow unsteadiness can be accommodated automatically. An efficient domain decomposition technique, where special care was taken to balance the memory requirement across processors, was developed as part of the work. The calculation was conducted in parallel mode on 128 CPUs of an SGI Origin 3000. Ten vibration cycles were obtained using over 2.2 CPU years, though the elapsed time was a week only. Steady-state flow measurements and predictions were found to be in good agreement. A comparison of the averaged unsteady flow and the steady-state flow revealed some discrepancies. It was concluded that, in due course, the methodology would be adopted by industry to perform routine numerical simulations of the unsteady flow through entire compressor assemblies with vibrating blades not only to minimise engine and rig tests but also to improve performance predictions.
NASA Technical Reports Server (NTRS)
Ashford, Gregory A.; Powell, Kenneth G.
1995-01-01
A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.
NASA Astrophysics Data System (ADS)
Nelson, Daniel A.; Jacobs, Gustaaf B.; Kopriva, David A.
2016-08-01
The effect of curved-boundary representation on the physics of the separated flow over a NACA 65(1)-412 airfoil is thoroughly investigated. A method is presented to approximate curved boundaries with a high-order discontinuous-Galerkin spectral element method for the solution of the Navier-Stokes equations. Multiblock quadrilateral element meshes are constructed with the grid generation software GridPro. The boundary of a NACA 65(1)-412 airfoil, defined by a cubic natural spline, is piecewise-approximated by isoparametric polynomial interpolants that represent the edges of boundary-fitted elements. Direct numerical simulation of the airfoil is performed on a coarse mesh and fine mesh with polynomial orders ranging from four to twelve. The accuracy of the curve fitting is investigated by comparing the flows computed on curved-sided meshes with those given by straight-sided meshes. Straight-sided meshes yield irregular wakes, whereas curved-sided meshes produce a regular Karman street wake. Straight-sided meshes also produce lower lift and higher viscous drag as compared with curved-sided meshes. When the mesh is refined by reducing the sizes of the elements, the lift decrease and viscous drag increase are less pronounced. The differences in the aerodynamic performance between the straight-sided meshes and the curved-sided meshes are concluded to be the result of artificial surface roughness introduced by the piecewise-linear boundary approximation provided by the straight-sided meshes.
NASA Astrophysics Data System (ADS)
Ashford, Gregory A.; Powell, Kenneth G.
1995-10-01
A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.
NASA Astrophysics Data System (ADS)
Tsuzuki, Yutaka
2015-09-01
This paper is concerned with a system of heat equations with hysteresis and Navier-Stokes equations. In Tsuzuki (J Math Anal Appl 423:877-897, 2015) an existence result is obtained for the problem in a 2-dimensional domain with the Navier-Stokes equation in a weak sense. However the result does not include uniqueness for the problem due to the low regularity for solutions. This paper establishes existence and uniqueness in 2- and 3-dimensional domains with the Navier-Stokes equation in a stronger sense. Moreover this work decides required height of regularity for the initial data by introducing the fractional power of the Stokes operator.
NASA Astrophysics Data System (ADS)
Young, D. L.; Tsai, C. H.; Wu, C. S.
2015-11-01
An alternative vector potential formulation is used to solve the Navier-Stokes (N-S) equations in 3D incompressible viscous flow problems with and without through-flow boundaries. Difficulties of the vector potential formulation include the implementation of boundary conditions for through-flow boundaries and the numerical treatment of fourth-order partial differential equations. The advantages on the other hand are the automatic satisfaction of the continuity equation; and pressure is decoupled from the velocity. The objective of this paper is to introduce the appropriate gauge and boundary conditions on the vector potential formulation by a localized meshless method. To handle the divergence-free property, a Coulomb gauge condition is enforced on the vector potential to ensure its existence and uniqueness mathematically. We further improve the algorithm to through-flow problems for the boundary conditions of vector potential by introducing the concept of Stokes' theorem. Based on this innovation, there is no need to include an additional variable to tackle the through-flow fields. This process will greatly simplify the imposition of boundary conditions by the vector potential approach. Under certain conditions, the coupled fourth-order partial differential equations can be easily solved by using this meshless local differential quadrature (LDQ) method. Due to the LDQ capability to deal with the high order differential equations, this algorithm is very attractive to solve this fourth-order vector potential formulation for the N-S equations as comparing to the conventional numerical schemes such as finite element or finite difference methods. The proposed vector potential formulation is simpler and has improved accuracy and efficiency compared to other pressure-free or pressure-coupled algorithms. This investigation can be regarded as the first complete study to obtain the N-S solutions by vector potential formulation through a LDQ method. Two classic 3D benchmark
Coupling Boltzmann and Navier-Stokes equations by friction
Bourgat, J.F.; Le Tallec, P. |; Tidriri, M.D.
1996-09-01
The aim of this paper is to introduce and validate a coupled Navier-Stokes Boltzman approach for the calculation of hypersonic rarefied flows around maneuvering vehicles. The proposed strategy uses locally a kinetic model in the boundary layer coupled through wall friction forces to a global Navier-Stokes solver. Different numerical experiments illustrate the potentialities of the method. 29 refs., 24 figs.
Ge, Liang; Sotiropoulos, Fotis
2007-08-01
A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [1]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow through a mechanical, bileaflet heart valve mounted in a model straight aorta with an anatomical-like triple sinus. PMID:19194533
Ge, Liang; Sotiropoulos, Fotis
2008-01-01
A novel numerical method is developed that integrates boundary-conforming grids with a sharp interface, immersed boundary methodology. The method is intended for simulating internal flows containing complex, moving immersed boundaries such as those encountered in several cardiovascular applications. The background domain (e.g the empty aorta) is discretized efficiently with a curvilinear boundary-fitted mesh while the complex moving immersed boundary (say a prosthetic heart valve) is treated with the sharp-interface, hybrid Cartesian/immersed-boundary approach of Gilmanov and Sotiropoulos [1]. To facilitate the implementation of this novel modeling paradigm in complex flow simulations, an accurate and efficient numerical method is developed for solving the unsteady, incompressible Navier-Stokes equations in generalized curvilinear coordinates. The method employs a novel, fully-curvilinear staggered grid discretization approach, which does not require either the explicit evaluation of the Christoffel symbols or the discretization of all three momentum equations at cell interfaces as done in previous formulations. The equations are integrated in time using an efficient, second-order accurate fractional step methodology coupled with a Jacobian-free, Newton-Krylov solver for the momentum equations and a GMRES solver enhanced with multigrid as preconditioner for the Poisson equation. Several numerical experiments are carried out on fine computational meshes to demonstrate the accuracy and efficiency of the proposed method for standard benchmark problems as well as for unsteady, pulsatile flow through a curved, pipe bend. To demonstrate the ability of the method to simulate flows with complex, moving immersed boundaries we apply it to calculate pulsatile, physiological flow through a mechanical, bileaflet heart valve mounted in a model straight aorta with an anatomical-like triple sinus. PMID:19194533
NASA Astrophysics Data System (ADS)
Ren, Dandan; Ou, Yaobin
2016-08-01
In this paper, we prove the incompressible limit of all-time strong solutions to the three-dimensional full compressible Navier-Stokes equations. Here the velocity field and temperature satisfy the Dirichlet boundary condition and convective boundary condition, respectively. The uniform estimates in both the Mach number {ɛin(0,overline{ɛ}]} and time {tin[0,∞)} are established by deriving a differential inequality with decay property, where {overline{ɛ} in(0,1]} is a constant. Based on these uniform estimates, the global solution of full compressible Navier-Stokes equations with "well-prepared" initial conditions converges to the one of isentropic incompressible Navier-Stokes equations as the Mach number goes to zero.
Automatic differentiation and Navier-Stokes.
Bischof, C.; Hovland, P.; Mohammadi, B.
1997-12-17
We describe the use of automatic differentiation (AD) to enhance a compressible Navier-Stokes model. With the solver, AD is used to accelerate convergence by more than an order of magnitude. Outside the solver, AD is used to compute the derivatives needed for optimization. We emphasize the potential for performance gains if the programmer does not treat AD as a black box, but instead utilizes high-level knowledge about the nature of the application.
NASA Technical Reports Server (NTRS)
Kwak, D.
1994-01-01
INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far
NASA Technical Reports Server (NTRS)
Smith, Crawford F.; Podleski, Steve D.
1993-01-01
The proper use of a computational fluid dynamics code requires a good understanding of the particular code being applied. In this report the application of CFL3D, a thin-layer Navier-Stokes code, is compared with the results obtained from PARC3D, a full Navier-Stokes code. In order to gain an understanding of the use of this code, a simple problem was chosen in which several key features of the code could be exercised. The problem chosen is a cone in supersonic flow at an angle of attack. The issues of grid resolution, grid blocking, and multigridding with CFL3D are explored. The use of multigridding resulted in a significant reduction in the computational time required to solve the problem. Solutions obtained are compared with the results using the full Navier-Stokes equations solver PARC3D. The results obtained with the CFL3D code compared well with the PARC3D solutions.
Time-accurate Navier-Stokes calculations with multigrid acceleration
NASA Technical Reports Server (NTRS)
Melson, N. Duane; Atkins, Harold L.; Sanetrik, Mark D.
1993-01-01
A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.
NASA Astrophysics Data System (ADS)
DeJong, Andrew
Numerical models of fluid-structure interaction have grown in importance due to increasing interest in environmental energy harvesting, airfoil-gust interactions, and bio-inspired formation flying. Powered by increasingly powerful parallel computers, such models seek to explain the fundamental physics behind the complex, unsteady fluid-structure phenomena. To this end, a high-fidelity computational model based on the high-order spectral difference method on 3D unstructured, dynamic meshes has been developed. The spectral difference method constructs continuous solution fields within each element with a Riemann solver to compute the inviscid fluxes at the element interfaces and an averaging mechanism to compute the viscous fluxes. This method has shown promise in the past as a highly accurate, yet sufficiently fast method for solving unsteady viscous compressible flows. The solver is monolithically coupled to the equations of motion of an elastically mounted 3-degree of freedom rigid bluff body undergoing flow-induced lift, drag, and torque. The mesh is deformed using 4 methods: an analytic function, Laplace equation, biharmonic equation, and a bi-elliptic equation with variable diffusivity. This single system of equations -- fluid and structure -- is advanced through time using a 5-stage, 4th-order Runge-Kutta scheme. Message Passing Interface is used to run the coupled system in parallel on up to 240 processors. The solver is validated against previously published numerical and experimental data for an elastically mounted cylinder. The effect of adding an upstream body and inducing wake galloping is observed.
NASA Technical Reports Server (NTRS)
Cannizzaro, Frank E.; Ash, Robert L.
1992-01-01
A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, upwind numerical techniques, multigrid acceleration, and multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available: van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multi-grid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with multiple topologies. The results shown in this dissertation were chosen to validate the computer code and display its geometric flexibility, which is provided by the multi-block structure.
NASA Technical Reports Server (NTRS)
Rosen, Bruce S.
1991-01-01
An upwind three-dimensional volume Navier-Stokes code is modified to facilitate modeling of complex geometries and flow fields represented by proposed National Aerospace Plane concepts. Code enhancements include an equilibrium air model, a generalized equilibrium gas model and several schemes to simplify treatment of complex geometric configurations. The code is also restructured for inclusion of an arbitrary number of independent and dependent variables. This latter capability is intended for eventual use to incorporate nonequilibrium/chemistry gas models, more sophisticated turbulence and transition models, or other physical phenomena which will require inclusion of additional variables and/or governing equations. Comparisons of computed results with experimental data and results obtained using other methods are presented for code validation purposes. Good correlation is obtained for all of the test cases considered, indicating the success of the current effort.
Gehrer, A.; Jericha, H.
1999-07-01
External heat transfer predictions are performed for two-dimensional turbine blade cascades. The Reynolds-averaged Navier-Stokes equations with algebraic (Arnone and Pacciani, 1998), one-equation (Spalart and Allmaras, 1994), and two-equation (low-Re {kappa}-{epsilon}, Biswas and Fukuyama, 1994) turbulence closures are solved with a fully implicit time-marching finite volume method. Comparisons with measurements (Arts et al., 1990; Arts, 1994) for a highly loaded transonic turbine nozzle guide vane cascade show good agreement in some cases, but also reveal problems with transition prediction and turbulence modeling. Special attention has been focused on the low-Re {kappa}-{epsilon} model concerning the influence of the inlet boundary condition for the {epsilon}-equation and problems in the stagnation point region.
Reacting Multi-Species Gas Capability for USM3D Flow Solver
NASA Technical Reports Server (NTRS)
Frink, Neal T.; Schuster, David M.
2012-01-01
The USM3D Navier-Stokes flow solver contributed heavily to the NASA Constellation Project (CxP) as a highly productive computational tool for generating the aerodynamic databases for the Ares I and V launch vehicles and Orion launch abort vehicle (LAV). USM3D is currently limited to ideal-gas flows, which are not adequate for modeling the chemistry or temperature effects of hot-gas jet flows. This task was initiated to create an efficient implementation of multi-species gas and equilibrium chemistry into the USM3D code to improve its predictive capabilities for hot jet impingement effects. The goal of this NASA Engineering and Safety Center (NESC) assessment was to implement and validate a simulation capability to handle real-gas effects in the USM3D code. This document contains the outcome of the NESC assessment.
NASA Astrophysics Data System (ADS)
Desjouy, Cyril; Ollivier, Sébastien; Marsden, Olivier; Karzova, Maria; Blanc-Benon, Philippe
2016-02-01
The local interactions occurring between incident and reflected shock waves in the vicinity of rigid surfaces are investigated. Both regular and irregular — also called von Neumann — regimes of reflection are studied, via experimental and numerical simulations. Shock waves are produced experimentally with a 20 kV electrical spark source which allows the generation of spherically diverging acoustic shocks. The behaviour of the resulting weak acoustic shocks near rigid boundaries is visualized with a Schlieren optical technique which allows the spatial structure of the shocks to be studied. In particular, the evolution of the Mach stem forming above a flat surface is examined, and its height is observed to be directly linked to the angle of incidence and the pressure amplitude of the incident shock. The propagation of an acoustic shock between two parallel rigid boundaries is also studied. It is shown that the strong interactions between the Mach stems emerging from the two boundaries can lead to a drastic modification of the morphology of the acoustic field in the waveguide. Experimental results are compared to numerical results obtained from high-order finite-difference based simulations of the 2D Navier-Stokes equations. The good agreement between the experimental distribution of the acoustic field and numerical results suggests that numerical simulations are promising as a predictive tool to study nonlinear acoustic propagation of acoustic waves in complex geometrical configurations with rigid boundaries.
A multidimensional flux function with applications to the Euler and Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.; Van Leer, Bram; Roe, Philip L.
1993-01-01
In the present grid-independent approximate Riemann solver for 2D and 3D flows that are governed by the Euler or Navier-Stokes equations, fluxes on grid faces are obtained by wave decomposition; the assumption of information-propagation in the velocity-difference directions leads to a more accurate resolution of shear and shock waves, when these are are oblique to the grid. The model, which yields significantly greater accuracy in both supersonic and subsonic first-order spatially accurate computations, describes the difference in states at each grid interface by the action of five waves.
Navier-Stokes simulations of WECS airfoil flowfields
Homicz, G.F.
1994-06-01
Sandia National Laboratories has initiated an effort to apply Computational Fluid Dynamics (CFD) to the study of WECS aerodynamics. Preliminary calculations are presented for the flow past a SAND 0018/50 airfoil. The flow solver used is F3D, an implicitly, finite-difference code which solves the Thin-Layer Navier-airfoil. The flow solver used is F3D, an implicit, finite-difference code which solves the Thin-Layer Navier-Stokes equations. 2D steady-state calculations are presented at various angles of attack, {alpha}. Sectional lift and drag coefficient, as well as surface pressure distributions, are compared with wind tunnel data, and exhibit reasonable agreement at low to moderate angles of attack. At high {alpha}, where the airfoil is stalled, a converged solution to the steady-state equations could not be obtained. The flowfield continued to change with successive iterations, which is consistent with the fact that the actual flow is inherently transient, and requires the solution of the full unsteady form of the equations.
NASA Technical Reports Server (NTRS)
Marvriplis, D. J.; Venkatakrishnan, V.
1995-01-01
An agglomeration multigrid strategy is developed and implemented for the solution of three-dimensional steady viscous flows. The method enables convergence acceleration with minimal additional memory overheads, and is completely automated, in that it can deal with grids of arbitrary construction. The multigrid technique is validated by comparing the delivered convergence rates with those obtained by a previously developed overset-mesh multigrid approach, and by demonstrating grid independent convergence rates for aerodynamic problems on very large grids. Prospects for further increases in multigrid efficiency for high-Reynolds number viscous flows on highly stretched meshes are discussed.
A dual potential formulation of the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Gegg, S. G.; Pletcher, R. H.; Steger, J. L.
1989-01-01
A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.
Progress Toward Overset-Grid Moving Body Capability for USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Pandyna, Mohagna J.; Frink, Neal T.; Noack, Ralph W.
2005-01-01
A static and dynamic Chimera overset-grid capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. Modifications to the solver primarily consist of a few strategic calls to the Donor interpolation Receptor Transaction library (DiRTlib) to facilitate communication of solution information between various grids. The assembly of multiple overlapping grids into a single-zone composite grid is performed by the Structured, Unstructured and Generalized Grid AssembleR (SUGGAR) code. Several test cases are presented to verify the implementation, assess overset-grid solution accuracy and convergence relative to single-grid solutions, and demonstrate the prescribed relative grid motion capability.
Time-Accurate, Unstructured-Mesh Navier-Stokes Computations with the Space-Time CESE Method
NASA Technical Reports Server (NTRS)
Chang, Chau-Lyan
2006-01-01
Application of the newly emerged space-time conservation element solution element (CESE) method to compressible Navier-Stokes equations is studied. In contrast to Euler equations solvers, several issues such as boundary conditions, numerical dissipation, and grid stiffness warrant systematic investigations and validations. Non-reflecting boundary conditions applied at the truncated boundary are also investigated from the stand point of acoustic wave propagation. Validations of the numerical solutions are performed by comparing with exact solutions for steady-state as well as time-accurate viscous flow problems. The test cases cover a broad speed regime for problems ranging from acoustic wave propagation to 3D hypersonic configurations. Model problems pertinent to hypersonic configurations demonstrate the effectiveness of the CESE method in treating flows with shocks, unsteady waves, and separations. Good agreement with exact solutions suggests that the space-time CESE method provides a viable alternative for time-accurate Navier-Stokes calculations of a broad range of problems.
Navier Stokes Theorem in Hydrology
NASA Astrophysics Data System (ADS)
Narayanan, M.
2005-12-01
In a paper presented at the 2004 AGU International Conference, the author outlined and stressed the importance of studying and teaching certain important mathematical techniques while developing a course in Hydrology and Fluid Mechanics. The Navier-Stokes equations are the foundation of fluid mechanics, and Stokes' theorem is used in nearly every branch of mechanics as well as electromagnetics. Stokes' Theorem also plays a vital role in many secondary theorems such as those pertaining to vorticity and circulation. Mathematically expressed, Stokes' theorem can be expressed by considering a surface S having a bounding curve C. Here, V is any sufficiently smooth vector field defined on the surface and its bounding curve C. In an article entitled "Corrections to Fluid Dynamics" R. F. Streater, (Open Systems and Information Dynamics, 10, 3-30, 2003.) proposes a kinetic model of a fluid in which five macroscopic fields, the mass, energy, and three components of momentum, are conserved. The dynamics is constructed using the methods of statistical dynamics, and results in a non-linear discrete-time Markov chain for random fields on a lattice. In the continuum limit he obtains a non-linear coupled parabolic system of field equations, showing a correction to the Navier-Stokes equations. In 2001, David Hoff published an article in Journees Equations aux derivees partielles. (Art. No. 7, 9 p.). His paper is entitled : Dynamics of Singularity Surfaces for Compressible Navier-Stokes Flows in Two Space Dimensions. In his paper, David Hoff proves the global existence of solutions of the Navier-Stokes equations of compressible, barotropic flow in two space dimensions with piecewise smooth initial data. These solutions remain piecewise smooth for all time, retaining simple jump discontinuities in the density and in the divergence of the velocity across a smooth curve, which is convected with the flow. The strengths of these discontinuities are shown to decay exponentially in time
Turbomachinery blade optimization using the Navier-Stokes equations
Chand, K.K.; Lee, K.D.
1997-12-01
A method is presented to perform aerodynamic design optimization of turbomachinery blades. The method couples a Navier-Stokes flow solver with a grid generator and numerical optimization algorithm to seek improved designs for transonic turbine blades. A fast and efficient multigrid, finite-volume flow solver provides accurate performance evaluations of potential designs. Design variables consist of smooth perturbations to the blade surface. A unique elliptic-hyperbolic grid generation method is used to regenerate a Navier-Stokes grid after perturbations have been added to the geometry. Designs are sought which improve a design objective while remaining within specified constraints. The method is demonstrated with two transonic turbine blades with different types and numbers of design variables.
NASA Technical Reports Server (NTRS)
Anderson, B. H.; Reddy, D. R.; Kapoor, K.
1993-01-01
A three-dimensional implicit Full Navier-Stokes (FNS) analysis and a 3D Reduced Navier-Stokes (RNS) initial value space marching solution technique has been applied to a class of separate flow problems within a diffusing S-duct configuration characterized as vortex-liftoff. Both Full Navier-Stokes and Reduced Navier-Stokes solution techniques were able to capture the overall flow physics of vortex lift-off, however more consideration must be given to the development of turbulence models for the prediction of the locations of separation and reattachment. This accounts for some of the discrepancies in the prediction of the relevant inlet distortion descriptors, particularly circumferential distortion. The 3D RNS solution technique adequately described the topological structure of flow separation associated with vortex lift-off.
NASA Astrophysics Data System (ADS)
Sahin, Mehmet
2005-11-01
A new semi-staggered finite volume method is presented for the solution of the incompressible Navier-Stokes equations on all-quadrilateral (2D)/hexahedral (3D) meshes. The velocity components are defined at element node points while the pressure term is defined at element centroids. The continuity equation is satisfied exactly within each elements. The checkerboard pressure oscillations are prevented using a special filtering matrix as a preconditioner for the saddle-point problem resulting from second-order discretization of the incompressible Navier-Stokes equations. The preconditioned saddle-point problem is solved using block preconditioners with GMRES solver. In order to achieve higher performance FORTRAN source code is based on highly efficient PETSc and HYPRE libraries. As test cases the 2D/3D lid-driven cavity flow problem and the 3D flow past array of circular cylinders are solved in order to verify the accuracy of the proposed method.
A Continuation and Bifurcation Technique for Navier-Stokes Flows
NASA Astrophysics Data System (ADS)
Sanchez, J.; Marques, F.; Lopez, J. M.
2002-07-01
An efficient numerical bifurcation and continuation method for the Navier-Stokes equations in cylindrical geometries is presented and applied to a nontrivial fluid dynamics problem, the flow in a cylindrical container driven by differential rotation. The large systems that result from discretizing the Navier-Stokes equations, especially in regimes where inertia is important, necessitate the use of iterative solvers which in turn need preconditioners. We use incomplete lower-upper decomposition (ILU) as an effective preconditioner for such systems and show the significant gain in efficiency when an incomplete LU of the full Jacobian is used instead of using only the Stokes operator. The computational cost, in terms of CPU time, grows with the size of the system (i.e., spatial resolution) according to a power law with exponent around 1.7, which is very modest compared to direct methods, indicating the appropriateness of the schemes for large nonlinear partial differential equation problems.
NASA Astrophysics Data System (ADS)
Kim, Bong-Sik
Three dimensional (3D) Navier-Stokes-alpha equations are considered for uniformly rotating geophysical fluid flows (large Coriolis parameter f = 2O). The Navier-Stokes-alpha equations are a nonlinear dispersive regularization of usual Navier-Stokes equations obtained by Lagrangian averaging. The focus is on the existence and global regularity of solutions of the 3D rotating Navier-Stokes-alpha equations and the uniform convergence of these solutions to those of the original 3D rotating Navier-Stokes equations for large Coriolis parameters f as alpha → 0. Methods are based on fast singular oscillating limits and results are obtained for periodic boundary conditions for all domain aspect ratios, including the case of three wave resonances which yields nonlinear "2½-dimensional" limit resonant equations for f → 0. The existence and global regularity of solutions of limit resonant equations is established, uniformly in alpha. Bootstrapping from global regularity of the limit equations, the existence of a regular solution of the full 3D rotating Navier-Stokes-alpha equations for large f for an infinite time is established. Then, the uniform convergence of a regular solution of the 3D rotating Navier-Stokes-alpha equations (alpha ≠ 0) to the one of the original 3D rotating NavierStokes equations (alpha = 0) for f large but fixed as alpha → 0 follows; this implies "shadowing" of trajectories of the limit dynamical systems by those of the perturbed alpha-dynamical systems. All the estimates are uniform in alpha, in contrast with previous estimates in the literature which blow up as alpha → 0. Finally, the existence of global attractors as well as exponential attractors is established for large f and the estimates are uniform in alpha.
Transonic airfoil and wing design using Navier-Stokes codes
NASA Technical Reports Server (NTRS)
Yu, N. J.; Campbell, R. L.
1992-01-01
An iterative design method has been implemented into 2D and 3D Navier-Stokes codes for the design of airfoils or wings with given target pressure distributions. The method begins with the analysis of an initial geometry, and obtains the analysis pressure distributions of that geometry. The differences between analysis pressures and target pressures are used to drive geometry changes through the use of a streamline curvature method. This paper describes the procedure that makes the iterative design method work for Navier-Stokes codes. Examples of 2D airfoil design, and 3D wing design are included. It is demonstrated that the method is highly effective for airfoil or wing design at flow conditions where no substantial separation occurs. Problems encountered in the airfoil design with shock induced flow separations are discussed.
NASA Technical Reports Server (NTRS)
Anderson, B. H.; Reddy, D. R.; Kapoor, K.
1993-01-01
A three-dimensional implicit Full Navier-Stokes (FNS) analysis and a 3D Reduced Navier Stokes (RNS) initial value space marching solution technique has been applied to a class of separated flow problems within a diffusing S-duct configuration characterized by vortex-liftoff. Both the FNS and the RNS solution technique were able to capture the overall flow physics of vortex lift-off, and gave remarkably similar results which agreed reasonably well with the experimental measured averaged performance parameters of engine face total pressure recovery and distortion. However, the Full Navier-Stokes and Reduced Navier-Stokes also consistently predicted separation further downstream in the M2129 inlet S-duct than was indicated by experimental data, thus compensating errors were present in the two Navier-Stokes analyses. The difficulties encountered in the Navier-Stokes separations analyses of the M2129 inlet S-duct center primarily on turbulence model issues, and these focused on two distinct but different phenomena, namely, (1) characterization of low skin friction adverse pressure gradient flows, and (2) description of the near wall behavior of flows characterized by vortex lift-off.
Preconditioned conjugate gradient methods for the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Ajmani, Kumud; Ng, Wing-Fai; Liou, Meng-Sing
1994-01-01
A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.
Finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.
1986-01-01
The space shuttle main engine (SSME) has extremely complex internal flow structure. The geometry of the flow domain is three-dimensional with complicated topology. The flow is compressible, viscous, and turbulent with large gradients in flow quantities and regions of recirculations. The analysis of the flow field in SSME involves several tedious steps. One is the geometrical modeling of the particular zone of the SSME being studied. Accessing the geometry definition, digitalizing it, and developing surface interpolations suitable for an interior grid generator require considerable amount of manual labor. There are several types of grid generators available with some general-purpose finite element programs. An efficient and robust computational scheme for solving 3D Navier-Stokes equations has to be implemented. Post processing software has to be adapted to visualize and analyze the computed 3D flow field. The progress made in a project to develop software for the analysis of the flow is discussed. The technical approach to the development of the finite element scheme and the relaxation procedure are discussed. The three dimensional finite element code for the compressible Navier-Stokes equations is listed.
A Navier-Stokes boundary element solver
NASA Technical Reports Server (NTRS)
Reddy, D. R.; Lafe, O.; Cheng, A. H-D.
1995-01-01
Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.
Navier-Stokes Computations on Commodity Computers
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Faulkner, Thomas R.
1998-01-01
In this paper we discuss and demonstrate the feasibility of solving high-fidelity, nonlinear computational fluid dynamics (CFD) problems of practical interest on commodity machines, namely Pentium Pro PC's. Such calculations have now become possible due to the progress in computational power and memory of the off-the-shelf commodity computers, along with the growth in bandwidth and communication speeds of networks. A widely used CFD code known as TLNS3D, which was developed originally on large shared memory computers was selected for this effort. This code has recently been ported to massively parallel processor (MPP) type machines, where natural partitioning along grid blocks is adopted in which one or more blocks are distributed to each of the available processors. In this paper, a similar approach is adapted to port this code to a cluster of Pentium Pro computers. The message passing among the processors is accomplished through the use of standard message passing interface (MPI) libraries. Scaling studies indicate fairly high level of parallelism on such clusters of commodity machines, thus making solutions to Navier-Stokes equations for practical problems more affordable.
Preconditioned conjugate gradient methods for the Navier-Stokes equations
Ajmani, K.; Ng, Wing Fai ); Liou, Meng Sing )
1994-01-01
A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulations. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU times on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow cases are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner. 34 refs., 15 figs.
NASA Astrophysics Data System (ADS)
Kandula, M.; Pearce, D. G.
1991-06-01
A steady incompressible three-dimensional viscous flow analysis has been conducted for the Space Shuttle external tank/orbiter propellant feed line disconnect flapper valves with upstream elbows. The Navier-Stokes code, INS3D, is modified to handle interior obstacles and a simple turbulence model. The flow solver is tested for stability and convergence in the presence of interior flappers. An under-relaxation scheme has been incorporated to improve the solution stability. Important flow characteristics such as secondary flows, recirculation, vortex and wake regions, and separated flows are observed. Computed values for forces, moments, and pressure drop are in satisfactory agreement with water flow test data covering a maximum tube Reynolds number of 3.5 million. The predicted hydrodynamical stability of the flappers correlates well with the measurements.
Transonic Navier-Stokes calculations about a 65 deg delta wing
NASA Technical Reports Server (NTRS)
Londenberg, W. Kelly
1994-01-01
A computational study has been conducted in which the CFL3D Navier-Stokes solver coupled with an algebraic and a one-equation nonequilibrium turbulence model has been used to predict the flow over a 65 degree delta wing at transonic conditions for Reynolds numbers ranging from 6 x 10(exp 6) to 120 x 10(exp 6) based on mean aerodynamic chord. Solutions obtained indicated that the computational method when used with the one-equation turbulence model predicts results that compare well with experiment for attached flow conditions. Comparisons with experimental pressure at separated conditions show that the computational method, even though primary flow-field features are predicted well, does not predict secondary flow features.
NASA Astrophysics Data System (ADS)
Stück, Arthur
2015-11-01
Inconsistent discrete expressions in the boundary treatment of Navier-Stokes solvers and in the definition of force objective functionals can lead to discrete-adjoint boundary treatments that are not a valid representation of the boundary conditions to the corresponding adjoint partial differential equations. The underlying problem is studied for an elementary 1D advection-diffusion problem first using a node-centred finite-volume discretisation. The defect of the boundary operators in the inconsistently defined discrete-adjoint problem leads to oscillations and becomes evident with the additional insight of the continuous-adjoint approach. A homogenisation of the discretisations for the primal boundary treatment and the force objective functional yields second-order functional accuracy and eliminates the defect in the discrete-adjoint boundary treatment. Subsequently, the issue is studied for aerodynamic Reynolds-averaged Navier-Stokes problems in conjunction with a standard finite-volume discretisation on median-dual grids and a strong implementation of noslip walls, found in many unstructured general-purpose flow solvers. Going out from a base-line discretisation of force objective functionals which is independent of the boundary treatment in the flow solver, two improved flux-consistent schemes are presented; based on either body wall-defined or farfield-defined control-volumes they resolve the dual inconsistency. The behaviour of the schemes is investigated on a sequence of grids in 2D and 3D.
About the Regularized Navier Stokes Equations
NASA Astrophysics Data System (ADS)
Cannone, Marco; Karch, Grzegorz
2005-03-01
The first goal of this paper is to study the large time behavior of solutions to the Cauchy problem for the 3-dimensional incompressible Navier Stokes system. The Marcinkiewicz space L3,∞ is used to prove some asymptotic stability results for solutions with infinite energy. Next, this approach is applied to the analysis of two classical “regularized” Navier Stokes systems. The first one was introduced by J. Leray and consists in “mollifying” the nonlinearity. The second one was proposed by J.-L. Lions, who added the artificial hyper-viscosity (-Δ)ℓ/ 2, ℓ > 2 to the model. It is shown in the present paper that, in the whole space, solutions to those modified models converge as t → ∞ toward solutions of the original Navier Stokes system.
A gas-kinetic BGK scheme for the compressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Xu, Kun
2000-01-01
This paper presents an improved gas-kinetic scheme based on the Bhatnagar-Gross-Krook (BGK) model for the compressible Navier-Stokes equations. The current method extends the previous gas-kinetic Navier-Stokes solver developed by Xu and Prendergast by implementing a general nonequilibrium state to represent the gas distribution function at the beginning of each time step. As a result, the requirement in the previous scheme, such as the particle collision time being less than the time step for the validity of the BGK Navier-Stokes solution, is removed. Therefore, the applicable regime of the current method is much enlarged and the Navier-Stokes solution can be obtained accurately regardless of the ratio between the collision time and the time step. The gas-kinetic Navier-Stokes solver developed by Chou and Baganoff is the limiting case of the current method, and it is valid only under such a limiting condition. Also, in this paper, the appropriate implementation of boundary condition for the kinetic scheme, different kinetic limiting cases, and the Prandtl number fix are presented. The connection among artificial dissipative central schemes, Godunov-type schemes, and the gas-kinetic BGK method is discussed. Many numerical tests are included to validate the current method.
Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Abdol-Hamid, Khaled S.; Campbell, Richard L.; Frink, Neal T.
2006-01-01
A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-epsilon two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user-specified region. The latter approach is applied for the solutions obtained using other one-and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.
Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Abdol-Harrid, Khaled S.; Campbell, Richard L.; Frink, Neal T.
2006-01-01
A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-e two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user specified region. The latter approach is applied for the solutions obtained using other one- and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.
LINFLUX-AE: A Turbomachinery Aeroelastic Code Based on a 3-D Linearized Euler Solver
NASA Technical Reports Server (NTRS)
Reddy, T. S. R.; Bakhle, M. A.; Trudell, J. J.; Mehmed, O.; Stefko, G. L.
2004-01-01
This report describes the development and validation of LINFLUX-AE, a turbomachinery aeroelastic code based on the linearized unsteady 3-D Euler solver, LINFLUX. A helical fan with flat plate geometry is selected as the test case for numerical validation. The steady solution required by LINFLUX is obtained from the nonlinear Euler/Navier Stokes solver TURBO-AE. The report briefly describes the salient features of LINFLUX and the details of the aeroelastic extension. The aeroelastic formulation is based on a modal approach. An eigenvalue formulation is used for flutter analysis. The unsteady aerodynamic forces required for flutter are obtained by running LINFLUX for each mode, interblade phase angle and frequency of interest. The unsteady aerodynamic forces for forced response analysis are obtained from LINFLUX for the prescribed excitation, interblade phase angle, and frequency. The forced response amplitude is calculated from the modal summation of the generalized displacements. The unsteady pressures, work done per cycle, eigenvalues and forced response amplitudes obtained from LINFLUX are compared with those obtained from LINSUB, TURBO-AE, ASTROP2, and ANSYS.
Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces
NASA Technical Reports Server (NTRS)
Wang, Chi R.
2005-01-01
This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one
Navier-Stokes analysis of radial turbine rotor performance
NASA Technical Reports Server (NTRS)
Larosiliere, L. M.
1993-01-01
An analysis of flow through a radial turbine rotor using the three-dimensional, thin-layer Navier-Stokes code RVC3D is described. The rotor is a solid version of an air-cooled metallic radial turbine having thick trailing edges, shroud clearance, and scalloped-backface clearance. Results are presented at the nominal operating condition using both a zero-clearance model and a model simulating the effects of the shroud and scalloped-backface clearance flows. A comparison with the available test data is made and details of the internal flow physics are discussed, allowing a better understanding of the complex flow distribution within the rotor.
Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Abdol-Hamid, Khaled S.; Frink, Neal T.
2009-01-01
Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.
An efficient non-linear multigrid procedure for the incompressible Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Sivaloganathan, S.; Shaw, G. J.
An efficient Full Approximation multigrid scheme for finite volume discretizations of the Navier-Stokes equations is presented. The algorithm is applied to the driven cavity test problem. Numerical results are presented and a comparison made with PACE, a Rolls-Royce industrial code, which uses the SIMPLE pressure correction method as an iterative solver.
NASA Technical Reports Server (NTRS)
Bui, Trong T.
1993-01-01
New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.
NASA Technical Reports Server (NTRS)
Bui, Trong T.
1993-01-01
New turbulence modeling options recently implemented for the 3D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Good agreements are obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(+) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. Test cases show that the highly optimized one- and two-equation turbulence models can be used in routine 3D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.
NASA Technical Reports Server (NTRS)
Mizukami, M.; Saunders, J. D.
1995-01-01
The supersonic diffuser of a Mach 2.68 bifurcated, rectangular, mixed-compression inlet was analyzed using a two-dimensional (2D) Navier-Stokes flow solver. Parametric studies were performed on turbulence models, computational grids and bleed models. The computer flowfield was substantially different from the original inviscid design, due to interactions of shocks, boundary layers, and bleed. Good agreement with experimental data was obtained in many aspects. Many of the discrepancies were thought to originate primarily from 3D effects. Therefore, a balance should be struck between expending resources on a high fidelity 2D simulation, and the inherent limitations of 2D analysis. The solutions were fairly insensitive to turbulence models, grids and bleed models. Overall, the k-e turbulence model, and the bleed models based on unchoked bleed hole discharge coefficients or uniform velocity are recommended. The 2D Navier-Stokes methods appear to be a useful tool for the design and analysis of supersonic inlets, by providing a higher fidelity simulation of the inlet flowfield than inviscid methods, in a reasonable turnaround time.
Navier-Stokes analysis of turbomachinery blade external heat transfer
NASA Astrophysics Data System (ADS)
Ameri, A. A.; Sockol, P. M.; Gorla, R. S. R.
1992-04-01
The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q - omega low Reynolds number, two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.
Navier-Stokes analysis of turbomachinery blade external heat transfer
NASA Technical Reports Server (NTRS)
Gorla, Rama S. R.
1991-01-01
The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q-omega low Reynolds number two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.
Towards Optimal Multigrid Efficiency for the Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Swanson, R. C.
2001-01-01
A fast multigrid solver for the steady incompressible Navier-Stokes equations is presented. Unlike time-marching schemes, this approach uses relaxation of the steady equations. Application of this method results in a discretization that correctly distinguishes between the advection and elliptic parts of the operator, allowing efficient smoothers to be constructed. Numerical solutions are shown for flow over a flat plate and a Karman-Trefftz airfoil. Using collective Gauss-Seidel line relaxation in both the vertical and horizontal directions, multigrid convergence behavior approaching that of O(N) methods is achieved. The computational efficiency of the numerical scheme is compared with that of a Runge-Kutta based multigrid method.
A finite element solver for 3-D compressible viscous flows
NASA Technical Reports Server (NTRS)
Reddy, K. C.; Reddy, J. N.; Nayani, S.
1990-01-01
Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.
NASA Technical Reports Server (NTRS)
Yoon, Seokkwan; Chang, Leon; Kwak, Dochan
1989-01-01
A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.
Evaluation of a Multigrid Scheme for the Incompressible Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Swanson, R. C.
2004-01-01
A fast multigrid solver for the steady, incompressible Navier-Stokes equations is presented. The multigrid solver is based upon a factorizable discrete scheme for the velocity-pressure form of the Navier-Stokes equations. This scheme correctly distinguishes between the advection-diffusion and elliptic parts of the operator, allowing efficient smoothers to be constructed. To evaluate the multigrid algorithm, solutions are computed for flow over a flat plate, parabola, and a Karman-Trefftz airfoil. Both nonlifting and lifting airfoil flows are considered, with a Reynolds number range of 200 to 800. Convergence and accuracy of the algorithm are discussed. Using Gauss-Seidel line relaxation in alternating directions, multigrid convergence behavior approaching that of O(N) methods is achieved. The computational efficiency of the numerical scheme is compared with that of Runge-Kutta and implicit upwind based multigrid methods.
NASA Technical Reports Server (NTRS)
Chen, Y. K.; Henline, W. D.
1993-01-01
The general boundary conditions including mass and energy balances of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver for predictions of the flow field, surface temperature, and surface ablation rates over re-entry space vehicles with ablating Thermal Protection Systems (TPS). The Navier-Stokes solver with general surface thermochemistry boundary conditions can predict more realistic solutions and provide useful information for the design of TPS. A test case with a proposed hypersonic test vehicle configuration and associated free stream conditions was developed. Solutions with various surface boundary conditions were obtained, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined. The solutions of the GASP code with complete ablating surface conditions were compared with those of the ASC code. The direction of future work is also discussed.
NASA Technical Reports Server (NTRS)
Mcmillin, S. Naomi; Thomas, James L.; Murman, Earll M.
1990-01-01
An Euler flow solver and a thin layer Navier-Stokes flow solver were used to numerically simulate the supersonic leeside flow fields over delta wings which were observed experimentally. Three delta wings with 75, 67.5, and 60 deg leading edge sweeps were computed over an angle-of-attack range of 4 to 20 deg at a Mach number 2.8. The Euler code and Navier-Stokes code predict equally well the primary flow structure where the flow is expected to be separated or attached at the leading edge based on the Stanbrook-Squire boundary. The Navier-Stokes code is capable of predicting both the primary and the secondary flow features for the parameter range investigated. For those flow conditions where the Euler code did not predict the correct type of primary flow structure, the Navier-Stokes code illustrated that the flow structure is sensitive to boundary layer model. In general, the laminar Navier-Stokes solutions agreed better with the experimental data, especially for the lower sweep delta wings. The computational results and a detailed re-examination of the experimental data resulted in a refinement of the flow classifications. This refinement in the flow classification results in the separation bubble with the shock flow type as the intermediate flow pattern between separated and attached flows.
Automated Euler and Navier-Stokes Database Generation for a Glide-Back Booster
NASA Technical Reports Server (NTRS)
Chaderjian, Neal M.; Rogers, Stuart E.; Aftosmis, Mike J.; Pandya, Shishir A.; Ahmad, Jasim U.; Tejnil, Edward
2004-01-01
The past two decades have seen a sustained increase in the use of high fidelity Computational Fluid Dynamics (CFD) in basic research, aircraft design, and the analysis of post-design issues. As the fidelity of a CFD method increases, the number of cases that can be readily and affordably computed greatly diminishes. However, computer speeds now exceed 2 GHz, hundreds of processors are currently available and more affordable, and advances in parallel CFD algorithms scale more readily with large numbers of processors. All of these factors make it feasible to compute thousands of high fidelity cases. However, there still remains the overwhelming task of monitoring the solution process. This paper presents an approach to automate the CFD solution process. A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment, the NASA Information Power Grid (IPG), using 13 computers located at 4 different geographical sites. Process automation and web-based access to a MySql database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The AeroDB framework is shown. The user submits/deletes jobs, monitors AeroDB's progress, and retrieves data and plots via a web portal. Once a job is in the database, a job launcher uses an IPG resource broker to decide which computers are best suited to run the job. Job/code requirements, the number of CPUs free on a remote system, and queue lengths are some of the parameters the broker takes into account. The Globus software provides secure services for user authentication, remote shell execution, and secure file transfers over an open network. AeroDB automatically decides when a job is completed. Currently, the Cart3D unstructured flow solver is used for the Euler equations, and the Overflow structured overset flow solver is used for the
NASA Technical Reports Server (NTRS)
Kiris, Cetin
1995-01-01
Development of an incompressible Navier-Stokes solution procedure was performed for the analysis of a liquid rocket engine pump components and for the mechanical heart assist devices. The solution procedure for the propulsion systems is applicable to incompressible Navier-Stokes flows in a steadily rotating frame of reference for any general complex configurations. The computer codes were tested on different complex configurations such as liquid rocket engine inducer and impellers. As a spin-off technology from the turbopump component simulations, the flow analysis for an axial heart pump was conducted. The baseline Left Ventricular Assist Device (LVAD) design was improved by adding an inducer geometry by adapting from the liquid rocket engine pump. The time-accurate mode of the incompressible Navier-Stokes code was validated with flapping foil experiment by using different domain decomposition methods. In the flapping foil experiment, two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions for a downstream larger stationary foil. Fairly good agreement was obtained between unsteady experimental data and numerical results from two different moving boundary procedures. Incompressible Navier-Stokes code (INS3D) has been extended for heat transfer applications. The temperature equation was written for both forced and natural convection phenomena. Flow in a square duct case was used for the validation of the code in both natural and forced convection.
The Proteus Navier-Stokes code
NASA Technical Reports Server (NTRS)
Towne, Charles E.; Bui, Trong T.; Cavicchi, Richard H.; Conley, Julianne M.; Molls, Frank B.; Schwab, John R.
1992-01-01
An effort is currently underway at NASA Lewis to develop two- and three-dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. The emphasis in the development of Proteus is not algorithm development or research on numerical methods, but rather the development of the code itself. The objective is to develop codes that are user-oriented, easily-modified, and well-documented. Well-proven, state-of-the-art solution algorithms are being used. Code readability, documentation (both internal and external), and validation are being emphasized. This paper is a status report on the Proteus development effort. The analysis and solution procedure are described briefly, and the various features in the code are summarized. The results from some of the validation cases that have been run are presented for both the two- and three-dimensional codes.
Reynolds-Averaged Navier-Stokes Simulation of a 2D Circulation Control Wind Tunnel Experiment
NASA Technical Reports Server (NTRS)
Allan, Brian G.; Jones, Greg; Lin, John C.
2011-01-01
Numerical simulations are performed using a Reynolds-averaged Navier-Stokes (RANS) flow solver for a circulation control airfoil. 2D and 3D simulation results are compared to a circulation control wind tunnel test conducted at the NASA Langley Basic Aerodynamics Research Tunnel (BART). The RANS simulations are compared to a low blowing case with a jet momentum coefficient, C(sub u), of 0:047 and a higher blowing case of 0.115. Three dimensional simulations of the model and tunnel walls show wall effects on the lift and airfoil surface pressures. These wall effects include a 4% decrease of the midspan sectional lift for the C(sub u) 0.115 blowing condition. Simulations comparing the performance of the Spalart Allmaras (SA) and Shear Stress Transport (SST) turbulence models are also made, showing the SST model compares best to the experimental data. A Rotational/Curvature Correction (RCC) to the turbulence model is also evaluated demonstrating an improvement in the CFD predictions.
Prediction of Business Jet Airloads Using The Overflow Navier-Stokes Code
NASA Technical Reports Server (NTRS)
Bounajem, Elias; Buning, Pieter G.
2001-01-01
The objective of this work is to evaluate the application of Navier-Stokes computational fluid dynamics technology, for the purpose of predicting off-design condition airloads on a business jet configuration in the transonic regime. The NASA Navier-Stokes flow solver OVERFLOW with Chimera overset grid capability, availability of several numerical schemes and convergence acceleration techniques was selected for this work. A set of scripts which have been compiled to reduce the time required for the grid generation process are described. Several turbulence models are evaluated in the presence of separated flow regions on the wing. Computed results are compared to available wind tunnel data for two Mach numbers and a range of angles-of-attack. Comparisons of wing surface pressure from numerical simulation and wind tunnel measurements show good agreement up to fairly high angles-of-attack.
NASA Technical Reports Server (NTRS)
Fletcher, Michael J.; Won, Mark J.; Cosentino, Gary B.; Te, Alexander
1993-01-01
Subsonic inlet ducts for advanced, high-performance aircraft are evolving towards complex three-dimensional shapes for reasons of overall integration and weight. These factors lead to diffuser geometries that may sacrifice inlet performance, unless careful attention to design details and boundary layer management techniques are employed. The ability of viscous computational fluid dynamic (CFD) analysis of such geometries to aid the aircraft configurator in this complex design problem is herein examined. The RANS-3D Reynolds-Averaged Navier-Stokes solver is applied to model the complex flowfield occurring in a representative diffuser geometry and the solutions are compared to experimental results from a static test of the inlet duct. The computational results are shown to compare very favorably with experimental results over a range of mass flow rates, including those involving large amounts of separation in the diffuser. In addition, a novel grid topology is presented, and two turbulence models are evaluated in this study as part of the RANS-3D code.
Hypersonic Navier Stokes Comparisons to Orbiter Flight Data
NASA Technical Reports Server (NTRS)
Campbell, Charles H.; Nompelis, Ioannis; Candler, Graham; Barnhart, Michael; Yoon, Seokkwan
2009-01-01
Hypersonic chemical nonequilibrium simulations of low earth orbit entry flow fields are becoming increasingly commonplace as software and computational capabilities become more capable. However, development of robust and accurate software to model these environments will always encounter a significant barrier in developing a suite of high quality calibration cases. The US3D hypersonic nonequilibrium Navier Stokes analysis capability has been favorably compared to a number of wind tunnel test cases. Extension of the calibration basis for this software to Orbiter flight conditions will provide an incremental increase in confidence. As part of the Orbiter Boundary Layer Transition Flight Experiment and the Hypersonic Thermodynamic Infrared Measurements project, NASA is performing entry flight testing on the Orbiter to provide valuable aerothermodynamic heating data. An increase in interest related to orbiter entry environments is resulting from this activity. With the advent of this new data, comparisons of the US3D software to the new flight testing data is warranted. This paper will provide information regarding the framework of analyses that will be applied with the US3D analysis tool. In addition, comparisons will be made to entry flight testing data provided by the Orbiter BLT Flight Experiment and HYTHIRM projects. If data from digital scans of the Orbiter windward surface become available, simulations will also be performed to characterize the difference in surface heating between the CAD reference OML and the digitized surface provided by the surface scans.
A three-dimensional structured/unstructured hybrid Navier-Stokes method for turbine blade rows
NASA Astrophysics Data System (ADS)
Tsung, F.-L.; Loellbach, J.; Kwon, O.; Hah, C.
1994-12-01
A three-dimensional viscous structured/unstructured hybrid scheme has been developed for numerical computation of high Reynolds number turbomachinery flows. The procedure allows an efficient structured solver to be employed in the densely clustered, high aspect-ratio grid around the viscous regions near solid surfaces, while employing an unstructured solver elsewhere in the flow domain to add flexibility in mesh generation. Test results for an inviscid flow over an external transonic wing and a Navier-Stokes flow for an internal annular cascade are presented.
Chaos Synchronization in Navier-Stokes Turbulence
NASA Astrophysics Data System (ADS)
Lalescu, Cristian; Meneveau, Charles; Eyink, Gregory
2013-03-01
Chaos synchronization (CS) has been studied for some time now (Pecora & Carroll 1990), for systems with only a few degrees of freedom as well as for systems described by partial differential equations (Boccaletti et al 2002). CS in general is said to be present in coupled dynamical systems when a specific property of each system has the same time evolution for all, even though the evolution itself is chaotic. The Navier-Stokes (NS) equations describe the velocity for a wide range of fluids, and their solutions are usually called turbulent if fluctuation amplitudes decrease as a power of their wavenumber. There have been some studies of CS for continuous systems (Kocarev et al 1997), but CS for NS turbulence seems not to have been investigated so far. We focus on the synchronization of the small scales of a turbulent flow for which the time history of large scales is prescribed. Our DNS results show that high-wavenumbers in turbulence are fully slaved to modes with wavenumbers up to a critical fraction of the Kolmogorov dissipation wavenumber. The motivation for our work is to study deeply sub-Kolmogorov scales in fully developed turbulence (Schumacher 2007), which we found to be recoverable even at very high Reynolds number from simulations with moderate resolutions. This work is supported by the National Science Foundation's CDI-II program, project CMMI-0941530
Chaos Synchronization in Navier-Stokes Turbulence
NASA Astrophysics Data System (ADS)
Lalescu, Cristian C.; Meneveau, Charles; Eyink, Gregory L.
2012-11-01
Chaos synchronization (CS) has been studied for some time now (Pecora & Carroll 1990), for systems with only a few degrees of freedom as well as for systems described by partial differential equations (Boccaletti et al. 2002). CS in general is said to be present in a pair of coupled dynamical systems when a specific property of each system has the same time evolution for both, even though the evolution itself is chaotic. There have been some studies of CS for systems with an infinite number of degrees of freedom (Kocarev et al. 1997), but CS for Navier-Stokes (NS) turbulence seems not to have been investigated so far. We focus on the synchronization of the small scales of a turbulent flow for which the time history of large scales is prescribed. We present DNS results which show that high-wavenumbers in turbulence are fully slaved to modes with wavenumbers up to a critical fraction of the Kolmogorov dissipation wavenumber. We compare our results with related ideas of ``approximate inertial manifolds.'' The motivation for our work is to study deeply sub-Kolmogorov scales in fully developed turbulence (Schumacher 2007), which we show are recoverable even at very high Reynolds number from simulations that only resolve down to about the Kolmogorov scale. This work is supported by the National Science Foundation's CDI-II program, project CMMI-0941530.
Scaling Navier-Stokes equation in nanotubes
NASA Astrophysics Data System (ADS)
Gǎrǎjeu, Mihail; Gouin, Henri; Saccomandi, Giuseppe
2013-08-01
On one hand, classical Monte Carlo and molecular dynamics simulations have been very useful in the study of liquids in nanotubes, enabling a wide variety of properties to be calculated in intuitive agreement with experiments. On the other hand, recent studies indicate that the theory of continuum breaks down only at the nanometer level; consequently flows through nanotubes still can be investigated with Navier-Stokes equations if we take suitable boundary conditions into account. The aim of this paper is to study the statics and dynamics of liquids in nanotubes by using methods of nonlinear continuum mechanics. We assume that the nanotube is filled with only a liquid phase; by using a second gradient theory the static profile of the liquid density in the tube is analytically obtained and compared with the profile issued from molecular dynamics simulation. Inside the tube there are two domains: a thin layer near the solid wall where the liquid density is non-uniform and a central core where the liquid density is uniform. In the dynamic case a closed form analytic solution seems to be no more possible, but by a scaling argument it is shown that, in the tube, two distinct domains connected at their frontiers still exist. The thin inhomogeneous layer near the solid wall can be interpreted in relation with the Navier length when the liquid slips on the boundary as it is expected by experiments and molecular dynamics calculations.
An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Korte, John J.
1991-01-01
An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required
A non-conforming 3D spherical harmonic transport solver
Van Criekingen, S.
2006-07-01
A new 3D transport solver for the time-independent Boltzmann transport equation has been developed. This solver is based on the second-order even-parity form of the transport equation. The angular discretization is performed through the expansion of the angular neutron flux in spherical harmonics (PN method). The novelty of this solver is the use of non-conforming finite elements for the spatial discretization. Such elements lead to a discontinuous flux approximation. This interface continuity requirement relaxation property is shared with mixed-dual formulations such as the ones based on Raviart-Thomas finite elements. Encouraging numerical results are presented. (authors)
Navier-Stokes computations for circulation control airfoils
NASA Technical Reports Server (NTRS)
Pulliam, Thomas H.; Jespersen, Dennis C.; Barth, Timothy J.
1987-01-01
Navier-Stokes computations of subsonic to transonic flow past airfoils with augmented lift due to rearward jet blowing over a curved trailing edge are presented. The approach uses a spiral grid topology. Solutions are obtained using a Navier-Stokes code which employs an implicit finite difference method, an algebraic turbulence model, and developments which improve stability, convergence, and accuracy. Results are compared against experiments for no jet blowing and moderate jet pressures and demonstrate the capability to compute these complicated flows.
Navier-Stokes computations for circulation controlled airfoils
NASA Technical Reports Server (NTRS)
Pulliam, T. H.; Jesperen, D. C.; Barth, T. J.
1986-01-01
Navier-Stokes computations of subsonic to transonic flow past airfoils with augmented lift due to rearward jet blowing over a curved trailing edge are presented. The approach uses a spiral grid topology. Solutions are obtained using a Navier-Stokes code which employs an implicit finite difference method, an algebraic turbulence model, and developments which improve stability, convergence, and accuracy. Results are compared against experiments for no jet blowing and moderate jet pressures and demonstrate the capability to compute these complicated flows.
On relaxation times in the Navier-Stokes-Voigt model
NASA Astrophysics Data System (ADS)
Layton, William J.; Rebholz, Leo G.
2013-03-01
We study analytically and numerically the relaxation time of flow evolution governed by the Navier-Stokes-Voigt (NSV) model. We first show that for the Taylor-Green vortex decay problem, NSV admits an exact solution which evolves slower than true fluid flow. Secondly, we show numerically for a channel flow test problem using standard discretisation methods that although NSV provides more regular solutions compared to usual Navier-Stokes solutions, NSV approximations take significantly longer to reach the steady state.
Compressible Navier-Stokes Equations with Revised Maxwell's Law
NASA Astrophysics Data System (ADS)
Hu, Yuxi; Racke, Reinhard
2016-05-01
We investigate the compressible Navier-Stokes equations where the constitutive law for the stress tensor given by Maxwell's law is revised to a system of relaxation equations for two parts of the tensor. The global well-posedness is proved as well as the compatibility with the classical compressible Navier-Stokes system in the sense that, for vanishing relaxation parameters, the solutions to the Maxwell system are shown to converge to solutions of the classical system.
Renumbering Methods to Unleash Multi-Threaded Approaches for a General Navier-Stokes Implementation
NASA Astrophysics Data System (ADS)
Vezolle, Pascal; Fournier, Yvan; Tallet, Nicolas; Heymans, Jerrold; D'Amora, Bruce
2010-09-01
Our investigation leverages the general industrial Navier-Stokes open-source Computational Fluid Dynamics (CFD) application, Code_Saturne, developed by Électricité de France (EDF). We deal with how to take advantage of the emerging processor features such as many-cores, Simultaneous Multi-Threading (SMT) and Thread Level Speculation (TLS), through a mixed MPI/multithreads approach. We focus here on the per-node performance improvements and present the constraints for a multithreads implementation to solve the general 3D Navier-Stokes equations using a finite volume discretization into polyhedral cells. We describe a simple and efficient mesh numbering scheme allowing us to introduce OpenMP and Thread Level Speculation implementations with minimal impact to overall code structure.
Three-Dimensional Navier-Stokes Calculations Using the Modified Space-Time CESE Method
NASA Technical Reports Server (NTRS)
Chang, Chau-lyan
2007-01-01
The space-time conservation element solution element (CESE) method is modified to address the robustness issues of high-aspect-ratio, viscous, near-wall meshes. In this new approach, the dependent variable gradients are evaluated using element edges and the corresponding neighboring solution elements while keeping the original flux integration procedure intact. As such, the excellent flux conservation property is retained and the new edge-based gradients evaluation significantly improves the robustness for high-aspect ratio meshes frequently encountered in three-dimensional, Navier-Stokes calculations. The order of accuracy of the proposed method is demonstrated for oblique acoustic wave propagation, shock-wave interaction, and hypersonic flows over a blunt body. The confirmed second-order convergence along with the enhanced robustness in handling hypersonic blunt body flow calculations makes the proposed approach a very competitive CFD framework for 3D Navier-Stokes simulations.
A concurrent hybrid Navier-Stokes/Euler approach to fluid dynamic computations
NASA Technical Reports Server (NTRS)
Tavella, Domingo A.; Djomehri, M. J.; Kislitzin, Katherine T.; Blake, Matthew W.; Erickson, Larry L.
1993-01-01
We present a methodology for the numerical simulation of flow fields by the simultaneous application of two distinct approaches to computational aerodynamics. We compute the three dimensional flow field of a missile at moderate angle of attack by dividing the flow field into two regions: a region near the surface where we use a structured grid and a Navier Stokes solver, and a region farther away from the surface where we utilize an unstructured grid and an Euler solver. The two solvers execute as independent UNIX processes either on the same machine or on two machines. The solvers communicate data across their common interfaces within the same machine or over the network. The computations indicate that extensively separated flow fields can be computed without significant distortion by combining viscous and inviscid solvers.
NASA Astrophysics Data System (ADS)
Tavelli, Maurizio; Dumbser, Michael
2016-08-01
unstructured meshes allows to discretize even complex physical domains with very coarse grids in both, space and time. The proposed method is verified for approximation polynomials of degree up to four in space and time by solving a series of typical 3D test problems and by comparing the obtained numerical results with available exact analytical solutions, or with other numerical or experimental reference data. To the knowledge of the authors, this is the first time that a space-time discontinuous Galerkin finite element method is presented for the three-dimensional incompressible Navier-Stokes equations on staggered unstructured tetrahedral grids.
NASA Astrophysics Data System (ADS)
Nan, Zhijie; Zheng, Xiaoxin
2016-09-01
We study Cauchy problem of the 3D Navier-Stokes equations with hyper-dissipation. By using the Fourier localization technique, we prove that the system has a unique global solution for large initial data in a critical Fourier-Herz space. More importantly, the energy of this solution is infinite.
CFL3D User's Manual (Version 5.0)
NASA Technical Reports Server (NTRS)
Krist, Sherrie L.; Biedron, Robert T.; Rumsey, Christopher L.
1998-01-01
This document is the User's Manual for the CFL3D computer code, a thin-layer Reynolds-averaged Navier-Stokes flow solver for structured multiple-zone grids. Descriptions of the code's input parameters, non-dimensionalizations, file formats, boundary conditions, and equations are included. Sample 2-D and 3-D test cases are also described, and many helpful hints for using the code are provided.
Wang, Z J
2012-12-06
The overriding objective for this project is to develop an efficient and accurate method for capturing strong discontinuities and fine smooth flow structures of disparate length scales with unstructured grids, and demonstrate its potentials for problems relevant to DOE. More specifically, we plan to achieve the following objectives: 1. Extend the SV method to three dimensions, and develop a fourth-order accurate SV scheme for tetrahedral grids. Optimize the SV partition by minimizing a form of the Lebesgue constant. Verify the order of accuracy using the scalar conservation laws with an analytical solution; 2. Extend the SV method to Navier-Stokes equations for the simulation of viscous flow problems. Two promising approaches to compute the viscous fluxes will be tested and analyzed; 3. Parallelize the 3D viscous SV flow solver using domain decomposition and message passing. Optimize the cache performance of the flow solver by designing data structures minimizing data access times; 4. Demonstrate the SV method with a wide range of flow problems including both discontinuities and complex smooth structures. The objectives remain the same as those outlines in the original proposal. We anticipate no technical obstacles in meeting these objectives.
The novel high-performance 3-D MT inverse solver
NASA Astrophysics Data System (ADS)
Kruglyakov, Mikhail; Geraskin, Alexey; Kuvshinov, Alexey
2016-04-01
We present novel, robust, scalable, and fast 3-D magnetotelluric (MT) inverse solver. The solver is written in multi-language paradigm to make it as efficient, readable and maintainable as possible. Separation of concerns and single responsibility concepts go through implementation of the solver. As a forward modelling engine a modern scalable solver extrEMe, based on contracting integral equation approach, is used. Iterative gradient-type (quasi-Newton) optimization scheme is invoked to search for (regularized) inverse problem solution, and adjoint source approach is used to calculate efficiently the gradient of the misfit. The inverse solver is able to deal with highly detailed and contrasting models, allows for working (separately or jointly) with any type of MT responses, and supports massive parallelization. Moreover, different parallelization strategies implemented in the code allow optimal usage of available computational resources for a given problem statement. To parameterize an inverse domain the so-called mask parameterization is implemented, which means that one can merge any subset of forward modelling cells in order to account for (usually) irregular distribution of observation sites. We report results of 3-D numerical experiments aimed at analysing the robustness, performance and scalability of the code. In particular, our computational experiments carried out at different platforms ranging from modern laptops to HPC Piz Daint (6th supercomputer in the world) demonstrate practically linear scalability of the code up to thousands of nodes.
A locally stabilized immersed boundary method for the compressible Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Brehm, C.; Hader, C.; Fasel, H. F.
2015-08-01
A higher-order immersed boundary method for solving the compressible Navier-Stokes equations is presented. The distinguishing feature of this new immersed boundary method is that the coefficients of the irregular finite-difference stencils in the vicinity of the immersed boundary are optimized to obtain improved numerical stability. This basic idea was introduced in a previous publication by the authors for the advection step in the projection method used to solve the incompressible Navier-Stokes equations. This paper extends the original approach to the compressible Navier-Stokes equations considering flux vector splitting schemes and viscous wall boundary conditions at the immersed geometry. In addition to the stencil optimization procedure for the convective terms, this paper discusses other key aspects of the method, such as imposing flux boundary conditions at the immersed boundary and the discretization of the viscous flux in the vicinity of the boundary. Extensive linear stability investigations of the immersed scheme confirm that a linearly stable method is obtained. The method of manufactured solutions is used to validate the expected higher-order accuracy and to study the error convergence properties of this new method. Steady and unsteady, 2D and 3D canonical test cases are used for validation of the immersed boundary approach. Finally, the method is employed to simulate the laminar to turbulent transition process of a hypersonic Mach 6 boundary layer flow over a porous wall and subsonic boundary layer flow over a three-dimensional spherical roughness element.
Fischer, P.F.
1996-12-31
Efficient solution of the Navier-Stokes equations in complex domains is dependent upon the availability of fast solvers for sparse linear systems. For unsteady incompressible flows, the pressure operator is the leading contributor to stiffness, as the characteristic propagation speed is infinite. In the context of operator splitting formulations, it is the pressure solve which is the most computationally challenging, despite its elliptic origins. We seek to improve existing spectral element iterative methods for the pressure solve in order to overcome the slow convergence frequently observed in the presence of highly refined grids or high-aspect ratio elements.
Explicit and implicit solution of the Navier-Stokes equations on a massively parallel computer
NASA Technical Reports Server (NTRS)
Levit, Creon; Jespersen, Dennis
1988-01-01
The design, implementation, and performance of a two-dimensional time-accurate Navier-Stokes solver for the CM2 supercomputer are described. The program uses a single processor for each grid point. Two different time-stepping methods have so far been implemented: an explicit third-order Runge-Kutta method and an implicit approximation-factorization method. The CM2 results are checked against those of a mature well-vectorized Cray 2 program, both for correctness and performance. The code is found to be correct, and the performance in some cases is up to several times that of the Cray 2.
Hypersonic nonequilibrium Navier-Stokes solutions over an ablating graphite nosetip
Chen, Y.K.; Henline, W.D.
1994-09-01
The general boundary conditions, including mass and energy balances, of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver GASP (General Aerodynamics Simulation Program). A test case with a proposed hypersonic test-vehicle configuration and associated freestream conditions was developed. The solutions of the GASP code with various surface boundary conditions were obtained and compared with those of the ASCC (ABRES Shape Change) code, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined. 22 refs.
Hypersonic nonequilibrium Navier-Stokes solutions over an ablating graphite nosetip
NASA Astrophysics Data System (ADS)
Chen, Y.-K.; Henline, W. D.
1994-09-01
The general boundary conditions, including mass and energy balances, of chemically equilibrated or nonequilibrated gas adjacent to ablating surfaces have been derived. A computer procedure based on these conditions was developed and interfaced with the Navier-Stokes solver GASP (General Aerodynamics Simulation Program). A test case with a proposed hypersonic test-vehicle configuration and associated freestream conditions was developed. The solutions of the GASP code with various surface boundary conditions were obtained and compared with those of the ASCC (ABRES Shape Change) code, and the effect of nonequilibrium gas as well as surface chemistry on surface heating and ablation rate were examined.
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Turkel, Eli
2006-01-01
We apply an unsteady Reynolds-averaged Navier-Stokes (URANS) solver for the simulation of a synthetic jet created by a single diaphragm piezoelectric actuator in quiescent air. This configuration was designated as Case 1 for the CFDVAL2004 workshop held at Williamsburg, Virginia, in March 2004. Time-averaged and instantaneous data for this case were obtained at NASA Langley Research Center, using multiple measurement techniques. Computational results for this case using one-equation Spalart-Allmaras and two-equation Menter's turbulence models are presented along with the experimental data. The effect of grid refinement, preconditioning and time-step variation are also examined in this paper.
Aerodynamic Design Optimization on Unstructured Meshes Using the Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Anderson, W. Kyle
1998-01-01
A discrete adjoint method is developed and demonstrated for aerodynamic design optimization on unstructured grids. The governing equations are the three-dimensional Reynolds-averaged Navier-Stokes equations coupled with a one-equation turbulence model. A discussion of the numerical implementation of the flow and adjoint equations is presented. Both compressible and incompressible solvers are differentiated and the accuracy of the sensitivity derivatives is verified by comparing with gradients obtained using finite differences. Several simplifying approximations to the complete linearization of the residual are also presented, and the resulting accuracy of the derivatives is examined. Demonstration optimizations for both compressible and incompressible flows are given.
Navier-Stokes computations past a prolate spheroid at incidence. II - High incidence case
NASA Astrophysics Data System (ADS)
Piquet, J.; Queutey, P.
1993-01-01
A fully elliptic numerical method for the solution of Reynolds-averaged Navier-Stokes equations (RANSEs) is developed and applied to the lifting flow past a prolate spheroid. The method uses a system of numerically generated curvilinear coordinates; it retains the Cartesian velocity components, a nonstaggered grid, and a segregated approach in which a pressure solver couples velocity and pressure fields. It is shown that the RANSEs could describe a strongly concentrated vortical flow past a 6:1 prolate spheroid at high incidence (30 deg) and a Reynolds number of 7.2 x 10 exp 6.
Numerical solutions of the complete Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Hassan, H. A.
1986-01-01
Using ideas from the kinetic theory, the Navier-Stokes equations are modified in such a way that they can be cast as a set of first order hyperbolic equations. This is achieved by incorporating time dependent terms into the definition of the stress tensor and the heat flux vectors. The boundary conditions are then determined from the theory of characteristics. Because the resulting equations reduce to the traditional Navier-Stokes equations when the steady state is reached, the present approach provides a straightforward scheme for the determination of inflow and outflow boundary conditions. The method is validated by comparing its predictions with known exact solutions of the steady Navier-Stokes equations.
On multigrid methods for the Navier-Stokes Computer
NASA Technical Reports Server (NTRS)
Nosenchuck, D. M.; Krist, S. E.; Zang, T. A.
1988-01-01
The overall architecture of the multipurpose parallel-processing Navier-Stokes Computer (NSC) being developed by Princeton and NASA Langley (Nosenchuck et al., 1986) is described and illustrated with extensive diagrams, and the NSC implementation of an elementary multigrid algorithm for simulating isotropic turbulence (based on solution of the incompressible time-dependent Navier-Stokes equations with constant viscosity) is characterized in detail. The present NSC design concept calls for 64 nodes, each with the performance of a class VI supercomputer, linked together by a fiber-optic hypercube network and joined to a front-end computer by a global bus. In this configuration, the NSC would have a storage capacity of over 32 Gword and a peak speed of over 40 Gflops. The multigrid Navier-Stokes code discussed would give sustained operation rates of about 25 Gflops.
What do the Navier-Stokes equations mean?
NASA Astrophysics Data System (ADS)
Schneiderbauer, Simon; Krieger, Michael
2014-01-01
The Navier-Stokes equations are nonlinear partial differential equations describing the motion of fluids. Due to their complicated mathematical form they are not part of secondary school education. A detailed discussion of fundamental physics—the conservation of mass and Newton’s second law—may, however, increase the understanding of the behaviour of fluids. Based on these principles the Navier-Stokes equations can be derived. This article attempts to make these equations available to a wider readership, especially teachers and undergraduate students. Therefore, in this article a derivation restricted to simple differential calculus is presented. Finally, we try to give answers to the questions ‘what is a fluid?’ and ‘what do the Navier-Stokes equations mean?’.
Lattice-gas automata for the Navier-Stokes equation
NASA Astrophysics Data System (ADS)
Frisch, U.; Hasslacher, B.; Pomeau, Y.
1986-04-01
It is shown that a class of deterministic lattice gases with discrete Boolean elements simulates the Navier-Stokes equations, and can be used to design simple, massively parallel computing machines. A hexagonal lattice gas (HLG) model consisting of a triangular lattice with hexagonal symmetry is developed, and is shown to lead to the two-dimensional Navier-Stokes equations. The three-dimensional formulation is obtained by a splitting method in which the nonlinear term in the three-dimensional Navier-Stokes equation is recasts as the sum of two terms, each containing spurious elements and each realizable on a different lattice. Freed slip and rigid boundary conditions are easily implemented. It is noted that lattice-gas models must be run at moderate Mach numbers to remain incompressible, and to avoid spurious high-order nonlinear terms. The model gives a concrete hydrodynamical example of how cellular automata can be used to simulate classical nonlinear fields.
A Nonlinear Modal Aeroelastic Solver for FUN3D
NASA Technical Reports Server (NTRS)
Goldman, Benjamin D.; Bartels, Robert E.; Biedron, Robert T.; Scott, Robert C.
2016-01-01
A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.
Exact solutions of the generalized Navier- Stokes equations for benchmarking
NASA Astrophysics Data System (ADS)
Bourchtein, Andrei
2002-08-01
The generalized Navier- Stokes equations for incompressible viscous flows through isotropic granular porous medium are studied. Some analytical classic solutions of the Navier- Stokes equations are generalized to the case of the considered equations. Obtained solutions of generalized equations reduce to classic ones as porosity effect disappears. Average velocity of generalized solutions is calculated and evaluated in two limiting regimes of flow. In the shallow conduit, the generalized flow rate approximates the free (without porous medium) flow rate and in the case of removed boundaries this approaches Darcy's law. The use of the derived exact solutions for benchmarking purposes is described. Copyright
Factorization of the Compressible Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Roberts, Thomas W.
2005-01-01
The Navier-Stokes equations for a Newtonian ideal gas are examined to determine the factorizable form of the equations relevant to the construction of a factorizable relaxation scheme. The principal linearization of the equations is found by examining the relative magnitude of the terms for short-wavelength errors. The principal part of the operator is then found. Comparison of the factors of the Navier-Stokes and Euler equations differ qualitatively because of the coupling of entropy and pressure through thermal diffusion. Special cases of the factorization are considered.
Some recent applications of Navier-Stokes codes to rotorcraft
NASA Technical Reports Server (NTRS)
Mccroskey, W. J.
1992-01-01
Many operational limitations of helicopters and other rotary-wing aircraft are due to nonlinear aerodynamic phenomena incuding unsteady, three-dimensional transonic and separated flow near the surfaces and highly vortical flow in the wakes of rotating blades. Modern computational fluid dynamics (CFD) technology offers new tools to study and simulate these complex flows. However, existing Euler and Navier-Stokes codes have to be modified significantly for rotorcraft applications, and the enormous computational requirements presently limit their use in routine design applications. Nevertheless, the Euler/Navier-Stokes technology is progressing in anticipation of future supercomputers that will enable meaningful calculations to be made for complete rotorcraft configurations.
Algorithm implementation on the Navier-Stokes computer
NASA Technical Reports Server (NTRS)
Krist, Steven E.; Zang, Thomas A.
1987-01-01
The Navier-Stokes Computer is a multi-purpose parallel-processing supercomputer which is currently under development at Princeton University. It consists of multiple local memory parallel processors, called Nodes, which are interconnected in a hypercube network. Details of the procedures involved in implementing an algorithm on the Navier-Stokes computer are presented. The particular finite difference algorithm considered in this analysis was developed for simulation of laminar-turbulent transition in wall bounded shear flows. Projected timing results for implementing this algorithm indicate that operation rates in excess of 42 GFLOPS are feasible on a 128 Node machine.
Pseudo-time algorithms for the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Swanson, R. C.; Turkel, E.
1986-01-01
A pseudo-time method is introduced to integrate the compressible Navier-Stokes equations to a steady state. This method is a generalization of a method used by Crocco and also by Allen and Cheng. We show that for a simple heat equation that this is just a renormalization of the time. For a convection-diffusion equation the renormalization is dependent only on the viscous terms. We implement the method for the Navier-Stokes equations using a Runge-Kutta type algorithm. This permits the time step to be chosen based on the inviscid model only. We also discuss the use of residual smoothing when viscous terms are present.
Recent Enhancements To The FUN3D Flow Solver For Moving-Mesh Applications
NASA Technical Reports Server (NTRS)
Biedron, Robert T,; Thomas, James L.
2009-01-01
An unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids has been extended to handle general mesh movement involving rigid, deforming, and overset meshes. Mesh deformation is achieved through analogy to elastic media by solving the linear elasticity equations. A general method for specifying the motion of moving bodies within the mesh has been implemented that allows for inherited motion through parent-child relationships, enabling simulations involving multiple moving bodies. Several example calculations are shown to illustrate the range of potential applications. For problems in which an isolated body is rotating with a fixed rate, a noninertial reference-frame formulation is available. An example calculation for a tilt-wing rotor is used to demonstrate that the time-dependent moving grid and noninertial formulations produce the same results in the limit of zero time-step size.
NASA Technical Reports Server (NTRS)
Goodrich, John W.
1991-01-01
An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.
NASA Technical Reports Server (NTRS)
Goodrich, John W.
1991-01-01
An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.
Analysis of regularized Navier-Stokes equations, 2
NASA Technical Reports Server (NTRS)
Ou, Yuh-Roung; Sritharan, S. S.
1989-01-01
A practically important regularization of the Navier-Stokes equations was analyzed. As a continuation of the previous work, the structure of the attractors characterizing the solutins was studied. Local as well as global invariant manifolds were found. Regularity properties of these manifolds are analyzed.
Symmetric approximations of the Navier-Stokes equations
Kobel'kov, G M
2002-08-31
A new method for the symmetric approximation of the non-stationary Navier-Stokes equations by a Cauchy-Kovalevskaya-type system is proposed. Properties of the modified problem are studied. In particular, the convergence as {epsilon}{yields}0 of the solutions of the modified problem to the solutions of the original problem on an infinite interval is established.
Symmetric approximations of the Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Kobel'kov, G. M.
2002-08-01
A new method for the symmetric approximation of the non-stationary Navier-Stokes equations by a Cauchy-Kovalevskaya-type system is proposed. Properties of the modified problem are studied. In particular, the convergence as \\varepsilon\\to0 of the solutions of the modified problem to the solutions of the original problem on an infinite interval is established.
From Petrov-Einstein to Navier-Stokes
NASA Astrophysics Data System (ADS)
Lysov, Vyacheslav
The fluid/gravity correspondence relates solutions of the incompressible Navier-Stokes equation to metrics which solve the Einstein equations. We propose propose two possible approaches to establish this correspondence: perturbative expansion for shear modes and large mean curvature expansion for algebraically special metrics. We show by explicit construction that for every solution of the incompressible Navier-Stokes equation in p+1 dimensions, there is an associated "dual" solution of the vacuum Einstein equations in p+2 dimensions. The dual geometry has an intrinsically flat time-like boundary segment whose extrinsic curvature is given by the stress tensor of the Navier-Stokes fluid. We consider a "near-horizon" limit in which hypersurface becomes highly accelerated. The near-horizon expansion in gravity is shown to be mathematically equivalent to the hydrodynamic expansion in fluid dynamics, and the Einstein equation reduces to the incompressible Navier-Stokes equation. It is shown that imposing a Petrov type I condition on the hypersurface geometry reduces the degrees of freedom in the extrinsic curvature to those of a fluid. Moreover, expanding around a limit in which the mean curvature of the embedding diverges, the leading-order Einstein constraint equations on hypersurface are shown to reduce to the non-linear incompressible Navier-Stokes equation for a fluid moving in hypersurface. We extend the fluid/gravity correspondence to include the magnetohydrodynamics/gravity correspondence, which translates solutions of the equations of magnetohydrodynamics (describing charged fluids) into geometries that satisfy the Einstein-Maxwell equations. We present an explicit example of this new correspondence in the context of flat Minkowski space. We show that a perturbative deformation of the Rindler wedge satisfies the Einstein-Maxwell equations provided that the parameters appearing in the expansion, which we interpret as fluid fields, satisfy the
The development of an intelligent interface to a computational fluid dynamics flow-solver code
NASA Technical Reports Server (NTRS)
Williams, Anthony D.
1988-01-01
Researchers at NASA Lewis are currently developing an 'intelligent' interface to aid in the development and use of large, computational fluid dynamics flow-solver codes for studying the internal fluid behavior of aerospace propulsion systems. This paper discusses the requirements, design, and implementation of an intelligent interface to Proteus, a general purpose, 3-D, Navier-Stokes flow solver. The interface is called PROTAIS to denote its introduction of artificial intelligence (AI) concepts to the Proteus code.
NASA Astrophysics Data System (ADS)
Miao, Sha; Hendrickson, Kelli; Liu, Yuming; Subramani, Hariprasad
2015-11-01
This work presents a novel and efficient Cartesian-grid based simulation capability for the study of an incompressible, turbulent gas layer over a liquid flow with disparate Reynolds numbers in two phases. This capability couples a turbulent gas-flow solver and a liquid-layer based on a second-order accurate Boundary Data Immersion Method (BDIM) at the deformable interface. The turbulent gas flow solver solves the incompressible Navier-Stokes equations via direct numerical simulation or through turbulence closure (unsteady Reynolds-Averaged Navier-Stokes Models) for Reynolds numbers O(106). In this application, a laminar liquid layer solution is obtained from depth-integrated Navier-Stokes equations utilizing shallow water wave assumptions. The immersed boundary method (BDIM) enforces the coupling at the deformable interface, the boundary conditions to turbulence closure equations and defines the domain geometry on the Cartesian grid. Validations are made for the turbulent gas channel flow over high-viscosity liquid. This simulation capability can be applied to problems in the oil and industrial sector such as channel and pipe flows with heavy oils as well as wind wave generation in shallow waters. Sponsored by the Chevron Energy Technology Company.
NASA Technical Reports Server (NTRS)
Atkins, Harold
1991-01-01
A multiple block multigrid method for the solution of the three dimensional Euler and Navier-Stokes equations is presented. The basic flow solver is a cell vertex method which employs central difference spatial approximations and Runge-Kutta time stepping. The use of local time stepping, implicit residual smoothing, multigrid techniques and variable coefficient numerical dissipation results in an efficient and robust scheme is discussed. The multiblock strategy places the block loop within the Runge-Kutta Loop such that accuracy and convergence are not affected by block boundaries. This has been verified by comparing the results of one and two block calculations in which the two block grid is generated by splitting the one block grid. Results are presented for both Euler and Navier-Stokes computations of wing/fuselage combinations.
NASA Astrophysics Data System (ADS)
Dutta, Vimala
1993-07-01
An implicit finite volume nodal point scheme has been developed for solving the two-dimensional compressible Navier-Stokes equations. The numerical scheme is evolved by efficiently combining the basic ideas of the implicit finite-difference scheme of Beam and Warming (1978) with those of nodal point schemes due to Hall (1985) and Ni (1982). The 2-D Navier-Stokes solver is implemented for steady, laminar/turbulent flows past airfoils by using C-type grids. Turbulence closure is achieved by employing the algebraic eddy-viscosity model of Baldwin and Lomax (1978). Results are presented for the NACA-0012 and RAE-2822 airfoil sections. Comparison of the aerodynamic coefficients with experimental results for the different test cases presented here establishes the validity and efficiency of the method.
Implementation and Validation of the Chien k-epsilon Turbulence Model in the Wind Navier-Stokes Code
NASA Technical Reports Server (NTRS)
Yoder, Dennis A.; Georgiadis, Nicholas J.
1999-01-01
The two equation k-epsilon turbulence model of Chien has been implemented in the WIND Navier-Stokes flow solver. Details of the numerical solution algorithm, initialization procedure, and stability enhancements are described. Results obtained with this version of the model are compared with those from the Chien k-epsilon model in the NPARC Navier-Stokes code and from the WIND SST model for three validation cases: the incompressible flow over a smooth flat plate, the incompressible flow over a backward facing step, and the shock-induced flow separation inside a transonic diffuser. The k-epsilon model results indicate that the WIND model functions very similarly to that in NPARC, though the WIND code appears to he slightly more accurate in the treatment of the near-wall region. Comparisons of the k-epsilon model results with those from the SST model were less definitive, as each model exhibited strengths and weaknesses for each particular case.
NASA Technical Reports Server (NTRS)
Ghosh, Amrit Raj
1996-01-01
The viscous, Navier-Stokes solver for turbomachinery applications, MSUTC has been modified to include the rotating frame formulation. The three-dimensional thin-layer Navier-Stokes equations have been cast in a rotating Cartesian frame enabling the freezing of grid motion. This also allows the flow-field associated with an isolated rotor to be viewed as a steady-state problem. Consequently, local time stepping can be used to accelerate convergence. The formulation is validated by running NASA's Rotor 67 as the test case. results are compared between the rotating frame code and the absolute frame code. The use of the rotating frame approach greatly enhances the performance of the code with respect to savings in computing time, without degradation of the solution.
Prediction of Unsteady Transitional Layers in Turbomachinery Using Navier Stokes Equations
NASA Technical Reports Server (NTRS)
Lakshminarayana, B.; Chernobrovkin, A.; Kang, D. J.
1998-01-01
The objective of the research reported in this presentation is to develop computational techniques for the prediction of unsteady transitional flows associated with the rotor stator interaction in turbomachinery. Three low-Reynolds number turbulence models are incorporated in two unsteady Navier-Stokes codes (one is pressure based and the other is time marching with Runge-Kutta time stepping) and evaluated for accuracy in predicting the onset and the end of unsteady transitional patches due to wake passing. The best model is then used for modification and improvement for the leading edge effect. An existing steady Navier-Stokes code was modified to include pseudo-time stepping, which provided acceleration from 5 to 25 times that of the original code. A systematic validation procedure was implemented to assess the effects of the grid, artificial dissipation, physical, and the pseudo-time step for an accurate prediction of transitional flows resulting from the rotor-stator interaction. The ability of the Navier-Stokes code to predict the unsteady transitional flow on a turbomachinery blade is demonstrated. The unsteady pressure and velocity fields are in good agreement with the experimental data and the prediction from the Euler/boundary layer approach. The numerical solver was able to capture all zones (wake induced transitional strip, wake induced turbulent strip, calmed region, etc.) associated with wake induced transition in a compressor cascade. Another significant step is the assessment of k-epsilon turbulence models, including the leading edge modifications. Best results were obtained from the FLB model. The LB model predicted earlier inception of the transition and shorter transition length. Modification of the k-epsilon model was found to be essential for an accurate prediction of the unsteady transitional flow in a compressor cascade. The CH model failed to predict the unsteady transitional flow. Predicted boundary layer was turbulent from the leading edge
Fast non-symmetric iterations and efficient preconditioning for Navier-Stokes equations
Silvester, D.; Elman, H.
1994-12-31
Discretisation of the steady-state Navier-Stokes equations: (u.grad)u-{nu}{del}{sup 2}u + grad p = 0; div u = 0 [1]. in some flow domain {Omega} {contained_in} IR{sup d}, (d = 2 or 3), gives a system of non-linear algebraic equations for discretised variables u (the velocity), and p (the pressure). The authors assume that appropriate boundary conditions are imposed. The non-linear equation system can be linearised using a fixed-point (Picard) iteration to give a matrix system which must be solved at every iteration. Part of this matrix is block diagonal, and consists of d convection-diffusion operators, one for each component of velocity. Two difficulties arise when solving this matrix equation. Firstly, the block diagonal part is not symmetric, although under certain conditions the symmetric part is positive definite. Secondly, the overall system is indefinite. This makes the design of fast and efficient iterative solvers for discretised Navier-Stokes operators an extremely challenging task.
Tetrahedral Finite-Volume Solutions to the Navier-Stokes Equations on Complex Configurations
NASA Technical Reports Server (NTRS)
Frink, Neal T.; Pirzadeh, Shahyar Z.
1998-01-01
A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the U.S. for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.
Solution of Reynolds-averaged Navier-Stokes equations by discontinuous Galerkin method
NASA Astrophysics Data System (ADS)
Kang, Sungwoo; Yoo, Jung Yul
2007-11-01
Discontinuous Galerkin method is a finite element method that allows discontinuities at inter-element boundaries. The discontinuities in the method are treated by approximate Riemann solvers. One important feature of the method is that it obtains high-order accuracy for unstructured mesh with no difficulty. Due to this feature, it can be useful for various practical applications to turbulence and aeroacoustics, but there are few problems to be solved before the method is applicable to practical flow problems. Due to discontinuous approximations in discontinuous Galerkin method, the treatments of viscous terms are complicated and expensive. Moreover, careful treatments of source terms in turbulence model equations are necessary for Reynolds-averaged Navier-Stokes equations to prevent blow-up of high-order-accurate simulations. In this study, we compare high-order accurate discontinuous Galerkin method with different viscous treatments and stabilization of source terms for compressible Reynold-averaged Navier-Stokes equations. Spalart-Allmaras or k-φ model is used for turbulence model. To compare the implemented formulations, steady turbulent flow over a flat plate and unsteady turbulent flow over cavity are solved.
Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme
NASA Technical Reports Server (NTRS)
Coirier, William J.; Powell, Kenneth G.
1995-01-01
A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.
Tetrahedral finite-volume solutions to the Navier-Stokes equations on complex configurations
NASA Astrophysics Data System (ADS)
Frink, N. T.; Pirzadeh, S. Z.
1999-09-01
A review of the algorithmic features and capabilities of the unstructured-grid flow solver USM3Dns is presented. This code, along with the tetrahedral grid generator, VGRIDns, is being extensively used throughout the USA for solving the Euler and Navier-Stokes equations on complex aerodynamic problems. Spatial discretization is accomplished by a tetrahedral cell-centered finite-volume formulation using Roe's upwind flux difference splitting. The fluxes are limited by either a Superbee or MinMod limiter. Solution reconstruction within the tetrahedral cells is accomplished with a simple, but novel, multidimensional analytical formula. Time is advanced by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the near-wall region of the boundary layer. The issues of accuracy and robustness of USM3Dns Navier-Stokes capabilities are addressed for a flat-plate boundary layer, and a full F-16 aircraft with external stores at transonic speed.
NASA Astrophysics Data System (ADS)
Jee, SolKeun; Moser, Robert D.
2012-08-01
This study provides a simple moving-grid scheme which is based on a modified conservative form of the incompressible Navier-Stokes equations for flow around a moving rigid body. The modified integral form is conservative and seeks the solution of the absolute velocity. This approach is different from previous conservative differential forms [1-3] whose reference frame is not inertial. Keeping the reference frame being inertial results in simpler mathematical derivation to the governing equation which includes one dyadic product of velocity vectors in the convective term, whereas the previous [2,3] needs to obtain the time derivative with respect to non-inertial frames causing an additional dyadic product in the convective term. The scheme is implemented in a second-order accurate Navier-Stokes solver and maintains the order of the accuracy. After this verification, the scheme is validated for a pitching airfoil with very high frequencies. The simulation results match very well with the experimental results [4,5], including vorticity fields and a net thrust force. This airfoil simulation also provides detailed vortical structures near the trailing edge and time-evolving aerodynamic forces that are used to investigate the mechanism of the thrust force generation and the effects of the trailing edge shape. The developed moving-grid scheme demonstrates its validity for a rapid oscillating motion.
Navier-Stokes simulation of the flow around an airfoil in Darrieus motion
Tchon, K.F.; Paraschivoiu, I. . Dept. of Mechanical Engineering)
1994-12-01
In order to study the dynamic stall phenomenon on a Darrieus wind turbine, the incompressible flow field around a moving airfoil is simulated using a noninertial stream function-vorticity formulation of the two-dimensional unsteady navier-Stokes equations. Spatial discretization is achieved by the streamline upwind Petrov-Galerkin finite element method on a hybrid mesh composed of a structured region of quadrilateral elements in the vicinity of solid boundaries, an unstructured region of triangular elements elsewhere, and a layer of infinite elements surrounding the domain and projecting the external boundary to infinity. Temporal discretization is achieved by an implicit second order finite difference scheme. At each time step, a nonlinear algebraic system is solved by a Newton method. To accelerate computations, the generalized minimum residual method with an incomplete triangular factorization preconditioning is used to solve the linearized Newton systems. The solver is applied to simulate the flow around a NACA 0015 airfoil in Darrieus motion and the results are compared to experimental observations. To the authors' knowledge, it is the first time that the simulation of such a motion has been performed using the Navier-Stokes equations.
Numerical algorithms for steady and unsteady incompressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Hafez, Mohammed; Dacles, Jennifer
1989-01-01
The numerical analysis of the incompressible Navier-Stokes equations are becoming important tools in the understanding of some fluid flow problems which are encountered in research as well as in industry. With the advent of the supercomputers, more realistic problems can be studied with a wider choice of numerical algorithms. An alternative formulation is presented for viscous incompressible flows. The incompressible Navier-Stokes equations are cast in a velocity/vorticity formulation. This formulation consists of solving the Poisson equations for the velocity components and the vorticity transport equation. Two numerical algorithms for the steady two-dimensional laminar flows are presented. The first method is based on the actual partial differential equations. This uses a finite-difference approximation of the governing equations on a staggered grid. The second method uses a finite element discretization with the vorticity transport equation approximated using a Galerkin approximation and the Poisson equations are obtained using a least squares method. The equations are solved efficiently using Newton's method and a banded direct matrix solver (LINPACK). The method is extended to steady three-dimensional laminar flows and applied to a cubic driven cavity using finite difference schemes and a staggered grid arrangement on a Cartesian mesh. The equations are solved iteratively using a plane zebra relaxation scheme. Currently, a two-dimensional, unsteady algorithm is being developed using a generalized coordinate system. The equations are discretized using a finite-volume approach. This work will then be extended to three-dimensional flows.
Navier-Stokes turbine heat transfer predictions using two-equation turbulence
NASA Technical Reports Server (NTRS)
Ameri, Ali A.; Arnone, Andrea
1992-01-01
Navier-Stokes calculations were carried out in order to predict the heat transfer rates on turbine blades. The calculations were performed using TRAF2D which is a two-dimensional, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using q-omega and k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a non-periodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the Rocketdyne Space Shuttle Main Engine (SSME) high pressure fuel turbine. The comparison with the experimental results, although generally favorable, serves to highlight the weaknesses of the turbulence models and the possible areas of improving these models for use in turbomachinery heat transfer calculations.
Cornetti, G.M.
1995-12-31
The 3D Navier-Stokes equations are solved via the Characteristic-Galerkin method extended to free boundary problems. A temporal discretization procedure is proposed for the case where a preferential direction to move mesh point exists, as in thin domains. Using a single layer of finite elements, the numerical results cover the so-called shallow water 2D approximation, showing the same wave propagation speed.
Use of Navier-Stokes analysis in section design
NASA Astrophysics Data System (ADS)
Nguyen, Phuc N.
1990-12-01
The Navier-Stokes analysis method and a design technique based on conformal mapping are combined to develop new 2-D thick sections. The Eppler-Somers design technique allows for fast design of arbitrary section shape. The well-validated David Taylor Navier-Stokes code is used to optimize the thickness of the section. From previous experimental results, the turbulence characteristics in the near-wake region correlate with the pressure spectra on the trailing edge of a 2-D lifting surface. Therefore, the turbulent kinetic energy, and the Reynolds shear stress are used as design parameters to develop new 2-D sections. Minimizing these parameters is assumed to provide desirable boundary layer and wake characteristics. The characteristics of one new section are compared with those of a baseline section to demonstrate the new foil design method.
Partial Averaged Navier-Stokes approach for cavitating flow
NASA Astrophysics Data System (ADS)
Zhang, L.; Zhang, Y. N.
2015-01-01
Partial Averaged Navier Stokes (PANS) is a numerical approach developed for studying practical engineering problems (e.g. cavitating flow inside hydroturbines) with a resonance cost and accuracy. One of the advantages of PANS is that it is suitable for any filter width, leading a bridging method from traditional Reynolds Averaged Navier-Stokes (RANS) to direct numerical simulations by choosing appropriate parameters. Comparing with RANS, the PANS model will inherit many physical nature from parent RANS but further resolve more scales of motion in great details, leading to PANS superior to RANS. As an important step for PANS approach, one need to identify appropriate physical filter-width control parameters e.g. ratios of unresolved-to-total kinetic energy and dissipation. In present paper, recent studies of cavitating flow based on PANS approach are introduced with a focus on the influences of filter-width control parameters on the simulation results.
Optimal control of thermally coupled Navier Stokes equations
NASA Technical Reports Server (NTRS)
Ito, Kazufumi; Scroggs, Jeffrey S.; Tran, Hien T.
1994-01-01
The optimal boundary temperature control of the stationary thermally coupled incompressible Navier-Stokes equation is considered. Well-posedness and existence of the optimal control and a necessary optimality condition are obtained. Optimization algorithms based on the augmented Lagrangian method with second order update are discussed. A test example motivated by control of transport process in the high pressure vapor transport (HVPT) reactor is presented to demonstrate the applicability of our theoretical results and proposed algorithm.
Compressible Navier Stokes Model with Inflow-Outflow Boundary Conditions
NASA Astrophysics Data System (ADS)
Novo, Sébastien
2005-11-01
In the paper [7], author gives a definition of weak solution to the nonsteady Navier Stokes system of equations which describes compressible and isentropic flows in some bounded region Ω with influx of fluid through a part of the boundary ∂Ω. Here, we present a way for proving existence of such solutions in the same situation as in [7] under the sole hypothesis γ > 3/2 for the adiabatic constant.
Finite element methods and Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Cuvelier, C.; Segal, A.; van Steenhoven, A. A.
This book is devoted to two and three-dimensional FEM analysis of the Navier-Stokes (NS) equations describing one flow of a viscous incompressible fluid. Three different approaches to the NS equations are described: a direct method, a penalty method, and a method that constructs discrete solenoidal vector fields. Subjects of current research which are important from the industrial/technological viewpoint are considered, including capillary-free boundaries, nonisothermal flows, turbulence, and non-Newtonian fluids.
A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations
NASA Technical Reports Server (NTRS)
Guruswamy, Guru P.; Lee, Henry
2014-01-01
The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no
NASA Technical Reports Server (NTRS)
Coirier, William J.; Powell, Kenneth G.
1994-01-01
A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: a gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.
NASA Technical Reports Server (NTRS)
Coirier, William J.; Powell, Kenneth G.
1995-01-01
A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.
Richter, Christiane; Kotz, Frederik; Giselbrecht, Stefan; Helmer, Dorothea; Rapp, Bastian E
2016-06-01
The fluid mechanics of microfluidics is distinctively simpler than the fluid mechanics of macroscopic systems. In macroscopic systems effects such as non-laminar flow, convection, gravity etc. need to be accounted for all of which can usually be neglected in microfluidic systems. Still, there exists only a very limited selection of channel cross-sections for which the Navier-Stokes equation for pressure-driven Poiseuille flow can be solved analytically. From these equations, velocity profiles as well as flow rates can be calculated. However, whenever a cross-section is not highly symmetric (rectangular, elliptical or circular) the Navier-Stokes equation can usually not be solved analytically. In all of these cases, numerical methods are required. However, in many instances it is not necessary to turn to complex numerical solver packages for deriving, e.g., the velocity profile of a more complex microfluidic channel cross-section. In this paper, a simple spreadsheet analysis tool (here: Microsoft Excel) will be used to implement a simple numerical scheme which allows solving the Navier-Stokes equation for arbitrary channel cross-sections. PMID:27233665
Attractors of three-dimensional fast-rotating Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Trahe, Markus
The three-dimensional (3-D) rotating Navier-Stokes equations describe the dynamics of rotating, incompressible, viscous fluids. In this work, they are considered with smooth, time-independent forces and the original statements implied by the classical "Taylor-Proudman Theorem" of geophysics are rigorously proved. It is shown that fully developed turbulence of 3-D fast-rotating fluids is essentially characterized by turbulence of two-dimensional (2-D) fluids in terms of numbers of degrees of freedom. In this context, the 3-D nonlinear "resonant limit equations", which arise in a non-linear averaging process as the rotation frequency O → infinity, are studied and optimal (2-D-type) upper bounds for fractal box and Hausdorff dimensions of the global attractor as well as upper bounds for box dimensions of exponential attractors are determined. Then, the convergence of exponential attractors for the full 3-D rotating Navier-Stokes equations to exponential attractors for the resonant limit equations as O → infinity in the sense of full Hausdorff-metric distances is established. This provides upper and lower semi-continuity of exponential attractors with respect to the rotation frequency and implies that the number of degrees of freedom (attractor dimension) of 3-D fast-rotating fluids is close to that of 2-D fluids. Finally, the algebraic-geometric structure of the Poincare curves, which control the resonances and small divisor estimates for partial differential equations, is further investigated; the 3-D nonlinear limit resonant operators are characterized by three-wave interactions governed by these curves. A new canonical transformation between those curves is constructed; with far-reaching consequences on the density of the latter.
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Singer, Bart A.
2003-01-01
We evaluate the applicability of a production computational fluid dynamics code for conducting detached eddy simulation for unsteady flows. A second-order accurate Navier-Stokes code developed at NASA Langley Research Center, known as TLNS3D, is used for these simulations. We focus our attention on high Reynolds number flow (Re = 5 x 10(sup 4) - 1.4 x 10(sup 5)) past a circular cylinder to simulate flows with large-scale separations. We consider two types of flow situations: one in which the flow at the separation point is laminar, and the other in which the flow is already turbulent when it detaches from the surface of the cylinder. Solutions are presented for two- and three-dimensional calculations using both the unsteady Reynolds-averaged Navier-Stokes paradigm and the detached eddy simulation treatment. All calculations use the standard Spalart-Allmaras turbulence model as the base model.
Shared Memory Parallelism for 3D Cartesian Discrete Ordinates Solver
NASA Astrophysics Data System (ADS)
Moustafa, Salli; Dutka-Malen, Ivan; Plagne, Laurent; Ponçot, Angélique; Ramet, Pierre
2014-06-01
This paper describes the design and the performance of DOMINO, a 3D Cartesian SN solver that implements two nested levels of parallelism (multicore+SIMD) on shared memory computation nodes. DOMINO is written in C++, a multi-paradigm programming language that enables the use of powerful and generic parallel programming tools such as Intel TBB and Eigen. These two libraries allow us to combine multi-thread parallelism with vector operations in an efficient and yet portable way. As a result, DOMINO can exploit the full power of modern multi-core processors and is able to tackle very large simulations, that usually require large HPC clusters, using a single computing node. For example, DOMINO solves a 3D full core PWR eigenvalue problem involving 26 energy groups, 288 angular directions (S16), 46 × 106 spatial cells and 1 × 1012 DoFs within 11 hours on a single 32-core SMP node. This represents a sustained performance of 235 GFlops and 40:74% of the SMP node peak performance for the DOMINO sweep implementation. The very high Flops/Watt ratio of DOMINO makes it a very interesting building block for a future many-nodes nuclear simulation tool.
An Evaluation of Parameters Influencing Jet Mixing Using the WIND Navier-stokes Code
NASA Technical Reports Server (NTRS)
Dembowski, Mary Ann; Georgiadis, Nicholas J.
2002-01-01
The WIND code, a Reynolds-averaged Navier-Stokes solver used for a variety of aerospace flow simulations, was investigated for a Mach 2 nozzle at a series of nozzle stagnation temperatures. Comparisons of WIND calculations are made to experimental measurements of axial velocity, Mach number, and stagnation temperature along the jet centerline. The primary objective was to investigate the capabilities of the two-equation turbulence models available in WIND, version 4.0, for the analysis of heated supersonic nozzle flows. The models examined were the Menter Shear Stress Transport (SST) model and the Chien k-epsilon model, with and without the compressibility correction due to Sarkar. It was observed that all of the turbulence models investigated produced solutions that did not agree well with the experimental measurements. The effects of freestream Mach number and turbulent Prandtl number specifications were also investigated.
NASA Technical Reports Server (NTRS)
Vatsa, Veer N.; Turkel, Eli L.
2004-01-01
We report research experience in applying an Unsteady Reynolds-Averaged Navier-Stokes (URANS) solver for the prediction of time-dependent flows in the presence of an active flow control device. The configuration under consideration is a synthetic jet created by a single diaphragm piezoelectric actuator in quiescent air. Time-averaged and instantaneous data for this case were obtained at Langley Research Center, using multiple measurement techniques. Computational results for this case using one-equation Spalart-Allmaras and two-equation Menter s turbulence models are presented here along with comparisons with the experimental data. The effect of grid refinement, preconditioning and time-step variation are also examined.
An efficient upwind/relaxation algorithm for the Euler and Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Zha, Ge-Cheng; Liu, Dao-Zhi; Ma, Tie-You
1989-05-01
An efficient Euler and full Navier-Stokes solver based on a flux splitting scheme is presented. The original Van Leer flux vector splitting form is extended to arbitrary body-fitted co-ordinates in the physical domain so that it can be used with a finite volume scheme. The block matrix is inverted by Gauss-Seidel iteration. It is verified that the often used reflection boundary condition will produce incorrect flux crossing the wall and cause too large numerical dissipation if flux vector splitting is used. To remove such errors, an appropriate tretment of wall boundary conditions is suggested. Inviscid and viscous steady transonic internal flows are analyzed, including the case of shock-induced boundary layer separation.
Analysis of a High-Lift Multi-Element Airfoil using a Navier-Stokes Code
NASA Technical Reports Server (NTRS)
Whitlock, Mark E.
1995-01-01
A thin-layer Navier-Stokes code, CFL3D, was utilized to compute the flow over a high-lift multi-element airfoil. This study was conducted to improve the prediction of high-lift flowfields using various turbulence models and improved glidding techniques. An overset Chimera grid system is used to model the three element airfoil geometry. The effects of wind tunnel wall modeling, changes to the grid density and distribution, and embedded grids are discussed. Computed pressure and lift coefficients using Spalart-Allmaras, Baldwin-Barth, and Menter's kappa-omega - Shear Stress Transport (SST) turbulence models are compared with experimental data. The ability of CFL3D to predict the effects on lift coefficient due to changes in Reynolds number changes is also discussed.
NASA Astrophysics Data System (ADS)
Li, Zhilin; Xiao, Li; Cai, Qin; Zhao, Hongkai; Luo, Ray
2015-09-01
In this paper, a new Navier-Stokes solver based on a finite difference approximation is proposed to solve incompressible flows on irregular domains with open, traction, and free boundary conditions, which can be applied to simulations of fluid structure interaction, implicit solvent model for biomolecular applications and other free boundary or interface problems. For some problems of this type, the projection method and the augmented immersed interface method (IIM) do not work well or does not work at all. The proposed new Navier-Stokes solver is based on the local pressure boundary method, and a semi-implicit augmented IIM. A fast Poisson solver can be used in our algorithm which gives us the potential for developing fast overall solvers in the future. The time discretization is based on a second order multi-step method. Numerical tests with exact solutions are presented to validate the accuracy of the method. Application to fluid structure interaction between an incompressible fluid and a compressible gas bubble is also presented.
Flux Based Surface Boundary Conditions for Navier-Stokes Simulations
NASA Astrophysics Data System (ADS)
Fertig, M.; Auweter-Kurtz, M.
2005-02-01
During re-entry high thermal combined with mechanical loads arise at the TPS surface of a re-entry vehicle. Due to low gas density, high Knudsen Numbers arise, which indicate rarefaction effects such as thermo-chemical non-equilibrium as well as temperature and velocity slip. With increasing altitude, local Knudsen Numbers predict the failure of continuum equations starting in the bow shock and at the surface. While local failure of the equations in the shock can be neglected for the determination of surface loads, local failure at the surface is not negligible. The validity of continuum models can be extended by emploing surface boundary equations accounting for temperature and velocity slip. A new flux based model has been developed originating on the Boltzmann Equation. Making use of the Enskog Method perturbed partition functions for a multi-component gas are determined from the Boltzmann Equation. By introduction of the moments of Boltzmann's Equation, Maxwell's Transport Equation can be obtained. Particles approaching the surface are distinguished from particles leaving the surface depending on their molecular velocities. Hence, mass, momentum and energy fluxes to the surface can be determined employing the collisional invariants. Reactive as well as scattering models can be easily introduced in order to compute the fluxes from the surface. Finally, flux differences are balanced with the continuum fluxes from the Navier-Stokes equations. Hence, the model is able to predict temperature and velocity slip at the surface of a re-entry vehicle under rarefied conditions. Moreover, it is valid in the continuum regime as well. The boundary equations are solved fully implicit and fully coupled with the non-equilibrium Navier-Stokes Code URANUS. Results are compared to DSMC simulations for the re-entry of the US Space Shuttle orbiter at high altitudes. Key words: Navier-Stokes; re-entry; slip; non-equilibrium.
Towards an ideal preconditioner for linearized Navier-Stokes problems
Murphy, M.F.
1996-12-31
Discretizing certain linearizations of the steady-state Navier-Stokes equations gives rise to nonsymmetric linear systems with indefinite symmetric part. We show that for such systems there exists a block diagonal preconditioner which gives convergence in three GMRES steps, independent of the mesh size and viscosity parameter (Reynolds number). While this {open_quotes}ideal{close_quotes} preconditioner is too expensive to be used in practice, it provides a useful insight into the problem. We then consider various approximations to the ideal preconditioner, and describe the eigenvalues of the preconditioned systems. Finally, we compare these preconditioners numerically, and present our conclusions.
Time-accurate Navier-Stokes calculations with multigrid acceleration
NASA Technical Reports Server (NTRS)
Melson, N. D.; Sanetrik, Mark D.; Atkins, Harold L.
1993-01-01
An efficient method for calculating unsteady flows is presented, with emphasis on a modified version of the thin-layer Navier-Stokes equations. Fourier stability analysis is used to illustrate the effect of treating the source term implicitly instead of explicity, as well as to illustrate other algorithmic choices. A 2D circular cylinder (with a Reynolds number of 1200 and a Mach number of 0.3) is calculated. The present scheme requires only about 10 percent of the computer time required by global minimum time stepping.
Smooth solutions of the Navier-Stokes equations
Pokhozhaev, S I
2014-02-28
We consider smooth solutions of the Cauchy problem for the Navier-Stokes equations on the scale of smooth functions which are periodic with respect to x∈R{sup 3}. We obtain existence theorems for global (with respect to t>0) and local solutions of the Cauchy problem. The statements of these depend on the smoothness and the norm of the initial vector function. Upper bounds for the behaviour of solutions in both classes, which depend on t, are also obtained. Bibliography: 10 titles.
A visual programming environment for the Navier-Stokes computer
NASA Technical Reports Server (NTRS)
Tomboulian, Sherryl; Crockett, Thomas W.; Middleton, David
1988-01-01
The Navier-Stokes computer is a high-performance, reconfigurable, pipelined machine designed to solve large computational fluid dynamics problems. Due to the complexity of the architecture, development of effective, high-level language compilers for the system appears to be a very difficult task. Consequently, a visual programming methodology has been developed which allows users to program the system at an architectural level by constructing diagrams of the pipeline configuration. These schematic program representations can then be checked for validity and automatically translated into machine code. The visual environment is illustrated by using a prototype graphical editor to program an example problem.
Airfoil design method using the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Malone, J. B.; Narramore, J. C.; Sankar, L. N.
1991-01-01
An airfoil design procedure is described that was incorporated into an existing 2-D Navier-Stokes airfoil analysis method. The resulting design method, an iterative procedure based on a residual-correction algorithm, permits the automated design of airfoil sections with prescribed surface pressure distributions. The inverse design method and the technique used to specify target pressure distributions are described. It presents several example problems to demonstrate application of the design procedure. It shows that this inverse design method develops useful airfoil configurations with a reasonable expenditure of computer resources.
Computation of 2D Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Chakrabartty, Sunil Kumar
Two schemes for computing two-dimensional Navier-Stokes equations are described and applied to laminar flow over a flat plate and viscous flow over a NACA0012 airfoil. The variation of local skin-friction coefficient with local Reynolds number is compared with the Blasius solution and that of Swanson and Turkel (1985). The effect of free-stream Mach number on the temperature profile is shown, and a comparison is made of velocity profile at M(infinity) = 0.50 and Re = 500, with no artificial viscosity used for stability. Pressure distributions, local skin friction distributions, and velocity profiles on the airfoil and wake are presented.
Turbulent solution of the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Deissler, R. G.
1980-01-01
The unaveraged Navier-Stokes equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three dimensional cosine velocity fluctuations and periodic boundary conditions are used. No mean gradients are present. The three components of the mean square velocity fluctuations are equal for the initial conditions chosen. The resulting solution shows characteristics of turbulence, such as the nonlinear excitation of small scale fluctuations. For the higher Reynolds numbers the initially nonrandom flow develops into an apparently random turbulence.
Navier-Stokes Analysis of a High Wing Transport High-Lift Configuration with Externally Blown Flaps
NASA Technical Reports Server (NTRS)
Slotnick, Jeffrey P.; An, Michael Y.; Mysko, Stephen J.; Yeh, David T.; Rogers, Stuart E.; Roth, Karlin; Baker, M.David; Nash, S.
2000-01-01
Insights and lessons learned from the aerodynamic analysis of the High Wing Transport (HWT) high-lift configuration are presented. Three-dimensional Navier-Stokes CFD simulations using the OVERFLOW flow solver are compared with high Reynolds test data obtained in the NASA Ames 12 Foot Pressure Wind Tunnel (PWT) facility. Computational analysis of the baseline HWT high-lift configuration with and without Externally Blown Flap (EBF) jet effects is highlighted. Several additional aerodynamic investigations, such as nacelle strake effectiveness and wake vortex studies, are presented. Technical capabilities and shortcomings of the computational method are discussed and summarized.
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Ghaffari, Farhad
2011-01-01
Numerical predictions of the longitudinal aerodynamic characteristics for the Ares I class of vehicles, along with the associated error estimate derived from an iterative convergence grid refinement, are presented. Computational results are based on the unstructured grid, Reynolds-averaged Navier-Stokes flow solver USM3D, with an assumption that the flow is fully turbulent over the entire vehicle. This effort was designed to complement the prior computational activities conducted over the past five years in support of the Ares I Project with the emphasis on the vehicle s last design cycle designated as the A106 configuration. Due to a lack of flight data for this particular design s outer mold line, the initial vehicle s aerodynamic predictions and the associated error estimates were first assessed and validated against the available experimental data at representative wind tunnel flow conditions pertinent to the ascent phase of the trajectory without including any propulsion effects. Subsequently, the established procedures were then applied to obtain the longitudinal aerodynamic predictions at the selected flight flow conditions. Sample computed results and the correlations with the experimental measurements are presented. In addition, the present analysis includes the relevant data to highlight the balance between the prediction accuracy against the grid size and, thus, the corresponding computer resource requirements for the computations at both wind tunnel and flight flow conditions. NOTE: Some details have been removed from selected plots and figures in compliance with the sensitive but unclassified (SBU) restrictions. However, the content still conveys the merits of the technical approach and the relevant results.
Modeling Vortex Generators in a Navier-Stokes Code
NASA Technical Reports Server (NTRS)
Dudek, Julianne C.
2011-01-01
A source-term model that simulates the effects of vortex generators was implemented into the Wind-US Navier-Stokes code. The source term added to the Navier-Stokes equations simulates the lift force that would result from a vane-type vortex generator in the flowfield. The implementation is user-friendly, requiring the user to specify only three quantities for each desired vortex generator: the range of grid points over which the force is to be applied and the planform area and angle of incidence of the physical vane. The model behavior was evaluated for subsonic flow in a rectangular duct with a single vane vortex generator, subsonic flow in an S-duct with 22 corotating vortex generators, and supersonic flow in a rectangular duct with a counter-rotating vortex-generator pair. The model was also used to successfully simulate microramps in supersonic flow by treating each microramp as a pair of vanes with opposite angles of incidence. The validation results indicate that the source-term vortex-generator model provides a useful tool for screening vortex-generator configurations and gives comparable results to solutions computed using gridded vanes.
Reliability enhancement of Navier-Stokes codes through convergence acceleration
NASA Technical Reports Server (NTRS)
Merkle, Charles L.; Dulikravich, George S.
1995-01-01
Methods for enhancing the reliability of Navier-Stokes computer codes through improving convergence characteristics are presented. The improving of these characteristics decreases the likelihood of code unreliability and user interventions in a design environment. The problem referred to as a 'stiffness' in the governing equations for propulsion-related flowfields is investigated, particularly in regard to common sources of equation stiffness that lead to convergence degradation of CFD algorithms. Von Neumann stability theory is employed as a tool to study the convergence difficulties involved. Based on the stability results, improved algorithms are devised to ensure efficient convergence in different situations. A number of test cases are considered to confirm a correlation between stability theory and numerical convergence. The examples of turbulent and reacting flow are presented, and a generalized form of the preconditioning matrix is derived to handle these problems, i.e., the problems involving additional differential equations for describing the transport of turbulent kinetic energy, dissipation rate and chemical species. Algorithms for unsteady computations are considered. The extension of the preconditioning techniques and algorithms derived for Navier-Stokes computations to three-dimensional flow problems is discussed. New methods to accelerate the convergence of iterative schemes for the numerical integration of systems of partial differential equtions are developed, with a special emphasis on the acceleration of convergence on highly clustered grids.
Navier-Stokes analysis of muzzle-blast-type waves
NASA Astrophysics Data System (ADS)
Baysal, O.
1986-05-01
A Navier-Stokes solution is presented as a mathematical model to muzzle-blast-type waves. The study has two novel features. First, it is a combined internal/external analysis relating barrel flow parameters to muzzle environment parameters. Second, the dissipative and dispersive effects of viscosity on the propagation phenomenon are captured. The investigation also serves as a numerical analysis of axisymmetric, high-pressure waves in an unsteady, viscous flow. Conservation-form Navier-Stokes equations are integrated by a two-step, explicit finite-difference scheme. The shocks are captured and treated by the inclusion of artificial dissipative terms. Turbulence is accounted for by an algebraic eddy-viscosity model. The internal flow is solved by a predictor-corrector method of characteristics with the shock fitted in; its results compare very well with the experimental data available. The numerical results obtained simulate the muzzle blast waves and show the effects of viscosity. Comparison with the classical spherical blast wave theory shows the deviation in propagation patterns of the axisymmetric and spherical waves.