Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall modeling”

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 145 records · Page 8

Higher order turbulence closure models

Theoretical models are developed and numerical studies conducted on various types of flows including both elliptic and parabolic. The purpose of this study is to find better higher order closure models for the computations of complex flows. This report summarizes three new achievements: (1) completion of the Reynolds-stress closure by developing a new pressure-strain correlation; (2) development of a parabolic code to compute jets and wakes; and, (3) application to a flow through a 180 deg turnaround duct by adopting a boundary fitted coordinate system. In the above mentioned models near-wall models are developed for pressure-strain correlation and third-moment, and incorporated into the transport equations. This addition improved the results considerably and is recommended for future computations. A new parabolic code to solve shear flows without coordinate tranformations is developed and incorporated in this study. This code uses the structure of the finite volume method to solve the governing equations implicitly. The code was validated with the experimental results available in the literature.

Amano, Ryoichi S.↗

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗

Tree-Based Ensemble Learning Models for Wall Temperature Predictions in Post-Critical Heat Flux Flow Regimes at Subcooled and Low-Quality Conditions

Accurately predicting post-critical heat flux (CHF) heat transfer is an important but challenging task in water-cooled reactor design and safety analysis. Although numerous heat transfer correlations have been developed to predict post-CHF heat transfer, these correlations are only applicable to relatively narrow ranges of flow conditions due to the complex physical nature of the post-CHF heat transfer regimes. In this paper, a large quantity of experimental data is collected and summarized from the literature for steady-state subcooled and low-quality film boiling regimes with water as the working fluid in vertical tubular test sections. In addition, a low-quality water film boiling (LWFB) database is consolidated with a total of 22,813 experimental data points, which cover a wide flow range of the system pressure from 0.1 to 9.0 MPa, mass flux from 25 to 2750 kg/m 2 s, and inlet subcooling from 1 to 70 °C. Two machine learning (ML) models, based on random forest (RF) and gradient boosted decision tree (GBDT), are trained and validated to predict wall temperatures in post-CHF flow regimes. The trained ML models demonstrate significantly improved accuracies compared to conventional empirical correlations. To further evaluate the performance of these two ML models from a statistical perspective, three criteria are investigated and three metrics are calculated to quantitatively assess the accuracy of these two ML models. For the full LWFB database, the root-mean-square errors between the measured and predicted wall temperatures by the GBDT and RF models are 5.7% and 6.2%, respectively, confirming the accuracy of the two ML models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗

Large-Eddy Simulation of Axisymmetric Compression Corner Flow

The Wall-Modeled Large Eddy Simulation (WMLES) approach is used to study the interaction of a shock wave with a high Reynolds number turbulent boundary layer. Since the near wall region is modeled, high Reynolds number turbulent flows can be simulated at a moderate computational cost. The case considered is that of an axisymmetric Mach 2.85 turbulent boundary layer over a 30 compression corner. The Reynolds number of the boundary layer upstream of the interaction based on momentum thickness (Re theta = u sub infinity theta/v sub infinity) is ~12,000. The geometry and flow conditions match the experiments of Dunagan et al. (NASA TM 88227, 1986). The simulations were performed using equilibrium and non-equilibrium wall models. The agreement with experiment is encouraging for the finest grid with respect to the separation bubble length, unsteady shock structure and wall pressure distribution. Sensitivity ofWMLES results to grid, wall model, and blockage effects in the tunnel are reported.

Iyer, Prahladh S.↗

Effect of Reynolds number variation on aerodynamics of a hydrogen-fueled transport concept at Mach 6

Two separate tests have been made on the same blended wing-body hydrogen-fueled transport model at a Mach number of about 6 and a range of Reynolds number (based on theoretical body length) of 1.577 to 55.36 X 10 to the 6th power. The results of these tests, made in a conventional hypersonic blowdown tunnel and a hypersonic shock tunnel, are presented through a range of angle of attack from -1 to 8 deg, with an extended study at a constant angle of attack of 3 deg. The model boundary layer flow appeared to be predominately turbulent except for the low Reynolds number shock tunnel tests. Model wall temperatures varied considerably; the blowdown tunnel varied from about 255 F to 340 F, whereas the shock tunnel had a constant 70 F model wall temperature. The experimental normal-force coefficients were essentially independent of Reynolds number. A current theoretical computer program was used to study the effect of Reynolds number. Theoretical predictions of normal-force coefficients were good, particularly at anticipated cruise angles of attack, that is 2 to 5 deg. Axial-force coefficients were generally underestimated for the turbulent skin friction conditions, and pitching-moment coefficients could not be predicted reliably.

Penland, Jim A.↗

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $ \$ $6B annually (at $ \$ $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year), and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE) missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

42 ENGINEERING↗

Surface Roughness Effects from Additive Manufacturing in High Efficiency Gas Turbine Combustion Systems

In the past, clean combustion systems were characterized by high system development costs due to growing complexity and escalating manufacturing costs from conventional manufacturing processes. As a result, high efficiency concepts are difficult to design and more difficult to be cost-effectively manufactured. In recent years, with the introduction of additive manufacturing (AM) technologies, the rapid prototyping and mass production processes of clean combustion systems are promising to be significantly simplified with significant reduction in terms of engine manufacturing cost and engine energy cost. Compared with conventional manufacturing built parts, the AM process enabled a simpler design to be adopted for the nozzle, reducing the number of required braze and weld joints from twenty-five to just five. The resulting nozzle was 25% lighter and five times more durable and contributed to a 15% reduction in fuel burn in comparison with the previous model produced. Applying these improvements to a fleet of 5,000 turbofan engines at 150 kN thrust would result in fuel cost savings of $6B annually (at $5/gallon Jet-A fuel price, 700 gallon/hour fuel consumption rate, and 2,300 operational hours per year),and reduce CO2 emissions by more than 25 million tons per year. These advances support core Department of Energy (DOE)missions in energy efficiency, improving productivity, and environmental sustainability. Maximizing the benefit of these new capabilities will require high prediction capability of high speed turbulent flow with wall modeled Large Eddy Simulation (LES). GE Aviation maintains that advanced simulation technology and supercomputing is required in order to provide the appropriate boundary conditions to realize the maximum potential of AM and Ceramic Matrix Composite to reduce cooling flow, a first order penalty on the Brayton cycle. The algorithm developed using the ANSYS/Fluent software would provide the foundation for entirely new avenues of research and development with potential multi-billion dollar impact to the US economy, and significant reduction in carbon based emissions across the aerospace and power generation industries. The formulation of this iWLES (integral wall model for LES) is generic and allows to capture the changes in flow dynamics that have a significant impact on the wall bounded flow characteristics, such as swirler effective area, bulk swirl number, local pressure distribution, exit velocity profile, and turbulence kinetic energy profile. A periodic channel flow with rough flat plate is simulated using LES (Wall-Adapting Local Eddy-viscosity) model to verify the implementation of iWLES in the Fluent User Defined Function (UDF). As the flow fields are highly sensitive due to surface roughness of the wall bounded flows in the combustion systems, there is significant potential to conduct further LES studies focusing on the turbulence boundary layer interaction. Accurately capturing near wall physics will elucidate the impact of rough surfaces on combustor flow and aero-thermal interactions.

99 GENERAL AND MISCELLANEOUS↗

High-Fidelity Simulations of a Mach 9 Shock-Turbulent Boundary Layer Interaction

High-fidelity large-eddy simulations (LES) and wall-modeled large-eddy simulations (WMLES) are performed on the interaction of a Mach 9 turbulent boundary layer with a 34-degree compression ramp. Our primary interest is the sensitivity of mean and unsteady wall quantities to grid density. Additionally within the context of WMLES, we seek to explore solution sensitivity to the wall-model exchange location.

42 ENGINEERING↗

Comparison of interference-free numerical results with sample experimental data for the AEDC wall-interference model at transonic and subsonic flow conditions

Numerical results obtained from two computer programs recently developed with NASA support and now available for use by others are compared with some sample experimental data taken on a rectangular-wing configuration in the AEDC 16-Foot Transonic Tunnel at transonic and subsonic flow conditions. This data was used in an AEDC investigation as reference data to deduce the tunnel-wall interference effects for corresponding data taken in a smaller tunnel. The comparisons were originally intended to see how well a current state-of-the-art transonic flow calculation for a simple 3-D wing agreed with data which was felt by experimentalists to be relatively interference-free. As a result of the discrepancies between the experimental data and computational results at the quoted angle of attack, it was then deduced from an approximate stress analysis that the sting had deflected appreciably. Thus, the comparisons themselves are not so meaningful, since the calculations must be repeated at the proper angle of attack. Of more importance, however, is a demonstration of the utility of currently available computational tools in the analysis and correlation of transonic experimental data.

Newman, P. A.↗

Expanded capabilities of the CarMa code in modeling resistive wall mode dynamics with 3-D conductors

In this work, an improved version of the CarMa code is presented, called CarMa-D, for the analysis of resistive wall modes (RWMs) in fusion devices, simultaneously considering the effects of volumetric three-dimensional conducting structures, plasma dynamics, toroidal rotation or drift-kinetic damping. Unlike static CarMa, CarMa-D does not rely on the simplifying assumptions such as neglecting the plasma mass, toroidal rotation and kinetic damping. The new coupling strategy is based on matrix-based Padé rational functions approximation of plasma a response. The arising mathematical model is formally equivalent to the original CarMa model, but with a higher number of degrees of freedom to model the dynamics of the plasma. CarMa-D overcomes the main limitations of the original CarMa, in particular: (i) the massless assumption for the plasma is removed, allowing modeling of global modes growing both on ideal kink time scales and in the typical RWM growth rate regime, with a suitable treatment of the model; (ii) the effects of toroidal plasma flow and drift kinetic damping can be included into the new model, providing a powerful tool to study macroscopic phenomena where both plasma dynamics and 3-D conducting structures play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A two-layer multiple-time-scale turbulence model and grid independence study

A two-layer multiple-time-scale turbulence model is presented. The near-wall model is based on the classical Kolmogorov-Prandtl turbulence hypothesis and the semi-empirical logarithmic law of the wall. In the two-layer model presented, the computational domain of the conservation of mass equation and the mean momentum equation penetrated up to the wall, where no slip boundary condition has been prescribed; and the near wall boundary of the turbulence equations has been located at the fully turbulent region, yet very close to the wall, where the standard wall function method has been applied. Thus, the conservation of mass constraint can be satisfied more rigorously in the two-layer model than in the standard wall function method. In most of the two-layer turbulence models, the number of grid points to be used inside the near-wall layer posed the issue of computational efficiency. The present finite element computational results showed that the grid independent solutions were obtained with as small as two grid points, i.e., one quadratic element, inside the near wall layer. Comparison of the computational results obtained by using the two-layer model and those obtained by using the wall function method is also presented.

Kim, S.-W.↗

Comments on the present state of second-order closure models for incompressible flows

Second-order closure models account for history and nonlocal effects of the mean velocity gradients on the Reynolds stress tensor. Turbulent flows involving body forces or curvature, Reynolds stress relaxational effects, and counter-gradient transport are usually better described. The topics are presented in viewgraph form and include: (1) the Reynolds stress transport equation; (2) issues in second-order closure modeling; and (3) near wall models.

Speziale, Charles G.↗

Heat Transfer in a Superelliptic Transition Duct

Local heat transfer measurements were experimentally mapped using a transient liquid-crystal heat transfer technique on the surface of a circular-to-rectangular transition duct. The transition duct had a length-to-diameter ratio of 1.5 and an exit-plane aspect ratio of 3. The crosssectional geometry was defined by the equation of a superellipse. The cross-sectional area was the same at the inlet and exit but varied up to 15 percent higher through the transition. The duct was preheated to a uniform temperature (nominally 64 C) before allowing room temperature air to be suddenly drawn through it. As the surface cooled, the resulting isothermal contours on the duct surface were revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record was made of the surface temperature and time data for all points on the duct surfaces during each test. Using this surface temperature-time data together with the temperature of the air flowing through the model and the initial temperature of the model wall, the heat transfer coefficient was calculated by employing the classic one-dimensional, semi-infinite wall heat transfer conduction model. Test results are reported for inlet diameter-based Reynolds numbers ranging from 0.4x106 to 2.4x106 and two grid-generated freestream turbulence intensities of about 1 percent, which is typical of wind tunnels, and up to 16 percent, which may be more typical of real engine conditions.

Poinsatte, Philip↗

A multiblock analysis for shuttle orbiter re-entry heating from Mach 24 to Mach 12

A multiblock, laminar heating analysis for the shuttle orbiter at three trajectory points ranging from Mach 24.3 to Mach 12.86 on re-entry is described. The analysis is performed using the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) with both a seven species chemical nonequilibrium model and an equilibrium model. A finite-catalytic-wall model appropriate for shuttle tiles at a radiative equilibrium wall temperature is applied. Computed heating levels are generally in good agreement with the flight data though a few rather large discrepancies remain unexplained. The multiblock relaxation strategy partitions the flowfield into manageable blocks requiring a fraction of the computational resources (time and memory) required by a full domain approach. In hot, the computational cost for a solution at even a single trajectory point would be prohibitively expensive at the given resolution without the multiblock approach. Converged blocks are reassembled to enable a fully coupled converged solution over the entire vehicle, starting from a nearly converged initial condition.

Gnoffo, Peter A.↗