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.
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.
Reynolds-Averaged Navier-Stokes Simulations of Two Partial-Span Flap Wing Experiments
NASA Technical Reports Server (NTRS)
Takalluk, M. A.; Laflin, Kelly R.
1998-01-01
Structured Reynolds Averaged Navier-Stokes simulations of two partial-span flap wing experiments were performed. The high-lift aerodynamic and aeroacoustic wind-tunnel experiments were conducted at both the NASA Ames 7-by 10-Foot Wind Tunnel and at the NASA Langley Quiet Flow Facility. The purpose of these tests was to accurately document the acoustic and aerodynamic characteristics associated with the principle airframe noise sources, including flap side-edge noise. Specific measurements were taken that can be used to validate analytic and computational models of the noise sources and associated aerodynamic for configurations and conditions approximating flight for transport aircraft. The numerical results are used to both calibrate a widely used CFD code, CFL3D, and to obtain details of flap side-edge flow features not discernible from experimental observations. Both experimental set-ups were numerically modeled by using multiple block structured grids. Various turbulence models, grid block-interface interaction methods and grid topologies were implemented. Numerical results of both simulations are in excellent agreement with experimental measurements and flow visualization observations. The flow field in the flap-edge region was adequately resolved to discern some crucial information about the flow physics and to substantiate the merger of the two vortical structures. As a result of these investigations, airframe noise modelers have proposed various simplified models which use the results obtained from the steady-state computations as input.
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.
Reynolds-Averaged Navier-Stokes Analysis of Zero Efflux Flow Control Over a Hump Model
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.
2006-01-01
The unsteady flow over a hump model with zero efflux oscillatory flow control is modeled computationally using the unsteady Reynolds-averaged Navier-Stokes equations. Three different turbulence models produce similar results, and do a reasonably good job predicting the general character of the unsteady surface pressure coefficients during the forced cycle. However, the turbulent shear stresses are underpredicted in magnitude inside the separation bubble, and the computed results predict too large a (mean) separation bubble compared with experiment. These missed predictions are consistent with earlier steady-state results using no-flow-control and steady suction, from a 2004 CFD validation workshop for synthetic jets.
Reynolds-Averaged Navier-Stokes Analysis of Zero Efflux Flow Control over a Hump Model
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.
2006-01-01
The unsteady flow over a hump model with zero efflux oscillatory flow control is modeled computationally using the unsteady Reynolds-averaged Navier-Stokes equations. Three different turbulence models produce similar results, and do a reasonably good job predicting the general character of the unsteady surface pressure coefficients during the forced cycle. However, the turbulent shear stresses are underpredicted in magnitude inside the separation bubble, and the computed results predict too large a (mean) separation bubble compared with experiment. These missed predictions are consistent with earlier steady-state results using no-flow-control and steady suction, from a 2004 CFD validation workshop for synthetic jets.
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.
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.
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.
Large Eddy And Reynolds-Averaged Navier-Stokes Simulations Of An Underexpanded Sonic Jet
NASA Astrophysics Data System (ADS)
Gorle, Catherine; Iaccarino, Gianluca
2011-05-01
The complex phenomena involved in scramjet propulsion can be investigated using Reynolds-averaged Navier- Stokes (RANS) simulations, but the computations can not fully represent all the physics involved. The mixing of fuel and air inside supersonic combustion chambers is one of the critical processes for successful operation of scramjet engines that requires modeling. The goal of the present work is the identification and characterization of the uncertainties in RANS mixing models, to support the development of an error model. In the first part of this paper simulations of an underexpanded hydrogen jet in quiescent air are described and compared to Schlieren images for verification of the jet injection modeling. In the second part RANS simulations of an underexpanded ethylene jet in a supersonic cross flow of Mach 2 are presented to identify the uncertainties in the mixing model, and quantify the sensitivity of the model outcome to these uncertainties.
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.
Baseline Validation of Unstructured Grid Reynolds-Averaged Navier-Stokes Toward Flow Control
NASA Technical Reports Server (NTRS)
Joslin, Ronald D.; Viken, Sally A.
2001-01-01
The value of the use of the Reynolds-averaged Navier-Stokes methodology for active flow control applications is assessed. An experimental flow control database exists for a NACA0015 airfoil modified at the leading edge to implement a fluidic actuator; hence, this configuration is used. Computational results are documented for the baseline wing configuration (no control) with the experimental results and assumes two-dimensional flow. The baseline wing configuration has discontinuities at the leading edge, trailing edge, and aft of midchord on the upper surface. A limited number of active flow control applications have been tested in the laboratory and in flight. These applications include dynamic stall control using a deformable leading edge, separation control for takeoff and landing flight conditions using piezoelectric devices, pulsed vortex generators, zero-net-mass oscillations, and thrust vectoring with zero-net-mass piezoelectric-driven oscillatory actuation. As yet, there is no definitive comparison with experimental data that indicates current computational capabilities can quantitatively predict the large aerodynamic performance gains achieved with active flow control in the laboratory. However, one study using the Reynolds-averaged Navier-Stokes (RANS) methodology has shown good quantitative agreement with experimental results for an isolated zero-net-mass actuator. In addition, some recent studies have used RANS to demonstrate qualitative performance gains compared with the experimental data for separation control on an airfoil. Those quantitative comparisons for both baseline and flow control cases indicated that computational results were in poor quantitative agreement with the experiments. The current research thrust will investigate the potential use of an unstructured grid RANS approach to predict aerodynamic performance for active flow control applications building on the early studies. First the computational results must quantitatively match
Reynolds-Averaged Navier-Stokes Solutions to Flat Plate Film Cooling Scenarios
NASA Technical Reports Server (NTRS)
Johnson, Perry L.; Shyam, Vikram; Hah, Chunill
2011-01-01
The predictions of several Reynolds-Averaged Navier-Stokes solutions for a baseline film cooling geometry are analyzed and compared with experimental data. The Fluent finite volume code was used to perform the computations with the realizable k-epsilon turbulence model. The film hole was angled at 35 to the crossflow with a Reynolds number of 17,400. Multiple length-to-diameter ratios (1.75 and 3.5) as well as momentum flux ratios (0.125 and 0.5) were simulated with various domains, boundary conditions, and grid refinements. The coolant to mainstream density ratio was maintained at 2.0 for all scenarios. Computational domain and boundary condition variations show the ability to reduce the computational cost as compared to previous studies. A number of grid refinement and coarsening variations are compared for further insights into the reduction of computational cost. Liberal refinement in the near hole region is valuable, especially for higher momentum jets that tend to lift-off and create a recirculating flow. A lack of proper refinement in the near hole region can severely diminish the accuracy of the solution, even in the far region. The effects of momentum ratio and hole length-to-diameter ratio are also discussed.
NASA Astrophysics Data System (ADS)
Ling, J.; Templeton, J.
2015-08-01
Reynolds Averaged Navier Stokes (RANS) models are widely used in industry to predict fluid flows, despite their acknowledged deficiencies. Not only do RANS models often produce inaccurate flow predictions, but there are very limited diagnostics available to assess RANS accuracy for a given flow configuration. If experimental or higher fidelity simulation results are not available for RANS validation, there is no reliable method to evaluate RANS accuracy. This paper explores the potential of utilizing machine learning algorithms to identify regions of high RANS uncertainty. Three different machine learning algorithms were evaluated: support vector machines, Adaboost decision trees, and random forests. The algorithms were trained on a database of canonical flow configurations for which validated direct numerical simulation or large eddy simulation results were available, and were used to classify RANS results on a point-by-point basis as having either high or low uncertainty, based on the breakdown of specific RANS modeling assumptions. Classifiers were developed for three different basic RANS eddy viscosity model assumptions: the isotropy of the eddy viscosity, the linearity of the Boussinesq hypothesis, and the non-negativity of the eddy viscosity. It is shown that these classifiers are able to generalize to flows substantially different from those on which they were trained. Feature selection techniques, model evaluation, and extrapolation detection are discussed in the context of turbulence modeling applications.
Towards Exploratory Aerodynamic Design using the Reynolds-Averaged Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Koo, David Tai Shun
The aerodynamic optimization framework Jetstream is applied to problems involving lift-constrained drag minimization using the Reynolds-averaged Navier-Stokes equations. A parallel Newton-Krylov algorithm is used to solve the governing equations on multiblock structured meshes; gradients are computed using the discrete-adjoint method. Geometry parameterization and mesh movement are integrated using B-spline control volumes. Drag minimization studies from past works are revisited and strategies are devised to improve optimization convergence. These strategies include linear constraints for geometric feasibility, robust flow solver parameters, and meshing with an O-O topology. The single-point and multi-point optimization of the NASA Common Research Model (CRM) wing geometry is presented. A rectangular NACA0012 wing is optimized with planform design variables, enabling significant changes in span, sweep, taper, and airfoil section. To demonstrate Jetstream's flexibility, a wing based on the B737-900 is optimized with nonplanar winglets, split-tip, and wingtip fence configurations. Finally, the box-wing optimization in subsonic flow is revisited.
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.
Ling, Julia; Templeton, Jeremy Alan
2015-08-04
Reynolds Averaged Navier Stokes (RANS) models are widely used in industry to predict fluid flows, despite their acknowledged deficiencies. Not only do RANS models often produce inaccurate flow predictions, but there are very limited diagnostics available to assess RANS accuracy for a given flow configuration. If experimental or higher fidelity simulation results are not available for RANS validation, there is no reliable method to evaluate RANS accuracy. This paper explores the potential of utilizing machine learning algorithms to identify regions of high RANS uncertainty. Three different machine learning algorithms were evaluated: support vector machines, Adaboost decision trees, and random forests. The algorithms were trained on a database of canonical flow configurations for which validated direct numerical simulation or large eddy simulation results were available, and were used to classify RANS results on a point-by-point basis as having either high or low uncertainty, based on the breakdown of specific RANS modeling assumptions. Classifiers were developed for three different basic RANS eddy viscosity model assumptions: the isotropy of the eddy viscosity, the linearity of the Boussinesq hypothesis, and the non-negativity of the eddy viscosity. It is shown that these classifiers are able to generalize to flows substantially different from those on which they were trained. As a result, feature selection techniques, model evaluation, and extrapolation detection are discussed in the context of turbulence modeling applications.
Ling, Julia; Templeton, Jeremy Alan
2015-08-04
Reynolds Averaged Navier Stokes (RANS) models are widely used in industry to predict fluid flows, despite their acknowledged deficiencies. Not only do RANS models often produce inaccurate flow predictions, but there are very limited diagnostics available to assess RANS accuracy for a given flow configuration. If experimental or higher fidelity simulation results are not available for RANS validation, there is no reliable method to evaluate RANS accuracy. This paper explores the potential of utilizing machine learning algorithms to identify regions of high RANS uncertainty. Three different machine learning algorithms were evaluated: support vector machines, Adaboost decision trees, and random forests.more » The algorithms were trained on a database of canonical flow configurations for which validated direct numerical simulation or large eddy simulation results were available, and were used to classify RANS results on a point-by-point basis as having either high or low uncertainty, based on the breakdown of specific RANS modeling assumptions. Classifiers were developed for three different basic RANS eddy viscosity model assumptions: the isotropy of the eddy viscosity, the linearity of the Boussinesq hypothesis, and the non-negativity of the eddy viscosity. It is shown that these classifiers are able to generalize to flows substantially different from those on which they were trained. As a result, feature selection techniques, model evaluation, and extrapolation detection are discussed in the context of turbulence modeling applications.« less
A Parallel Newton-Krylov-Schur Algorithm for the Reynolds-Averaged Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Osusky, Michal
be readily extended to high order, in combination with the Newton-Krylov-Schur iterative method to produce a powerful parallel algorithm for the numerical solution of the Reynolds-averaged Navier-Stokes equations. The algorithm can efficiently solve the flow over a range of clean geometries, making it suitable for use at the core of an optimization algorithm.
Large Eddy/Reynolds-Averaged Navier-Stokes Simulations of CUBRC Base Heating Experiments
NASA Technical Reports Server (NTRS)
Salazar, Giovanni; Edwards, Jack R.; Amar, Adam J.
2012-01-01
ven with great advances in computational techniques and computing power during recent decades, the modeling of unsteady separated flows, such as those encountered in the wake of a re-entry vehicle, continues to be one of the most challenging problems in CFD. Of most interest to the aerothermodynamics community is accurately predicting transient heating loads on the base of a blunt body, which would result in reduced uncertainties and safety margins when designing a re-entry vehicle. However, the prediction of heat transfer can vary widely depending on the turbulence model employed. Therefore, selecting a turbulence model which realistically captures as much of the flow physics as possible will result in improved results. Reynolds Averaged Navier Stokes (RANS) models have become increasingly popular due to their good performance with attached flows, and the relatively quick turnaround time to obtain results. However, RANS methods cannot accurately simulate unsteady separated wake flows, and running direct numerical simulation (DNS) on such complex flows is currently too computationally expensive. Large Eddy Simulation (LES) techniques allow for the computation of the large eddies, which contain most of the Reynolds stress, while modeling the smaller (subgrid) eddies. This results in models which are more computationally expensive than RANS methods, but not as prohibitive as DNS. By complimenting an LES approach with a RANS model, a hybrid LES/RANS method resolves the larger turbulent scales away from surfaces with LES, and switches to a RANS model inside boundary layers. As pointed out by Bertin et al., this type of hybrid approach has shown a lot of promise for predicting turbulent flows, but work is needed to verify that these models work well in hypersonic flows. The very limited amounts of flight and experimental data available presents an additional challenge for researchers. Recently, a joint study by NASA and CUBRC has focused on collecting heat transfer data
NASA Technical Reports Server (NTRS)
Mehta, Unmeel; Lomax, Harvard
1981-01-01
During the past five years, numerous pioneering archival publications have appeared that have presented computer solutions of the mass-weighted, time-averaged Navier-Stokes equations for transonic problems pertinent to the aircraft industry. These solutions have been pathfinders of developments that could evolve into a major new technological capability, namely the computational Navier-Stokes technology, for the aircraft industry. So far these simulations have demonstrated that computational techniques, and computer capabilities have advanced to the point where it is possible to solve forms of the Navier-Stokes equations for transonic research problems. At present there are two major shortcomings of the technology: limited computer speed and memory, and difficulties in turbulence modelling and in computation of complex three-dimensional geometries. These limitations and difficulties are the pacing items of the continuing developments, although the one item that will most likely turn out to be the most crucial to the progress of this technology is turbulence modelling. The objective of this presentation is to discuss the state of the art of this technology and suggest possible future areas of research. We now discuss some of the flow conditions for which the Navier-Stokes equations appear to be required. On an airfoil there are four different types of interaction of a shock wave with a boundary layer: (1) shock-boundary-layer interaction with no separation, (2) shock-induced turbulent separation with immediate reattachment (we refer to this as a shock-induced separation bubble), (3) shock-induced turbulent separation without reattachment, and (4) shock-induced separation bubble with trailing edge separation.
NASA Technical Reports Server (NTRS)
Cappelli, Daniele; Mansour, Nagi N.
2012-01-01
Separation can be seen in most aerodynamic flows, but accurate prediction of separated flows is still a challenging problem for computational fluid dynamics (CFD) tools. The behavior of several Reynolds Averaged Navier-Stokes (RANS) models in predicting the separated ow over a wall-mounted hump is studied. The strengths and weaknesses of the most popular RANS models (Spalart-Allmaras, k-epsilon, k-omega, k-omega-SST) are evaluated using the open source software OpenFOAM. The hump ow modeled in this work has been documented in the 2004 CFD Validation Workshop on Synthetic Jets and Turbulent Separation Control. Only the baseline case is treated; the slot flow control cases are not considered in this paper. Particular attention is given to predicting the size of the recirculation bubble, the position of the reattachment point, and the velocity profiles downstream of the hump.
Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor
NASA Technical Reports Server (NTRS)
Fulton, Jesse A.; Edwards, Jack R.; Hassan, Hassan A.; Rockwell, Robert; Goyne, Christopher; McDaniel, James; Smith, Chad; Cutler, Andrew; Johansen, Craig; Danehy, Paul M.; Kouchi, Toshinori
2012-01-01
Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused
F-16XL Hybrid Reynolds-Averaged Navier-Stokes/Large Eddy Simulation on Unstructured Grids
NASA Technical Reports Server (NTRS)
Park, Michael A.; Abdol-Hamid, Khaled S.; Elmiligui, Alaa
2015-01-01
This study continues the Cranked Arrow Wing Aerodynamics Program, International (CAWAPI) investigation with the FUN3D and USM3D flow solvers. CAWAPI was established to study the F-16XL, because it provides a unique opportunity to fuse fight test, wind tunnel test, and simulation to understand the aerodynamic features of swept wings. The high-lift performance of the cranked-arrow wing planform is critical for recent and past supersonic transport design concepts. Simulations of the low speed high angle of attack Flight Condition 25 are compared: Detached Eddy Simulation (DES), Modi ed Delayed Detached Eddy Simulation (MDDES), and the Spalart-Allmaras (SA) RANS model. Iso- surfaces of Q criterion show the development of coherent primary and secondary vortices on the upper surface of the wing that spiral, burst, and commingle. SA produces higher pressure peaks nearer to the leading-edge of the wing than flight test measurements. Mean DES and MDDES pressures better predict the flight test measurements, especially on the outer wing section. Vorticies and vortex-vortex interaction impact unsteady surface pressures. USM3D showed many sharp tones in volume points spectra near the wing apex with low broadband noise and FUN3D showed more broadband noise with weaker tones. Spectra of the volume points near the outer wing leading-edge was primarily broadband for both codes. Without unsteady flight measurements, the flight pressure environment can not be used to validate the simulations containing tonal or broadband spectra. Mean forces and moment are very similar between FUN3D models and between USM3D models. Spectra of the unsteady forces and moment are broadband with a few sharp peaks for USM3D.
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.; Greenblatt, David
2007-01-01
This is an expanded version of a limited-length paper that appeared at the 5th International Symposium on Turbulence and Shear Flow Phenomena by the same authors. A computational study was performed for steady and oscillatory flow control over a hump model with flow separation to assess how well the steady and unsteady Reynolds-averaged Navier-Stokes equations predict trends due to Reynolds number, control magnitude, and control frequency. As demonstrated in earlier studies, the hump model case is useful because it clearly demonstrates a failing in all known turbulence models: they under-predict the turbulent shear stress in the separated region and consequently reattachment occurs too far downstream. In spite of this known failing, three different turbulence models were employed to determine if trends can be captured even though absolute levels are not. Overall the three turbulence models showed very similar trends as experiment for steady suction, but only agreed qualitatively with some of the trends for oscillatory control.
NASA Astrophysics Data System (ADS)
Yang, Yang; Kær, Søren Knudsen
2012-04-01
The flow structure of one isothermal swirling case in the Sydney swirl flame database was studied using two numerical methods. Results from the Reynolds-averaged Navier-Stokes (RANS) approach and large eddy simulation (LES) were compared with experimental measurements. The simulations were applied in two different Cartesian grids which were investigated by a grid independence study for RANS and a post-estimator for LES. The RNG k-ɛ turbulence model was used in RANS and dynamic Smagorinsky-Lilly model was used as the sub-grid scale model in LES. A validation study and cross comparison of ensemble average and root mean square (RMS) results showed LES outperforms RANS statistic results. Flow field results indicated that both approaches could capture dominant flow structures, like vortex breakdown (VB), and precessing vortex core (PVC). Streamlines indicate that the formation mechanisms of VB deducted from the two methods were different. The vorticity field was also studied using a velocity gradient based method. This research gained in-depth understanding of isothermal swirling flow.
NASA Astrophysics Data System (ADS)
Schoenawa, Stefan; Hartmann, Ralf
2014-04-01
In this article we consider the development of Discontinuous Galerkin (DG) methods for the numerical approximation of the Reynolds-averaged Navier-Stokes (RANS) equations with the shear-stress transport (SST) model by Menter. This turbulence model is based on a blending of the Wilcox k-ω model used near the wall and the k-ɛ model used in the rest of the domain where the blending functions depend on the distance to the nearest wall. For the computation of the distance of each quadrature point in the domain to the nearest of the curved, piecewise polynomial wall boundaries, we propose a stabilized continuous finite element (FE) discretization of the eikonal equation. Furthermore, we propose a new wall boundary condition for the dissipation rate ω based on the projection of the analytic near-wall behavior of ω onto the discrete ansatz space of the DG discretization. Finally, we introduce an artificial viscosity to the discretization of the turbulence kinetic energy (k-)equation to suppress oscillations of k near the underresolved boundary layer edge. The wall distance computation based on the continuous FE discretization of the eikonal equation is demonstrated for an internal and three external/aerodynamic flow geometries including a three-element high-lift configuration. The DG discretization of the RANS equations with the SST model is demonstrated for turbulent flows past a flat plate and the RAE2822 airfoil (Cases 9 and 10). The results are compared to the underlying k-ω model and experimental data.
NASA Astrophysics Data System (ADS)
Hillewaere, J.; Dooms, D.; Van Quekelberghe, B.; Degroote, J.; Vierendeels, J.; De Roeck, G.; Lombaert, G.; Degrande, G.
2012-04-01
During a storm in October 2002, wind induced ovalling vibrations were observed on several empty silos of a closely spaced group (pitch-to-diameter ratio of 1.05) consisting of 8 by 5 silos in the port of Antwerp (Belgium). Numerical simulations of the turbulent wind flow are performed to clarify the occurrence of the observed ovalling vibrations near the lee side corner of the group by studying the dynamic wind pressures on the silo surfaces and linking to the dynamic properties of the silo structures. As the orientation of the group largely affects the pressure distribution around the cylinders of the group, the influence of the angle of incidence of the wind flow on these ovalling vibrations is examined while other parameters, such as spacing ratio and Reynolds number are unchanged. To achieve results within a reasonable computation time, 2D unsteady Reynolds averaged Navier-Stokes (URANS) equations using Menter's shear stress transport turbulence model were performed. In order to elucidate the influence of the applied turbulence model and to qualitatively validate the spatial and temporal discretization of the 2D highly turbulent post-critical (Re=1.24×107) flow simulations for the silo group, single cylinder simulations were used. The geometric resemblance of the group arrangement with rectangular cylinders on the one hand and of the interstitial spaces with tube arrays (e.g. heat exchangers) on the other hand is used to qualitatively compare the observed flow phenomena. The simulations show that the silo group can be treated neither as a tube array nor as a solid bluff body. Subsequent linking of dynamic wind pressures to dynamic properties of the silo structures reveals strong narrow band frequency peaks in the turbulent pressure coefficient spectra of the silos near the lee side corners of the group that match the structural natural frequencies of the third and fourth ovalling mode shape of the silos. This match indicates a forced, resonant response which
NASA Astrophysics Data System (ADS)
Gorlé, C.; Iaccarino, G.
2013-05-01
Reynolds-averaged Navier-Stokes (RANS) simulations are a practical approach for solving complex multi-physics turbulent flows, but the underlying assumptions of the turbulence models introduce errors and uncertainties in the simulation outcome. The flow in scramjet combustors is an example of such a complex flow and the accurate characterization of safety and operability limits of these engines using RANS simulations requires an assessment of the model uncertainty. The objective of this paper is to present a framework for the epistemic uncertainty quantification of turbulence and mixing models in RANS simulations. The capabilities of the methodology will be demonstrated by performing simulations of the mixing of an underexpanded jet in a supersonic cross flow, which involves many flow features observed in scramjet engines. The fundamental sources of uncertainty in the RANS simulations are the models used for the Reynolds stresses in the momentum equations and the turbulent scalar fluxes in the scalar transport equations. The methodology consists in directly perturbing the modeled quantities in the equations, thereby establishing a method that is completely independent of the initial model form to overcome the limitations of traditional sensitivity studies. The perturbations are defined in terms of the decomposed Reynolds stress tensor, i.e., the tensor magnitude and the eigenvalues and eigenvectors of the normalized anisotropy tensor. The turbulent scalar fluxes are perturbed by using the perturbed Reynolds stresses in a generalized gradient diffusion model formulation and by changing the model constant. The perturbations were parameterized based on a comparison between the Reynolds stresses obtained from a baseline RANS simulation and those obtained from a large-eddy simulation database. Subsequently an optimization problem was solved, varying the parameters in the perturbation functions to maximize a quantity of interest that quantifies the downstream mixing. The
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.
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.
NASA Astrophysics Data System (ADS)
Schilling, Oleg
2012-11-01
A multicomponent, weighted essentially nonoscillatory implementation of several four-equation K- ɛ and K- L based Reynolds-averaged Navier-Stokes models is used to simulate reshocked Richtmyer-Meshkov turbulent mixing at various Mach and Atwood numbers. One class of models is based on mechanical turbulence coupled to scalar variance and its dissipation rate, and the other is based on mechanical turbulence coupled to mass flux and the density-specific volume correlation. The predicted evolution of the mixing layer, molecular mixing and other quantities obtained from these models are systematically intercompared, as well as compared to experimental shock tube data. The relative advantages and disadvantages of the various models are discussed. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Mackay, Kyle K.; Schilling, Oleg
2012-11-01
An implicit-in-time multicomponent, weighted essentially nonoscillatory implementation of a four-equation K- ɛ based Reynolds-averaged Navier-Stokes model is used to simulate Rayleigh-Taylor turbulent mixing at Atwood numbers ranging from 0 . 05 - 0 . 9 . The mechanical turbulence equations are coupled to modeled transport equations for the scalar variance and its dissipation rate. The predicted evolution of the mixing layer, molecular mixing and other quantities are compared to available experimental data, as well as to analytical self-similar solutions. The predictive capability of the model is evaluated, and several parametric studies are also presented. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.
NASA Astrophysics Data System (ADS)
Kovalev, R. V.; Kudryavtsev, V. V.; Churakov, D. A.
2015-06-01
This paper deals with different engineering correlations for prediction of transitional hypersonic flows both for transition onset prediction and transition zone description. In spite of recent progress in simulation of stability equations and direct methods, these approaches remain main tool in everyday design practice. Very close approach is using correlations incorporated in the Reynolds-averaged Navier-Stokes (RANS) equations, e. g., γ-Reθ transition model. But this model has some limitations when applied to supersonic flows and some approaches for subduing these limitations are considered. Also, the effect of local heating/cooling on laminar-turbulent transition (LTT) on sharp cone is analyzed on the basis of both algebraic correlations and RANS models.
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.
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.
NASA Astrophysics Data System (ADS)
Schilling, Oleg; Finn, Thomas
2013-11-01
A multicomponent, four-equation Reynolds-averaged Navier-Stokes model based on the K- ɛ model and coupled to density variance and density variance dissipation rate equations is used to simulate Rayleigh-Taylor turbulent mixing experiments previously performed in a gas channel facility at Texas A&M University at Atwood numbers 0 . 04 and 0 . 6 [A. Banerjee et al., J. Fluid Mech. 659, 127 (2010)]. The model was previously used to simulate a large set of Richtmyer-Meshkov instability experiments, and is implemented in a weighted essentially nonoscillatory code. Model predictions of the Rayleigh-Taylor mixing layer width growth, volume fraction, mixing layer parameters, molecular mixing parameter, density-velocity correlation, and turbulent kinetic energy are compared to experimental data. The sensitivity of the predicted quantities to variations in the principal model coefficients and initial conditions is discussed. This work was performed under the auspices of the U.S. Department of Energy by LLNL under Contract DE-AC52-07NA27344, and by Texas A&M University under the Department of Energy National Nuclear Security Administration research grant DE-FG52-09NA29462.
NASA Astrophysics Data System (ADS)
Xiao, Heng; Wu, Jinlong; Wang, Jianxun; Sun, Rui; Roy, Christopher J.
2015-11-01
For many practical flows, the turbulence models are the most important source of uncertainty in Reynolds-Averaged Navier-Stokes (RANS) predictions. In this work, we develop an open-box, physics-informed Bayesian framework for quantifying the model-form uncertainties in RANS simulations. Uncertainties are introduced directly to the Reynolds stresses and are represented with compact parameterization accounting for empirical prior knowledge and physical constraints (e.g., realizability, smoothness, and symmetry). An iterative ensemble Kalman method is used to incorporate the prior information with available observation data in a Bayesian framework to posterior distributions of the Reynolds stresses and other quantities of interest. Two representative cases, the flow over periodic hills and the flow in a square duct, are used to evaluate the performance of the proposed framework. Simulation results suggest that the obtained posterior mean has significantly better agreement with the benchmark data compared to the baseline simulation, even with very sparse observations. At most locations, the posterior distribution adequately represents the model-form uncertainties.
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.
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.
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.
NASA Astrophysics Data System (ADS)
Morán-López, J. T.; Schilling, O.
2014-05-01
Turbulent mixing generated by shock-driven acceleration of a perturbed interface is simulated using a new multi-component Reynolds-averaged Navier-Stokes (RANS) model closed with a two-equation - model. The model is implemented in a hydrodynamics code using a third-order weighted essentially non-oscillatory finite-difference method for the advection terms and a second-order central difference method for the gradients in the source and diffusion terms. In the present reshocked Richtmyer-Meshkov instability and mixing study, an incident shock with Mach number is generated in air and progresses into a sulfur hexafluoride test section. The time evolution of the predicted mixing layer widths corresponding to six shock tube test section lengths are compared with experimental measurements and three-dimensional multi-mode numerical simulations. The mixing layer widths are also compared with the analytical self-similar power-law solution of the simplified model equations prior to reshock. A set of model coefficients and initial conditions specific to these six experiments is established, for which the widths before and after reshock agree very well with experimental and numerical simulation data. A second set of general coefficients that accommodates a broader range of incident shock Mach numbers, Atwood numbers, and test section lengths is also established by incorporating additional experimental data for , , and with and with and previous RANS modeling. Terms in the budgets of the turbulent kinetic energy and dissipation rate equations are examined to evaluate the relative importance of turbulence production, dissipation and diffusion mechanisms during mixing. Convergence results for the mixing layer widths, mean fields, and turbulent fields under grid refinement are presented for each of the cases.
NASA Astrophysics Data System (ADS)
Spotts, Nathan
As modern trends in commercial aircraft design move toward high-bypass-ratio fan systems of increasing diameter with shorter, nonaxisymmetric nacelle geometries, inlet distortion is becoming common in all operating regimes. The distortion may induce aerodynamic instabilities within the fan system, leading to catastrophic damage to fan blades, should the surge margin be exceeded. Even in the absence of system instability, the heterogeneity of the flow affects aerodynamic performance significantly. Therefore, an understanding of fan-distortion interaction is critical to aircraft engine system design. This thesis research elucidates the complex fluid dynamics and fan-distortion interaction by means of computational fluid dynamics (CFD) modeling of a complete engine fan system; including rotor, stator, spinner, nacelle and nozzle; under conditions typical of those encountered by commercial aircraft. The CFD simulations, based on a Reynolds-averaged Navier-Stokes (RANS) approach, were unsteady, three-dimensional, and of a full-annulus geometry. A thorough, systematic validation has been performed for configurations from a single passage of a rotor to a full-annulus system by comparing the predicted flow characteristics and aerodynamic performance to those found in literature. The original contributions of this research include the integration of a complete engine fan system, based on the NASA rotor 67 transonic stage and representative of the propulsion systems in commercial aircraft, and a benchmark case for unsteady RANS simulations of distorted flow in such a geometry under realistic operating conditions. This study is unique in that the complex flow dynamics, resulting from fan-distortion interaction, were illustrated in a practical geometry under realistic operating conditions. For example, the compressive stage is shown to influence upstream static pressure distributions and thus suppress separation of flow on the nacelle. Knowledge of such flow physics is
NASA Astrophysics Data System (ADS)
Arko, Bryan M.
Design trends for the low-pressure turbine (LPT) section of modern gas turbine engines include increasing the loading per airfoil, which promises a decreased airfoil count resulting in reduced manufacturing and operating costs. Accurate Reynolds-Averaged Navier-Stokes predictions of separated boundary layers and transition to turbulence are needed, as the lack of an economical and reliable computational model has contributed to this high-lift concept not reaching its full potential. Presented here for what is believed to be the first time applied to low-Re computations of high-lift linear cascade simulations is the Abe-Kondoh-Nagano (AKN) linear low-Re two-equation turbulence model which utilizes the Kolmogorov velocity scale for improved predictions of separated boundary layers. A second turbulence model investigated is the Kato-Launder modified version of the AKN, denoted MPAKN, which damps turbulent production in highly strained regions of flow. Fully Laminar solutions have also been calculated in an effort to elucidate the transitional quality of the turbulence model solutions. Time accurate simulations of three modern high-lift blades at a Reynolds number of 25,000 are compared to experimental data and higher-order computations in order to judge the accuracy of the results, where it is shown that the RANS simulations with highly refined grids can produce both quantitatively and qualitatively similar separation behavior as found in experiments. In particular, the MPAKN model is shown to predict the correct boundary layer behavior for all three blades, and evidence of transition is found through inspection of the components of the Reynolds Stress Tensor, spectral analysis, and the turbulence production parameter. Unfortunately, definitively stating that transition is occurring becomes an uncertain task, as similar evidence of the transition process is found in the Laminar predictions. This reveals that boundary layer reattachment may be a result of laminar
NASA Technical Reports Server (NTRS)
Brown, James L.
2014-01-01
Examined is sensitivity of separation extent, wall pressure and heating to variation of primary input flow parameters, such as Mach and Reynolds numbers and shock strength, for 2D and Axisymmetric Hypersonic Shock Wave Turbulent Boundary Layer interactions obtained by Navier-Stokes methods using the SST turbulence model. Baseline parametric sensitivity response is provided in part by comparison with vetted experiments, and in part through updated correlations based on free interaction theory concepts. A recent database compilation of hypersonic 2D shock-wave/turbulent boundary layer experiments extensively used in a prior related uncertainty analysis provides the foundation for this updated correlation approach, as well as for more conventional validation. The primary CFD method for this work is DPLR, one of NASA's real-gas aerothermodynamic production RANS codes. Comparisons are also made with CFL3D, one of NASA's mature perfect-gas RANS codes. Deficiencies in predicted separation response of RANS/SST solutions to parametric variations of test conditions are summarized, along with recommendations as to future turbulence approach.
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.
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.
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.
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.
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.
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.
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)
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.
NASA Technical Reports Server (NTRS)
Chaussee, Denny S.
1988-01-01
A computational fluid dynamics tool has been developed capable of analyzing the viscous supersonic/hypersonic flow about realistic configurations. This techniques can predict the flow in regions of canopies, wings, and canards in addition to the usual simple symmetric configurations. It also allows for interactions between aerodynamic surfaces such as the vortex interaction between canards and wings.
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.
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.
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.
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.
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.
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)
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.
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)
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)
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)
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 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.
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.
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 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)
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
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)
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.
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.
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 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.
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.
Three-dimensional Navier-Stokes analysis of turbine passage heat transfer
NASA Technical Reports Server (NTRS)
Ameri, Ali A.; Arnone, Andrea
1991-01-01
The three-dimensional Reynolds-averaged Navier-Stokes equations are numerically solved to obtain the pressure distribution and heat transfer rates on the endwalls and the blades of two linear turbine cascades. The TRAF3D code which has recently been developed in a joint project between researchers from the University of Florence and NASA Lewis Research Center is used. The effect of turbulence is taken into account by using the eddy viscosity hypothesis and the two-layer mixing length model of Baldwin and Lomax. Predictions of surface heat transfer are made for Langston's cascade and compared with the data obtained for that cascade by Graziani. The comparison was found to be favorable. The code is also applied to a linear transonic rotor cascade to predict the pressure distributions and heat transfer rates.
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
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
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.
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 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.
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.
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
Numerical solutions of the complete Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Hassan, H. A.
1988-01-01
The physical phenomena within supersonic flows that sustain chemical reactions are investigated. An earlier study to develop accurate physical models for supersonic reacting flowfields focused on 2-D laminar shear layers. The objective is to examine the mixing and subsequent combustion within turbulent reacting shear layers. To conduct this study, a computer program has been written to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method uses a cell-centered finite volume approach and a Runge Kutta time stepping scheme. The Reynolds averaged equations are closed using the eddy viscosity concept. Several zero-equation models have been tested by making calculations for an H2-air nonreacting coaxial jet flow. Comparisons made with experimental data show that Cohen's eddy viscosity model provides best agreement. The finite rate chemistry model used in the study of 2-D laminar shear layers is incorporated into the computer program and data is compared from a recent experiment performed at NASA Langley.
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
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)
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.
Partially-Averaged Navier Stokes Model for Turbulence: Implementation and Validation
NASA Technical Reports Server (NTRS)
Girimaji, Sharath S.; Abdol-Hamid, Khaled S.
2005-01-01
Partially-averaged Navier Stokes (PANS) is a suite of turbulence closure models of various modeled-to-resolved scale ratios ranging from Reynolds-averaged Navier Stokes (RANS) to Navier-Stokes (direct numerical simulations). The objective of PANS, like hybrid models, is to resolve large scale structures at reasonable computational expense. The modeled-to-resolved scale ratio or the level of physical resolution in PANS is quantified by two parameters: the unresolved-to-total ratios of kinetic energy (f(sub k)) and dissipation (f(sub epsilon)). The unresolved-scale stress is modeled with the Boussinesq approximation and modeled transport equations are solved for the unresolved kinetic energy and dissipation. In this paper, we first present a brief discussion of the PANS philosophy followed by a description of the implementation procedure and finally perform preliminary evaluation in benchmark problems.
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.
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.
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)
Rumsey, Christopher L.
2009-01-01
In current practice, it is often difficult to draw firm conclusions about turbulence model accuracy when performing multi-code CFD studies ostensibly using the same model because of inconsistencies in model formulation or implementation in different codes. This paper describes an effort to improve the consistency, verification, and validation of turbulence models within the aerospace community through a website database of verification and validation cases. Some of the variants of two widely-used turbulence models are described, and two independent computer codes (one structured and one unstructured) are used in conjunction with two specific versions of these models to demonstrate consistency with grid refinement for several representative problems. Naming conventions, implementation consistency, and thorough grid resolution studies are key factors necessary for success.
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.
Aerodynamics of thrust vectoring by Navier-Stokes solutions
NASA Technical Reports Server (NTRS)
Tseng, Jing-Biau; Lan, C. Edward
1991-01-01
Induced aerodynamics from thrust vectoring are investigated by a computational fluid dynamic method. A thin-layer Reynolds-averaged Navier-Stokes code with multiblock capability is used. Jet properties are specified on the nozzle exit plane to simulate the jet momentum. Results for a rectangular jet in a cross flow are compared with data to verify the code. Further verification of the calculation is made by comparing the numerical results with transonic data for a wing-body combination. Additional calculations were performed to elucidate the following thrust vectoring effects: the thrust vectoring effect on shock and expansion waves, induced effects on nearby surfaces, and the thrust vectoring effect on the leading edge vortex.
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.
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.
Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Swanson, Roy C., Jr.; Turkel, Eli; Rossow, C.-C.
2007-01-01
The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.
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.
Properties of the Residual Stress of the Temporally Filtered Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Pruett, C. D.; Gatski, T. B.; Grosch, C. E.; Thacker, W. D.
2002-01-01
The development of a unifying framework among direct numerical simulations, large-eddy simulations, and statistically averaged formulations of the Navier-Stokes equations, is of current interest. Toward that goal, the properties of the residual (subgrid-scale) stress of the temporally filtered Navier-Stokes equations are carefully examined. Causal time-domain filters, parameterized by a temporal filter width 0 less than Delta less than infinity, are considered. For several reasons, the differential forms of such filters are preferred to their corresponding integral forms; among these, storage requirements for differential forms are typically much less than for integral forms and, for some filters, are independent of Delta. The behavior of the residual stress in the limits of both vanishing and in infinite filter widths is examined. It is shown analytically that, in the limit Delta to 0, the residual stress vanishes, in which case the Navier-Stokes equations are recovered from the temporally filtered equations. Alternately, in the limit Delta to infinity, the residual stress is equivalent to the long-time averaged stress, and the Reynolds-averaged Navier-Stokes equations are recovered from the temporally filtered equations. The predicted behavior at the asymptotic limits of filter width is further validated by numerical simulations of the temporally filtered forced, viscous Burger's equation. Finally, finite filter widths are also considered, and a priori analyses of temporal similarity and temporal approximate deconvolution models of the residual stress are conducted.
Partially-averaged Navier-Stokes method for turbulent thermal plume
NASA Astrophysics Data System (ADS)
Kumar, Rajesh; Dewan, Anupam
2015-12-01
In this paper, the partially-averaged Navier-Stokes (PANS) simulation is performed for a turbulent thermal plume. The aim of the paper is to assess the PANS method for modeling buoyancy-driven flows at a reasonable computational cost. PANS is a turbulence closure model which is developed to be used as a bridging model ranging from the direct numerical simulation to the Reynolds-averaged Navier-Stokes simulation by varying the level of resolution. The PANS computations are performed for various values of the filter-width to evaluate the sensitivity of the filter-widths to the computed flow statistics. The present simulations have been carried out employing a source code buoyantPimpleFOAM based on the OpenFOAM platform. In order to capture the effect of buoyancy on turbulence, the generalized gradient diffusion hypothesis is employed to model the production of turbulence due to buoyancy. A detailed comparison of the time-averaged and turbulent statistics obtained from the PANS simulations with the experimental data and LES results reported in the literature has been presented. The present results have also been compared with the results of the unsteady Reynolds-averaged Navier-Stokes solutions. The PANS model is shown to enhance the computing capability significantly in predicting buoyancy-driven flows compared with those of URANS model. Finally, various important unsteady flow structures of turbulent thermal plume have been visualized from the instantaneous flow statistics obtained using the PANS simulations.
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.
A fully vectorized numerical solution of the incompressible Navier-Stokes equations. Ph.D. Thesis
NASA Technical Reports Server (NTRS)
Patel, N.
1983-01-01
A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.
Transonic Navier-Stokes solutions of three-dimensional afterbody flows
NASA Technical Reports Server (NTRS)
Compton, William B., III; Thomas, James L.; Abeyounis, William K.; Mason, Mary L.
1989-01-01
The performance of a three-dimensional Navier-Stokes solution technique in predicting the transonic flow past a nonaxisymmetric nozzle was investigated. The investigation was conducted at free-stream Mach numbers ranging from 0.60 to 0.94 and an angle of attack of 0 degrees. The numerical solution procedure employs the three-dimensional, unsteady, Reynolds-averaged Navier-Stokes equations written in strong conservation form, a thin layer assumption, and the Baldwin-Lomax turbulence model. The equations are solved by using the finite-volume principle in conjunction with an approximately factored upwind-biased numerical algorithm. In the numerical procedure, the jet exhaust is represented by a solid sting. Wind-tunnel data with the jet exhaust simulated by high pressure air were also obtained to compare with the numerical calculations.
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.
Navier-Stokes simulation of transonic wing flow fields using a zonal grid approach
NASA Technical Reports Server (NTRS)
Chaderjian, Neal M.
1988-01-01
The transonic Navier-Stokes code was used to simulate flow fields about isolated wings for workshop wind-tunnel and free-air cases using the thin-layer Reynolds-averaged Navier-Stokes equations. An implicit finite-difference scheme based on a diagonal version of the Beam-Warming algorithm was used to integrate the governing equations. A zonal grid approach was used to allow efficient grid refinement near the wing surface. The flow field was sensitive to the turbulent transition model, and flow unsteadiness was observed for a wind-tunnel case but not for the corresponding free-air case. The specification of experimental pressure at the wind-tunnel exit plane is the primary reason for the difference of these two numerical 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.
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.
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.
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.
Navier-Stokes simulations of unsteady transonic flow phenomena
NASA Technical Reports Server (NTRS)
Atwood, C. A.
1992-01-01
Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical
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.
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.
Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description
NASA Technical Reports Server (NTRS)
Towne, Charles E.; Schwab, John R.; Bui, Trong T.
1993-01-01
A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.
Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide
NASA Technical Reports Server (NTRS)
Towne, Charles E.; Schwab, John R.; Bui, Trong T.
1993-01-01
A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This User's Guide describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.
Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference
NASA Technical Reports Server (NTRS)
Towne, Charles E.; Schwab, John R.; Bui, Trong T.
1993-01-01
A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.
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.
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.
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.
Multigrid calculations of 3-D turbulent viscous flows
NASA Technical Reports Server (NTRS)
Yokota, Jeffrey W.
1989-01-01
Convergence properties of a multigrid algorithm, developed to calculate compressible viscous flows, are analyzed by a vector sequence eigenvalue estimate. The full 3-D Reynolds-averaged Navier-Stokes equations are integrated by an implicit multigrid scheme while a k-epsilon turbulence model is solved, uncoupled from the flow equations. Estimates of the eigenvalue structure for both single and multigrid calculations are compared in an attempt to analyze the process as well as the results of the multigrid technique. The flow through an annular turbine is used to illustrate the scheme's ability to calculate complex 3-D flows.
NASA Technical Reports Server (NTRS)
Dorney, Daniel J.
1996-01-01
Experimental data from jet-engine tests have indicated that unsteady blade-row interaction effects can have a significant impact on the efficiency of low-pressure turbine stages. Measured turbine efficiencies at takeoff can be as much as two points higher than those at cruise conditions. Preliminary studies indicate that Reynolds number effects may contribute to the lower efficiencies at cruise conditions. In the current study, numerical experiments have been performed to quantify the Reynolds number dependence of unsteady wake/separation bubble interaction on the performance of a low-pressure turbine.
Laskowski, Gregory Michael
2005-12-01
Flows with strong curvature present a challenge for turbulence models, specifically eddy viscosity type models which assume isotropy and a linear and instantaneous equilibrium relation between stress and strain. Results obtained from three different codes and two different linear eddy viscosity turbulence models are compared to a DNS simulation in order to gain some perspective on the turbulence modeling capability of SIERRA/Fuego. The Fuego v2f results are superior to the more common two-layer k-e model results obtained with both a commercial and research code in terms of the concave near wall behavior predictions. However, near the convex wall, including the separated region, little improvement is gained using the v2f model and in general the turbulent kinetic energy prediction is fair at best.
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.
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.
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.
A new nonlinear turbulence model based on Partially-Averaged Navier-Stokes Equations
NASA Astrophysics Data System (ADS)
Liu, J. T.; Wu, Y. L.; Cai, C.; Liu, S. H.; Wang, L. Q.
2013-12-01
Partially-averaged Navier-Stokes (PANS) Model was recognized as a Reynolds-averaged Navier-Stokes (RANS) to direct numerical simulation (DNS) bridging method. PANS model was purported for any filter width-from RANS to DNS. PANS method also shared some similarities with the currently popular URANS (unsteady RANS) method. In this paper, a new PANS model was proposed, which was based on RNG k-ε turbulence model. The Standard and RNG k-ε turbulence model were both isotropic models, as well as PANS models. The sheer stress in those PANS models was solved by linear equation. The linear hypothesis was not accurate in the simulation of complex flow, such as stall phenomenon. The sheer stress here was solved by nonlinear method proposed by Ehrhard. Then, the nonlinear PANS model was set up. The pressure coefficient of the suction side of the NACA0015 hydrofoil was predicted. The result of pressure coefficient agrees well with experimental result, which proves that the nonlinear PANS model can capture the high pressure gradient flow. A low specific centrifugal pump was used to verify the capacity of the nonlinear PANS model. The comparison between the simulation results of the centrifugal pump and Particle Image Velocimetry (PIV) results proves that the nonlinear PANS model can be used in the prediction of complex flow field.
The Proteus Navier-Stokes code. [two and three dimensional computational fluid dynamics
NASA Technical Reports Server (NTRS)
Towne, Charles E.; 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. Proteus solves the Reynolds-averaged, unsteady, compressible Navier-Stokes equations in strong conservation law form. Turbulence is modeled using a Baldwin-Lomax based algebraic eddy viscosity model. In addition, options are available to solve thin layer or Euler equations, and to eliminate the energy equation by assuming constant stagnation enthalpy. An extensive series of validation cases have been run, primarily using the two dimensional planar/axisymmetric version of the code. Several flows were computed that have exact solution such as: fully developed channel and pipe flow; Couette flow with and without pressure gradients; unsteady Couette flow formation; flow near a suddenly accelerated flat plate; flow between concentric rotating cylinders; and flow near a rotating disk. The two dimensional version of the Proteus code has been released, and the three dimensional code is scheduled for release in late 1991.
NASA Astrophysics Data System (ADS)
Luo, Dahai; Yan, Chao; Wang, Xiaoyong
2015-02-01
Separation commonly exists in the flows around flight vehicles and also in the internal combustor flows. Simulation of high-speed turbulent-separated flows using a reliable computational design tool is crucial for the development of supersonic and hypersonic vehicles. In this paper, we present the computational results of supersonic base and ramped-cavity flows at high Reynolds numbers using the partially averaged Navier-Stokes (PANS) method. The current PANS models are based on the Menter SST turbulence model and also the Wilcox k-ω model. Results from PANS simulations are compared in detail with the available experimental data. The effect of the resolution control parameter fk (the ratio of unresolved-to-total kinetic energy) relevant to the PANS method is investigated. More turbulent flow structures are resolved as expected with decreasing fk, but it does not mean better results can be obtained. Spatially varying and dynamically updated fk in PANS simulations has been performed. Results from variable fk PANS simulations show good agreement with the experiment and great improvement when compared to Reynolds averaged Navier-Stokes (RANS) computation and constant fk PANS simulations.
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.
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.
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.
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.
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.
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 vertex-based finite-volume algorithm for the Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Chakrabartty, S. K.; Dhanalakshmi, K.
1993-07-01
A vertex-based, finite-volume algorithm has been developed to solve the Reynolds-averaged Navier-Stokes equations without thin-layer approximation. An explicit, five-stage Runge-Kutta, time-stepping scheme has been used for time integration along with different acceleration techniques to reach the steady state. A code employing multi-block grid structure has been developed. This code can accept any type of grid topology. As test cases, the turbulent flow past RAE-2822 and NACA-0012 airfoils, and the laminar flow past a cropped delta wing at ten degrees angle of attack have been computed and the results compared with available numerical and experimental results. The Baldwin-Lomax turbulence model has been used in the case of turbulent flows.
Navier-Stokes and Comprehensive Analysis Performance Predictions of the NREL Phase VI Experiment
NASA Technical Reports Server (NTRS)
Duque, Earl P. N.; Burklund, Michael D.; Johnson, Wayne
2003-01-01
A vortex lattice code, CAMRAD II, and a Reynolds-Averaged Navier-Stoke code, OVERFLOW-D2, were used to predict the aerodynamic performance of a two-bladed horizontal axis wind turbine. All computations were compared with experimental data that was collected at the NASA Ames Research Center 80- by 120-Foot Wind Tunnel. Computations were performed for both axial as well as yawed operating conditions. Various stall delay models and dynamics stall models were used by the CAMRAD II code. Comparisons between the experimental data and computed aerodynamic loads show that the OVERFLOW-D2 code can accurately predict the power and spanwise loading of a wind turbine rotor.
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.
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.
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.
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)
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.
Summary of EASM Turbulence Models in CFL3D With Validation Test Cases
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.; Gatski, Thomas B.
2003-01-01
This paper summarizes the Explicit Algebraic Stress Model in k-omega form (EASM-ko) and in k-epsilon form (EASM-ke) in the Reynolds-averaged Navier-Stokes code CFL3D. These models have been actively used over the last several years in CFL3D, and have undergone some minor modifications during that time. Details of the equations and method for coding the latest versions of the models are given, and numerous validation cases are presented. This paper serves as a validation archive for these models.
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.
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.
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 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.
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.
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.
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.
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.
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
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.
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.
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.
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
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.
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.
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.
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.
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.
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.
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.
Multitasking a three-dimensional Navier-Stokes algorithm on the Cray-2
NASA Technical Reports Server (NTRS)
Swisshelm, Julie M.
1989-01-01
A three-dimensional computational aerodynamics algorithm has been multitasked for efficient parallel execution on the Cray-2. It provides a means for examining the multitasking performance of a complete CFD application code. An embedded zonal multigrid scheme is used to solve the Reynolds-averaged Navier-Stokes equations for an internal flow model problem. The explicit nature of each component of the method allows a spatial partitioning of the computational domain to achieve a well-balanced task load for MIMD computers with vector-processing capability. Experiments have been conducted with both two- and three-dimensional multitasked cases. The best speedup attained by an individual task group was 3.54 on four processors of the Cray-2, while the entire solver yielded a speedup of 2.67 on four processors for the three-dimensional case. The multiprocessing efficiency of various types of computational tasks is examined, performance on two Cray-2s with different memory access speeds is compared, and extrapolation to larger problems is discussed.
Navier-Stokes computation of compressible turbulent flows with a second order closure, part 1
NASA Technical Reports Server (NTRS)
Haminh, Hieu; Kollmann, Wolfgang; Vandromme, Dany
1990-01-01
A second order closure turbulence model for compressible flows is developed and implemented in a 2D Reynolds-averaged Navier-Stokes solver. From the beginning where a kappa-epsilon turbulence model was implemented in the bidiagonal implicit method of MACCORMACK (referred to as the MAC3 code) to the final stage of implementing a full second order closure in the efficient line Gauss-Seidel algorithm, numerous work was done, individually and collectively. Besides the collaboration itself, the final product of this work is a second order closure derived from the Launder, Reece, and Rodi model to account for near wall effects, which has been called FRAME model, which stands for FRench-AMerican-Effort. During the reporting period, two different problems were worked out. The first was to provide Ames researchers with a reliable compressible boundary layer code including a wide collection of turbulence models for quick testing of new terms, both in two equations and in second order closure (LRR and FRAME). The second topic was to complete the implementation of the FRAME model in the MAC5 code. The work related to these two different contributions is reported. dilatation in presence of stron shocks. This work, which has been conducted during a work at the Center for Turbulence Research with Zeman aimed also to cros-check earlier assumptions by Rubesin and Vandromme.
Navier-Stokes Simulation of a Heavy Lift Slowed-Rotor Compound Helicopter Configuration
NASA Technical Reports Server (NTRS)
Allan, Brian G.; Jenkins, Luther N.; Yao, Chung-Sheng; Bartram, Scott M.; Hallissy, Jim B.; Harris, Jerome; Noonan, Kevin W.; Wong, Oliver D.; Jones, Henry E.; Malovrh, Brendon D.; reis, Deane G.; Mace, W. Derry
2009-01-01
Time accurate numerical simulations were performed using the Reynolds-averaged Navier-Stokes (RANS) flow solver OVERFLOW for a heavy lift, slowed-rotor, compound helicopter configuration, tested at the NASA Langley 14- by 22-Foot Subsonic Tunnel. The primary purpose of these simulations is to provide support for the development of a large field of view Particle Imaging Velocimetry (PIV) flow measurement technique supported by the Subsonic Rotary Wing (SRW) project under the NASA Fundamental Aeronautics program. These simulations provide a better understanding of the rotor and body wake flows and helped to define PIV measurement locations as well as requirements for validation of flow solver codes. The large field PIV system can measure the three-dimensional velocity flow field in a 0.914m by 1.83m plane. PIV measurements were performed upstream and downstream of the vertical tail section and are compared to simulation results. The simulations are also used to better understand the tunnel wall and body/rotor support effects by comparing simulations with and without tunnel floor/ceiling walls and supports. Comparisons are also made to the experimental force and moment data for the body and rotor.
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.
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.
PAB3D: Its History in the Use of Turbulence Models in the Simulation of Jet and Nozzle Flows
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Pao, S. Paul; Hunter, Craig A.; Deere, Karen A.; Massey, Steven J.; Elmiligui, Alaa
2006-01-01
This is a review paper for PAB3D s history in the implementation of turbulence models for simulating jet and nozzle flows. We describe different turbulence models used in the simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations use modified linear or nonlinear two-equation models to account for supersonic flow as well as high temperature mixing. Two multiscale-type turbulence models are used for unsteady flow simulations. These models require modifications to the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k-epsilon) model with a RANS/LES 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 (PANS) formulation. All of these models are implemented in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational methods, code implementation, computed results for a wide range of nozzle configurations at various operating conditions, 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.
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.
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
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.
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.
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.
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.
Advancing the theoretical foundation of the partially-averaged Navier-Stokes approach
NASA Astrophysics Data System (ADS)
Reyes, Dasia Ann
The goal of this dissertation is to consolidate the theoretical foundation of variable-resolution (VR) methods in general and the partially-averaged Navier-Stokes (PANS) approach in particular. The accurate simulation of complex turbulent flows remains an outstanding challenge in modern computational fluid dynamics. High-fidelity approaches such as direct numerical simulations (DNS) and large-eddy simulation (LES) are not typically feasible for complex engineering simulations with current computational technologies. Low-fidelity approaches such as Reynolds-averaged Navier-Stokes (RANS), although widely used, are inherently inadequate for turbulent flows with complex flow features. VR bridging methods fill the gap between DNS and RANS by allowing a tunable degree of resolution ranging from RANS to DNS. While the utility of VR methods is well established, the mathematical foundations and physical characterization require further development. This dissertation focuses on the physical attributes of fluctuations in partially-resolved simulations of turbulence. The specific objectives are to: (i) establish a framework for assessing the physical fidelity of VR methods to examine PANS fluctuations; (ii) investigate PANS simulations subject to multiple resolution changes; (iii) examine turbulent transport closure modeling for partially-resolved fields; (iv) examine the effect of filter control parameters in the limit of spectral cut-off in the dissipative region; and (v) validate low-Reynolds number corrections with RANS for eventual implementation with PANS. While the validation methods are carried out in the context of PANS, they are considered appropriate for all VR bridging methods. The key findings of this dissertation are summarized as follows. The Kolmogorov hypotheses are suitably adapted to describe fluctuations of partially-resolved turbulence fields, and the PANS partially-resolved field is physically consistent with the adapted Kolmogorov hypotheses. PANS
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.
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.
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.
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.
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.
NASA Astrophysics Data System (ADS)
Pan, Xinghong
2016-06-01
Consider an axisymmetric suitable weak solution of 3D incompressible Navier-Stokes equations with nontrivial swirl, v =vrer +vθeθ +vzez. Let z denote the axis of symmetry and r be the distance to the z-axis. If the solution satisfies a slightly supercritical assumption (that is, | v | ≤ C(ln | ln r |)/α r for α ∈ [ 0 , 0.028 ] when r is small), then we prove that v is regular. This extends the results in [6,16,18] where regularities under critical assumptions, such as | v | ≤Cr, were proven. As a useful tool in the proof of our main result, an upper-bound estimate to the fundamental solution of the parabolic equation with a critical drift term will be given in the last part of this paper.
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.
An Evaluation of Architectural Platforms for Parallel Navier-Stokes Computations
NASA Technical Reports Server (NTRS)
Jayasimha, D. N.; Hayder, M. E.; Pillay, S. K.
1996-01-01
We study the computational, communication, and scalability characteristics of a computational fluid dynamics application, which solves the time accurate flow field of a jet using the compressible Navier-Stokes equations, on a variety of parallel architecture platforms. The platforms chosen for this study are a cluster of workstations (the LACE experimental testbed at NASA Lewis), a shared memory multiprocessor (the Cray YMP), and distributed memory multiprocessors with different topologies - the IBM SP and the Cray T3D. We investigate the impact of various networks connecting the cluster of workstations on the performance of the application and the overheads induced by popular message passing libraries used for parallelization. The work also highlights the importance of matching the memory bandwidth to the processor speed for good single processor performance. By studying the performance of an application on a variety of architectures, we are able to point out the strengths and weaknesses of each of the example computing platforms.
Parallelizing Navier-Stokes Computations on a Variety of Architectural Platforms
NASA Technical Reports Server (NTRS)
Jayasimha, D. N.; Hayder, M. E.; Pillay, S. K.
1997-01-01
We study the computational, communication, and scalability characteristics of a Computational Fluid Dynamics application, which solves the time accurate flow field of a jet using the compressible Navier-Stokes equations, on a variety of parallel architectural platforms. The platforms chosen for this study are a cluster of workstations (the LACE experimental testbed at NASA Lewis), a shared memory multiprocessor (the Cray YMP), distributed memory multiprocessors with different topologies-the IBM SP and the Cray T3D. We investigate the impact of various networks, connecting the cluster of workstations, on the performance of the application and the overheads induced by popular message passing libraries used for parallelization. The work also highlights the importance of matching the memory bandwidth to the processor speed for good single processor performance. By studying the performance of an application on a variety of architectures, we are able to point out the strengths and weaknesses of each of the example computing platforms.
Recent Analytical and Numerical Results for The Navier-Stokes-Voigt Model and Related Models
NASA Astrophysics Data System (ADS)
Larios, Adam; Titi, Edriss; Petersen, Mark; Wingate, Beth
2010-11-01
The equations which govern the motions of fluids are notoriously difficult to handle both mathematically and computationally. Recently, a new approach to these equations, known as the Voigt-regularization, has been investigated as both a numerical and analytical regularization for the 3D Navier-Stokes equations, the Euler equations, and related fluid models. This inviscid regularization is related to the alpha-models of turbulent flow; however, it overcomes many of the problems present in those models. I will discuss recent work on the Voigt-regularization, as well as a new criterion for the finite-time blow-up of the Euler equations based on their Voigt-regularization. Time permitting, I will discuss some numerical results, as well as applications of this technique to the Magnetohydrodynamic (MHD) equations and various equations of ocean dynamics.
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.
Investigation of Navier-Stokes Code Verification and Design Optimization
NASA Technical Reports Server (NTRS)
Vaidyanathan, Rajkumar
2004-01-01
With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization
Investigation of Navier-Stokes code verification and design optimization
NASA Astrophysics Data System (ADS)
Vaidyanathan, Rajkumar
With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a finer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-epsilonturbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi
Unsteady three-dimensional thin-layer Navier-Stokes solutions for turbomachinery in transonic flow
NASA Astrophysics Data System (ADS)
Chen, Jen Ping
1991-02-01
This dissertation concerns the development of a numerical scheme to simulate the flow field generated by rotating machinery. The implicit Euler solver 'TURBO' developed at MSU and in use at NASA Lewis Research Center for analyzing compressible flows in general rotating machinery was augmented by including viscous effects. This was accomplished by making the thin-layer approximation with a Baldwin-Lomax turbulence model to simplify the Reynolds-averaged Navier-Stokes equations. This code is an implicit finite volume scheme with flux Jacobians evaluated by flux-vector-splitting and residual fluxes by Roe's flux-difference-splitting. Diffusive flux terms were treated explicitly in order to preserve the existing code structure. A modified two pass matrix solver based on the Gauss-Seidel iterative method was used in place of the standard LU approximate factorization to enhance the code's stability. Newton subiterations were used for unsteady flow computations. A localized grid distortion technique was applied for data communication between blocks in relative motion. An inlet boundary condition for preserving total reservoir conditions and a new approach for computing wall shear stress for nonorthogonal grids are also presented. Three engineering problems have been examined by this computational technique with fair results. Hamilton-Standard's advanced propfan design SR7 served as a test case for external flows. NASA's Rotor 67 with a rotating blade now and NASA's Stage 67 with both rotating and static blade rows with a different number of blades in each row test the code's ability for internal flow. Although some improvements are still needed for this code to become a real engineering design tool, the decent results obtained from applying this code to the above three real-world engineering problems suggest that it is moving in the right direction.
Navier-Stokes Simulation of the Canard-Wing-Body Longitudinal Dynamic Stability Characteristics
NASA Technical Reports Server (NTRS)
Tu, Eugene L.; VanDalsem, William R. (Technical Monitor)
1996-01-01
Many modern aircraft are canard-configured for aircraft control and improved aerodynamic performance. Canards can often enhance aircraft cruise performance, maneuverability and agility. For close-coupled canard configurations, the aerodynamic interaction between the canard and wing significantly changes the flow characteristics of the wing. In unsteady flow, such changes in the flow structure and performance of wings can be quite pronounced. Accurate modeling of the unsteady aerodynamics is essential for potential CFD design and analysis of such configurations. A time-accurate numerical simulation is performed to study the unsteady aerodynamic interaction between a canard and wing with emphasis on the effects of the canard on the configuration's dynamic response characteristics. The thin-layer Reynolds-averaged Navier-Stokes Equations with various turbulence models are used in this study. Computations are made on a generic, analytically-defined, close-coupled canard-wing-body configuration which has been the subject of numerous previously published experimental studies during the 1970's to mid-80's. More recently, a series of steady-flow simulations has been performed and published by the author. In the current study, the configuration is given prescribed ramp and oscillatory motions in order to predict characteristics such as the damping-in-pitch and oscillatory longitudinal stability parameters. The current computations are made at high-subsonic and transonic Mach numbers, moderate angles-of- attack from -4 to 20 degrees, and at various pitch rates and reduced frequencies. Comparisons of pressures and integrated force quantities (e.g. lift, drag, pitching moment and selected dynamic parameters) are made with other published computational results and available experimental data. Results showing the unsteady effects of the canard on surface pressures, integrated forces, canard-wing vortex interaction and vortex breakdown will be presented.
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.
SSME thrust chamber simulation using Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Przekwas, A. J.; Singhal, A. K.; Tam, L. T.
1984-01-01
The capability of the PHOENICS fluid dynamics code in predicting two-dimensional, compressible, and reacting flow in the combustion chamber and nozzle of the space shuttle main engine (SSME) was evaluated. A non-orthogonal body fitted coordinate system was used to represent the nozzle geometry. The Navier-Stokes equations were solved for the entire nozzle with a turbulence model. The wall boundary conditions were calculated based on the wall functions which account for pressure gradients. Results of the demonstration test case reveal all expected features of the transonic nozzle flows. Of particular interest are the locations of normal and barrel shocks, and regions of highest temperature gradients. Calculated performance (global) parameters such as thrust chamber flow rate, thrust, and specific impulse are also in good agreement with available data.
Application of Navier-Stokes analysis to stall flutter
NASA Technical Reports Server (NTRS)
Wu, J. C.; Srivastava, R.; Sankar, L. N.
1988-01-01
A solution procedure was developed to investigate the two-dimensional, one- or two-dimensional flutter characteristics of arbitrary airfoils. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. The flutter characteristics of airfoils in subsonic speed at high angles of attack and airfoils in high subsonic and transonic speeds at low angles of attack are investigated. The stall flutter characteristics are also predicted using the same procedure.
Convergence analysis of WLS based solution of Navier Stokes equation
NASA Astrophysics Data System (ADS)
Kosec, G.
2016-06-01
A numerical solution of a Navier-Stokes problem based on the Weighted Least Squares (WLS) approximation of velocity and pressure fields is presented in this paper. The approximation function is constructed over the local support, i.e., a sub cluster of computational nodes. Besides local approximation of the fields also the pressure-velocity algorithm is constructed locally. The presented solution procedure is demonstrated on two classical fluid-flow benchmark tests, i.e., lid-driven cavity and backward-facing step problem. The method is validated through comparison against already published data on regular nodal distributions and convergence analyses. In addition the method is also tested on irregular nodal distributions. Results are presented in terms of cross-section velocity profiles and convergence plots.
The Navier-Stokes Equations in Nonendpoint Borderline Lorentz Spaces
NASA Astrophysics Data System (ADS)
Phuc, Nguyen Cong
2015-12-01
It is shown both locally and globally that {L_t^{∞}(L_x^{3,q})} solutions to the three-dimensional Navier-Stokes equations are regular provided {q≠∞}. Here {L_x^{3,q}}, {0 < q ≤∞}, is an increasing scale of Lorentz spaces containing {L^3_x}. Thus the result provides an improvement of a result by Escauriaza et al. (Uspekhi Mat Nauk 58:3-44, 2003; translation in Russ Math Surv 58, 211-250, 2003), which treated the case q = 3. A new local energy bound and a new {ɛ}-regularity criterion are combined with the backward uniqueness theory of parabolic equations to obtain the result. A weak-strong uniqueness of Leray-Hopf weak solutions in {L_t^{∞}(L_x^{3,q})}, {q≠∞}, is also obtained as a consequence.
Perturbation of eigenvalues of preconditioned Navier-Stokes operators
Elman, H.C.
1996-12-31
We study the sensitivity of algebraic eigenvalue problems associated with matrices arising from linearization and discretization of the steady-state Navier-Stokes equations. In particular, for several choices of preconditioners applied to the system of discrete equations, we derive upper bounds on perturbations of eigenvalues as functions of the viscosity and discretization mesh size. The bounds suggest that the sensitivity of the eigenvalues is at worst linear in the inverse of the viscosity and quadratic in the inverse of the mesh size, and that scaling can be used to decrease the sensitivity in some cases. Experimental results supplement these results and confirm the relatively mild dependence on viscosity. They also indicate a dependence on the mesh size of magnitude smaller than the analysis suggests.
Navier-Stokes analysis of turbine blade heat transfer
NASA Technical Reports Server (NTRS)
Boyle, R. J.
1990-01-01
Comparisons with experimental heat transfer and surface pressures were made for seven turbine vane and blade geometries using a quasi-three-dimensional thin-layer Navier-Stokes analysis. Comparisons are made for cases with both separated and unseparated flow over a range of Reynolds numbers and freestream turbulence intensities. The analysis used a modified Baldwin-Lomax turbulent eddy viscosity mode. Modifications were made to account for the effects of: (1) freestream turbulence on both transition and leading edge heat transfer; (2) strong favorable pressure gradients on relaminarization; and (3) variable turbulent Prandtl number heat transfer. In addition, the effect of heat transfer on the near wall model of Deissler is compared with the Van Driest model.
Navier-Stokes Dynamics by a Discrete Boltzmann Model
NASA Technical Reports Server (NTRS)
Rubinstein, Robet
2010-01-01
This work investigates the possibility of particle-based algorithms for the Navier-Stokes equations and higher order continuum approximations of the Boltzmann equation; such algorithms would generalize the well-known Pullin scheme for the Euler equations. One such method is proposed in the context of a discrete velocity model of the Boltzmann equation. Preliminary results on shock structure are consistent with the expectation that the shock should be much broader than the near discontinuity predicted by the Pullin scheme, yet narrower than the prediction of the Boltzmann equation. We discuss the extension of this essentially deterministic method to a stochastic particle method that, like DSMC, samples the distribution function rather than resolving it completely.
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.
Iterative methods for compressible Navier-Stokes and Euler equations
Tang, W.P.; Forsyth, P.A.
1996-12-31
This workshop will focus on methods for solution of compressible Navier-Stokes and Euler equations. In particular, attention will be focused on the interaction between the methods used to solve the non-linear algebraic equations (e.g. full Newton or first order Jacobian) and the resulting large sparse systems. Various types of block and incomplete LU factorization will be discussed, as well as stability issues, and the use of Newton-Krylov methods. These techniques will be demonstrated on a variety of model transonic and supersonic airfoil problems. Applications to industrial CFD problems will also be presented. Experience with the use of C++ for solution of large scale problems will also be discussed. The format for this workshop will be four fifteen minute talks, followed by a roundtable discussion.
Time Integration Schemes for the Unsteady Navier-stokes Equations
NASA Technical Reports Server (NTRS)
Bijl, Hester; Carpenter, Mark H.; Vatsa, Veer N.
2001-01-01
The efficiency and accuracy of several time integration schemes are investigated for the unsteady Navier-Stokes equations. This study focuses on the efficiency of higher-order Runge-Kutta schemes in comparison with the popular Backward Differencing Formulations. For this comparison an unsteady two-dimensional laminar flow problem is chosen, i.e., flow around a circular cylinder at Re = 1200. It is concluded that for realistic error tolerances (smaller than 10(exp -1)) fourth-and fifth-order Runge-Kutta schemes are the most efficient. For reasons of robustness and computer storage, the fourth-order Runge-Kutta method is recommended. The efficiency of the fourth-order Runge-Kutta scheme exceeds that of second-order Backward Difference Formula by a factor of 2.5 at engineering error tolerance levels (10(exp -1) to 10(exp -2)). Efficiency gains are more dramatic at smaller tolerances.
Incompressible Navier-Stokes computations of rotating flows
NASA Astrophysics Data System (ADS)
Kiris, Cetin; Chang, Leon; Kwak, Dochan; Rogers, Stuart
1993-01-01
Flow through pump components, such as an inducer and an impeller, is efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. The resulting computer code is applied to the flow analysis inside a generic rocket engine pump inducer, a fuel pump impeller, and SSME high-pressure fuel turbopump impeller. Numerical results of inducer flow are compared with experimental measurements. Flow analyses at 80-, 100-, and 120-percent of design conditions are presented.
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 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 Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP
NASA Technical Reports Server (NTRS)
Vadyak, Joseph; Shrewsbury, George E.; Narramore, Jim C.; Montry, Gary; Holst, Terry; Kwak, Dochan (Technical Monitor)
1995-01-01
The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time
Numerical study on comparison of Navier-Stokes and Burgers equations
NASA Astrophysics Data System (ADS)
Ohkitani, Koji; Dowker, Mark
2012-05-01
We compare freely decaying evolution of the Navier-Stokes equations with that of the 3D Burgers equations with the same kinematic viscosity and the same incompressible initial data by using direct numerical simulations. The Burgers equations are well-known to be regular by a maximum principle [A. A. Kiselev and O. A. Ladyzenskaya, "On existence and uniqueness of the solutions of the nonstationary problem for a viscous incompressible fluid," Izv. Akad. Nauk SSSR Ser. Mat. 21, 655 (1957); A. A. Kiselev and O. A. Ladyzenskaya, Am. Math. Soc. Transl. 24, 79 (1957)] unlike the Navier-Stokes equations. It is found in the Burgers equations that the potential part of velocity becomes large in comparison with the solenoidal part which decays more quickly. The probability distribution of the nonlocal term -{u}\\cdot nabla p, which spoils the maximum principle, in the local energy budget is studied in detail. It is basically symmetric, i.e., it can be either positive or negative with fluctuations. Its joint probability density functions with 1/2|{u}|^2 and with 1/2|{ω }|^2 are also found to be symmetric, fluctuating at the same times as the probability density function of -{u}\\cdot nabla p. A power-law relationship is found in the mathematical bound for the enstrophy growth dfrac{dQ}{dt} + 2 ν P ∝ left(Q^a P^bright)^α , where Q and P denote the enstrophy and the palinstrophy, respectively, and the exponents a and b are determined by calculus inequalities. We propose to quantify nonlinearity depletion by the exponent α on this basis.
NASA Technical Reports Server (NTRS)
Bonhaus, Daryl L.; Wornom, Stephen F.
1991-01-01
Two codes which solve the 3-D Thin Layer Navier-Stokes (TLNS) equations are used to compute the steady state flow for two test cases representing typical finite wings at transonic conditions. Several grids of C-O topology and varying point densities are used to determine the effects of grid refinement. After a description of each code and test case, standards for determining code efficiency and accuracy are defined and applied to determine the relative performance of the two codes in predicting turbulent transonic wing flows. Comparisons of computed surface pressure distributions with experimental data are made.
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.
NASA Technical Reports Server (NTRS)
Iannetti, Anthony C.; Moder, Jeffery P.
2010-01-01
Developing physics-based tools to aid in reducing harmful combustion emissions, like Nitrogen Oxides (NOx), Carbon Monoxide (CO), Unburnt Hydrocarbons (UHC s), and Sulfur Dioxides (SOx), is an important goal of aeronautics research at NASA. As part of that effort, NASA Glenn Research Center is performing a detailed assessment and validation of an in-house combustion CFD code known as the National Combustion Code (NCC) for turbulent reacting flows. To assess the current capabilities of NCC for simulating turbulent reacting flows with liquid jet fuel injection, a set of Single Swirler Lean Direct Injection (LDI) experiments performed at the University of Cincinnati was chosen as an initial validation data set. This Jet-A/air combustion experiment operates at a lean equivalence ratio of 0.75 at atmospheric pressure and has a 4 percent static pressure drop across the swirler. Detailed comparisons of NCC predictions for gas temperature and gaseous emissions (CO and NOx) against this experiment are considered in a previous work. The current paper is focused on detailed comparisons of the spray characteristics (radial profiles of drop size distribution and at several radial rakes) from NCC simulations against the experimental data. Comparisons against experimental data show that the use of the correlation for primary spray break-up implemented by Raju in the NCC produces most realistic results, but this result needs to be improved. Given the single or ten step chemical kinetics models, use of a spray size correlation gives similar, acceptable results
NASA Technical Reports Server (NTRS)
Shyam, Vikram
2010-01-01
A preprocessor for the Computational Fluid Dynamics (CFD) code TURBO has been developed and tested. The preprocessor converts grids produced by GridPro (Program Development Company (PDC)) into a format readable by TURBO and generates the necessary input files associated with the grid. The preprocessor also generates information that enables the user to decide how to allocate the computational load in a multiple block per processor scenario.
Reynolds averaged theory of turbulent shear flows over undulating beds and formation of sand waves.
Bose, Sujit K; Dey, Subhasish
2009-09-01
Based on the Reynolds averaged Navier-Stokes (RANS) equations and the time-averaged continuity equation, a theory of turbulent shear flow over an undulating sand bed is developed addressing the instability criterion of plane sand beds in free-surface flows leading to the formation of sand waves. In the analysis, the integration of RANS equations leads to generalized Saint Venant equations, in which the time-averaged streamwise velocity is characterized by a power law obtained from turbulence closure, treating the curvilinear streamlines by the Boussinesq approximation. As a consequence, the modified pressure distribution has a departure from the traditionally linear hydrostatic pressure distribution. The instability analysis of a plane sand bed yields the curves of the Froude number versus nondimensional wave number, determining an instability zone for which at lower Froude numbers (less than 0.8), the plane bed becomes unstable with the formation of dunes; whereas at higher Froude numbers, the plane bed becomes unstable with the formation of standing waves and antidunes. For higher Froude numbers, the experimental data for antidunes lie within the unstable zone; while for lower Froude numbers, the same is found for dunes with some experimental scatter. PMID:19905209
Stable boundary conditions and difference schemes for Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Dutt, P.
1985-01-01
The Navier-Stokes equations can be viewed as an incompletely elliptic perturbation of the Euler equations. By using the entropy function for the Euler equations as a measure of energy for the Navier-Stokes equations, it was possible to obtain nonlinear energy estimates for the mixed initial boundary value problem. These estimates are used to derive boundary conditions which guarantee L2 boundedness even when the Reynolds number tends to infinity. Finally, a new difference scheme for modelling the Navier-Stokes equations in multidimensions for which it is possible to obtain discrete energy estimates exactly analogous to those we obtained for the differential equation was proposed.
Mathematical analysis of the Navier-Stokes equations with non standard boundary conditions
NASA Technical Reports Server (NTRS)
Tidriri, M. D.
1995-01-01
One of the major applications of the domain decomposition time marching algorithm is the coupling of the Navier-Stokes systems with Boltzmann equations in order to compute transitional flows. Another important application is the coupling of a global Navier-Stokes problem with a local one in order to use different modelizations and/or discretizations. Both of these applications involve a global Navier-Stokes system with nonstandard boundary conditions. The purpose of this work is to prove, using the classical Leray-Schauder theory, that these boundary conditions are admissible and lead to a well posed problem.
Incompressible Navier-Stokes calculations in pump flows
NASA Astrophysics Data System (ADS)
Kiris, Cetin; Chang, Leon; Kwak, Dochan
1993-07-01
Flow through pump components, such as the SSME-HPFTP Impeller and an advanced rocket pump impeller, is efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. The resulting computer code is applied to the flow analysis inside an 11-inch SSME High Pressure Fuel Turbopump impeller, and an advanced rocket pump impeller. Numerical results of SSME-HPFTP impeller flow are compared with experimental measurements. In the advanced pump impeller, the effects of exit and shroud cavities are investigated. Flow analyses at design conditions will be presented.
Incompressible Navier-Stokes Calculations in Pump Flows
NASA Technical Reports Server (NTRS)
Kiris, Cetin; Chang, Leon; Kwak, Dochan
1993-01-01
Flow through pump components, such as the SSME-HPFTP Impeller and an advanced rocket pump impeller, is efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. The resulting computer code is applied to the flow analysis inside an 11-inch SSME High Pressure Fuel Turbopump impeller, and an advanced rocket pump impeller. Numerical results of SSME-HPFTP impeller flow are compared with experimental measurements. In the advanced pump impeller, the effects of exit and shroud cavities are investigated. Flow analyses at design conditions will be presented.
A Data Parallel Multizone Navier-Stokes Code
NASA Technical Reports Server (NTRS)
Jespersen, Dennis C.; Levit, Creon; Kwak, Dochan (Technical Monitor)
1995-01-01
We have developed a data parallel multizone compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the "chimera" approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. The design choices can be summarized as: 1. finite differences on structured grids; 2. implicit time-stepping with either distributed solves or data motion and local solves; 3. sequential stepping through multiple zones with interzone data transfer via a distributed data structure. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran (HPF). One interesting feature is the issue of turbulence modeling, where the architecture of a parallel machine makes the use of an algebraic turbulence model awkward, whereas models based on transport equations are more natural. We will present some performance figures for the code on the CM-5, and consider the issues involved in transitioning the code to HPF for portability to other parallel platforms.
Navier-Stokes computations of aft end flow fields
NASA Astrophysics Data System (ADS)
Weinberg, B. C.; McDonald, H.; Shamroth, S. J.
1982-05-01
A Navier-Stokes code to solve the aft end flow field of missile type configurations is presented. The consistently split linearized block implicit method of McDonald and Briley is employed in modified form to handle L-shaped domains with sharp reentrant corners. Appropriate boundary conditions are applied for the supersonic flow in particular at the outer boundary so that waves generated within the flow field are allowed to pass out of the computational domain without reflecting back into it. An adaptive grid option has been incorporated into the code and has been exercised by following the shear layer in a model backstep problem. Results are presented for the supersonic turbulent flow over a nozzle boattail configuration with and without jet exhaust and the results are compared with experiment. Calculations of the 2-D turbulent supersonic flow over a right angle back step with shear layer reattachment on a 20 deg ramp are also shown, and compared with experiments. The computation shows the qualitative physical behavior of the flows and there is generally good agreement with the experimental velocity profiles through most of the free shear layer and the ramp reattachment zone.
Design efficiency evaluation for transonic airfoil optimization - A case for Navier-Stokes design
NASA Technical Reports Server (NTRS)
Hager, J. O.; Eyi, S.; Lee, K. D.
1993-01-01
A constrained-optimization design method which improves the aerodynamic performance of transonic airfoils is evaluated from a design-quality and design-efficiency viewpoint. Design efficiency is a measure of the performance improvement and the design time (CPU time). Total-airfoil design and upper-surface design are performed using the Euler and Navier-Stokes equations with several grids, and are evaluated using the Navier-Stokes equations to determine the anticipated physical design response. Even though the cost of the Euler design is lower than Navier-Stokes design, the Navier-Stokes evaluation indicates that the Euler design does not necessarily improve the aerodynamic performance. Therefore, the design optimization should be based on an accurate flow simulation to achieve an actual performance improvement, and the design time is a secondary concern.
Commercial turbofan engine exhaust nozzle flow analyses using PAB3D
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.; Uenishi, K.; Carlson, John R.; Keith, B. D.
1992-01-01
Recent developments of a three-dimensional (PAB3D) code have paved the way for a computational investigation of complex aircraft aerodynamic components. The PAB3D code was developed for solving the simplified Reynolds Averaged Navier-Stokes equations in a three-dimensional multiblock/multizone structured mesh domain. The present analysis was applied to commercial turbofan exhaust flow systems. Solution sensitivity to grid density is presented. Laminar flow solutions were developed for all grids and two-equation k-epsilon solutions were developed for selected grids. Static pressure distributions, mass flow and thrust quantities were calculated for on-design engine operating conditions. Good agreement between predicted surface static pressures and experimental data was observed at different locations. Mass flow was predicted within 0.2 percent of experimental data. Thrust forces were typically within 0.4 percent of experimental data.
CFL3D Contribution to the AIAA Supersonic Shock Boundary Layer Interaction Workshop
NASA Technical Reports Server (NTRS)
Rumsey, Christopher L.
2010-01-01
This paper documents the CFL3D contribution to the AIAA Supersonic Shock Boundary Layer Interaction Workshop, held in Orlando, Florida in January 2010. CFL3D is a Reynolds-averaged Navier-Stokes code. Four shock boundary layer interaction cases are computed using a one-equation turbulence model widely used for other aerodynamic problems of interest. Two of the cases have experimental data available at the workshop, and two of the cases do not. The effect of grid, flux scheme, and thin-layer approximation are investigated. Comparisons are made to the available experimental data. All four cases exhibit strong three-dimensional behavior in and near the interaction regions, resulting from influences of the tunnel side-walls.
Improved Convergence and Robustness of USM3D Solutions on Mixed Element Grids (Invited)
NASA Technical Reports Server (NTRS)
Pandya, Mohagna J.; Diskin, Boris; Thomas, James L.; Frink, Neal T.
2015-01-01
Several improvements to the mixed-element USM3D discretization and defect-correction schemes have been made. A new methodology for nonlinear iterations, called the Hierarchical Adaptive Nonlinear Iteration Scheme (HANIS), has been developed and implemented. It provides two additional hierarchies around a simple and approximate preconditioner of USM3D. The hierarchies are a matrix-free linear solver for the exact linearization of Reynolds-averaged Navier Stokes (RANS) equations and a nonlinear control of the solution update. Two variants of the new methodology are assessed on four benchmark cases, namely, a zero-pressure gradient flat plate, a bump-in-channel configuration, the NACA 0012 airfoil, and a NASA Common Research Model configuration. The new methodology provides a convergence acceleration factor of 1.4 to 13 over the baseline solver technology.
NASA Technical Reports Server (NTRS)
Liu, C. H.; Wong, T. C.; Kandil, O. A.
1988-01-01
The two-dimensional flow over a blunt leading-edge plate is simulated on the basis of an Euler/Navier-Stokes zonal scheme. The scheme uses an implicit upwind finite-volume scheme, which is based on the van Leer flux-vector splitting. It is shown that the Euler/Navier-Stokes zonal scheme with downstream boundary-layer compatibility conditions is accurate and efficient.
Multi-zonal Navier-Stokes code with the LU-SGS scheme
NASA Technical Reports Server (NTRS)
Klopfer, G. H.; Yoon, S.
1993-01-01
The LU-SGS (lower upper symmetric Gauss Seidel) algorithm has been implemented into the Compressible Navier-Stokes, Finite Volume (CNSFV) code and validated with a multizonal Navier-Stokes simulation of a transonic turbulent flow around an Onera M6 transport wing. The convergence rate and robustness of the code have been improved and the computational cost has been reduced by at least a factor of 2 over the diagonal Beam-Warming scheme.
Single-grid spectral collocation for the Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Bernardi, Christine; Canuto, Claudio; Maday, Yvon; Metivet, Brigitte
1988-01-01
The aim of the paper is to study a collocation spectral method to approximate the Navier-Stokes equations: only one grid is used, which is built from the nodes of a Gauss-Lobatto quadrature formula, either of Legendre or of Chebyshev type. The convergence is proven for the Stokes problem provided with inhomogeneous Dirichlet conditions, then thoroughly analyzed for the Navier-Stokes equations. The practical implementation algorithm is presented, together with numerical results.
Navier-Stokes and potential theory solutions for ahelicopter fuselage and comparison with experiment
NASA Technical Reports Server (NTRS)
Chaffin, Mark S.; Berry, John D.
1994-01-01
A thin-layer Navier-Stokes code and a panel method code are used to predict the flow over a generic helicopter fuselage. The computational results are compared with pressure data at four experimental conditions. Both methods produce results that agree with the experimental pressure data. However, separation patterns and other viscous flow features from the Navier-Stokes code solution are shown that cannot be easily modeled with the panel method.
Navier-Stokes solution for steady two-dimensional transonic cascade flows
Kwon, O.K.
1987-01-01
A robust, time-marching Navier-Stokes solution procedure based on the explicit hopscotch method is presented for solution of steady, two-dimensional, transonic turbine cascade flows. The method is applied to the strong conservation form of the unsteady Navier-Stokes equations written in arbitrary curvilinear coordinates. Cascade flow solutions are obtained on an orthogonal, body-conforming ''O'' grid with the standard k-epsilon turbulence model. Computed results are presented and compared with experimental data.
Adapting a Navier-Stokes code to the ICL-DAP
NASA Technical Reports Server (NTRS)
Grosch, C. E.
1985-01-01
The results of an experiment are reported, i.c., to adapt a Navier-Stokes code, originally developed on a serial computer, to concurrent processing on the CL Distributed Array Processor (DAP). The algorithm used in solving the Navier-Stokes equations is briefly described. The architecture of the DAP and DAP FORTRAN are also described. The modifications of the algorithm so as to fit the DAP are given and discussed. Finally, performance results are given and conclusions are drawn.