Search NASA⌕ Search

SEARCH · Search NASA

Results for “BOUNDARY LAYER - EFFECT OF COMPRESSIBILITY”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Near-wall model for compressible turbulent boundary layers based on an inverse velocity transformation

In this work, a near-wall model, which couples the inverse of a recently developed compressible velocity transformation (Griffin et al., Proc. Natl Acad. Sci., vol. 118, 2021, p. 34) and an algebraic temperature–velocity relation, is developed for high-speed turbulent boundary layers. As input, the model requires the mean flow state at one wall-normal height in the inner layer of the boundary layer and at the boundary-layer edge. As output, the model can predict mean temperature and velocity profiles across the entire inner layer, as well as the wall shear stress and heat flux. The model is tested in an a priori sense using a wide database of direct numerical simulation high-Mach-number turbulent channel flows, pipe flows and boundary layers (48 cases, with edge Mach numbers in the range 0.77–11, and semi-local friction Reynolds numbers in the range 170–5700). The present model is significantly more accurate than the classical ordinary differential equation (ODE) model for all cases tested. The model is deployed as a wall model for large-eddy simulations in channel flows with bulk Mach numbers in the range 0.7–4 and friction Reynolds numbers in the range 320–1800. When compared to the classical framework, in the a posteriori sense, the present method greatly improves the predicted heat flux, wall stress, and temperature and velocity profiles, especially in cases with strong heat transfer. In addition, the present model solves one ODE instead of two, and has a computational cost and implementation complexity similar to that of the commonly used ODE model.

42 ENGINEERING↗
DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗
DOE OSTI2020

Lift Equivalence and Cancellation for Airfoil Surge–Pitch–Plunge Oscillations

A NACA 0018 airfoil in freestream velocity is oscillated in longitudinal, transverse, and angle-of-attack directions with respect to the freestream velocity, known as surge, plunge, and pitch. The lift-based equivalence method introduces phase shifts between these three motions to construct in-phase sinusoidal components for maximum lift, waveform construction. Lift cancellation is also determined with the exact negative pitch and plunge motion amplitudes found from the equivalence method to achieve out-of-phase wave destruction. Lift cancellation occurs when a combination of these motions is sought to obtain a constant lift magnitude throughout the oscillation cycle. To achieve both equivalence and cancellation of lift, a prescribed pure pitch amplitude through the Theodorsen theory equates the corresponding equivalent plunge amplitude and pitch–plunge phase shift. These Theodorsen, linear superposition findings of pitch–plunge are leveraged toward the Greenberg theory to determine a closed-form, surge–pitch–plunge solution through the addition of a surge–plunge phase shift and optimal surge amplitude for lift cancellation. The lift cancellation surge–pitch–plunge amplitudes define the equivalence amplitude investigated here and theoretically limit the experiment to combinations of the first lift harmonic of the Greenberg theory. The analytical results are then compared with experimental lift force measurements and dye visualization. Finally, the normalized lift differences due to unsteady wake and boundary-layer behavior are examined to explore the extents of the Greenberg theory for these cases of lift-based equivalence and cancellation.

42 ENGINEERING↗
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗
DOE OSTI2025

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↗
DOE OSTI2022

Compressive properties of silicone Bouligand structures

Abstract This article presents an examination of silicone, Bouligand lattices in compression. Appearing frequently in biological organisms and manufacturing design, Bouligand structures comprise layers of parallel strands or fibers oriented in a helical fashion. They can exhibit exceptional fracture resistance when composed of rigid or composite materials. The behavior of elastomeric Bouligand structures, however, is less well understood. Additively manufactured (AM) elastomeric lattices have applications in stress mitigation, medical devices, and soft robotics. This article demonstrates that Bouligand structures are a useful addition to the design space of AM elastomers. By adjusting the layer-rotation parameters, lattice stress can increase by more than 300% without altering the porosity. Additionally, we introduce path length metrics that help explain the observed relationship between layer rotation and compression response. Impact statement Additive manufacturing (AM) continues to push the boundary of manufacturable structures and enhance the ability to robustly design for specific properties and behaviors. The more we understand the design space of a novel AM microstructure, the greater its application range. In this article, we describe the mechanical behavior of helicoidal, elastomeric lattices and introduce path length metrics to help explain their stress response. We show that these structures can exhibit a large range of mechanical behaviors in compression, making them well suited for applications such as stress mitigation and impact absorption. Additionally, the path length metrics could become useful design tools and may be applicable to a larger set of cellular structures. These findings expand our ability to rapidly design materials with highly specific and customizable properties to meet the needs of modern engineering challenges. Graphical abstract

Van Meerbeek, I. M. (ORCID:000000021258132X)↗
DOE OSTI2024

Generalized grain boundary constitutive description implemented in a strain-gradient large-strain FFT-based formulation: Application to nano-metallic laminates

Herein this paper presents a general treatment of grain boundary constitutive behavior in the context of strain-gradient (SG) plasticity, and its numerical implementation in a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. Two novel grain boundary constitutive equations are proposed, allowing for more accurate description of the Burgers vector flow at the grain boundary. The capabilities of the generalized SG-LS-EVPFFT formulation are illustrated for the case of kink-band formation during layer-parallel compression of nano-metallic laminates (NMLs), requiring consideration of the interaction between dislocations and interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stress corrosion cracking mechanism of cold spray coating on a galvanically similar substrate

Here, the chloride-induced stress corrosion cracking (CISCC) mechanism of cold spray (CS) coating on a galvanically similar substrate is investigated. Arc welded 304L stainless steel (SS) specimens are loaded into four-point bend fixtures, cold sprayed with 304L SS, then immersed in boiling MgCl 2 . Interconnected porosity forms through crevice corrosion along CS splat boundaries, allowing corrosive species to penetrate through the CS layer. Nevertheless, the substrate is resistant to CISCC likely because of compressive stress introduced by peening during CS particle impacts. These findings underscore the importance of residual stress in the environmental degradation of CS coatings or repairs of engineering structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Gas permeation properties of amorphous zeolitic imidazolate framework membranes made by atomic/molecular layer deposition

Amorphous metal-organic framework (MOF) membranes are desirable because they may retain some of the molecular sieving properties of their crystalline counterparts while being free of grain boundary defects, which often hinder the consistent achievement of high membrane performance. However, current methods, like melting and compression, for fabricating amorphous MOF membranes involve multi-step processes that require the formation of a crystalline membrane first, that is then amorphized, and therefore, could be challenging to scale. Here, we utilize atomic/molecular layer deposition (ALD/MLD) of diethylzinc (DEZ) and 2-methylimidazole (2mIm) to directly synthesize ultrathin amorphous zeolitic imidazolate framework (aZIF) deposits on γ-alumina-coated α-alumina supports. As the number of ALD/MLD cycles increased from 10 to 300, gas permeances decreased while ideal selectivities increased. Mixture separation factors for C3H6/C3H8, CO2/N2, and H2/C3H8 as high as 4, 37, and 194, respectively, were obtained. At 200 °C and 2.5 bar equimolar feed of H2 and C3H8, a H2/C3H8 mixture separation factor of 185 is obtained with an H2 permeance of ca. 4.15x10-8 mol/m2-s-Pa (124 GPU). Additionally, the membrane achieves a CO2/N2 mixture separation factor of 37 at 25 °C and 1 bar with a CO2 permeance of ca. 3.96x10-8 mol/m2-s-Pa (118 GPU). Considering the vast array of compositionally distinct aZIFs that can be potentially deposited by this approach, an enormously large parameter space for membrane design is emerging to be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Progress report on analytic and numerical studies of x-ray-induced impulse

In a variety of high-energy-density (HED) systems, x-rays of a given energy are used to generate shockwaves, bulk motion, and impulse in materials. As radiation energy is deposited within the material, the material is heated. The number of photons of a certain wavelength absorbed is determined by the spectral intensity of the radiation and the material’s opacity evaluated at that wavelength. This heating results in a pressure increase dictated by the material’s equation of state. Depending on the intensity of the absorbed radiation, it may also cause the material to change phase into a liquid, gas, or plasma. The increased pressure drives the heated surface layer to blow off, imparting impulse and sending a compressive wave into the bulk of the material. Additionally, the compression wave interacts with the solid boundary of the material, resulting in a tensile wave that may cause the material to spall. The impulse generated by the deposition of x-ray energy within the sample can be modeled using purely analytical methods, e.g. the Bethe, Bade, Averell, and Yost (BBAY) model. However, the blow-off process is rather complicated, and proper modeling efforts must account for material ejected by spallation, vaporization, jetting, and plasma ablation. For this reason, analytical models have an unclosed term describing the final energy of the blown-off material Ef(z). Prior modeling efforts have arbitrarily fixed this at some value or modeled it using a limiting set of thermodynamic assumptions. The work we are currently performing uses validated simulations using sophisticated photon transport, equation of state, and strength models/data to provide a fit for Ef(z) that is useful for predictive calculation of impulse. We will apply our methodology and show results for different materials and x-ray sources. This work is particularly useful for the design of experiments studying x-ray impulse generation. The present report is outlined as follows. Section 1.2 describes a series of HED experiments investigating x-ray-generated impulse in materials. Section 1.3 describes the computational approach we employ in this study, and presents validation results against the aforementioned experiments. Section 1.4 introduces analytical models for impulse generation, as well as our method for utilizing impulse from simulations to close the models. We discuss the concept of impulse-spectrum sensitivity, it’s application to uncertainty quantification, and derive a very useful analytical expression for it in 1.5. We then summarize recent progress in this project and discuss future work in section 1.6.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Development and Validation of Predictive Models for In-Cylinder Radiation and Wall Heat Transfer

The lack of accurate submodels for in-cylinder radiation and heat transfer has been identified as a key shortcoming in developing truly predictive, physics-based computational fluid dynamics (CFD) models that can be used to develop and design combustion systems for advanced high-efficiency, low-emissions engines. Measurements of wall boundary layers in engines show discrepancies of up to 100% with respect to standard CFD boundary layer models. And analysis of in-cylinder radiation based on the most recent spectral property databases and high-fidelity radiative transfer equation solvers has shown that at operating pressures and exhaust-gas recirculation levels typical of modern compression-ignition truck engines, molecular gas radiation (mainly carbon dioxide and water vapor) can be more important than soot particle radiation, and that a significant fraction of the emitted radiation (50% or more) can be reabsorbed before reaching the walls, thereby acting to redistribute energy in the combustion chamber. In this project, an integrated experimental and simulation approach has been used to generate new physical insight into in-cylinder heat-transfer processes in engines, and to develop computationally practicable CFD-based models that accurately represent the physics. New experimental data have been taken in an optically accessible single-cylinder research engine, and in a multi-cylinder production engine. Those data have provided new insight into the radiative environment in engines. New CFD models have been developed and validated for in-cylinder heat transfer, focusing mainly on radiative heat transfer and coupling between turbulent boundary layer wall heat transfer and radiative heat transfer. The CFD models have been developed and implemented in an open-source software framework that can be shared with other researchers. The models are compatible both with unsteady Reynolds-averaged Navier Stokes (URANS) formulations and with large-eddy simulations (LES) of in-cylinder processes in engines.

02 PETROLEUM↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Nucleation of Frank Dislocation during the Squeeze-Out Process in Boundary Lubrication: A Molecular Dynamics Study

Liquid–vapor molecular dynamics (LVMD) simulations are performed to reinvestigate the phase transition and solvation force oscillation behavior of a simple argon liquid film confined between two solid surfaces. Our simulations present a novel scenario in which the n → n – 1 layering transitions are accompanied by the formation, climb, and annihilation of Frank partial dislocations during the squeeze-out process under compression. This is indicated by the splitting of the repulsive peaks in the solvation force profile. The detailed analysis reveals that the formation–climb–annihilation mechanism of Frank dislocation occurs during approach and disappears during receding, which would result in force hysteresis. In combination with our recent works, this study provides new insights into the physical property of nanoconfined lubricant films in boundary lubrication.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Kink mechanism in Cu/Nb nanolaminates explored by $\mathcal{in}$ $\mathcal{situ}$ pillar compression

We report Nano metallic laminates (NMLs) exhibit different failure modes depending on the loading conditions due to their mechanical anisotropies. Kinking is a typical failure mode in many NMLs compressed along a layer-parallel direction. However, a detailed description of the microstructure evolution during kink band (KB) formation and an in-depth understanding of the formation mechanisms are lacking. In this work, the KB process is investigated in Cu/Nb NMLs by in situ micro pillar compression in the scanning electron microscope (SEM) along a layer-parallel direction. Post-mortem S/TEM and transmission Kikuchi diffraction (TKD) analyses show that kink banding leads to significant microstructure changes characterized by an accumulation of geometrically necessary dislocations (GNDs) and of tilt geometrically necessary boundaries (GNBs) near KB boundaries (KBBs). The distinct microstructure evolution implies that KB formation is facilitated by the inhomogeneous microstructures resulting in constrained deformation modes. Specifically, dislocations active on slip planes nearly parallel to the interfaces make a major contribution to kink evolution after the onset of kinking. Once layer-parallel slip systems are activated, preexisting lattice dislocations and dislocations nucleating from interfaces will accumulate as GNDs near KBBs via the stochastic storage of lattice dislocations that have certain Burgers vectors. GNDs can further transform into GNBs via cross-slip and climb driven processes near the KBB. Furthermore, GNBs near KBBs can grow by incorporating more GNDs or by coalescence to accommodate the KB evolution. We further hypothesize that microstructural perturbations and their ensuing stresses can initiate KB formation in Cu/Nb NMLs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Deformation behavior of crystalline/amorphous Al-Si nanocomposites with nanolaminate or nanofibrous microstructures

Deformation mechanisms in sputter-deposited crystalline Al/amorphous Si nanocomposites with nanolaminate or nanofibrous morphology are characterized by nanoindentation, micropillar compression testing and transmission electron microscopy (TEM). The nanofibrous composite having crystalline Al nanofibers with ~40–50 nm in length and 15–20 nm in diameter embedded in amorphous Si exhibits strain hardening to a maximum flow stress of 2.9 GPa and no shear band (SB) formation in compression up to plastic strain exceeding 24%. On the other hand, nanolaminate composite that is composed of 80 nm crystalline Al layers and 20 nm amorphous Si layers exhibits catastrophic SBs starting at plastic strains in the range of 5–10%. Cross-sectional TEM of the deformed samples reveals a high density of stacking faults and twin boundaries in Al nanofibers and no microshear bands in the nanofibrous composite, suggesting plastic deformation in amorphous Si phase and crystalline Al nanofibers. Further, molecular dynamics simulations revealed that the plastic deformation in amorphous Si phase in the cosputtered films could be favored by the decrease in flow strength of amorphous Si with increasing Al solute concentration trapped in Si.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗