Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “total variation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

Computational Aerothermodynamic Simulation Issues on Unstructured Grids

The synthesis of physical models for gas chemistry and turbulence from the structured grid codes LAURA and VULCAN into the unstructured grid code FUN3D is described. A directionally Symmetric, Total Variation Diminishing (STVD) algorithm and an entropy fix (eigenvalue limiter) keyed to local cell Reynolds number are introduced to improve solution quality for hypersonic aeroheating applications. A simple grid-adaptation procedure is incorporated within the flow solver. Simulations of flow over an ellipsoid (perfect gas, inviscid), Shuttle Orbiter (viscous, chemical nonequilibrium) and comparisons to the structured grid solvers LAURA (cylinder, Shuttle Orbiter) and VULCAN (flat plate) are presented to show current capabilities. The quality of heating in 3D stagnation regions is very sensitive to algorithm options in general, high aspect ratio tetrahedral elements complicate the simulation of high Reynolds number, viscous flow as compared to locally structured meshes aligned with the flow.

Gnoffo, Peter A.↗

Computations of Flow over a Hump Model Using Higher Order Method with Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, P.↗

Transition in a Supersonic Boundary Layer Due to Acoustic Disturbances

The boundary layer receptivity process due to the interaction of three-dimensional slow and fast acoustic disturbances with a blunted flat plate is numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 10(exp 6)/inch. The computations are performed with and without two-dimensional isolated roughness element located near the leading edge. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The simulations showed that the linear instability waves are generated very close to the leading edge. The wavelength of the disturbances inside the boundary layer first increases gradually and becomes longer than the wavelength for the instability waves within a short distance from the leading edge. The wavelength then decreases gradually and merges with the wavelength for the Tollmien-Schlichting wave. The initial amplitudes of the instability waves near the neutral points, the receptivity coefficients, are about 1.20 and 0.07 times the amplitude of the free-stream disturbances for the slow and the fast waves respectively. It was also revealed that small isolated roughness element does not enhance the receptivity process for the given nose bluntness.

Balakumar, P.↗

Computations of Flow Over a Hump Model Using Higher Order Method With Turbulence Modeling

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k-omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as they are observed in the experiments. The computed separation regions are much longer than that are observed. However, the percentage reduction in the separation region in the steady suction case is closer to that was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than that are measured and it points towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Balakumar, Ponnampalam↗

Transition in a Supersonic Boundary Layer due to Acoustic Disturbances

The boundary layer receptivity process due to the interaction of three-dimensional slow and fast acoustic disturbances with a blunted flat plate is numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 106/inch. The computations are performed with and without two-dimensional isolated roughness element located near the leading edge. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The simulations showed that the linear instability waves are generated very close to the leading edge. The wavelength of the disturbances inside the boundary layer first increases gradually and becomes longer than the wavelength for the instability waves within a short distance from the leading edge. The wavelength then decreases gradually and merges with the wavelength for the Tollmien_Schlichting wave. The initial amplitudes of the instability waves near the neutral points, the receptivity coefficients, are about 1.20 and 0.07 times the amplitude of the free-stream disturbances for the slow and the fast waves respectively. It was also revealed that small isolated roughness element does not enhance the receptivity process for the given nose bluntness.

Balakumar, Ponnampalam↗

Stability of Supersonic Boundary Layers Over Blunt Wedges

Receptivity and stability of supersonic boundary layers over blunt flat plates and wedges are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. Computations are performed for a flat plate with leading edge thicknesses of 0.0001, 0.001, 0.005 and 0.01 inches that give Reynolds numbers based on the leading edge thickness ranging from 1000 to 10000. Calculations are also performed for a wedge of 10 degrees half angle with different leading edge radii 0.001 and 0.01 inches. The linear stability results showed that the bluntness has a strong stabilizing effect on the stability of two-dimensional boundary layers. The transition Reynolds number for a flat plate with a leading edge thickness of 0.01 inches is about 3.5 times larger than it is for the Blasius boundary layer. It was also revealed that boundary layers on blunt wedges are far more stable than on blunt flat plates.

Balakumar, Ponnampalam↗

Receptivity of Hypersonic Boundary Layers Due to Acoustic Disturbances over Blunt Cone

The transition process induced by the interaction of acoustic disturbances in the free-stream with boundary layers over a 5-degree straight cone and a wedge with blunt tips is numerically investigated at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves the Navier-Stokes equations are solved in axisymmetric coordinates. The governing equations are solved using the 5th -order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, acoustic disturbances are introduced at the outer boundary of the computational domain and unsteady simulations are performed. Generation and evolution of instability waves and the receptivity of boundary layer to slow and fast acoustic waves are investigated. The mean flow data are compared with the experimental results. The results show that the instability waves are generated near the leading edge and the non-parallel effects are stronger near the nose region for the flow over the cone than that over a wedge. It is also found that the boundary layer is much more receptive to slow acoustic wave (by almost a factor of 67) as compared to the fast wave.

Kara, K.↗

Simulation of Stagnation Region Heating in Hypersonic Flow on Tetrahedral Grids

Hypersonic flow simulations using the node based, unstructured grid code FUN3D are presented. Applications include simple (cylinder) and complex (towed ballute) configurations. Emphasis throughout is on computation of stagnation region heating in hypersonic flow on tetrahedral grids. Hypersonic flow over a cylinder provides a simple test problem for exposing any flaws in a simulation algorithm with regard to its ability to compute accurate heating on such grids. Such flaws predominantly derive from the quality of the captured shock. The importance of pure tetrahedral formulations are discussed. Algorithm adjustments for the baseline Roe / Symmetric, Total-Variation-Diminishing (STVD) formulation to deal with simulation accuracy are presented. Formulations of surface normal gradients to compute heating and diffusion to the surface as needed for a radiative equilibrium wall boundary condition and finite catalytic wall boundary in the node-based unstructured environment are developed. A satisfactory resolution of the heating problem on tetrahedral grids is not realized here; however, a definition of a test problem, and discussion of observed algorithm behaviors to date are presented in order to promote further research on this important problem.

Gnoffo, Peter A.↗

Effects of Nose Bluntness on Stability of Hypersonic Boundary Layers over Blunt Cone

Receptivity and stability of hypersonic boundary layers are numerically investigated for boundary layer flows over a 5-degree straight cone at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves, we solve the Navier-Stokes equations in axisymmetric coordinates. The governing equations are solved using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. Generation of instability waves from leading edge region and receptivity of boundary layer to slow acoustic waves are investigated. Computations are performed for a cone with nose radii of 0.001, 0.05 and 0.10 inches that give Reynolds numbers based on the nose radii ranging from 650 to 130,000. The linear stability results showed that the bluntness has a strong stabilizing effect on the stability of axisymmetric boundary layers. The transition Reynolds number for a cone with the nose Reynolds number of 65,000 is increased by a factor of 1.82 compared to that for a sharp cone. The receptivity coefficient for a sharp cone is about 4.23 and it is very small, approx.10(exp -3), for large bluntness.

Kara, K.↗

Computations of Flow Over the Hump Model Using Higher-Order Method With Turbulence Modeling

The flow over the two-dimensional hump model is computed by solving the RANS equations with kappa-omega (SST) model. The governing equations, the flow equations and the turbulent equations, are solved using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using explicit third order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The WENO and the TVD methods and the formulas are explained in [1] and the application of ENO method to N-S equations is given in [2]. The solution method implemented in this computation is described in detail in [3].

Balakumar, P.↗

Receptivity of Supersonic Boundary Layers Due To Acoustic Disturbances Over Blunt Cones

Receptivity and stability of supersonic boundary layers over a 5-degree straight cone with a blunt tip are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 106/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The linear stability results showed that bluntness has less stabilizing effects on the stability of boundary layers over cones than on flat plates and wedges. The unsteady simulations of the interaction of plane threedimensional acoustic waves with the cone showed that the modulation of wavelength and the generation of instability waves first occurred near the leading edge in the plane where the constant acoustic phase lines are perpendicular to the cone axis. Further downstream, this instability region spreads in the azimuthal direction from this plane.

Balakumar, P.↗

High-Order Energy Stable WENO Schemes

A new third-order Energy Stable Weighted Essentially NonOscillatory (ESWENO) finite difference scheme for scalar and vector linear hyperbolic equations with piecewise continuous initial conditions is developed. The new scheme is proven to be stable in the energy norm for both continuous and discontinuous solutions. In contrast to the existing high-resolution shock-capturing schemes, no assumption that the reconstruction should be total variation bounded (TVB) is explicitly required to prove stability of the new scheme. A rigorous truncation error analysis is presented showing that the accuracy of the 3rd-order ESWENO scheme is drastically improved if the tuning parameters of the weight functions satisfy certain criteria. Numerical results show that the new ESWENO scheme is stable and significantly outperforms the conventional third-order WENO finite difference scheme of Jiang and Shu in terms of accuracy, while providing essentially nonoscillatory solutions near strong discontinuities.

Yamaleev, Nail K.↗

Boundary Layer Receptivity due to Roughness and Freestream Sound for Supersonic Flows over Axisymmetric Cones

The receptivity of supersonic boundary layers due to the interaction of a two-dimensional acoustic wave with a three-dimensional roughness is numerically investigated over a 5-degree sharp tipped straight cone at a free stream Mach number of 3.5 and a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. The simulations are performed with different roughness heights and roughness locations. The steady computations with the roughness reveal that the perturbations induced by the roughness remain very close to the wall and do not penetrate to the outer part of the boundary layer. The effects of the roughness persist for about 10 boundary layer thicknesses downstream of each element. The unsteady simulation shows the generation of unstable three-dimensional first mode disturbances within a short distance from the roughness. However, the amplitudes of the instability waves are very small near the neutral point, on the order of 10(exp -4) compared to the free stream disturbances. The receptivity coefficients due to the roughness elements located very close to the tip are one order of magnitude smaller than those due to the roughness elements located closer to the neutral point.

Balakumar, P.↗

Effects of Wall Cooling on Hypersonic Boundary Layer Receptivity Over a Cone

Effects of wall cooling on the receptivity process induced by the interaction of slow acoustic disturbances in the free-stream are numerically investigated for a boundary layer flow over a 5-degrees straight cone. The free-stream Mach number is 6.0 and the Reynolds number is 7.8x10(exp 6)/ft. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using 3rd-order total variation diminishing (T VD) Runge-K utta scheme for time integration. Computations are performed for a cone with nose radius of 0.001 inch for adiabatic wall temperature (T(sub aw)), 0.75*T(sub aw), 0.5*T(sub aw), 0.40*T(sub aw), 0.30*T(sub aw), and 0.20*T(sub aw). Once the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. Generation of instability waves from leading edge region and receptivity of boundary layer to slow acoustic waves are investigated. Computations showed that wall cooling has strong stabilization effect on the first mode disturbances as was observed in the experiments. T ransition location moved to upstream when wall cooling was applied It is also found that the boundary layer is much more receptive to fast acoustic wave (by almost a factor of 50). When simulations performed using the same forcing frequency growth of the second mode disturbances are delayed with wall cooling and they attained values two times higher than that of adiabatic case. In 0.20*T(sub aw) case the transition Reynolds number is doubled compared to adiabatic conditions. The receptivity coefficient for adiabatic wall case (804 R) is 1.5225 and for highly cooled cones (241, and 161 R); they are in the order of 10(exp -3).

Kara, K.↗

Stability of Supersonic Boundary Layers on a Cone at an Angle of Attack

The stability and receptivity of three-dimensional supersonic boundary layers over a 7deg sharp tipped straight cone at an angle of attack of 4.2deg is numerically investigated at a free stream Mach number of 3.5 and at two high Reynolds numbers, 0.25 and 0.50x10(exp 6)/inch. The generation and evolution of stationary crossflow vortices are also investigated by performing simulations with three-dimensional roughness elements located on the surface of the cone. The flow fields with and without the roughness elements are obtained by solving the full Navier-Stokes equations in cylindrical coordinates using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using the third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. Stability computations reveal that the azimuthal wavenumbers are in the range of m approx. 25-50 for the most amplified traveling disturbances and in the range of m approx. 40-70 for the stationary disturbances. The N-Factor computations predicted that transition would occur further forward in the middle of the cone compared to the transition fronts near the windward and the leeward planes. The simulations revealed that the crossflow vortices originating from the nose region propagate towards the leeward plane. No perturbations were observed in the lower part of the cone.

Balakumar, Ponnampalam↗

Receptivity of Hypersonic Boundary Layers over Straight and Flared Cones

The effects of adverse pressure gradients on the receptivity and stability of hypersonic boundary layers were numerically investigated. Simulations were performed for boundary layer flows over a straight cone and two flared cones. The steady and the unsteady flow fields were obtained by solving the two-dimensional Navier-Stokes equations in axi-symmetric coordinates using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The mean boundary layer profiles were analyzed using local stability and non-local parabolized stability equations (PSE) methods. After the most amplified disturbances were identified, two-dimensional plane acoustic waves were introduced at the outer boundary of the computational domain and time accurate simulations were performed. The adverse pressure gradient was found to affect the boundary layer stability in two important ways. Firstly, the frequency of the most amplified second-mode disturbance was increased relative to the zero pressure gradient case. Secondly, the amplification of first- and second-mode disturbances was increased. Although an adverse pressure gradient enhances instability wave growth rates, small nose-tip bluntness was found to delay transition due to the low receptivity coefficient and the resulting weak initial amplitude of the instability waves. The computed and measured amplitude-frequency spectrums in all three cases agree very well in terms of frequency and the shape except for the amplitude.

Balakumar, Ponnampalam↗

Stability of Hypersonic Boundary Layers on a Cone at an Angle of Attack

The stability and receptivity of a three-dimensional hypersonic boundary layer over a 7deg half-angle straight cone at an angle of attack of 6deg is numerically investigated at a freestream Mach number of 6.0 and a Reynolds number of 10.4x10(exp 6)/m. The generation and evolution of stationary crossflow vortices are also investigated by performing simulations with three-dimensional roughness elements located on the surface of the cone. The flow fields with and without the roughness elements are obtained by solving the full Navier- Stokes equations in cylindrical coordinates using a fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and a third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. Stability computations produced azimuthal wavenumbers in the range of m approx. 20-50 for the most amplified traveling disturbances and in the range of m approx.30-70 for the stationary disturbances. The frequency of the unstable second-mode ranges from 400 kHz to 900 kHz along the windward ray. The N-Factor computations predicted transition would occur more forward on the sides of the cone as compared to the transition fronts near the windward and the leeward rays. The simulations also show the crossflow vortices originating from the nose region propagate towards the leeward ray. No perturbations were observed toward the windward half of the cone.

Balakumar, P.↗

Convective Weather Forecast Accuracy Analysis at Center and Sector Levels

This paper presents a detailed convective forecast accuracy analysis at center and sector levels. The study is aimed to provide more meaningful forecast verification measures to aviation community, as well as to obtain useful information leading to the improvements in the weather translation capacity models. In general, the vast majority of forecast verification efforts over past decades have been on the calculation of traditional standard verification measure scores over forecast and observation data analyses onto grids. These verification measures based on the binary classification have been applied in quality assurance of weather forecast products at the national level for many years. Our research focuses on the forecast at the center and sector levels. We calculate the standard forecast verification measure scores for en-route air traffic centers and sectors first, followed by conducting the forecast validation analysis and related verification measures for weather intensities and locations at centers and sectors levels. An approach to improve the prediction of sector weather coverage by multiple sector forecasts is then developed. The weather severe intensity assessment was carried out by using the correlations between forecast and actual weather observation airspace coverage. The weather forecast accuracy on horizontal location was assessed by examining the forecast errors. The improvement in prediction of weather coverage was determined by the correlation between actual sector weather coverage and prediction. observed and forecasted Convective Weather Avoidance Model (CWAM) data collected from June to September in 2007. CWAM zero-minute forecast data with aircraft avoidance probability of 60% and 80% are used as the actual weather observation. All forecast measurements are based on 30-minute, 60- minute, 90-minute, and 120-minute forecasts with the same avoidance probabilities. The forecast accuracy analysis for times under one-hour showed that the errors in intensity and location for center forecast are relatively low. For example, 1-hour forecast intensity and horizontal location errors for ZDC center were about 0.12 and 0.13. However, the correlation between sector 1-hour forecast and actual weather coverage was weak, for sector ZDC32, about 32% of the total variation of observation weather intensity was unexplained by forecast; the sector horizontal location error was about 0.10. The paper also introduces an approach to estimate the sector three-dimensional actual weather coverage by using multiple sector forecasts, which turned out to produce better predictions. Using Multiple Linear Regression (MLR) model for this approach, the correlations between actual observation and the multiple sector forecast model prediction improved by several percents at 95% confidence level in comparison with single sector forecast.

Wang, Yao↗