Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Wall Model”

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 19 records

Wall-modeled LES of Turbulent Flow Over a Two-dimensional Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a nominally two-dimensional Gaussian-shaped bump geometry to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on high-fidelity numerical simulation of Uzun & Malik (AIAA Journal 2022). The oncoming flow Mach number is 0.2, with the bump length-based Reynolds number (ReL) of 2 million. In our previous study, while WMLES with the constant coefficient Vreman subgrid scale model and equilibrium wall model performed satisfactorily at lower ReL = 1 million, it failed to predict flow separation at the higher Re, contrary to the experimental observations. We investigate the sensitivity of WMLES to different subgrid scale models, wall models, and grid resolution and topology on flow separation by comparing with available data.

Computational Fluid Dynamics

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.

Bryson Frank

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the acceleration and separation regions of the flow. An unstructured finite-volume solver is used along with an equilibrium wall model. The Mach number of the oncoming flow selected for the simulation is 0.176, and two Reynolds numbers (Re) are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Spanwise periodic results for the lower Re case are compared with the available DNS data, while the full three-dimensional simulation results for the higher Re case are compared to available experimental data. Sensitivity of WMLES results are assessed for a number of factors including grid resolution, wall model exchange location, type of wall model, and unstructured grid topology. The WMLES results for the lower Re case agree well with available DNS data in terms of the wall pressure variation, velocity and turbulent stress profile comparisons. The skin-friction predictions show a reduced tendency to separate. This appears to be accentuated for the higher Re case, for the medium resolution grid used in this study, where WMLES does not predict any flow separation in the mean, which is in contrast to the large separation bubble observed in the experiments.

turbulence

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer

Transonic Lift and Drag Predictions using Wall Modelled Large Eddy Simulations

Wall-modelled Large Eddy Simulations (WMLES) of the NASA Common Research Model (CRM) at transonic conditions and various angles of attacks leading up to and including shock-induced flow separation are performed using the LAVA computational framework. The simulations are shown to accurately predict the lift curve slope and the onset of separation characterized by the break in the pitching moment. Furthermore, careful assessment of skin-friction drag at cruise condition and its subsequent decrease with increasing angles of attack are shown to be in agreement with viscous sublayer resolving Reynolds Averaged Navier Stokes(RANS) simulations. The small differences between WMLES and RANS appear to be of the same order as differences seen between two RANS models at the cruise-point. Some sensitivity is observed to the coefficient used in the subgrid scale model, although this can be reconciled by noting the low chord- and shock-incidence Reynolds numbers (Rex≈106) seen in the outboard regions of the wing, along with the uncertainties associated with tripping and the numerical transition that occurs near the leading edge. Among the primary drawbacks seen in the WMLES predictions is the subdued increase in wave-drag with increasing angle-of-attack and the shock intensity when compared with experimental data. Preliminary assessment suggests that additional span- and stream-wise grid resolution is likely needed in mid- and out-board portions of the wing to better resolve the shock-induced separated flow and to further investigate prediction accuracy of unsteady temporal characteristics of the problem.

ARMD

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on the ongoing CFD validation experiments of Slotnick [NATO STO-MP-AVT-307, 2019]. The oncoming flow Mach number is 0.176, and two Reynolds numbers are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Preliminary Reynolds-averaged Navier-Stokes simulations are first performed to assess the effects of Mach number, Reynolds number, tunnel top and sidewall effects. Finally, WMLES results with an equilibrium wall model will be presented at two Reynolds numbers to assess their performance for this flow by making detailed comparisons with available experimental and higher-fidelity numerical data.

turbulence, wall model, large eddy simulation

Predictions of LAGOON Nose Landing Gear Flow and Noise Using Wall-Modeled Large-Eddy Simulations

Wall-modeled large-eddy simulations (WMLESs) of the LAGOON nose landing gear are conducted with compressible Navier–Stokes equations and immersed boundary technique using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. The simulations are conducted using six different Cartesian octree meshes for the grid sensitivity analysis of the near-field and far-field numerical predictions, where the far-field noise results are computed with the Ffowcs Williams–Hawkings acoustic analogy. The effects of numerical tripping induced at the exact locations of the tripping devices in the experiments are also examined. In general, better comparison with the experimental results are shown for the the near-field results obtained with the simulations under the effects of numerical tripping. The effects of tripping are not significant on the far-field noise calculations and the results have reasonable comparison with the experimental data in the low and medium frequency ranges when an impermeable formulation of the acoustic analogy is used.

CST

Near-wall modelling of compressible turbulent flows

Work was carried out to extend the near-wall models formulated for the incompressible Reynolds stress equations to compressible flows. The idea of splitting the compressible dissipation function into a solenoidal part that is not sensitive to changes of compressibility indicators and a compressible part that is directly affected by these changes is adopted. This means that all models involving the dissipation rate could be expressed in terms of the solenoidal dissipation rate and an equation governing its transport could be formulated to close the set of compressible Reynolds stress equations. The near-wall modelling of the dissipation rate equation is investigated and its behavior near a wall is studied in detail using k-epsilon closure. It is found that all existing modelled equations give the wrong behavior for the dissipation rate near a wall. Improvements are suggested and the resultant behavior is found to be in good agreement with near-wall data. Furthermore, the present modified k-epsilon closure is used too calculate a flat plate boundary layer and the results are compared with four existing k-epsilon closures. These comparisons show that all closures tested give essentially the same flow properties, except in a region very close to the wall. In this region, the present k-epsilon closure calculations are in better agreement with measurements and direct simulation data; in particular, the behavior of the dissipation rate.

So, Ronald M. C.

Near-wall modeling of the dissipation rate equation

Near-wall modeling of the dissipation rate equation is investigated and its asymptotic behavior is studied in detail using a k-epsilon model. It is found that all existing modeled dissipation rate equations predict an incorrect behavior for the dissipation rate near a wall. An improvement is proposed and the resulting near-wall dissipation rate distribution is found to be similar to that given by numerical simulation data. To further validate the improved k-epsilon model, it is used to calculate flat-plate turbulent boundary-layer flows at high- as well as low-turbulence Reynolds numbers, and the results are compared with measurements, numerical simulation data, and the calculations of three different two-equation models. These comparisons show that all the models tested give essentially the same flow properties away from the wall; significant differences only occur in a region very close to the wall. In this region, the calculations of the improved k-epsilon model are in better agreement with measurements and numerical simulation data. In particular, the modeled distribution of the dissipation rate is significantly improved and a maximum is predicted at the wall instead of away from the wall. Furthermore, the improved k-epsilon model is found to be the most asymptotically consistent among the four different two-equation models examined.

So, R. M. C.

A Wall-Modeled LES Perspective for the High Lift Common Research Model Using LAVA

A new immersed boundary Wall-Modelled Large Eddy Simulation (WMLES) formulation is developed to study high-lift aerodynamics on the NASA High-Lift Common Research Model (HL-CRM). A sequence of Cartesian Octree grids with sizes ranging from 100 Million through 2.02 Billion grid points is utilized to systematically assess grid-sensitivity and convergence for the in-tunnel (QinetiQ) configuration of the model, and remarkable agreement between the immersed boundary and the curvilinear body-aligned WMLES formulations is reported on grids with comparable resolutions. In the free-air configuration of the model, consistent predictions between the Curvilinear Overset and the Cartesian Octree formulations are reported for angles of attack up to C(L,max) at a=19.57. However, some differences in the onset of stall are seen between the two methods for a>20°: while the curvilinear WMLES experiences wing-root separation with increasing angle of attack (Topology A), the Cartesian Octree formulation shows a different flow topology characterized by boundary layer weakness on the main element, emanating from the pylon-wing attachment (Topology B). In order to obtain further insight into the two-distinct topologies, carefully designed numerical experiments to isolate effects of the model standoff and the tunnel wall-boundary layers are conducted using the immersed boundary WMLES formulation. The increased incidence angle-of-attack on the inboard portion of the wing due to the standoff is shown to be sufficient for triggering a switch from Topology B to Topology A in Cartesian WMLES. The role of the floor boundary layer is further examined in detail by identification of additional corner-flow vorticity generated by the viscous juncture flow interactions between the floor boundary layer and the standoff leading to formation of a strong coherent and persistent vortex on the belly-side of the fuselage. The intensity of this vortex is shown to increase with the thickness of the floor boundary layer. A further increase in the incidence angle of attack near the leading-edge strake caused by the presence of this belly-side vortex is quantified for two-distinct floor boundary layers. Both of the floor boundary layers considered result in the onset of large scale wing-root separation at a=21.47in non-confined (free-air) configurations.

TTT

Near-wall modelling of compressible turbulent flows

Work was carried out to formulate near-wall models for the equations governing the transport of the temperature-variance and its dissipation rate. With these equations properly modeled, a foundation is laid for their extension together with the heat-flux equations to compressible flows. This extension is carried out in a manner similar to that used to extend the incompressible near-wall Reynolds-stress models to compressible flows. The methodology used to accomplish the extension of the near-wall Reynolds-stress models is examined and the actual extension of the models for the Reynolds-stress equations and the near-wall dissipation-rate equation to compressible flows is given. Then the formulation of the near-wall models for the equations governing the transport of the temperature variance and its dissipation rate is discussed. Finally, a sample calculation of a flat plate compressible turbulent boundary-layer flow with adiabatic wall boundary condition and a free-stream Mach number of 2.5 using a two-equation near-wall closure is presented. The results show that the near-wall two-equation closure formulated for compressible flows is quite valid and the calculated properties are in good agreement with measurements. Furthermore, the near-wall behavior of the turbulence statistics and structure parameters is consistent with that found in incompressible flows.

So, Ronald M. C.

Wall-modeled LES of the Three-dimensional Speed Bump Experiment

We evaluate the performance of wall-modeled large eddy simulation (WMLES) in predicting smooth-body turbulent flow separation for a three-dimensional Gaussian-shaped speed bump geometry. The Reynolds number based on the bump length is 2 million, and an unstructured compressible finite-volume solver is used with an equilibrium wall model and dynamic subgrid scale model. Spanwise periodic simulations of the centerline two-dimensional bump were used to assess grid resolution requirements, and resolving the thin internal layer in the accelerating region was found to be necessary to correctly capture the downstream separated flow region. Based on these insights, an optimized grid that was smaller by a factor of two provided comparable accuracy to a finer grid that has been used by us and other researchers in past studies. Using insights gained from the spanwise periodic simulation, an unstructured polyhderal grid with about 250 million cells was used for the three-dimensional (3D) configuration with inviscid tunnel side and top walls. Detailed comparisons of wall skin-friction coefficient, wall pressure, velocity and turbulent stresses with available experimental data indicated excellent agreement. While the low-Reynolds Number Spalart-Allmaras RANS model with rotation/curvature correction gave qualitatively good agreement with experiments, WMLES showed significantly more accurate quantitative predictions of the flowfield. The 3D WMLES showed good agreement with experiments in the separated region, and the centerplane results indicated a different separation topology compared to the spanwise periodic simulation.

Computational Fluid Dynamics

Development of a turbulence near-wall model and its application to separated and reattached flows

A numerical study is reported of flow and heat transfer in the separated and reattached flows created by an abrupt pipe expansion. In the study attention has been given primarily to the development of turbulence near-wall models based on the wall function. An efficient numerical method is also employed for the computation of high-velocity flows. The computed results are compared with experimental data obtained earlier. Generally, better results are obtained by employing the present near-wall models, and among them the three-layer model is superior to the two-layer one.

Amano, R. S.

Towards a Viscous Wall Model for Immersed Boundary Methods

Immersed boundary methods are frequently employed for simulating flows at low Reynolds numbers or for applications where viscous boundary layer effects can be neglected. The primary shortcoming of Cartesian mesh immersed boundary methods is the inability of efficiently resolving thin turbulent boundary layers in high-Reynolds number flow application. The inefficiency of resolving the thin boundary is associated with the use of constant aspect ratio Cartesian grid cells. Conventional CFD approaches can efficiently resolve the large wall normal gradients by utilizing large aspect ratio cells near the wall. This paper presents different approaches for immersed boundary methods to account for the viscous boundary layer interaction with the flow-field away from the walls. Different wall modeling approaches proposed in previous research studies are addressed and compared to a new integral boundary layer based approach. In contrast to common wall-modeling approaches that usually only utilize local flow information, the integral boundary layer based approach keeps the streamwise history of the boundary layer. This allows the method to remain effective at much larger y+ values than local wall modeling approaches. After a theoretical discussion of the different approaches, the method is applied to increasingly more challenging flow fields including fully attached, separated, and shock-induced separated (laminar and turbulent) flows.

Immersed

High-Order Wall-Modeled Large-Eddy Simulation of High-Lift Configuration

This paper presents the assessment of several recent enhancements for a high-order wall-modeled large-eddy simulation (WMLES) approach and demonstrates order independence with a fixed data exchange location in the wall model. The two enhancements include the use of isotropic tetrahedral elements to improve accuracy and an explicit subgrid-scale model, the Vreman model, to improve accuracy and robustness. The [Formula: see text] study focused on the high-lift Common Research Model (HL-CRM) at the angle of attack of 19.57 deg, a benchmark problem from the 4th AIAA High-Lift Prediction Workshop. Solution polynomial orders of [Formula: see text], and 5 were used in the study. The study demonstrated [Formula: see text] independence in integrated forces, pitch moment, velocity profile in the wall-normal direction, and surface flow topology. It also showed that a [Formula: see text] order of at least 3 ([Formula: see text]) was needed to correctly predict the external inviscid flow and the surface flow topology. Thereafter, [Formula: see text] simulations over several other angles of attack demonstrated that the high-order WMLES approach can correctly predict the maximum lift and flow separation regions for HL-CRM with about 40 million degrees of freedom (DOF) compared to at least 250 million DOF required by second-order methods.

Engineering

Turbulent transport modeling of shear flows around an aerodynamic wing. Development of turbulent near-wall model and its application to recirculating flows

Progress in implementing and refining two near-wall turbulence models in which the near-wall region is divided into either two or three zones is outlined. These models were successfully applied to the computation of recirculating flows. The research was further extended to obtaining experimental results of two different recirculating flow conditions in order to check the validity of the present models. Two different experimental apparatuses were set up: axisymmetric turbulent impinging jets on a flat plate, and turbulent flows in a circular pipe with a abrupt pipe expansion. It is shown that generally better results are obtained by using the present near-wall models, and among the models the three-zone model is superior to the two-zone model.

Amano, R. S.

Wall Modeled Large Eddy Simulations for NASA’s Jet Noise Consensus Database of Single-Flow, Round, Convergent Jets

A campaign of wall-modeled large-eddy simulations (WMLES) using structured curvilinear overlapping grids has been performed with the Launch Ascent and Vehicle Aerodynamics(LAVA) computational fluid dynamics (CFD) software to predict jet noise for single-stream axisymmetric round jets. The simulations address the new Prediction Uncertainty Reduction(PUR) technical challenge within the context of NASA’s Commercial Supersonic Technology(CST) project. The goal of PUR is to quantify and reduce uncertainties from scale-resolving simulations to assess noise characteristics of next generation quiet supersonic commercial jets during takeoff and landing conditions where the noise from the exhaust jet dominates. The focus of this effort is to generate a simulation database for single-stream axisymmetric round nozzles at several conditions both for static (no ambient co-flow), which is the focus of this article, and in-flight (flight stream co-flow) conditions, which are beyond the current scope. Nine different flow conditions ranging in jet exit Mach number from 0.38 to 1.0 with nozzle temperature ratios (NTR) from 0.84 to 2.7 have been conducted. Details of the structured overset grids, numerical discretization and wall-model are provided. Near-field comparisons to PIV data show great agreement for both velocity and normal stresses, however a systematic TKE overshoot at the nozzle exit is seen in the lip line shear-layer. A permeable Ffowcs Williams Hawkings (FWH) surface, enclosing the jet, is used to predict far-field noise from the simulated flow-field. Comparison of CFD predictions to microphone array measurements demonstrate excellent agreement within the resolved frequency range. A systematic under-prediction of far-aft observer angles larger than 150 degrees has been observed across all simulations. We achieved a cost reduction of an order of magnitude for these WMLES compared to an earlier study of this configuration due to algorithmic and software improvements. The accuracy of the results and short turnaround time demonstrate that WMLES within the LAVA framework is a cost-effective approach for jet noise predictions that could soon be incorporated into the design cycle of jet noise reduction technologies.

CST