In the present article, we study the possibilities of machine learning for estimation of seeing at the Maidanak Astronomical Observatory (38∘40′24′′N, 66∘53′47′′E) using only Era-5 re-analysis data. Seeing is usually associated with the integral of the turbulence strength Cn2(z) over the height z. Based on the seeing measurements accumulated over 13 years we created ensemble models of multi-layer neural networks under the machine learning framework, including training and validation. For the first time in the world, we have simulated optical turbulence (seeing variations) during night-time with deep neural networks trained on a 13-year database of astronomical seeing. A set of neural networks for simulations of night-time seeing variations have been obtained. For these neural networks the linear correlation coefficient ranges from 0.48 to 0.68. We show that modeled seeing with neural networks is well described through meteorological parameters, which include wind speed components, air temperature, humidity and turbulent surface stresses. One of the fundamental new results is that the structure of small-scale (optical) turbulence over the Maidanak Astronomical Observatory does not depend or depends negligibly on the vortex component of atmospheric flows.

This paper attempts to clarify an long-standing issue about the number of unknowns in the Reynolds-Averaged Navier-Stokes equations (RANS). This study shows that all perspectives regarding the numbers of unknowns in the RANS stem from the misinterpretation of the Reynolds stress tensor. The current literature consider that the Reynolds stress tensor has six unknown components; however, this study shows that the Reynolds stress tensor actually has only three unknown components, namely the three components of fluctuation velocity. This understanding might shed a light to understand the well-known closure problem of turbulence.

This paper attempts to clarify an issue regarding the lasting unsolved problem of turbulence, namely the closure problem. This study shows that all perspectives regarding the numbers of un- known quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor. The current literatures have a consensus that the Reynolds stress tensor has six unknowns; however, this study shows that the Reynolds stress tensor actually has only three ones, namely the three components of fluctuation velocity. With this new understanding, the closed turbulence equations for incompressible flows are proposed.

This paper shown that turbulence closure problem is not an issue at all. All mistakes in the literature regarding the numbers of unknown quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor, i.e., all literatures have stated that the symmetric Reynolds stress tensor has six unknowns; however, it actually has only three unknowns, i.e., the three components of fluctuation velocity. We shown the integral-differential equations of the Reynolds mean and fluctuation equations have exactly eight equations, which equal to the numbers of quantities in total, namely, three components of mean velocity, three components of fluctuation velocity, one mean pressure and one fluctuation pressure. That is why we claim in this paper, that the Reynolds Navier-Stokes turbulence equations of incompressible flow are closed rather than unclosed. This study may help to solve the puzzle that has eluded scientists and mathematicians for centuries.

The relations between the kinetic energy spectrum and the second order longitudinal structure function in two dimensions are derived, and several examples are considered. The forward conversion (from spectrum to structure function) is illustrated first with idealized power law spectra, representing turbulent inertial ranges. The forward conversion is also applied to the zonal kinetic energy spectrum of Nastrom and Gage (1985) and the result agrees well with the longitudinal structure function of Lindborg (1999). The inverse conversion (from structure function to spectrum) is tested with data from 2D turbulence simulations. When applied to the theoretical structure function (derived from the forward conversion of the spectrum), the result closely resembles the original spectrum, except at the largest wavenumbers. However the inverse conversion is much less successful when applied to the structure function obtained from pairs of particles in the flow. This is because the inverse conversion favors large pair separations, which are typically noisy with particle data. Fitting the structure function to a polynomial improves the result, but not sufficiently to distinguish the correct inertial range dependencies. Furthermore the inversion of non-local spectra is largely unsuccessful. Thus it appears that focusing on structure functions with Lagrangian data is preferable to estimating spectra.

Magnetic helicity plays a tremendously important role when it is different from zero on average. Most notably, it leads to the phenomenon of an inverse cascade. Here, we consider decaying magnetohydrodynamic turbulence as well as some less common examples of magnetic evolution under the Hall effect and ambipolar diffusion, as well as cases in which the magnetic field evolution is constrained by the presence of an asymmetry in the number density of chiral fermions, whose spin is systematically either aligned or anti-aligned with its momentum. In all those cases, there is a new conserved quantity: the Hosking integral. We present quantitative scaling results for the magnetic integral scale as well as the magnetic energy density and its spectrum. We also compare with cases were also a magnetic version of the Saffman integral is initially finite.

Capabilities of hovering unmanned aerial vehicles (UAVs) in low-altitude sensing of atmospheric turbulence with high spatial resolution are studied experimentally. The experiment was carried out at the Basic Experimental Observatory of the V.E. Zuev Institute of Atmospheric Optics SB RAS. UAV findings are compared with objective data on the state of the atmosphere that were measured by acoustic anemometers installed at weather towers at altitudes of 4, 10, and 27 m. Profiles of atmospheric turbulence were recorded with three UAVs hovering at the same heights at a distance of 5 m from the anemometers. The behavior of the longitudinal and lateral wind velocity components is studied in the frequency band 0-10 Hz before and after one-minute smoothing, which reveals long turbulent wind gusts. The discrepancy between the UAV and anemometer data is analyzed. It is found that after smoothing, the discrepancy does not exceed 0.5 m/s in 95% cases. This value meets the requirements of the World Meteorological Organization for some applications in wind velocity measurements. The spectral and correlation analysis of the UAV and anemometer measurements is carried out. The profiles of the longitudinal and lateral scales of turbulence are examined based on the UAV data and are found to be in agreement with the anemometer data to a good accuracy.

Atmospheric turbulence, which produces chaotic motions in the planetary boundary layer, can inhibit mixing between fast-reacting species produced or released at different locations. This segregation process modifies the effective rate at which reactions occur between these species, and are not appropriately accounted for in coarse resolution models since these models assume complete mixing of tracers within each gridbox. Here we present a few examples of LES-based simulations applied to chemically reactive species in a forested area with high emissions of biogenic hydrocarbons, an urban area rich in anthropogenic emissions; and a maritime area with high emissions of reduced sulfur species.

This study revisits the Reynolds-averaged Navier--Stokes equations (RANS) and finds that the existing literature is erroneous regarding the primary unknowns and the number of independent unknowns in the RANS. The literature claims that the Reynolds stress tensor has six independent unknowns, but in fact the six unknowns can be reduced to three that are functions of the three velocity fluctuation components, because the Reynolds stress tensor is simply an integration of a second-order dyadic tensor of flow velocity fluctuations rather than a general symmetric tensor. This difficult situation is resolved by returning to the time of Reynolds in 1895 and revisiting Reynolds' averaging formulation of turbulence. The study of turbulence modeling could focus on the velocity fluctuations instead of on the Reynolds stress. An advantage of modeling the velocity fluctuations is, from both physical and experimental perspectives, that the velocity fluctuation components are observable whereas the Reynolds stress tensor is not.

This paper proposed an explicit and simple representation of velocity fluctuation and the Reynolds stress tensor in terms of the mean velocity field. The proposed turbulence equations are closed. The proposed formulations reveal that the mean vorticity is the key source of producing turbulence. It is found that there are no velocity fluctuation and turbulence if there were no vorticity. As a natural consequence, the laminar- turbulence transition condition was obtained in a rational way.

Based on author's previous work [Sun, B. The Reynolds Navier-Stokes Turbulence Equations of Incompressible Flow Are Closed Rather Than Unclosed. Preprints 2018, 2018060461 (doi: 10.20944/preprints201806.0461.v1)], this paper proposed an explicit representation of velocity fluctuation and formulated the Reynolds stress tensor in terms of the mean velocity field. The proposed closed Reynolds Navier-Stokes turbulence formulations reveal that the mean vorticity is the key source of producing turbulence.

We consider the developed turbulence of capillary waves on shallow water. Analytic theory shows that an isotropic cascade spectrum is unstable with the respect to small angular perturbations, in particular, to spontaneous breakdown of the reflection symmetry and generation of nonzero momentum. By computer modeling we show that indeed a random pumping, generating on average zero momentum, produces turbulence with a nonzero total momentum. A strongly anisotropic large-scale pumping produces turbulence whose degree of anisotropy decreases along a cascade. It tends to saturation in the inertial interval and then further decreases in the dissipation interval. Surprisingly, neither the direction of the total momentum nor the direction of the compensated spectrum anisotropy is locked by our square box preferred directions (side or diagonal) but fluctuate.

In this study, an experimentally validated computational model was developed to investigate the hydrodynamics in a rotor-stator vortex RVR agglomeration reactor having a rotating disc at the centre with two shrouded outer plates. A numerical simulation was performed using a simplified form of the reactor geometry to compute the 3D flow field in batch mode operations. Thereafter, the model was validated using data from a 2D Particle Image Velocimetry (PIV) flow analysis performed during the design of the reactor. Using different operating speeds—70, 90, 110 and 130 rpm, the flow fields were computed numerically followed by a comprehensive data analysis. The simulation results showed separated boundary layers on the rotating disc and the stator. The flow field within the reactor is characterized by a rotational plane circular forced vortex flow in which the streamlines are concentric circles with a rotational vortex. Overall, the results of the numerical simulation demonstrate a fairly good agreement between the CFD model and the experimental data as well as the available theoretical predictions. The swirl ratio β was found to be approximately 0.4044, 0.4038, 0.4044 and 0.4043 for operating speeds of N=70, 90, 110 and 130 rpm respectively. In terms of the spatial distribution, the turbulence intensity and kinetic energy are concentrated on the outer region of the reactor while the axial velocity showed a decreasing intensity towards the shroud. However, a comparison of the CFD and experimental predictions of the tangential velocity and the vorticity amplitude profiles shows that these parameters were under-predicted by the experimental analysis which could be attributed to some of the experimental limitations rather than the robustness of the CFD model or numerical code.

The present work is a numerical simulation of velocity and mixture fraction fields in a turbulent non-reaction Propane jet flow issuing into parallel co-flowing air, in isothermal condition which has been experimentally described here: http://www.sandia.gov/TNF/DataArch/ProJet.html. The objective is a better understanding of the flow structure and mixing process, a situation in which there is no chemical interaction and heat-transfer. The two-equation Realizable k-ε eddy viscosity turbulence model has been used to simulate the turbulent flow field on a 2D plane (i.e., on a 5-degree sector of the experimental domain), because Realizable k-ε more accurately predicts the spreading rate of both planar and round jets and presents the best proficiency in comparison with all versions of the k-ε models. Afterward, axial and radial profiles of mean velocities, turbulence energy, mean mixture fraction, the mixture fraction half radius (Lf), and the mass flux diagram have been numerically elicited for grid independent mesh (mesh B) and compared with corresponding experimental data to assess the numerical model. To obtain turbulence kinetic energy, k, due to lack of w' in the experimental data, the assumption of w'=v' seems to be valid. Simulations have demonstrated that mean mixture fraction (at radial profiles at locations x/D=0, 4, 15, 30 and 50) and its half radius, Lf, which characterizes jet width expanse, are prominently well-captured; Moreover, mean vertical velocity fields (in X-direction: Umean) have revealed less accuracy but still conspicuously well-captured. However, mean vertical velocity fields (in Y-direction: Vmean) have disclosed less resolution; Likewise, turbulence kinetic energy, k, have manifested moderate accuracy (between Umean and Vmean). It should be noticed that although numerical results for absolute pressure, p, have been obtained on aforementioned sections, there were no experimental data to compare with them. Thus, the corresponding numerical data has not been demonstrated in this study.

Turbulent dissipation is a central issue in the star and galaxy formation process. Its fundamental property of space-time intermittency, well characterised in incompressible laboratory experiments, remains elusive in cosmic turbulence. Progress requires the combination of state-of-the-art modelling, numerical simulations and observations. The power of such a combination is illustrated here where the statistical method intended to locate extrema of velocity shears in a turbulent field is applied to numerical simulations of compressible magneto-hydrodynamical (MHD) turbulence dedicated to dissipation scales and to a nearby turbulent diffuse molecular cloud. We demonstrate that short-spacing increments of observables can detect strongly dissipative structures. In our simulations, we compute structure functions of various synthetic observables and show that they verify Extended Self-Similarity. This allows to compute their intermittency exponents and we show how they could help constraining some properties of the underlying three-dimensional turbulence. In observations of a turbulent cloud close to the Sun in our Galaxy, a remarkable coherent structure of velocity shear extremum is disclosed . At the location of the largest velocity shear, this coherent structure is found to foster 10 × more carbon monoxide molecules than standard diffuse molecular gas, an enrichment supported by models of non-equilibrium warm chemistry triggered by turbulent dissipation. The power of the combination between modelling and observations is also illustrated by observations of the CH+ cation that provide unique quantitative informations on the kinetic energy trail in the massive, multi-phase and turbulent circum-galactic medium of a galaxy group at redshift z = 2.8.

This paper presents the results of simulations of particle-laden air-solid jet flow in long straight tubes using CFD-DEM, along with an analysis of coarse-graining. Although previous studies have used CFD-DEM for similar flows, these have typically been in a dilute flow regime where uncoupled simulations are accurate. However, fully coupled simulations can introduce issues, necessitating validation studies to ensure that all coupling parameters are effectively used and that the physics is accurately represented. This paper validates the simulations against two different experimental studies, with fluid Reynolds numbers between 10,000 and 40,000 and Stokes numbers between 5.6 and 50. Interestingly, the profiles of mean particle velocity exhibit fewer discrepancies as the Stokes number increases but more discrepancies for the root mean square velocity compared to the experiments. Particle number flux is consistent with the experiments after the nozzle exit. Coarse-graining is also applied to the same simulations, achieving relatively accurate results. However, as expected, scaling of contact collision frequencies, forces, and stresses cannot be achieved, meaning that coarse-graining may be useful for comparing designs or operating parameters on an industrial scale but falls short when measuring the total energy dissipation of one experiment.

This brief review introduces some general properties of lower atmospheric motions that may be exploited by birds, linking historical findings and current knowledge from both observational and modeling perspectives. The relative importance of turbulent mixing resulting from wind shear and that from buoyant thermals is emphasized as well as the determinants of the geometry of the resulting vertical motions.

This paper presents an analysis of linear viscous stress Favre-Averaged turbulence models for computational fluid dynamics (CFD) of fully turbulent round jets with a long straight tube geometry in the near field. Although similar work has been performed in the past with very relevant solutions, considerations were not given for the issues and limitations involved with coupling between an Eulerian and Lagrangian phase, such as in fully two-way coupled CFD-DEM. These issues include limitations on solution domain, mesh cell size, wall modelling, and momentum coupling between the two phases in relation to the particles size. Therefore, within these considerations, solutions are provided to the Navier-Stokes equations with various turbulence models using a three-dimensional wedge long straight tube geometry for fully developed turbulence flow. Simulations are performed with a Reynolds number of 15000 and 50000 using two different tube diameters. It is found that a modified k−ε turbulence model achieved the most agreeable results for both the velocity and turbulent flow fields between these two flow regimes, while a modified k−ω SST/BSL also provided suitable results.

Sea fog event over the Eastern Yellow Sea on 15–16 April 2012 was reproduced in the Weather Research and Forecasting (WRF) simulation with high-resolution to investigate the roles of phys-ical processes and synoptic-scale flows on advection fog with sea surface warming. Initially, longwave radiative cooling (LRC) with negative sensible heat flux (SHF) due to the turbulence af-ter sunset triggered a formation of cloud at the surface under the moist advection with a southerly wind. This is a conventional type of advection fog. At night, continuous cooling due to longwave radiation and SHF near the surface modulated the change of the SHF from negative to positive, resulting in a drastic increase in the latent heat flux (LHF) that provided sufficient moisture at lower atmosphere (self-moistening). This is a favorable condition for advection fog with sea sur-face heating (ssH), and this transition represents advection fog with ssH. Enhanced turbulent mixing driven by a buoyancy force increased the depth of the sea fog with a gradual rise in the marine atmospheric boundary layer (MABL) height, even at nighttime. In addition, cold advec-tion with a prevailing northerly wind at the top of the MABL led to a drastic increase in turbulent mixing and the MABL height, which resulted in rapid growth of the height of sea fog due to ver-tical diffusion. After sunrise, shortwave radiative warming in the fog layers offsetting the LRC near the surface weakened thermal instability, which contributed to the reduction in the MABL height, even during the daytime. In addition, dry advection of northerly wind induced dissipa-tion of the fog via evaporation. An additional sensitivity test of sea surface salinity showed weaker and shallower sea fog than the control due to the decrease in both the LHF and local self-moistening.

The role of energy and entropy in decomposition of turbulent velocity flow-fields is to be shown in the paper. The decomposition methods based on energy concept are taken into account, namely the Proper Orthogonal Decomposition (POD) and its extension Bi-Orthogonal Decomposition (BOD). Entropy motivated view on the decomposed modes could offer new possibilities in the modes physical interpretation and in Reduced Order Modelling (ROM) strategy efficiency evaluation.

A method for solving the turbulence equations embedded in the sigma ocean general ocean circulation model is proposed. Like the general circulation model, the turbulence equations are solved using the splitting method by physical processes. The turbulence equations are split into two main stages describing transport-diffusion and generation-dissipation processes. Parameterization of turbulence in the framework of equations allows, at the generation-dissipation stage, to use both numerical and analytical solutions and to ensure high efficiency of the algorithm. The results of large-scale ocean dynamics simulation taking into account the parameterization of vertical turbulent exchange are considered. Numerical experiments were carried out using k-omega turbulence model embedded to the Institute of Numerical Mathematics Ocean general circulation Model (INMOM). Both the circulation and turbulence models are solved using the splitting method with respect to physical processes. The coupled model is used to simulate the hydrophysical fields of the North Atlantic and Arctic Oceans for 1948--2009. The model has a horizontal resolution of 0.25 degree and 40 sigma-levels along the vertical. The sensitivity of the solution to the changes in mixing parameterization is studied. Experiments demonstrate that taking into account the climatic annual mean buoyancy frequency improves the reproduction of large-scale ocean characteristics. There is a positive effect of Prandtl number variations for reproducing the upper mixed layer depth. The experiments also demonstrate the computational effectiveness of the proposed approach in solving the turbulence equations.

In this paper, we studied and analyzed the behavior of velocity profiles and the distribution of pressures in a pipe connected to an accessory. In this case, it was a 90º elbow fitting in which water flows as an incompressible flow. This study was developed numerically through the fluent program of ANSYS. It is of practical interest in engineering because it is based on results regarding the internal structure of a model of incompressible flows. In the same way, it provides a basis for the proper design and implementation of velocity profiles within an elbow.In the agricultural industry, suitable systems can be implemented for crop irrigation in this way with the results obtained to improve the design in the fluid distribution system.We used a feasible κ-ε turbulence model because this model causes less turbulence than the standard κ-ε model. A mesh independence study was developed, in which it was observed that combining different densities of meshes for the same geometry causes some variables, so a fine meshing was used to carry out the study so that the results are as close to reality as possible.

The study of the change in entropy during the formation of stars in collapsing molecular clouds is an important area of research in astrophysics. As the cloud decreases in size, the presence of gases leads to an increase in entropy. However, the extent of this change in entropy is influenced by the amount of space available for the gases. In this scientific investigation, I conduct a detailed analysis of entropy on a small scale within the contracting cloud, following the principles of the second law of thermodynamics. Initially, my focus is on a turbulent cloud composed of N particles with a mass of M, which collapses under the force of gravity. It becomes clear that in order for the collapse to occur, gravitational energy must overcome other pressure forces that oppose it. This process results in the rapid movement of particles within the cloud, ultimately leading to a change in entropy

Subgrid scale models for unresolved turbulent fluxes are investigated with a focus on combustion for space propulsion applications. An extension to the gradient model is proposed, introducing a dependency on the local burning regime. The dynamic behavior of the model's coefficients is investigated, and scaling laws are studied. The discussed models are validated using a DNS database of a high-pressure turbulent fuel-rich methane-oxygen diffusion flame. The operating point and turbulence characteristics are selected to resemble those of modern combustors for space propulsion applications to support the future usage of the devised model in this context.

In this paper, for the Reynolds-averaged Navier-Stokes equations, a self-closed turbulence model without any adjustable parameter is formulated and a simple model is proposed. The validation rule for self-closed turbulence model is rigorously derived from the Reynolds-averaged Navier-Stokes equation. The rule is not effected by turbulence modelling on the Reynolds stresses. As an application, two dimensional turbulent boundary layers, and plate turbulent boundary layers have been formulated.

Localized turbulence is common in geophysical flows, where the roles of rotation and stratification are paramount. In this study, we investigate the evolution of a stratified turbulent cloud under rotation. Recognizing that a turbulent cloud is composed of vortices of varying scales and shapes, we start our investigation with a single eddy. Compared to an eddy under pure rotation, the stratified eddy shows the physical manifestation of a known potential vorticity mode, appearing as a static stable vortex. In addition, the expected shift from inertial waves to inertial-gravity waves is observed. In our numerical simulations of the turbulent cloud, carried out at a constant Rossby number over a range of Froude numbers, stratification causes columnar structures to deviate from vertical alignment. This deviation increases with increasing stratification, slowing the expansion rate of the cloud. The observed characteristics of these columnar structures are consistent with the predictions of linear theory, particularly in their tilt angles and vertical growth rates, suggesting a significant influence of inertial-gravity waves. Using Lagrangian particle tracking, we have identified regions where wave activity dominates over turbulence. In scenarios of milder stratification, these inertial-gravity waves are responsible for a significant energy transfer away from the turbulent cloud, a phenomenon that attenuates with increasing stratification.

An aneurysm is a vascular malformation, it can be classified according to its location (cerebral, aortic); and according its shape (saccular, fusiform, and mycotic). The study of blood flow interaction with the aneurysms has gained attention from physicists an engineers in recent years. Shear stresses, oscillatory shear index (OSI), gradient oscillatory number (GON), and residence time had been used as potential variables to describe the hemodynamics as well as the origin and evolution of aneurysm. However, the causes of aneurysms and the hemodynamics conditions that promote the growth of these structures are still under debate. The present work presents numerical simulation on three aneurysms: 2 aortic and 1 cerebral. The results shows that the blood rheology is not relevant for aortic aneurysms, however, for cerebral aneurysm could play a relevant role in the hemodynamics. The turbulence models evaluated show equivalent results in both cases. Lastly, a simulation considering the fluid-structure interaction (FSI) shows that this phenomenon is the dominant factor for aneurysms simulation.

Tropical cyclones have long been known to be powered by turbulent enthalpy fluxes from the ocean surface and retarded by turbulent momentum fluxes into the surface. Here were review evidence that the development and structure of these storms are also partially controlled by turbulence in the outflow near the storm top. Finally, we present new research that shows that tropical cyclone-like, low-aspect-ratio vortices are most likely in systems in which the bottom heat flux is controlled by mechanical turbulence and the top boundary is insulating.

A multiscale modeling concept for numerical simulation of multiphysics turbulent flow utilizing map-based advection is described. The approach is outlined with emphasis on its theoretical foundations and physical interpretations in order to establish the context for subsequent presentation of the associated numerical algorithms and the results of validation studies. The model formulation is a synthesis of existing methods, modified and extended in order to obtain a qualitatively new capability. The salient feature of the approach is that time advancement of the flow is fully resolved both spatially and temporally, albeit with modeled advancement processes restricted to one spatial dimension. This one-dimensional advancement is the basis of a bottom-up modeling approach in which three-dimensional space is discretized into under-resolved mesh cells, each of which contains an instantiation of the modeled one-dimensional advancement. Filtering is done only to provide inputs to a pressure correction that enforces continuity and to obtain mesh-scale-filtered outputs if desired. The one-dimensional advancement, the pressure correction, and coupling of one-dimensional instantiations using a Lagrangian implementation of mesh-resolved volume fluxes is sufficient to advance the three-dimensional flow without time advancing coarse-grained equations, a feature that motivates the designation of the approach as autonomous microscale evolution (AME). In this sense, the one-dimensional treatment is not a closure because there are no unclosed terms to evaluate. However, the approach is additionally suitable for use as a subgrid-scale closure of existing large-eddy-simulation methods. The potential capabilities and limitations of both of these implementations of the approach are assessed conceptually and with reference to demonstrated capabilities of related methods.

With the advent of submesoscale O(1km) permitting basin-scale ocean simulations, the seasonality in the mesoscale O(50km) eddies with kinetic energies peaking in summer has been commonly attributed to submesoscale eddies feeding back onto the mesoscale via an inverse energy cascade under the constraint of stratification and Earth’s rotation. In contrast, by running a 101-member, seasonally forced, three-layer quasi-geostrophic (QG) ensemble configured to represent an idealized double-gyre system of the subtropical and subpolar basin, we find that the mesoscale kinetic energy shows a seasonality consistent with the summer peak without resolving the submesoscales; by definition, a QG model only resolves small Rossby number dynamics (O(Ro)≪1) while as submesoscale dynamics are associated with O(Ro)∼1. Here, by quantifying the Lorenz cycle of the mean and eddy energy, defined as the ensemble mean and fluctuations about the mean respectively, we propose a different mechanism from the inverse energy cascade by which the stabilization and strengthening of the western-boundary current during summer due to increased stratification leads to a shedding of stronger mesoscale eddies from the separated jet. Conversely, the opposite occurs during the winter; the separated jet destablizes and results in overall lower mean and eddy kinetic energies despite the domain being more susceptible to baroclinic instability from weaker stratification.

The perspective suggests an unorthodox interpretation of cancer dynamics in terms of strange attractors and fluid turbulence. The relevance of strange attractors in tumor growth-invasion dynamics and chemical pattern formation are reviewed. Chemical turbulence, fractals and chaos are discussed as the principal characteristics of tumor ecosystems and their complex cybernetics. To conclude, cancers are presented as complex systems exhibiting chaotic gene expression dynamics.

We have developed an unsteady and non-linear wind synopsis simulator called RIAM-COMPACT (Research Institute for Applied Mechanics, Kyushu University, COMputational Prediction of Airflow over Complex Terrain) in order to simulate the airflow on a microscale, i.e., a few tens of km or less. In RIAM-COMPACT, the large-eddy simulation (LES) has been adopted for turbulence modeling. LES is a technique in which the structures of relatively large eddies are directly simulated and smaller eddies are modeled using a sub-grid scale model. In the present study, we have conducted the numerical wind diagnoses for Taikoyama Wind Farm nacelle separation accident in Japan. The simulation results suggest that all six wind turbines at the Taikoyama Wind Farm are subject to significant influence from separated flow (terrain-induced turbulence) which is generated due to the topographic irregularities in the vicinity of the wind turbines. A proposal has been also made on reconstruction of the wind farm.

1) Background: the budget of non-thermal energy in galaxy clusters is not well constrained, owing to the observational and theoretical difficulties in studying these diluted plasmas on large scales. 2) Method: we use recent cosmological simulations with complex physics in order to connect the emergence of non-thermal energy to the underlying evolution of gas and dark matter. 3) Results: the impact of non-thermal energy (e.g. cosmic rays, magnetic fields and turbulent motions) is found to increase in the outer region of galaxy clusters. Within numerical and theoretical uncertainties, turbulent motions dominate the budget of non-thermal energy in most of the cosmic volume. 4) Conclusion: assessing the distribution non-thermal energy in galaxy clusters is crucial to perform high-precision cosmology in the future. Constraining the level of non-thermal energy in cluster outskirts will improve our understanding of the acceleration of relativistic particles by cosmic shocks and of the origin of extragalactic magnetic fields.

Over the last few years there has been much effort put into the development and validation of stochastic models of the trajectories of swarming insects. These models typically assume that the positions and velocities of swarming insects can be represented by continuous jointly Markovian processes. These models are first-order autoregressive processes. In more sophisticated models, second-order autoregressive processes the positions, velocities and accelerations of swarming insects are collectively Markovian. Although it is mathematically conceivable that this hierarchy of stochastic models could be extended to higher orders, here I show that such a procedure would not be well-based biologically because some terms in these models represent processes that have the potential destabilize insect flight dynamics. This prediction is supported by an analysis of pre-existing data for laboratory swarms of the non-biting midge Chironomus riparius. I suggest that the Reynolds number is a finely tuned property of swarming, as swarms may disintegrate at both sufficiently low and sufficiently high Reynolds numbers.

The shear stress generated by wall turbulence is the main cause of wall friction resistance in turbulent flow through pipes. This paper investigates the impact of inserting rest periods (regions of constant Reynolds number) within the pulsating operating cycle of velocity on the fully turbulent flow at large time-averaged Reynolds numbers, using the Large Eddy Simulation (LES) method. The study aims to explore the effect of increasing rest periods within pulsations on drag resistance. The dimensionless shear stress and drag reduction rates during different time periods of rest were analyzed. Numerical simulation results indicate that the pulsating velocity operation mode does not necessarily lead to drag reduction; it may even result in increased resistance. Inserting rest periods within the pulsation cycle can achieve drag reduction effects, with the maximum drag reduction rate reaching 21.8\%. Comparisons with experimental data and Direct Numerical Simulation (DNS) from the literature validate the feasibility of using the LES method for pipe pulsating operation modes.

Laser-induced fluorescence spectroscopy (LIFS) is actively used for remote sensing of atmospheric aerosols and currently is one of the most sensitive and selective techniques for determining small concentrations of substances inside particles. The use of high-power femtosecond laser sources for LIFS-based remote sensing of aerosol contributes to the development of new-generation fluorescence atmospheric lidars since it makes it possible overcoming the energy threshold for the nonlinear-optical effects of multiphoton absorption in particles and receiving the emission signal at long distances in the atmosphere. Our study is aimed to the development and experimental demonstration the technique of nonlinear laser-induced fluorescence spectroscopy (NLIFS) based on the remote excitation of aerosol fluorescent emission stimulated by a spatially-structured high-power femtosecond laser pulse. Importantly, for the first time to our knowledge, we demonstrate the advances of using stochastically-structured plasma-free intense light channels (postfilaments) specially formed by propagation of femtosecond laser radiation through a turbulent air layer for improving NLIFS efficiency. A multiple increase in the received signal of two-photon-excited fluorescence of polydisperse dyed aqueous aerosol by the structured postfilaments is reported.

An analysis of the features of measurements and correction of phase distortions in optical waves propagating in the atmosphere at various levels of turbulence was performed. It is shown that with increasing intensity fluctuations, the limiting capabilities of phase correction decrease, and the phase of an optical wave that has passed through a turbulence layer consists of two components: potential and vortex. It was found that even in the region of weak fluctuations there is an overlap of spectral filtering functions for intensity and phase fluctuations. Areas of turbulence inhomogeneities have been identified that will have mutual influence and negatively affect the operation of the phase meter. It is noted that correlation functions, both phase and intensity, are less susceptible to this compared to structural functions. The results of experimental studies of the reconstruction of the wavefront of laser radiation distorted by atmospheric turbulence using a Shack-Hartmann wavefront sensor during vignetting and central screening of the entrance pupil of the optical system are presented. Studies have been carried out on the propagation of laser radiation along a horizontal atmospheric path for various levels of turbulence. The results are analyzed in terms of Zernike polynomials.

This paper draws upon the theory of turbulence in physics to explore the similarities and differences between volatility in financial markets and turbulent phenomena at a statistical physics level. By analogizing the dynamics of financial markets with fluid turbulence, we have developed an innovative analytical framework aimed at deepening the understanding of the complexity inherent in financial markets. Our research methodology includes a comparative analysis of simulated time series data for the Standard & Poor's 500 Index and fully developed turbulent velocity time series. We focus on the similarities in key statistical physical properties, including probability distributions, correlation structures, and power spectral densities. Furthermore, we have established a model of capital flow in financial markets and proposed corresponding solutions. Through computational simulations and data analysis, we find that while financial market volatility shares some statistical characteristics with turbulent phenomena, significant differences exist in the shapes of probability distributions and the temporal scales of correlations. This suggests that although financial markets exhibit patterns akin to turbulence, the behavior of this complex, multi-variable driven system does not completely correspond to natural turbulence, highlighting the limitations of directly applying turbulence theory to financial market analysis. Our study offers a new perspective for the theories of financial markets and the field of econophysics, providing fresh insights into the complexity of financial markets and contributing to the prevention and management of financial risk.

In situ aircraft observations in typhoons have been scarce. This paper documents and analyses the aircraft and dropsonde data collected in Super Typhoon Saola (2023) over the northern part of the South China Sea. The wind and turbulence structures of the typhoon are investigated. The turbulence intensities are quantified in terms of turbulent kinetic energy (TKE) and eddy dissipate rate (EDR), and these data are compared with other available estimates of turbulence intensities, such as those based on weather radars, meteorological satellites, and numerical weather predic-tion (NWP) models. It is found that the TKE and EDR are closely correlated, and they are con-sistent with the weather radar/satellite observations as well as NWP-based outputs. Furthermore, the boundary layer inflow, vertical wind profiles, and atmospheric stability are analysed based on the dropsonde observations. The analysed results would advance the understanding of typhoon structures and offer references for the validation of remote sensing observations and NWP models.

Renewable energy research becomes increasingly necessary as it grows in importance. Wind turbines are an excellent method of harnessing wind energy, which is one of the most important renewable energy sources. In light of the importance of studying wind turbine performance, this study describes the background and principles of wind turbine operation. Computational fluid dynamics (CFD) was used to evaluate the hydrodynamic performance of two types of offshore and onshore wind turbines. Mechanical power and thrust force are calculated using the CFD model for both types of turbines in a steady state based on the number of blades and rotors. The results of three-dimensional simulations are compared with those of experimental testing. Based on the method and assumptions used, the offshore turbine's mechanical torque showed an average deviation of about 4%. Increased wind speed increased mechanical power and thrust, while turbulence caused irregular pressure distributions on surfaces as a result of turbulence intensity. Furthermore, as the free wind speed increases, the blade tip and hub speed will also increase, expanding the vortices created by that flow.

Quantum turbulence is characterized by many degrees of freedom interacting non-linearly to produce disordered states, both in space and time. The advances in trapping, cooling, and tuning the interparticle interactions in atomic Bose-Einstein condensates (BECs) make them excellent candidates for studying quantum turbulence. In this work, we investigate the decaying regime of quantum turbulence in a trapped BEC. Although much progress has been made in understanding quantum turbulence, other strategies are needed to overcome some intrinsic difficulties. We present an alternative way of investigating this phenomenon by defining and computing a characteristic length scale, which possesses relevant characteristics to study the establishment of the quantum turbulent regime. One intrinsic difficulty related to these systems is that absorption images of BECs are projected to a plane, thus eliminating some of the information present in the original momentum distribution. We overcome this difficulty by exploring the symmetry of the cloud, which allows us to reconstruct the three-dimensional momentum distributions with the inverse Abel transform. We present our analysis with both the two- and three-dimensional momentum distributions, discussing their similarities and differences. We argue that the characteristic length allows us to visualize the time evolution of the turbulent state intuitively.

The mean flow and turbulence statistics of the flow through a simplified urban environment, which is an active research area in order to improve the knowledge of turbulent flow in cities, is investigated. This is useful for civil engineering, pedestrian comfort and for health concerns caused by pollutant spreading. In this work, we provide analysis of the turbulence statistics obtained from well-resolved large-eddy simulations (LES). A detailed analysis of this database reveals the impact of the geometry of the urban array on the flow characteristics and provides for a good description of the turbulent features of the flow within a simplified urban environment. The most prominent features of this complex flow include coherent vortical structures such as the so-called arch vortex, the horseshoe vortex and the roof vortex. These structures of the flow have been identified by an analysis of the turbulence statistics. The influence of the geometry of the urban environment (and particularly the street width and the building height) on the overall flow behavior have also been studied. Finally, the well-resolved LES results were compared with the experimental database from Monnier et al. to discuss differences and similarities between the respective urban configurations.

With improved measurement and modelling technology, variability has emerged as an essential feature in non-equilibrium processes. While traditionally, mean values and variance have been heavily used, they are not appropriate in describing extreme events where a significant deviation from mean values often occurs. Furthermore, stationary Probability Density Functions (PDFs) miss crucial information about the dynamics associated with variability. It is thus critical to go beyond a traditional approach and deal with time-dependent PDFs. Here, we consider atmospheric data from the Whole Atmosphere Community Climate Model (WACCM) and calculate time-dependent PDFs and the information length from these PDFs, which is the total number of statistically different states that a system passes through in time. Time-dependent PDFs are shown to be non-Gaussian in general, and the information length calculated from these PDFs shed us a new perspective of understanding variabilities, correlation among different variables and regions. Specifically, we calculate time-dependent PDFs and information length and show that the information length tends to increase with the altitude albeit in a complex form. This tendency is more robust for flows/shears than temperature. Also, much similarity among flows and shears in the information length is found in comparison with the temperature. This means a stronger correlation among flows/shears because of a strong coupling through gravity waves in this particular WACCM model. We also find the increase of the information length with the latitude and interesting hemispheric asymmetry for flows/shears/temperature, a stronger anti-correlation (correlation) between flows/shears and temperature at a higher (low) latitude. These results also suggest the importance of high latitude/altitude in the information budge in the Earth's atmosphere, the spatial gradient of the information as a useful proxy for the transport of physical quantities.

Monin-Obukhov similarity theory (MOST) overestimates wind shear in some atmospheric stable conditions, i.e. Richardson number $R_f<0.25$. The overestimated wind shear that leads to an under-predicted friction wind speed and a lower ambient turbulence intensity for a given hub-height reference wind speed and a given roughness length, could influence wake modeling of a wind turbine. This work investigates the side effects of the breakdown of MOST on wake modeling under stable conditions and makes some modifications to the flow similarity functions to eliminate these side effects. Based on a field measurement in a wind farm, we firstly show that MOST predicts a larger wind shear for the atmospheric stability parameter $\zeta>0.1$ and proposes new flow similarity functions without constraining $R_f$ to limit the overestimated wind shear by MOST. Next, different turbulence models based on MOST and a modified one based on the new similarity functions are investigated through numerical simulations. These turbulence models are combined with the actuator disk model (AD) and Reynolds-averaged Navier–Stokes equations (RANS) to model wind turbine wakes under stable conditions. As compared to measurements, numerical results show that turbulence models based on MOST result in larger wake deficits and slower wake recovery rate with a square root of the mean-squared-error (RSME) of wake deficit in the range of 0.07-0.18. This overestimated wake effect is improved by applying the new similarity functions and the RSME of wake deficit is averagely reduced by 0.05. Finally, we check the role of the under-predicted turbulence intensity playing in the larger wake deficit predicted by models based MOST. Additional numerical simulations using the modified turbulence model are carried out, in which the roughness length is reduced to impose a hub-height ambient turbulence intensity equivalent to the MOST case. Simulation results show that reducing turbulence intensity enhances wake effects, however, it cannot reproduce the large wake deficit predicted by models based on MOST, which suggests that the overestimated wake effect by MOST could be also related to the overestimated wind shear.

The structure of turbulent flow over non-flat surfaces is a topic of major interest in practical applications in both engineering and geophysical settings. A lot of work has been done in the fully rough regime at high Reynolds numbers where the effect on the outer layer turbulence structure and the resulting friction drag is well documented. It turns out that surface topology plays a significant role on the flow drag especially in the transitional roughness regime and therefore, hard to characterize. Survey of literature shows that roughness function depends on the interaction of roughness height, flow Reynolds number and topology shape. In addition, if the surface topology contains large enough scales then it can impact the outer layer dynamics and in turn modulate the total frictional force. Therefore, it is important to understand the mechanisms underlying drag increase from systematically varied surface undulations in order to better interpret quantifications based on mean statistics such as roughness function. In this study, we explore the mechanisms that modulate the turbulence structure over a two-dimensional (2D) sinusoidal wavy surface with a fixed amplitude, but varying slope. To accomplish this, we model the turbulent flow between two infinitely wide 2D wavy plates at a friction Reynolds number, $Re_{\tau}=180$. We pursue two different but related flavors of analysis. The first one adopts a roughness characterization flavor of such wavy surfaces. The second one focuses on understanding the non-equilibrium near surface turbulence structure and their impact on roughness characterization. Analysis of the different statistical quantifications show strong dependence on wave slope for the roughness function indicating drag increase due to enhanced turbulent stresses resulting from increased production of vertical velocity variance from the surface undulations.

The Boltzmann-Gibbs (BG) entropy has been used in a wide variety of problems for more than a century. It is well known that BG entropy is extensive, but for certain systems such as those dictated by long-range interactions, the entropy must be non-extensive. Tsallis entropy possesses non-extensive characteristics, which is parametrized by a variable q (q = 1 being the classic BG limit), but unless q is determined from microscopic dynamics, the model remains but a phenomenological tool. To this date very few examples have emerged in which q can be computed from first principles. This paper shows that the space plasma environment, which is governed by long-range collective electromagnetic interaction, represents a perfect example for which the q parameter can be computed from micro-physics. By taking the electron velocity distribution function measured in the heliospheric environment into account, and considering them to be in quasi equilibrium state with electrostatic turbulence known as the quasi-thermal noise, it is shown that the value corresponding to q = 9/13 = 0.6923 may be deduced. This prediction is verified against observation made by spacecraft, and it is shown to be in excellent agreement.

Modern radio spectrometers make measurement of polarized intensity as a function of Faraday depth possible. I investigate the effect of depolarization along a model line of sight. I model sightlines with two components informed by observations: a diffuse interstellar medium with a lognormal electron density distribution and a narrow, denser component simulating a spiral arm or H~{\sc ii} region, all with synchrotron-emitting gas mixed in. I then calculate the polarized intensity from 300-1800 MHz and calculate the resulting Faraday depth spectrum. The idealized synthetic observations show far more Faraday complexity than is observed in Global Magneto-Ionic Medium Survey observations. In a model with a very nearby H~{\sc ii} region observed at low frequencies, most of the effects of a ``depolarization wall'' are evident: the H~{\sc ii} region depolarizes background emission and less (but not zero) information from beyond the H~{\sc ii} region reaches the observer. In other cases, the effects are not so clear, as significant amounts of information reach the observer even through significant depolarization, and it is not clear that low-frequency observations sample largely different volumes of the interstellar medium than high-frequency observations. The observed Faraday depth can be randomized such that it does not always have any correlation with the true Faraday depth.

In the present study, field observation wind data from the time of the wind turbine blade damage accident on Shiratakiyama Wind Farm were analyzed in detail. In parallel, high-resolution LES turbulence simulations were performed in order to examine the model’s ability to numerically reproduce terrain-induced turbulence. The comparison of the observed and simulated time series (1 second average values) from a 10 minute period from the time of the accident led to the conclusion that the settings of the horizontal grid resolution and time increment are important to numerically reproduce the terrain-induced turbulence that caused the wind turbine blade damage accident on Shiratakiyama Wind Farm. A spectral analysis of the same set of observed and simulated data revealed that the simulated data reproduced the energy cascade of the actual terrain-induced turbulence well.

At the Atsumi Wind Farm in Aichi Prefecture, Japan, damage to wind turbines occurred frequently due to terrain-induced turbulence. In the present study, numerical analyses of terrain-induced turbulence were conducted by reproducing the topography in the vicinity of the wind turbine sites in high resolution and using RIAM-COMPACT natural terrain version, which is based on large eddy simulation (LES). The results of the diagnoses indicated that, in the case of south-easterly wind, terrain-induced turbulence is generated at a small terrain feature located upstream of Wind Turbine (WT) #2, which serves as the origin of the turbulence. At the Atsumi Wind Farm, a combination of the series of wind diagnoses and on-site operation experience led to a decision to adopt an "automatic shutdown program" for WTs #1 and #2. Here, "automatic shutdown program" refers to the automatic suspension of wind turbine operation upon the wind speed and direction meeting the conditions associated with significant effects of terrain-induced turbulence at a wind turbine site. The adoption of the "automatic shutdown program" has successfully resulted in a large reduction in the number of occurrences of wind turbine damage, thus, creating major positive economic effects. 1) a reduction in the repair costs by 9.322 million yen per year per wind turbine, 2) an increase in the availability factor by 8.05%, and 3) an increase in the capacity factor by 1.7%.

Computational fluid dynamics is used to study the impact of the support structure of a tidal turbine on performance and the downstream wake characteristics. A high-fidelity computational model of a dual rotor, contra-rotating tidal turbine in a large channel domain is presented, with turbulence modelled using large eddy simulation. Actuator lines represent the turbine blades, permitting the analysis of transient flow features and turbine diagnostics. The following four cases are considered: the flow in an unexploited, empty channel; flow in a channel containing the rotors; flow in a channel containing the support structure; and flow in a channel with both rotors and support structure. The results indicate that the support structure contributes significantly to the behaviour of the turbine and to turbulence levels downstream, even when the rotors are upstream. This implies that inclusion of the turbine structure, or some parametrisation thereof, is a prerequisite for the realistic prediction of turbine performance and reliability, particularly for array layouts where wake effects become significant.

The recent demand for reducing pollutant and CO2 emissions from internal combustion engines has created an urgent need for the development of ultra-lean burn engines that ensure combus-tion stability and reducing knocking tendency. One of the most promising methods is the Pre-chamber Spark Ignition (PCSI) system, where a jet of high-energy reactive gases, produced by pilot combustion in a pre-chamber, ignites the main combustion event in the cylinder. Given the intricate phenomena inherent to PCSI systems, conducting 3D CFD studies is imperative for a comprehensive analysis and optimal design. Furthermore, the detailed CFD model, combined with the calibrated 0-D/1-D model, is anticipated to yield a significant amount of new data that would be difficult to obtain through experimental approaches, making it essential for advancing our understanding and optimization of these systems. However, recently published papers re-port that some state-of-the-art models developed for conventional SI operation may not be pre-dictive under the challenging conditions of TJI systems, particularly under lean conditions. This review elucidates the reasons behind the widespread adoption of CFD in the optimal design of PCSI engines. It presents significant examples and delves into the potential and challenges of employing CFD not just as a predictive tool but also as a design instrument for enhancing PCSI engine performance.

As the demand for fossil fuels continues to rise and their availability becomes increasingly limited, it is crucial to explore innovative approaches to optimize fuel combustion in existing systems. Among various engines, diesel engines are heavily reliant on fluid motion within the engine cylinder, which plays a critical role in fuel-air mixing. This abstract focuses on the potential of dynamic mathematical model simulation for multi-compartment rotor compressed combustion engines, aiming to revolutionize power generation by enhancing fuel burning efficiency. By leveraging advanced simulation techniques, this research seeks to unlock new insights into fluid dynamics and optimize fuel-air mixing, leading to more efficient and sustainable energy utilization. Through this exploration, we aim to pave the way for groundbreaking advancements in the field of combustion engines, ultimately contributing to a greener and more sustainable future. In the realm of internal combustion engines, the fluid dynamics within the engine cylinder greatly influence the combustion processes and heat transfer. For compression ignition (C.I) engines, both the bulk gas motion and turbulence characteristics play vital roles in engine performance. This study focuses on enhancing the combustibility of the fuel-air mixture in C.I engines by modifying the engine design to induce turbulence through squish and tumble flows. The modification involves creating multiple compartments on the rotor crown, comprising three small chambers spaced 120º apart. By implementing this design, the aim is to optimize the fuel-air mixing and improve the overall combustion efficiency. Through comprehensive analysis and experimentation, we seek to revolutionize C.I engine technology and pave the way for more efficient and environmentally-friendly power generation. In this study, the performance of a compression ignition (C.I) engine with a multiple compartment rotor has been analyzed under motoring conditions using computational fluid dynamics (CFD) simulations with FLUENT software. The obtained results have been compared with those of a base C.I engine. The analysis focused on the effects of the modified engine design on the tumble ratio and squish velocity, which are important parameters for fuel-air mixing and combustion efficiency. The results indicated a significant improvement with the modified engine, with a 35% increase in tumble ratio and a 31% increase in squish velocity compared to the base engine. These findings suggest that the modified engine design has the potential to enhance fuel-air mixing and combustion performance, leading to improved overall engine efficiency. This research presents promising insights into the optimization of C.I engine design and contributes to the development of more efficient and environmentally-friendly power generation technologies..

The flow around thick airfoils at low Reynolds numbers and the possibility of using energy methods to control flowfield, increase lift and reduce drag are considered. For calculations, a series of propulsion airfoils is used, obtained by a combination of solving the inverse problem of aerodynamics and stochastic methods of global optimization. The calculations take into account the transition from the laminar flow regime to the turbulent one, which has a significant effect on the lift and airfoil drag. The applicability of various turbulence models for airfoil modelling at low Reynolds numbers is discussed. The results of numerical simulation make it possible to determine the dependence of the lift coefficient on the angle of attack and the optimal flow rate of air taken from the airfoil surface for each angle of attack. The results of numerical simulation are compared with the data of a physical experiment, and conclusions are drawn about the accuracy of various turbulence models.

The dynamics of the topological charge of a vortex optical beam propagating in a turbulent air with accounting for the cubic nonlinearity is theoretically considered. On a number of examples, we show that the optical beam self-focusing manifests itself ambiguously depending on optical wave power. At near-critical values of beam power, self-focusing leads to enhanced spatial localization of optical vortices and substantial suppression of vortex walk-off relative to the beam axis caused by air turbulence. However, with increasing of optical intensity the modulation instability imposed by cubic nonlinearity becomes significant and contributes jointly with medium turbulence to faster divergence of vortex beams.

This article examines the wavelet transforms and the machine learning methods in accelerating large eddy simulations of turbulent flows. An overarching challenge of large eddy simulations is to accurately represent the cascade of kinetic energy across the cut-off length scale. Taylor hypothesis suggests that the energy cascade occurs through the process of vortex stretching. However, Richardson hypothesis assumes the self-amplification of the strain field and a hierarchical break down of large eddies lead the energy cascade. The large eddy simulations typically employs the self-amplification of the strain in formulating subgrid models. However, several studies also proposed subgrid models based on vortex stretching. The wavelet-based coherent vortex simulation of turbulence directly accounts for vortex stretching in overall forward scatter of energy, while allowing local backscatter. The wavelet-based large eddy simulation adapts the grid to capture the creation of small-scale eddies, while adopting subgrid models based on the self-amplification of strain. The advancement of artificial intelligence in turbulence modelling is currently evolving around accelerating the numerical simulations of turbulent flow. However, there is a clear connection between the application of wavelet transfroms and neural networks for directly solving the Navier-Stokes equation, indicating some potential benefits of wavelet methods over the neural networks.

In large eddy simulation (LES) of turbulent flows, the most critical dynamical processes to be considered by dynamic subgrid models to account for an average cascade of kinetic energy from the largest to the smallest scales of the flow is not fully clear. Furthermore, evidence of vortex stretching being the primary mechanism of the cascade is not out of the question. In this article, we study some essential statistical characteristics of vortex stretching and its role in dynamic approaches of modeling subgrid-scale turbulence. We have compared the interaction of subgrid stresses with the filtered quantities among four models using invariants of the velocity gradient tensor. This technique is a single unified approach to studying a wide range of length scales in the turbulent flow. In addition, it also provides a rational basis for the statistical characteristics a subgrid model must serve in physical space to ensure an appropriate cascade of kinetic energy. Results indicate that the stretching mechanism extracts energy from the large-scale straining motion and passes it onto small-scale stretched vortices.

The adoption of Machine Learning (ML) for building emulators for complex physical processes has seen an exponential rise in the recent years. While neural networks are good function approximators, optimizing the hyper-parameters of the network to reach a global minimum is not trivial, and often needs human knowl- edge and expertise. In this light, automatic ML or autoML methods have gained large interest as they automate the process of network hyper-parameter tuning. In addition, Neural Architecture Search (NAS) has shown promising outcomes for improving model performance. While autoML methods have grown in popularity for image, text and other applications, their effectiveness for high-dimensional, complex scientific datasets remains to be investigated. In this work, a data driven emulator for turbulence closure terms in the context of Large Eddy Simulation (LES) models is trained using Artificial Neural Networks and an autoML frame- work based on Bayesian Optimization, incorporating priors to jointly optimize the hyper-parameters as well as conduct a full neural network architecture search to converge to a global minima, is proposed. Additionally, we compare the effect of using different network weight initializations and optimizers such as ADAM, SGDM and RMSProp, to explore the best performing setting. Weight and function space similarities during the optimization trajectory are investigated, and critical differences in the learning process evolution are noted and compared to theory. We observe ADAM optimizer and Glorot initialization consistently performs better, while RMSProp outperforms SGDM as the latter appears to have been stuck at a local minima. Therefore, this autoML BayesOpt framework provides a means to choose the best hyper-parameter settings for a given dataset.

The problem of simulating wakes in a stratified oceanic environment with active background turbulence is considered. Anisotropic RANS turbulence models are tested against laboratory and eddy-resolving models of the problem. An important aspect of our work is to acknowledge that the environment is not quiescent; therefore, additional sources are included in the models to provide a non-zero background turbulence. The RANS models are found to reproduce some key features from the eddy-resolving and laboratory descriptions of the problem. Tests using the freestream sources show the intuitive result that background turbulence causes more rapid wake growth and decay.

The application of naturally ventilated pig buildings (NVPBs) with outdoor exercise yards is on the rise mainly due to animal welfare considerations, while the issue of emissions from the buildings to the surrounding environment is important. Since air pollutants are mainly transported by airflow, the knowledge on the airflow characteristics downwind the building is required. The objective of this research was to investigate airflow properties downwind of a NVPB with a roofed outdoor exercise yard for roof slopes of 5°, 15°, and 25°. Air velocities downwind a 1:50 scaled NVPB model were measured using a Laser Doppler Anemometer in a large boundary layer wind tunnel. A region with reduced mean air velocities was found along the downwind side of the building with a distance up to 0.5 m (i.e. 3.8 times building height), in which the emission concentration might be high. It was found that a smaller roof slope (i.e. 5° slope) resulted in a higher and shorter wake zone and thus a shorter air pollutant dispersion distance. It was concluded that a smaller roof slope could contribute to the dilution of air pollutants and a lower air pollutant concentration near the ground.

Wind tunnels are devices that enable researchers to study the flow over objects of interest, the forces acting on them and their interaction with the flow, which is nowadays playing an increasingly important role due to noise pollution. Since the first closed circuit wind tunnel with variable cross-section was built in G¨ottingen, its Prandtl configuration has little change. The wind tunnel with Prandtl configuration has four corners and vanes, more than 50% of the total pressure loss are caused by the corners and vanes. How to reduce the total pressure loss is a world class problem in the wind tunnel design. This study attempts to propose a novel configuration of wind tunnel, where the corners have been replaced by semi-circular tunnel. Sun wind tunnel 2 has only two corners and vanes, while Sun wind tunnel 1 has no corners and vanes at all. It is expected the new wind tunnel can reduce the total pressure loss from 50% to 10%.

We review recent progress in high-speed orbital angular momentum (OAM) multiplexed free-space optical communication systems. The outdoor atmospheric turbulence is emulated by an indoor turbulence emulator, which is based on split-step beam propagation method. Adaptive optics, channel coding, Huffman coding combined with LDPC coding, and spatial offset are used for turbulence mitigation; while OAM multiplexing and wavelength-division multiplexing (WDM) are applied to boost aggregate capacity.

In this study, Fluent software 2020R2 is used to simulate vertical axis turbine in two dimensions, aiming to analyze the influence of turbulence intensity on wake characteristics. The detailed flow field data are obtained by simulating the operation of turbine under different turbulence intensities. On this basis, the influence factor of turbulence intensity are creatively introduced into the wake mathematical formula of vertical axis turbine established by Lam, and some parameter calculations of the formula are supplemented and perfected. The results show that the theoretical model can well reflect the influence of turbulence intensity on the wake distribution of vertical axis turbines. The formula can provide a theoretical basis for the design and optimization of vertical axis turbines and the array arrangement, and can be used to predict the wake characteristics under different turbulence intensity operating conditions.

The amount of heat transfer in gas turbine blades depends on cooling techniques and various flow phenomena. The effect of eddies, passing shock waves and free flow turbulence has been noticed by many researchers since the beginning. The focus of the upcoming work is on the numerical investigation of the effects of turbulence intensity and Reynolds number on the amount of heat transfer from the blade surface with internal cooling. The SST model has been used to solve this problem, which has been compared with the experimental work to ensure the correctness of the numerical simulation. In this regard, by changing the free flow turbulence intensity with values of 1, 5, 7, and 10% and at three Reynolds numbers of 150,000, 350,000, and 750,000, the changes in heat transfer coefficients have been investigated. The results show that with the increase of the turbulence intensity at different Reynolds numbers, due to the positive effect of the flow turbulence on the suppression of separation and promotion of the boundary layer transition, the blade surface heat transferincreases and this increase is more evident at higher Reynolds numbers.

We develop a quantitative microscopic theory of decaying Turbulence by studying the dimensional reduction of the Navier-Stokes loop equation for the velocity circulation. We have found an infinite dimensional manifold of solutions of the Navier-Stokes loop equation\cite{M93, M23PR} for the Wilson loop in decaying Turbulence in arbitrary dimension $d >2$. This family of solutions corresponds to a fractal curve in complex space $\mathbb C^d$, described by an algebraic equation between consecutive positions plus a nonlinear periodicity condition. We derive the constrained SDE for the evolution of the fractal curve at a fixed moment of physical time as a function of an auxiliary stochastic time. We expect this stochastic process to cover our fixed manifold of the solutions of the decaying Turbulence. The energy density of the fluid decays as $\mathcal E_0/t$, where $\mathcal E_0$ is an initial dissipation rate. Presumably, we have found a new phase of extreme Turbulence yet to be observed in real or numerical experiments.

The newly developed friction-interlocking armor unit, called ‘Double cube’ (DC), has been designed to improve the performance of the concrete armor unit used in river/marine structures against currents/waves actions. DC unit, is an octagonal-shaped block, made up of two parts: an upper cube set on a lower base that is either square or octagonal in shape. The innovative design aims to provide a good performance in terms of stability, high tolerance placement with various contact points and ease of placement. DC’s shape and placement enhance the stability by bringing the center of gravity closer to the underlayer and providing a large contact surfaces with surrounding blocks that reduce the chance of extraction and limit movements (rocking, lifting) via "keystone" effect. The characteristics of this new unit provides a relatively high hydraulic stability number for the armor layer (NS=2.9), a favorable hydraulic performance due to energy dissipation from turbulence and aeration, as well as a high roughness coefficient (γf≈0.46) helping to reduce the overtopping.

We present large eddy simulations of a mid-latitude open ocean front using a modified state of the art computational fluid dynamics code. We investigate the energy and information fluxes at the submesoscale/small-scale range in the absence of any atmospheric forcing. We find submesoscale conditions (Ro∼1, Ri∼1) near surface, within baroclinic structures related with partially imbalanced frontogenetic activity. Near surface, the simulations show a significant scale coupling on scales larger than ∼103(m). This is manifested as a strong direct energy cascade and intense mutual communication between scales, where the latter was evaluated using an estimator based on Mutual Information Theory. At scales smaller than ∼103(m) the results show near-zero energy flux, however, at this scale range, the estimator of mutual communication still shows values corresponding with a significant level of communication between them. This fact motivates to investigate the nature of the self-organized turbulent motion at this scale range with weak energetic coupling but where communication between scales is still significant and to inquire into the existence of synchronization or functional relationships between scales with emphasis on eventual underlying non-local processes.

The gauge field equation for fluid mechanics established in Part 1 is developed into a first-order scattering theory in the simplified case of a two-dimensional incompressible flow over a flat plate. This is used to present a model for the origin of Tollmien-Schlichting (TS) waves based on scattering between fluid particles. As a result, analytical formulae for the maximum amplification factor and the transition point from laminar to turbulent flow in the boundary layer are obtained. The mathematical transformations from the stationary field equations in Part 1 to a scattering theory with time evolution along the flow axis using Wick rotation are elaborated in detail.

In this review we summarize existing trends of flow control used to improve the aerodynamic efficiency of wings. We first discuss active methods to control turbulence, starting with flat-plate geometries and building towards the more complicated flow around wings. Then, we discuss active approaches to control separation, a crucial aspect towards achieving high aerodynamic efficiency. Furthermore, we highlight methods relying on turbulence simulation, and discuss various levels of modelling. Finally, we thoroughly revise data-driven methods, their application to flow control, and focus on deep reinforcement learning (DRL). We conclude that this methodology has the potential to discover novel control strategies in complex turbulent flows of aerodynamic relevance.

Quantum turbulence deals with the phenomenon of turbulence in quantum fluids, such as superfluid helium and trapped Bose-Einstein condensates (BECs). Although much progress has been made in understanding quantum turbulence, several fundamental questions remain to be answered. In this work, we investigated the entropy of a trapped BEC in several regimes, including equilibrium, small excitations, the onset of turbulence, and a turbulent state. We considered the time evolution when the system is perturbed and let to evolve after the external excitation is turned off. We derived an expression for the entropy consistent with the accessible experimental data, that is, using the assumption that the momentum distribution is well-known. We related the excitation amplitude to different stages of the perturbed system, and we found distinct features of the entropy in each of them. In particular, we observed a sudden increase in the entropy following the establishment of a particle cascade. We argue that entropy and related quantities can be used to investigate and characterize quantum turbulence.

Therapeutic proteins are used to successfully treat hemophilia, Crohn’s Disease, diabetes, and cancer. Recent product recalls have occurred because of sub-visible particle formation resulting from the inherent instability of proteins. It has been suggested that particle formation is associated with late stage processing steps of filling, shipping, and delivery. Previous works demonstrated that cavitation might occur in therapeutic vials subjected to agitation or accidentally dropped, but that mitigation can be achieved with fluid property manipulation. The goal of this research was to (1) assess the risk of cavitation under common pharmaceutical manufacturing conditions (i.e., pipe contraction and pumps), (2) establish a simple threshold criterion for when particulate will form, and (3) suggest a series of mitigation techniques based on these thresholds. To accomplish these tasks, computational fluid dynamic simulations for a variety of pipe contraction and fluid properties were performed. The results of this research show that reducing the turbulence in a fluid system will reduce the likelihood of cavitation. Additionally, threshold bounds were created that establish a definitive transition at which cavitation will occur.

Cancers remain the lead cause of disease-related, pediatric death in North America. The emerging field of complex systems has redefined cancer networks as a computational system with intractable algorithmic complexity. Herein, a tumor and its heterogeneous phenotypes are discussed as dynamical systems having multiple, strange attractors. Machine learning, network science and algorithmic information dynamics are discussed as current tools for cancer network reconstruction. Deep Learning architectures and computational fluid models are proposed for better forecasting gene expression patterns in cancer ecosystems. Cancer cell decision-making is investigated within the framework of complex systems and complexity theory.

Cancers are complex, adaptive ecosystems. They remain the leading cause of disease-related death among children in North America. As we approach computational oncology and Deep Learning Healthcare, our mathematical models of cancer dynamics must be revised. Recent findings support the perspective that cancer-microenvironment interactions consist of turbulent flows. As such, cancer pattern formation, protein-folding and metastatic invasion are discussed herein as processes driven by chemical turbulence within the framework of complex systems theory. Current state-of-the-art quantitative approaches used in reconstructing cancer stem cell networks are reviewed. To conclude, cancer stem cells are presented as strange attractors of the Waddington landscape.

Three-phase interactions (metal-slag-argon) in ladle stirring operations have strong effects on metal-slag mass transfer processes. Specifically, the thickness of the slag controls the fluid turbulence to an extent that once trespassing a critical thickness, increases of stirring strength have not further effects on the flow. To analyze these conditions, a physical model considering the three phases was built to study liquid turbulence in the proximities of the metal-slag interface. A velocity probe placed close to the interface permitted the continuous monitoring and statistical analyses of turbulence. The slag-eye opening was found to be strongly dependent on the stirring conditions, and the mixing times decreased with thin slag thicknesses. Slag entrainment was enhanced with thick slag layers, and high flow rates of the gas phase. A multiphase model was developed to simulate these results finding a good agreement between experimental and numerical results.

While laminar flow heat transfer and mixing in microfluidic geometries has been investigated experimentally, as has the effect of geometry-induced turbulence in microfluidic flow (it is well documented that turbulence increases convective heat transfer in macrofluidic flow), little literature exists investigating the effect of electrokinetically-induced turbulence on heat transfer at the micro scale. Using recently observed experimental data, this work employed computational fluid dynamics coupled with electromagnetic simulations to determine if electrokinetically-forced, low-Reynolds number turbulence could be observed in a rectangular microchannel with using Newtonian fluids. Analysis of the results was done via comparison to the experimental criteria defined for turbulent flow. This work shows that, even with a simplified simulation setup, computational fluid dynamics (CFD) software can produce results comparable to experimental observations of low-Reynolds turbulence in microchannels using Newtonian fluids. In addition to comparing simulated velocities and turbulent energies to experimental data this work also presents initial data on the effects of electrokinetic forcing on microfluidic flow based on entropy generation rates.

Most of the structures in flowing water are a challenge to their stability and sustainable with different flow conditions. Recent, renewable energy research and development covers ocean and river energy platform in which flow of water drag considered in various conversion devices towards the offshore and onshore establishment. Various energy platforms have been suggested for offshore development. However, the stability of these platforms in water is a serious concern. To study the water interaction over circular and square cross-section cascade system under the water has been carried out. Water flow around the pillars or column of the energy platform are analyzed through simulation software. Very low velocity 0.5 m/s has been considered to analyze the system. Total fifteen numbers of cascade pillars having circular and square cross-section area were considered. K-ε turbulence model is adopted to calculate the flow interaction to the column. A velocity, pressure, and energy fields are found around the column.

To address the issue of low prediction accuracy for sub-grid scale turbulence-radiation interaction (SGS-TRI) characteristics, which impacts the precision of remote infrared signal calculations for gas jet, a term for correlation of radiation-participating component fluctuations with temperature fluctuation was introduced into the current sub-grid Snegirev equation; Simultaneously, we established a new method that is more convenient in complex projects for solving the dynamic sub-grid model coefficients, e.g., dynamic Smagorinsky model and temperature variance model. The validity of these algorithms was confirmed through direct numerical simulations (DNS) for isotropic turbulence and test data of high-temperature air jet flow. The Multi-Scale Multi-Group Wide-Band (MSMGWB) k-distribution model was used to calculate radiative transfer in strongly non-isothermal and inhomogeneous media such as the remote infrared sensing of jet plume scenario considering influences of TRI and SGS-TRI. The study for axisymmetric jet plume revealed that the TRI accounts for more than 1/3 of total remote infrared signal in the 3-5micron band, with the SGS-TRI effect contributing more than 10%. Additionally, this influence was found to further increase with increasing detection distances.

Many designs for natural convection heat sinks, and semi-empirical correlations have been proposed in the recent years, but they are only valid in a limited range of Elenbaas number $El$ and were mostly tested for laminar flows. To alleviate those limits, parametric studies with 2D and quasi-3D models were carried out with ranges of Grashof numbers up to $1.55 \cdot 10^{11}$ and Elenbaas numbers up to $3.42 \cdot 10^7$. Ansys Fluent’s laminar, transition-SST, SST k-$\omega$, and k-$\epsilon$ models were applied. In addition, when used in this valid range, \textit{i.e.} mean Elenbaas numbers, with the simplified quasi-3D model, the transition-SST model can predict better results, over-estimating the heat flux by $10$ to $15\%$ compared to semi-empirical correlations. The 2D model was not deemed satisfying regarding turbulence models. Consequently, a quasi-3D model was developed: it appeared to be an efficient trade-off between computational time and prediction accuracy, in particular for turbulence models. New grouping factors were also found, to ensure proper dimensioning of natural convection heat sinks. They correspond to non-dimensional parameters that dictate the physical behaviour of the heat sink with respect to the semi-empirical correlations. Typically, the ratio of the spacing to the optimal spacing predicted by Bar-Cohen's correlation turned out to be an appropriate grouping factor with a threshold of $1$, above which the fins could safely be considered as isolated, thus greatly simplifying all further calculations.

The objective of this study, which is divided into two parts, is twofold: to address long standing challenges in the sensing of atmospheric turbulence-induced wavefront aberrations in strong scintillation conditions via comparative analysis of several basic scintillation resistant wavefront sensing (SR-WFS) architectures and iterative phase retrieval (IPR) techniques (Part I, this paper), and to develop a framework for the potential integration of SR-WFS techniques into practical closed-loop non-astronomical atmospheric adaptive optics (AO) systems (Part II). In this paper we consider basic SR-WFS mathematical models and phase retrieval algorithms, tradeoffs in sensor design and phase retrieval technique implementation, and methodologies for WFS parameter optimization and performance assessment. The analysis is based on wave-optics numerical simulations imitating realistic turbulence-induced phase aberrations and intensity scintillations, and optical field propagation inside SR-WFSs. Several potential issues important for practical implementation of SR-WFS and IPR techniques such as requirements for phase retrieval computational grid resolution, tolerance with respect to optical element misalignments, and the impact of camera noise and input light non-monochromaticity are also considered. Results demonstrate that major wavefront sensing requirements desirable for AO operation in strong intensity scintillations can potentially be achieved by transitioning to novel SR-WFS architectures based on iterative phase retrieval techniques.

In this study, the unsteady Reynolds-averaged Navier-Stokes (URANS) equations are employed in conjunction with the Menter Shear Stress Transport (SST)-Scale-Adaptive Simulation (SAS) turbulence model in compressible flow, unstructured mesh and complex geometry. Whereas, other scale-resolving approaches in space and time, such as Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES), supply more comprehensive information about the turbulent energy spectra of the fluctuating component of the flow; however, previously mentioned avenues imply computationally intensive situations, usually performed over structured meshes and relatively simply geometries. In contrast, the SAS approach is designed according to “physically” prescribed length scales of the flow. More precisely, it operates by locally comparing the length scale of the modelled turbulence to the von Karman length scale (which is prescribed according to local fluid velocity derivatives and depends on the grid point distribution in the computational domain). This length scale ratio allows the flow to dynamically adjust the local eddy viscosity in order to better capture the large scale motions (LSM) in unsteady regions of URANS simulations. While SAS may be constrained to model only low flow frequencies or waivenumbers (i.e., LSM), its versatility and computational low-cost make it attractive to obtain a quick first insight of the flow physics; particularly, in those situations dominated by strong flow unsteadiness. The selected numerical application is the well-known M219 three-dimensional rectangular acoustic cavity from the literature at two different freestream Mach numbers, M∞ (0.85 and 1.35) and a length-to-depth ratio of 5:1. Thus, we consider the “deep configuration" in experiments by Henshaw. The SST-SAS model has demonstrated a satisfactory compromise between simplicity, accuracy and flow physics description.

It is imperative to examine the performance of wind turbines both in steady and unsteady states due to their significant role in generating renewable energy by utilizing wind power. The computational fluid dynamics method (CFD) is one of the most efficient tools in this field, and the unsteady state has a great deal of importance when predicting turbine performance. Using the CFD method, the hydrodynamic performance of two types of offshore and onshore horizontal-axis wind turbines in the unsteady state was examined. Offshore wind turbines were modeled using the SST turbulence model and onshore wind turbines were modeled using the large eddy simulation (LES). Three-dimensional results of the simulation were compared and validated with experimental results, and good agreement between the two results indicates that the assumptions and method were accurate. Vortices were also studied in relation to the blade tip and their formation and development. Increasing the number of blades results in a reduction in torque, as well as a rise in vortices caused by the tips of the rotor blades. Due to the intensity of turbulence, the distribution of pressure on the surfaces will become irregular as the wind speed increases.

This paper summaries the main findings from the full-scale CFD analyses conducted at SINTEF Ocean on the case of MV REGAL, which is one of the benchmark vessels studied in the ongoing joint industry project JoRes. The numerical approach is described in detail, and comparative results are presented regarding the propeller open water characteristics, ship towing resistance, and ship self-propulsion performance. The focus of numerical investigations is on the assessment of the existing simulation best practices applied to a ship scale case in a blind simulation exercise, and performance of thereof with different turbulence modelling methods. The results are compared directly with full-scale performance predictions based on model tests conducted at SINTEF Ocean and sea trails data obtained in the JoRes project.

The unsteady flow around the monopile under the influence of wave and vibration is studied by numerical method. Based on VOF model and SST Turbulence model, combined with Intermittent Transition model and dynamic mesh method, wave generation and pile vibration in numerical flume are simulated. The flow field around vibrating monopile in a calm water flume and a wave flume is simulated and the simulation results are compared. The calculation results show that the monopile vibration has a great disturbance to the unsteady flow around the pile, and the disturbance degree intensifies with the increase of the vibration amplitude of the monopile: Under the action of inertia, the vibration of the monopile leads to the reattachment of the shedding vortices; Under the influence of the reattachment vortices, the asymmetry of the vorticity on the side wall of the monopile intensifies; The intensification of the disturbance causes the reduction of the phase lead angle of the bed shear stress on the side of the single pile, and an additional high bed shear stress appears on the side of the monopile. In addition, the vortex structure formed alternately around the monopile during the vibration process leads to the periodic change of the high bed shear stress.

Few studies have examined the role of mediating mechanisms played by bricolage in the relationship between strategic flexibility and new product development, even though most of the studies view strategic flexibility as a capability to cope with environmental turbulence and promote product innovation. The bricolage approach embraces the resource-based understanding of the product development process as a concrete process of recombining resources in the process and exploring the mechanisms behind the bricolage approach. This study is the first of the kind that has examined the empirical effect of environmental turbulence on new product development through serial mediation of strategic flexibility and bricolage. Based on our results, environmental turbulence positively and significantly affects new product development through the serial mediation of strategic flexibility and bricolage. Bricolage is one of the mechanisms by which strategic flexibility may play a part in supporting product innovation. Our contribution to strategic flexibility research comes from examining the effect of strategic flexibility on bricolage and new product development from the point that environmental turbulence compels the high-tech industries to remain strategically flexible for their sustainability and growth. Besides providing new evidence on the relationship between environmental turbulence and product innovation, this research also provides new insight into the economic opportunity of new product development in transition economies such as China.

Existing Reynolds Averaged Navier-Stokes based transition models do not accurately predict separation induced transition for low pressure turbines. Therefore, in this study, a novel framework based on computational fluids dynamics driven machine learning coupled with multi-expression and multi-objective optimization is explored to develop models which can improve the transition prediction for the T106A low pressure turbine at an isentropic exit Reynolds number of Re2is=100,000. Model formulations are proposed for the transfer and laminar eddy viscosity terms of the laminar kinetic energy transition model using seven non-dimensional pi groups. The multi-objective optimization approach makes use of cost functions based on the suction-side wall-shear stress and the pressure coefficient. A family of solutions is thus developed, whose performance is assessed using Pareto analysis and in terms of physical characteristics of separated-flow transition. Two models are found which bring the wall-shear stress profile in the separated region at least two times closer to the reference high-fidelity data than the baseline transition model. As these models are able to accurately predict the flow coming off the blade trailing edge, they are also able to significantly enhance the wake-mixing prediction over the baseline model.

The study of turbulent heat transfer in liquid metal flows has gained interest because of applications in several industrial fields. The common assumption of similarity between the dynamical and thermal turbulence, namely the Reynolds analogy, has been proven to be not valid for these fluids. Many methods have been proposed in order to overcome the difficulties encountered in a proper definition of the turbulent heat flux, such as global or local correlations for the turbulent Prandtl number or four parameter turbulence models. In this work we assess a four parameter logarithmic turbulence model for liquid metals based on RANS approach. Several simulation results considering fluids with Pr = 0.01 and Pr = 0.025 are reported in order to show the validity of this approach. The Kays turbulence model is also assessed and compared with integral heat transfer correlations for a wide range of Peclet numbers.

In this paper we develop from first principles a unique law pertaining to the flow of fluids through closed conduits. This law, which we call “Quinn’s Law”, may be described as follows: When fluids are forced to flow through closed conduits under the driving force of a pressure gradient, there is a linear relationship between the normalized dimensionless pressure gradient, P_Q, and the normalized dimensionless fluid current, C_Q. The relationship is expressed mathematically as: P_Q = k₁ +k₂C_Q. This linear relationship remains the same whether the conduit is filled with or devoid of solid obstacles. The law differentiates, however, between a packed and an empty conduit by virtue of the tortuosity of the fluid path, which is seamlessly accommodated within the normalization framework of the law itself. When movement of the fluid is very close to being at rest, i.e., very slow, this relationship has the unique minimum constant value of k₁, and as the fluid acceleration increases, it varies with a slope of k₂ as a function of normalized fluid current. Quinn’s Law is validated herein by applying it to the data from published classical studies of measured permeability in both packed and empty conduits, as well as to the data generated by home grown experiments performed in the author’s own laboratory.

The heat transfer and fluid mechanics analysis in agitated vessels is usually performed using the Dittus-Boelter correlation. In this study, we propose the analysis of an agitated vessel with a draft tube through the analysis of dimensionless numbers correlations. The heat transfer and fluid mechanics simulation was conducted using ANSYS Fluent software. The simulation used the Sliding Mesh approach and the SST k-ω turbulent model to capture the vessel's fluid flow and heat transfer phenomena. A grid independence study was conducted by varying the number of time steps to obtain an appropriate time step size for the simulation. The convective heat transfer coefficient, as indicated by the mean values of the Nusselt number, was calculated. The local values along the radial coordinate of the heat transfer surface were also determined. These results provide insights into the distribution of heat transfer rates within the vessel, which can help optimize process efficiency and equipment design. The simulation results were analyzed by varying the c and m coefficients in the Nusselt correlation presented as Nu= c Rem Pr1/3. Furthermore, correlations describing the mean Nusselt number at the bottom and wall of the vessel are presented and compared with existing literature, contributing to advancing knowledge in this field.

Optical communications are described and analyzed by shaped beams; we challenge the ‎effects of parameters that impact the profile of a Flat Top Gaussian (FTG) beam. When ‎the laser beam propagates throughout the atmosphere, it can be influenced by different ‎optical phenomena including scattering, absorption, and turbulence due to changes in the ‎scintillation index and the forms of intensity that are displayed in the source and receiver ‎planes. In this project, the FTG laser beam that propagates through a weak turbulent ‎region is numerically investigated using open-source software. This simulation will be ‎performed according to a mathematical model based on the split-step beam propagation ‎method. Intensity distributions at the source plane and the received average intensity in ‎atmospheric turbulence are calculated, and additional contour is in the transducer plane. ‎The scintillation index, structure constant, source size, and other parameters, are applied ‎in the Rytov method to quantify the weak turbulent model. Moreover, these parameters ‎are analyzed in near-field propagation. Also, the effects of the beam’s scintillation and ‎beam wander are determined. All results simulated are discussed and compared with the ‎TEM00 Gaussian beam. Finally, these results are compared to measurements in the ‎experimental part of the work.

Atmospheric turbulence characteristics are essential in determining the quality of astronomical images and implementing adaptive optics systems. In this study, the vertical distributions of optical turbulence at the Peak Terskol observatory (43.27472 ° N 42.50083 ° E, 3127 m a.s.l.) using the Era-5 re-analysis and scintillation measurements are investigated. For the closest reanalysis grid node to the observatory, vertical profiles of the structural constant of the air refractive index turbulent fluctuations Cn2 were obtained. The calculated Cn2(z) vertical profiles are compared with the vertical distribution of turbulence intensity obtained from tomographic measurements with Shack-Hartmann sensor. The atmospheric coherence length at the location of Terskol Peak was estimated. Using combination of atmospheric models and scheme paramaterization of turbulence, Cn2(z) profiles at Mt. Kurapdag were obtained. The values of atmospheric coherence length at the Peak Terskol are compared with estimated values of this length at the ten astronomical sites including Ali, Lenghu and Daocheng.

The aeroelastic response of lightweight low-speed aircraft with slender wings under extreme flow turbulence intensity is not well-understood. Experiments on a commmercial unmanned aerial vehicle (UAV) with a 3-meter wingspan and aspect ratio of 13.6 were performed in a large open-return wind tunnel with extreme flow turbulence intensity of T∞≈10%. The structural dynamics of the wing in the bending mode and the flow beneath the wing to capture the effect of aileron deflection was measured using laser displacement sensors and tomographic particle image velocimetry (PIV) techniques. The unsteady lift produced by the wing was also measured using a high-capacity load cell at an angle of attack, α of 2 degrees for three freestream velocities U∞ of 13.4 m/s, 17.9 m/s, and 26.8 m/s representing the UAV’s stall through cruise speed. It was found that high flow turbulence intensity with large integral length scales relative to the wing chord plays a dominant role in the large unsteady lift and wing displacements measured. The power spectral density (PSD) of wing structural vibration shows that flow shedding from the wing and the integral length scales have a significant impact on the overall power inherent in the bending vibration of the wing. Computation of vorticity iso-surfaces in the flow measurement volume surrounding the aileron reveal a striking observation; aileron deflection, δa of 10∘becomes less effective in producing additional spanwise vorticity, which is proportional to circulation and lift, at U∞ of 26.8 m/s because the freestream already has elevated levels of vorticity. A paradigm shift in design is suggested for light aircraft structures with slender wings operating in highly turbulent flow, which is to employ multiple control surfaces in order to respond to this flow and mitigate large bending or torsion displacements and the probability of structural failure.

We study the Kelvinons: monopole ring solutions to the Euler equations, regularized as the Burgers vortex in the viscous core. There is finite anomalous dissipation in the inviscid limit. However, in the anomalous Hamiltonian, some terms are growing as logarithms of Reynolds number; these terms come from the core of the Burgers vortex. In our theory, the turbulent multifractal phenomenon is similar to asymptotic freedom in QCD, with these logarithmic terms summed up by an RG equation. The small effective coupling does not imply small velocity; on the contrary, velocity is large compared to its fluctuations, which opens the way for a quantitative theory. In the leading order in the perturbation theory in this effective coupling constant, we compute running multifractal dimensions for high moments of velocity circulation, in good agreement with the data for quantum Turbulence and available data for classical Turbulence. The logarithmic dependence of fractal dimensions on the loop size comes from the running coupling in anomalous dimensions. This slow logarithmic drift of fractal dimensions would be barely observable at Reynolds numbers achievable at modern DNS.

The year 2022 was characterized by several major events that impacted the business world in Indonesia and even globally after the COVID-19 pandemic devastated all business sectors in Indonesia and the world. The current business environment is experiencing accelerating turbulence characterized by increasing intensity of business competition, rapid changes in the market, and industrial connections climate. The industrial connection climate of an institution can be measured by the pattern of relationships that occur between workers and management. Industrial connections are fundamental to the sustainability of an organization and are understood as an important factor in influencing institutional performance. This research aims to look at the impact of the existence of industrial connections instruments, turbulence in the business environment on the industrial connections climate. Besides, how the industrial connections climate relates to firm performance. The research design is an inferential quantitative, using sampling and snowball sampling methods with total sample are 406. The data collected were analyzed using Structural Equation Model using LISREL version 8.70. This research used a survey to build an SEM model that tested the effect of organizational turbulence and industrial connections instruments on organizational climate directly, and tested mediation on organizational performance. All related indicators have outer loading more than 0.5 and T-value more than 1.96. Therefore, the conclusion is that all measurements are valid. Composite Reliability (CR) and VE for all indicators are >0.7 and >0.5 respectively. This research shows that industrial assosiation climate has a positive effect on business turbulence and industrial connections instruments. However, industrial connections instruments show a greater impact on the assosiation climate. This is because industrial connections instruments are the media used in the implementation of industrial connections in organizations. This research also confirms that industrial assosiation climate is positively related to organizational performance. The result of the indirect influence (IE) of 0.20 with t value is 8.47 means that the climate of Industrial connections is able to mediate the influence of Industrial connections Facilities on Organizational Performance.

The capabilities of a quadcopter in the hover mode for low-altitude sensing of atmospheric turbulence with high spatial resolution in urban areas characterized by complex orography are investigated. The studies were carried out in different seasons (winter, spring, summer, and fall), and the quadcopter hovered in the immediate vicinity of ultrasonic weather stations. The DJI Phantom 4 Pro quadcopter and AMK-03 ultrasonic weather stations installed in different places of the studied territory were used in the experiment. The smoothing procedure was used to main regularities in the behavior of the longitudinal and lateral spectra of turbulence in the inertial and energy production ranges. The longitudinal and lateral turbulence scales were estimated by the least-square fit method with the von Karman model as a regression curve. It is shown that the turbulence spectra obtained with DJI Phantom 4 Pro and AMK-03 generally coincide with minor differences observed in the high-frequency region of the spectrum. In the inertial range, the behavior of the turbulence spectra shows that they obey the Kolmogorov-Obukhov “5/3” law. In the energy production range, the longitudinal and lateral turbulence scales and their ratio measured by DJI Phantom 4 Pro and AMK-03 agree to a good accuracy. Discrepancies in the data obtained with the quadcopter and the ultrasonic weather stations at the territory with complex orography are explained by the partial correlation of the wind velocity series at different measurement points and the influence of the inhomogeneous surface.

Turbulent flow is the first and fundamental physical phenomena to evaluate when optimising cost and reducing emissions from an Anode Baking Furnace (ABF). Gas flow patterns, velocity field, pressure drop, shear stress, and turbulent dissipation rate variables are the main operational parameters to be optimised, considering a specific geometry. Computational Fluid Dynamics (CFD) allows simulating physical phenomena using numerical methods with computer resources. In particular, the finite element method is one of the most used methods to solve the flow equations. This method requires a discretisation of the geometry of the ABF, called mesh. Hence, mesh is the main input to the finite element method. A suitable mesh for applying a discretisation method determines whether the problem can be simulated or not. Generating an appropriate mesh remains a challenge to perform accurate simulations. In this work, a comparison between meshes generated using two mesh generation tools is presented. Results of different study cases are included.

The local climate effects of utility-scale photovoltaic (PV) plants have also attracted more and more attention with the rapid development of PV industry. Turbulent fluxes and intensity characteristics for the PV plant and the adjacent reference site are investigated in the gobi area in Wujianqu in Xinjiang. Results indicate that near surface boundary layer is more unstable during the daytime while is more stable at the night in the PV plant than the reference site. The average roughness is 0.089m and 0.041m for the PV plant and the reference site respectively. The turbulence intensity in the vertical direction (I_w) of the PV plant is higher than the reference site especially during the daytime. The turbulent kinetic energy of the PV plant during the daytime is lower than the reference site.The momentum in the PV plant is higher than the reference site, especially during the daytime. Compared to the reference site, the PV plant has higher sensible heat flux and less latent heat flux. Site specific turbulent intensity relationships were developed for the PV plant and the reference site by analyzing wind components as a function of stability parameter (z/L). The turbulent components of wind followed 1/3 power law in the unstable conditions and stable conditions in the PV plant and the reference site.

Based on the theory of turbulence, equations are derived for estimations of turbulent fluctuations of the longitudinal and transverse components of the wind velocity during ideal hovering of a quadcopter in a turbulent atmosphere. We present the results of experiments which were carried out on the territory of the Geophysical Observatory of Institute of Monitoring of Climatic and Ecological Systems, Siberian Branch, Russian Academy of Sciences, located in Tomsk Akademgorodok on the territory with complex orography, in a parkland zone with buildings of research institutes and motorways. Time series of turbulent fluctuations of the longitudinal and transverse components of wind velocity fluctuations were received with the use of an automated weather station, and time series of estimates of these components, from data of a DJI Phantom 4 Pro quadcopter during hovering. According to the automated weather station data, anisotropy was observed in one experiment during measurements in the atmosphere, but this phenomenon was not observed in the other experiment: the fluctuation spectra of all components of wind speed fluctuations coincide. The spectra of fluctuations of the longitudinal and transverse wind velocity components based on the automated weather station data and UAV telemetry are presented. The fluctuation spectra of these components for the automated weather station data and quadcopter generally coincide. The behavior of the spectra coincides with the spectrum which corresponds to Kolmogorov–Obukhov “–5/3” law within the inertial range. The turbulent spectra of the wind velocity fluctuations obtained with the use of the automatic weather station and with the unmanned aerial vehicle differ in the high-frequency spectral region.

As it turns out, both isotope scaling and density limits are phenomena closely linked to fluid closure. The necessity to include ion viscosity arises for both phenomena. Thus, we have added ion viscosity to our model. The experimental isotope scaling has been successfully recovered in our fluid model through parameter scans. Although ion viscosity typically exerts a small effect, the density limit is manifested by increasing the density by approximately tenfold from the typical experimental density. In our case, this increase originates from the density in the Cyclone base case. Notably, these phenomena would not manifest with a gyro-Landau fluid closure. The isotope scaling is nullified by the addition of a gyro-Landau term, while the density limit results from permitting ion viscosity to become comparable to the gyro-Landau term. The mechanism of zonal flows, demonstrated analytically for the Dimits upshift, yields insights into the isotope scaling observed in experiments. In our approach, ion viscosity is introduced in place of the Landau fluid resonances found in some fluid models. This implies that the mechanism of isotope scaling operates at the level of fluid closure in connection with the generation of zonal flows. The strength of zonal flows in our model has been verified, particularly in connection with the successful simulation of the nonlinear Dimits shift. Consequently, a role is played by our approach in the temperature perturbation part of the Reynolds stress.

Data obtained by Direct Numerical Simulations (DNS) of pressure-driven turbulent channel flow are studied in the range 180 ≤ Reτ≤ 10,000. Reynolds number effects on the Mean Velocity Profile (MVP) and second order statistics are analyzed with a view of finding logarithmic behavior in the overlap region or even further from the wall, well in the boundary layer’s outer region. Values of Kármán constant for the MVPs and the Townsend-Perry constants for the streamwise and spanwise fluctuation variances are determined for the Reynolds numbers considered. A data-driven model of the MVP, proposed and validated for zero pressure gradient flow over a flat plate, is employed for pressure-driven channel flow by appropriately adjusting Coles’ strength of the wake function parameter, Π. There is excellent agreement between the analytic model predictions of MVP and the DNS-computed MVP as well as of the Reynolds shear stress profile. The skin friction coefficient Cf is calculated analytically. The agreement between the analytical model predictions and the DNS-based computed discrete values of Cf is excellent.

In addition to their importance in statistical thermodynamics, probabilistic entropy measurements are crucial for understanding and analyzing complex systems, with diverse applications in time series and one-dimensional profiles. However, extending these methods to two- and three-dimensional data still requires further development. In this study, we present a new method to classify spatiotemporal processes based on entropy measurements. To test and validate the method, we selected four classes of similar processes related to the evolution of random patterns: dynamic colored noises ((i)~white and (ii)~red); (iii)~weak turbulence from reaction-diffusion; (iv)~hydrodynamic fully developed turbulence, and (v)~plasma turbulence from MHD. Considering seven possible ways to measure entropy from a matrix, we present the method as a parameter space composed of the two best separating measures of the five selected classes. The results highlight better combined performance of Shannon Permutation Entropy ($S^p_H$) and a new approach based on Tsallis Spectral Permutation Entropy ($S^s_q$). Notably, our observations reveal the segregation of reaction terms in this $S^p_H \times S^s_q$ space, a result that identifies specific sectors for each class of dynamic process, and can be used to train machine learning models for automatic classification of complex spatiotemporal patterns.

This is the second paper in a cycle investigating the exact solution of loop equations in decaying turbulence. We perform numerical simulations of the Euler ensemble, suggested in the previous work, as a solution to the loop equations. We designed novel algorithms for simulation, which take a small amount of computer RAM so that only the CPU time grows linearly with the size of the system and its statistics. This algorithm allows us to simulate systems with billions of discontinuity points on our fractal curve, dual to the decaying turbulence. For the vorticity correlation function, we obtain \textbf{quantum fractal laws} with regime changes between universal discrete values of indexes. The traditional description of turbulence with fractal or multifractal scaling laws does not apply here. The measured conditional probabilities of fluctuating variables are smooth functions of the logarithm of scale, with statistical errors being negligible at our large number of random samples $T = 167,772,160$ with $N = 200,000,000$ points on a fractal curve. The quantum jumps arise from the analytical distribution of scaling variable $X = \frac{\cot^2(\pi p/q)}{q^2}$ where $\frac{p}{q}$ is a random fraction. This distribution is related to the totient summatory function and is a discontinuous function of $X$. In particular, the energy spectrum is a devil's staircase with tilted uneven steps with slopes between $-1$ and $-2.03$. The logarithmic derivative of energy decay $n(t)$ as a function of time is jumping down the stairs of universal levels between $1.50$ and $1.53$. The quantitative verification of this quantization would require more precise experimental data.