Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW VELOCITY, HIGH - SUPERSONIC”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: recent advances

Basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of wave-kinetic theory. Recently, this framework has been re-examined from first principles, which has led to more accurate yet still tractable "improved" wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow "curvature" (the second radial derivative of the flow velocity) on dynamics and stability of drift waves and zonal flows. We overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

fusion plasmas↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Toward ultra-efficient high-fidelity predictions of wind turbine wakes: Augmenting the accuracy of engineering models with machine learning

This study proposes a novel machine learning (ML) methodology for the efficient and cost-effective prediction of high-fidelity three-dimensional velocity fields in the wake of utility-scale turbines. The model consists of an autoencoder convolutional neural network with U-Net skipped connections, fine-tuned using high-fidelity data from large-eddy simulations (LES). The trained model takes the low-fidelity velocity field cost-effectively generated from the analytical engineering wake model as input and produces the high-fidelity velocity fields. The accuracy of the proposed ML model is demonstrated in a utility-scale wind farm for which datasets of wake flow fields were previously generated using LES under various wind speeds, wind directions, and yaw angles. Comparing the ML model results with those of LES, the ML model was shown to reduce the error in the prediction from 20% obtained from the Gauss Curl hybrid (GCH) model to less than 5%. In addition, the ML model captured the non-symmetric wake deflection observed for opposing yaw angles for wake steering cases, demonstrating a greater accuracy than the GCH model. The computational cost of the ML model is on par with that of the analytical wake model while generating numerical outcomes nearly as accurate as those of the high-fidelity LES.

Mechanics↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Experimental characterization of extreme temperature granular flows for solar thermal energy transport and storage

High-temperature, dense granular flows along an inclined plane were considered for solar thermal energy transport and storage with sintered bauxite particles. A series of experiments was performed for particle inlet temperatures of ~ 200, 400, 600, and 800 °C to understand the mechanisms of granular flows at extreme temperatures. Mass flow rates were measured using a load cell and free-surface velocities were measured and computed using particle image velocimetry. Surface temperatures were measured using infrared cameras. A significant decrease in steady-state particle mass flow rate was observed with increasing temperature due to changing flow properties. A decrease in bulk particle free-surface velocities was observed at higher temperatures. Free-surface velocity measurement error between experiments were within 20% of the average. The particle surface temperatures decreased from inlet to outlet with larger gradients at higher temperatures observed due to increasing convection and radiative heat losses. Here, a decrease in temperature was observed along the side walls due to a decrease in particle velocities.

14 SOLAR ENERGY↗
DOE OSTI2023

Wintertime Supraglacial Lake Drainage Cascade Triggers Large‐Scale Ice Flow Response in Greenland

Surface melt forces summertime ice-flow accelerations on glaciers and ice sheets. Here, we show that large meltwater-forced accelerations also occur during wintertime in Greenland. We document supraglacial lakes (SGLs) draining in cascades at unusually high elevation, causing an expansive flow acceleration over a ~5,200 km 2 region during winter. The three-component interferometric surface velocity field and decomposition modeling reveal the underlying flood propagation with unprecedented detail as it traveled over 160 km from the drainage site to the margin, providing novel constraints on subglacial water pathways, drainage morphology, and links with basal sliding. The triggering SGLs continuously grew over 40 years and suddenly released decades of stored meltwater, demonstrating surface melting can impact dynamics well beyond melt production. We show these events are likely common and thus their cumulative impact on dynamics should be further evaluated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Turbulence generation by shock interaction with a highly nonuniform medium

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, $\sigma_\rho/\bar{\rho}\sim1$) and strong shocks (shock Mach numbers $\bar{M}_s\gtrsim10$). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effects of radiation pressure on the evaporative wind of HD 209458b

ABSTRACT The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet’s host star has been proposed as a mechanism to drive the escaping atmosphere into a ‘cometary’ tail and explain the high velocities observed in systems where mass-loss is occurring. In this paper, we present results from high-resolution 3D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outwards from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-α flux required to drive the photoevaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-α fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Directional Variations in Tidal Flow Multifractality and Intermittency

We explored distinct directional variations of multifractal and intermittent characteristics of ebb and flood flow velocities at nodule point, WA, tidal energy site, and complementary inspection on the East River, NY using scaling exponents of the structure function, distribution flatness, detrending moving average (DMA) analysis, multifractal detrended fluctuation analysis (MF-DFA), and high-order spectral moments. Our findings reveal that tidal flow presents higher long-range dependence (LRD) and intermittent levels for the ebb flow sections, whereas the flood flow sections exhibit a higher degree of multifractality and greater sensitivity to larger magnitude of turbulent fluctuations. We demonstrate that long-range dependence predominantly contributes to multifractal behavior in both ebb and flood flows, as evidenced by the significantly reduced multifractal spectrum width for temporally randomly permuted time series. Moreover, spectral kurtosis analysis uncovers a higher intermittent level across all frequency scales for ebb flow sections and reveals a distinct pattern of tidal flow intermittency differing from the monotonically increasing intermittent level observed in wall-bounded and grid turbulence. Finally, we demonstrate that DMA, MF-DFA, and high-order spectral moments provide more comprehensive insights than structure function scaling exponents and PDF flatness methods.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗
DOE OSTI2024

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2022

An improved pressure drop correlation for modeling localized effects in a pebble bed reactor

Advances in the development of pebble bed reactors (PBRs) has created a desire for accurate and cost-effective simulation tools for design scoping studies and safety analysis. The current state-of-the-art for these simulations is the use of porous media models, although these models rely on correlations to capture the effects of flow features that are not explicitly modeled. One of the areas where correlation accuracy is currently lacking is in the near-wall region of the bed. In this region, the presence of the wall causes the pebbles to pack more orderly, drastically changing the geometry and flow behavior in this region. This work presents a new generalized pressure drop correlation for PBRs based on the KTA equation. A high-to-low methodology is applied, where large eddy simulation (LES) is performed on two beds of 1568 and 1700 pebbles to generate a high-fidelity dataset. The flow fields are then averaged in time and separated into concentric rings of 0.05 D peb width. Average porosity, velocity, and pressure drop are extracted for each ring and the friction and form losses are calculated. The Reynolds number range for this study is 625–10,000, and thus the form losses are dominant over the friction losses. The form losses across the rings are investigated, and a correction term for the form loss calculation is determined and applied to the KTA equation to drastically improve the capability of modeling localized porosity effects in a porous media code. Finally, the improved correlation reduces near-wall velocity prediction error from over 50% with the KTA correlation to around 5%. Agreement in pressure drop prediction between LES and porous media simulations is also improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: Recent advances

The basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of the wave-kinetic theory. Recently, this framework has been reexamined from first principles, which has led to more accurate yet still tractable “improved” wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow “curvature” (the second radial derivative of the flow velocity) on the dynamics and stability of drift waves and zonal flows. Here, we overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurements of the Flow in the Vicinity of an Additively Manufactured Turbine Leading-Edge Using X-Ray Particle Tracking Velocimetry

X-ray particle tracking velocimetry (PTV) is performed, for the first time, to measure the velocity field inside a leading-edge of a turbine blade made by laser-additive-manufacturing (LAM) process. The traditional showerhead holes were replaced by a porous matrix in the leading-edge. The flow through such a leading-edge piece cannot be faithfully recreated by traditional prototype testing methods due to the surface roughness and imperfections caused by LAM process. Hence, direct measurement is the only option. However, it is difficult to measure flow inside such pieces with traditional velocimetry measurements due to the existence of metallic walls. Moreover, small internal size and high flow speeds call for a measurement technique with high spatial and temporal resolutions. To address these issues, we performed time-resolved X-ray PTV using the Advanced Photon Source (APS) synchrotron facility at the Argonne National Laboratory (ANL). A hydraulic system was constructed to run water, mixed with seeding particles, through the leading-edge piece. A high-speed camera captured the images of the seeding particles, which were later processed to create particle tracks. The time-averaged velocity field showed distinct pairs of vortices located in front of the porous outlet inside the leading-edge piece. The inlet channel showed reversed flow due to partial obstruction by the porous inlet of the test piece. Such knowledge of the flow field inside a leading-edge of a turbine blade will help us to design better cooling paths leading to higher cooling efficiency and increased life-span of a turbine blade.

42 ENGINEERING↗
DOE OSTI2024

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↗