Search NASA⌕ Search

SEARCH · Search NASA

Results for “K-W”

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.

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Implementation, Realization and an Effective Solver of Two-Equation Turbulence Models

Currently, when the Reynolds-Averaged Navier-Stokes (RANS) equations are solved using turbulence modeling, most often the one-equation model of Spalart and Allmaras is used. Then, it is only necessary to solve the RANS equations in conjunction with a single transport equation for modeling turbulence. For this model, considerable assessment and analysis has been performed, allowing the possibility of a reliable solution method for an eddy viscosity required to compute the Reynolds stresses in the RANS equations. Such evaluation along with analysis has not been achieved to realize similar performance with two-equation models of the k-w type. The primary objective of this paper is to present and discuss the components of an effective numerical algorithm for solving the RANS equations and the two transport equations of k-w type turbulence models. All the important details of the turbulence model as actually implemented are given, which is sometimes not done in various papers considering such modeling. The viability and effectiveness of this solution algorithm are demonstrated by solving both two-dimensional and three-dimensional aerodynamic flows. In all applications, a linear rate of convergence without oscillations or other evidence of unstable behavior is observed. This behavior is also particularly true when the proposed algorithm is applied to systematically renewed mesh sequences, which is generally not observed with algorithms solving more than one transport equation. Thus, numerical integration errors are systematically reduced, allowing for a significantly more reliable assessment of the effectiveness of the turbulence model. Additionally, in this paper, analysis of the solution algorithm, including linear stability, is also performed for a particular flow problem.

Turbulence↗
DOE OSTI2021

Review and Analysis of Heat Transfer Correlations For Horizontal Pseudocritical CO 2 Heat Exchanger Applications

There is significant interest in the development of supercritical carbon dioxide (sCO 2 ) power cycles because of the potential for smaller and more energy efficient systems than a steam Rankine cycle. Heat exchanger designs typically use empirical correlations, but the applicability of these correlations near the CO 2 critical point is a potential issue. Though numerous correlations have been proposed in the literature, there are some disagreements when it comes to the accuracy. The current work recognizes the role of thermophysical properties, and its impacts on the heat transfer correlations and cycle efficiency. Heat transfer correlations proposed for horizontal flow inside circular pipes were analyzed with the help of numerical simulations. Steady state RANS simulations were performed using SST k-w turbulence model to evaluate the Nusselt number empirical correlations. It was found that the most of correlations (except Yoon) produced a Nusselt number that differed significantly with the one predicted numerically. Some of the correlations were developed for pure forced convection regime and as mentioned in Lin et al. do not account for mixed convection or free convection effects. Based on the limited observation, it appears that Yoon et al. predictions match well with the numerically predicted Nusselt Numbers. However, further analysis is required understand the applicability of various correlations and is contingent on accurate measurements or predictions of wall temperature profiles in the axial and the circumferential directions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1982

Operation of the J-series thruster using inert gas

Electron bombardment ion thrusters using inert gases are candidates for large space systems. The J-Series 30 cm diameter thruster, designed for operation up to 3 k-W with mercury, is at a state of technology readiness. The characteristics of operation with xenon, krypton, and argon propellants in a J-Series thruster with that obtained with mercury are compared. The performance of the discharge chamber, ion optics, and neutralizer and the overall efficiency as functions of input power and specific impulse and thruster lifetime were evaluated. As expected, the discharge chamber performance with inert gases decreased with decreasing atomic mass. Aspects of the J-Series thruster design which would require modification to provide operation at high power with insert gases were identified.

Rawlin, V. K.↗
NASA NTRS1990

One- and two-dimensional Stirling machine simulation using experimentally generated flow turbulence models

Investigations of one- and two-dimensional (1- or 2-D) simulations of Stirling machines centered around experimental data generated by the U. of Minnesota Mechanical Engineering Test Rig (METR) are covered. This rig was used to investigate oscillating flows about a zero mean with emphasis on laminar/turbulent flow transitions in tubes. The Space Power Demonstrator Engine (SPDE) and in particular, its heater, were the subjects of the simulations. The heater was treated as a 1- or 2-D entity in an otherwise 1-D system. The 2-D flow effects impacted the transient flow predictions in the heater itself but did not have a major impact on overall system performance. Information propagation effects may be a significant issue in the simulation (if not the performance) of high-frequency, high-pressure Stirling machines. This was investigated further by comparing a simulation against an experimentally validated analytic solution for the fluid dynamics of a transmission line. The applicability of the pressure-linking algorithm for compressible flows may be limited by characteristic number (defined as flow path information traverses per cycle); this warrants further study. Lastly the METR was simulated in 1- and 2-D. A two-parameter k-w foldback function turbulence model was developed and tested against a limited set of METR experimental data.

Goldberg, Louis F.↗
NASA NTRS1995

The applicability of turbulence models to aerodynamic and propulsion flowfields at McDonnell-Douglas Aerospace

The objective of this viewgraph presentation is to evaluate turbulence models for integrated aircraft components such as the forebody, wing, inlet, diffuser, nozzle, and afterbody. The one-equation models have replaced the algebraic models as the baseline turbulence models. The Spalart-Allmaras one-equation model consistently performs better than the Baldwin-Barth model, particularly in the log-layer and free shear layers. Also, the Sparlart-Allmaras model is not grid dependent like the Baldwin-Barth model. No general turbulence model exists for all engineering applications. The Spalart-Allmaras one-equation model and the Chien k-epsilon models are the preferred turbulence models. Although the two-equation models often better predict the flow field, they may take from two to five times the CPU time. Future directions are in further benchmarking the Menter blended k-w/k-epsilon and algorithmic improvements to reduce CPU time of the two-equation model.

Kral, Linda D.↗
NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗
NASA NTRS2001

Two-Equation Turbulence Models for Prediction of Heat Transfer on a Transonic Turbine Blade

Two versions of the two-equation k-omega model and a shear stress transport (SST) model are used in a three-dimensional, multi-block, Navier-Stokes code to compare the detailed heat transfer measurements on a transonic turbine blade. It is found that the SST model resolves the passage vortex better on the suction side of the blade, thus yielding a better comparison with the experimental data than either of the k-w models. However, the comparison is still deficient on the suction side of the blade. Use of the SST model does require the computation of distance from a wall, which for a multiblock grid, such as in the present case, can be complicated. However, a relatively easy fix for this problem was devised. Also addressed are issues such as (1) computation of the production term in the turbulence equations for aerodynamic applications, and (2) the relation between the computational and experimental values for the turbulence length scale, and its influence on the passage vortex on the suction side of the turbine blade.

Garg, Vijay K.↗
NASA NTRS2009

Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation

This report gives an overview of physical modeling enhancements to the Wind-US flow solver which were made to improve the capabilities for simulation of hypersonic flows and the reliability of computations to complement hypersonic testing. The improvements include advanced turbulence models, a bypass transition model, a conjugate (or closely coupled to vehicle structure) conduction-convection heat transfer capability, and an upgraded high-speed combustion solver. A Mach 5 shock-wave boundary layer interaction problem is used to investigate the benefits of k- s and k-w based explicit algebraic stress turbulence models relative to linear two-equation models. The bypass transition model is validated using data from experiments for incompressible boundary layers and a Mach 7.9 cone flow. The conjugate heat transfer method is validated for a test case involving reacting H2-O2 rocket exhaust over cooled calorimeter panels. A dual-mode scramjet configuration is investigated using both a simplified 1-step kinetics mechanism and an 8-step mechanism. Additionally, variations in the turbulent Prandtl and Schmidt numbers are considered for this scramjet configuration.

Georgiadis, Nicholas J.↗
NASA NTRS2015

Effect of Interfacial Turbulence and Accommodation Coefficient on CFD Predictions of Pressurization and Pressure Control in Cryogenic Storage Tank

Laminar models agree closely with the pressure evolution and vapor phase temperature stratification but under-predict liquid temperatures. Turbulent SST k-w and k-e models under-predict the pressurization rate and extent of stratification in the vapor but represent liquid temperature distributions fairly well. These conclusions seem to equally apply to large cryogenic tank simulations as well as small scale simulant fluid pressurization cases. Appropriate turbulent models that represent both interfacial and bulk vapor phase turbulence with greater fidelity are needed. Application of LES models to the tank pressurization problem can serve as a starting point.

Condensation Evaporation↗
NASA NTRS2015

Uncertainty Quantification of Turbulence Model Closure Coefficients for Transonic Wall-Bounded Flows

The goal of this work was to quantify the uncertainty and sensitivity of commonly used turbulence models in Reynolds-Averaged Navier-Stokes codes due to uncertainty in the values of closure coefficients for transonic, wall-bounded flows and to rank the contribution of each coefficient to uncertainty in various output flow quantities of interest. Specifically, uncertainty quantification of turbulence model closure coefficients was performed for transonic flow over an axisymmetric bump at zero degrees angle of attack and the RAE 2822 transonic airfoil at a lift coefficient of 0.744. Three turbulence models were considered: the Spalart-Allmaras Model, Wilcox (2006) k-w Model, and the Menter Shear-Stress Trans- port Model. The FUN3D code developed by NASA Langley Research Center was used as the flow solver. The uncertainty quantification analysis employed stochastic expansions based on non-intrusive polynomial chaos as an efficient means of uncertainty propagation. Several integrated and point-quantities are considered as uncertain outputs for both CFD problems. All closure coefficients were treated as epistemic uncertain variables represented with intervals. Sobol indices were used to rank the relative contributions of each closure coefficient to the total uncertainty in the output quantities of interest. This study identified a number of closure coefficients for each turbulence model for which more information will reduce the amount of uncertainty in the output significantly for transonic, wall-bounded flows.

Schaefer, John↗
NASA NTRS2020

Acoustic Optimization for Anti-Phase Asymmetric Rotor

This investigation seeks to optimize the implementation of anti-phase alternating trailing edge (TE) patterns for rotor noise suppression. The design objective is to maximize reduction of noise perceived by the community while maintaining the aerodynamic thrust. Computations using a three-dimensional Unsteady-Reynolds-Averaged-Navier-Stokes (URANS) with k-w Shear Stress Transport (SST) turbulence model and Ffowcs-Williams and Hawkings (FW-H) formula are used to obtain aerodynamic thrust and far-field noise level. A parametric acoustic study of 13 configurations of KDE rotor with variable alternating trailing edge period, alternating trailing edge length, and trailing edge deflection angle is conducted. The best design candidate for the KDE rotor has a four-period TE waveform which results in a reduction in far-field noise level of 2.1 dB in the hover condition and a reduction of 1.1 dB in the forward flight condition at 9.7 m/s. A further parametric acoustic study is conducted for a different rotor manufactured by APC. Six APC rotor design candidates are simulated. The best design candidate 4H for the APC rotor results in a reduction in far-field noise level of 4.0 dB in the hover condition and a reduction of 1.3 dB in the forward flight condition at 9.7 m/s. A series of acoustic experiments in the Penn State University (PSU) anechoic chamber have been conducted. In the forward flight condition at 9.7 m/s, the APC anti-phase 4H rotor offers clear evidence of noise suppression capability across a wide range of the azimuthal angle. In the broadband frequency range of 2000-4000 Hz, the APC anti-phase 4H rotor produces as much as 6 dB noise reduction. The experimental results appear to confirm the noise suppression capability of the proposed anti-phase rotor design concepts.

Xiong, Juntao↗
NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗
NASA NTRSDate not supplied

On Generalizably Improving RANS Predictions of Flow Separation and Reattachment

This work presents progress in an ongoing data-driven turbulence modeling effort at NASA Langley Research Center to improve predictions for separated flows. Two augmentation strategies are presented – one to improve predictions for the point of reattachment after smooth-body separation, and another to improve predictions for flow separation due to prolonged adverse pressure gradients. The reattachment augmentation was inspired from lessons learned during attempts at using field inversion and machine learning (FIML); however, it is prescribed as a conventional analytic closed-form expression. Results show that, when introduced into the SST k-w model, it shortens the otherwise overpredicted separation bubble length for three cases characterized by different Reynolds numbers and flow geometries. The second augmentation, which addresses flow separation under prolonged adverse pressure gradients, is inferred via FIML and applied as a neural network within the Spalart-Allmaras model. It uses features that are chosen such that the presented results show an improved predictive accuracy for the training case without affecting other cases where the augmentation is not needed. Further testing is required to establish the generalizability of this augmentation. The rationale behind the chosen features that serve as inputs to both these augmentations is also presented for the benefit of the reader.

machine learning↗
NASA NTRSDate not supplied

Loci/STREAM Sharp Interface Model Validation of Cryogenic Propellant Tank Self-Pressurization

This paper documents a model validation effort for the prediction of cryogenic propellant tank selfpressurization physics using the Loci/STREAM sharp interface CFD model developed at NASA MSFC. Cryogenic fluid management (CFM) applications for long duration spaceflight missions require accurate modeling of heat transfer and thermodynamics. Therefore, detailed validation of those physics is important to ensure simulation accuracy. In the past, the MSFC Propulsion Fluid Dynamics Branch has validated and used the Loci/STREAM Volume of Fluid (VOF) module extensively for propellant slosh applications, which are dominated by fluid motion. For cases where the interface gas/liquid interface is essentially static, the Loci/STREAM sharp interface was developed for more rapid and efficient CFM analysis. In the present study, Loci/STREAM sharp interface simulations were compared with the NASA K-site self-pressurization experiment of a flight weight, partially full, LH2 tank in normal gravity. This study investigated several parametric effects and made a number of important observations on the performance of the Loci/STREAM sharp interface method: - First, a homogeneous model was derived to verify the CFD model. The homogeneous model assumes both gas and liquid phases are uniform and are at saturation temperature. The predicted pressure rise rates at 3.5 W/m2 and 2.0 W/m2 with a 49% fill level are all lower than the experimental measurement. This is expected due to the observed temperature stratification near the gas-liquid interface during testing. - The k-w and k-e turbulence models predict comparable pressurization rates, which are much lower than the experiment measurement. On the other hand, the pressure values and pressurization rates from both models are almost the same as those of the homogeneous model. The current turbulence model is too diffusive and is not capable of capturing the thermal stratification in the liquid and ullage. - When using the laminar flow model, the pressurization rate is found to be very close to the experimental value when a quasi-steady flow field is allowed to develop by continuously venting the tank prior to pressurization as was done during testing. On the other hand, without initial venting, the rate is not constant. Additionally, the development time to reach a quasi-steady flow field differed from the experiment due to the lack of turbulent mixing away from the gas-liquid interface that is expected in reality. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 15% of the experimental value for 49% fill level at 3.5 W/m 2 and is within 18% for 2.0 W/m 2 heat loading. - Using the laminar flow model and appropriate initial conditioning, the pressurization rate is within 18% of the experimental value for a 29% fill level and within 15% for an 83% fill level.

H. Q. Yang↗
NASA NTRSDate not supplied

Novel Turbulence Model Tailored for Complex Rocket Motor Exhaust Jets

This paper describes a new RANS turbulence model intended to give improved accuracy in turbulence that is affected by compressibility, rotation and curvature, and variable density, making it well-suited to complex rocket motor exhaust jets. Four key elements of the model are: the use of simple-averaged (not mass-averaged) turbulence decomposition; a second-order Taylor series for the equilibrated turbulent fluxes; empirical modeling that accounts for the combined effects of aforementioned turbulence modifiers; and model tuning on large-eddy simulation data for mixing layers for a wide range of these turbulence modifiers. The model is implemented as an extension of the k-w SST model. Results include comparisons to detached eddy simulation results for hot and cold co-flowing exhaust jets, and to experimental data for supersonic retropropulsion on the CobraMRV Mars-lander concept vehicle. The model is named TMRC for its constituents elements, the Taylor series with Mach number, rotation and curvature effects.

turbulence modeling↗
DOE OSTI2020

Materials Data on K3W by Materials Project

K3W is alpha bismuth trifluoride structured and crystallizes in the cubic Fm-3m space group. The structure is three-dimensional. there are two inequivalent K sites. In the first K site, K is bonded to four equivalent K and four equivalent W atoms to form a mixture of distorted corner, edge, and face-sharing KK4W4 tetrahedra. All K–K bond lengths are 3.91 Å. All K–W bond lengths are 3.91 Å. In the second K site, K is bonded in a body-centered cubic geometry to eight equivalent K atoms. W is bonded in a body-centered cubic geometry to eight equivalent K atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on K3W by Materials Project

K3W is alpha bismuth trifluoride structured and crystallizes in the tetragonal I4/mmm space group. The structure is three-dimensional. there are two inequivalent K sites. In the first K site, K is bonded in a 4-coordinate geometry to four equivalent K and four equivalent W atoms. All K–K bond lengths are 3.99 Å. All K–W bond lengths are 3.99 Å. In the second K site, K is bonded to eight equivalent K and four equivalent W atoms to form distorted KK8W4 cuboctahedra that share corners with four equivalent KK8W4 cuboctahedra, corners with eight equivalent WK12 cuboctahedra, edges with eight equivalent KK8W4 cuboctahedra, faces with four equivalent KK8W4 cuboctahedra, and faces with six equivalent WK12 cuboctahedra. All K–W bond lengths are 4.18 Å. W is bonded to twelve K atoms to form WK12 cuboctahedra that share corners with four equivalent WK12 cuboctahedra, corners with eight equivalent KK8W4 cuboctahedra, edges with eight equivalent WK12 cuboctahedra, faces with four equivalent WK12 cuboctahedra, and faces with six equivalent KK8W4 cuboctahedra.

36 MATERIALS SCIENCE↗