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.

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.

This paper brings to fruition recently found \cite{migdal2023exact} exact reduction of decaying turbulence in the Navier-Stokes equation in $3 + 1$ dimensions to a Number Theory problem of finding the statistical limit of the Euler ensemble. We reformulate the Euler ensemble as a Markov chain and show the equivalence of this formulation to the \textbf{quantum} statistical theory of $N$ fermions on a ring, with an external field related to the random fractions of $\pi$. We find the solution of this system in the turbulent limit $N\to \infty, \nu \to 0$ in terms of a complex trajectory (instanton) providing a saddle point to the path integral over the density of these fermions. This results in an analytic formula for the observable correlation function of vorticity in wavevector space. This is a full solution of decaying turbulence from the first principle without approximations or fitted parameters. The energy spectrum decays as $k^{-\frac{7}{2}}$. We compute resulting integrals in \Mathematica{} and present effective index $n(t) = -t\partial_t \log E $ for the energy decay as a function of time Fig.\ref{fig::NPlot}. The asymptotic value of the effective index in energy decay $n(\infty) = \frac{5}{4}$, but the universal function $n(t)$ is neither constant nor linear.

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.

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.

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..

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.

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the Kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the Kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously-emitted thermal fluctuations and the background turbulence leads to the formation of non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the Kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

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 review is a first attempt at bringing together various concepts from research on wall- and magnetically-bounded turbulent flows. Brief reviews of both fields are provided: The main similarities identified are coherent (turbulent) structures, flow generation and transport barriers. Examples are provided and discussed.

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.

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.

Computational rhinology is a specialized branch of biomechanics leveraging engineering techniques for mathematical modelling and simulation to complement the medical field of rhinology. Computational rhinology has already contributed significantly to advancing our understanding of the nasal function, including airflow patterns, mucosal cooling, particle deposition, and drug delivery, and is foreseen as a crucial element in e.g. the development of virtual surgery as a clinical, patient-specific decision support tool. The current paper delves into the field of computational rhinology from a nasal airflow perspective, highlighting the use of computational fluid dynamics to enhance diagnostics and treatment of breathing disorders. The paper consists of three distinct parts – an introduction to and review of the field of computational rhinology, a review of published literature on in vitro and in silico studies of nasal airflow, and the presentation and analysis of previously unpublished high fidelity CFD simulation data of in silico rhinomanometry. While the two first parts of the paper summarize the current status and challenges in the application of computational tools in rhinology, the last part addresses the gross disagreement commonly observed when comparing in silico and in vivo rhinomanometry results. It is concluded that this discrepancy cannot readily be explained by CFD model deficiencies caused by poor choice of turbulence model, insufficient spatial or temporal resolution, or neglecting transient effects. Hence, alternative explanations such as nasal cavity compliance or drag effects due to nasal hair should be investigated.

Unsupervised physics-informed deep learning can be used to solve computational physics problems by training neural networks to satisfy the underlying equations and boundary conditions without labeled data. Parameters such as network architecture and training method determine the training success. However, the best choice is unknown a priori as it is case specific. Here, we investigated network shapes, sizes, and types for unsupervised physics-informed deep learning of the two-dimensional Reynolds-averaged flow around cylinders. We trained mixed-variable networks and compared them to traditional models. Several network architectures with different shape factors and sizes were evaluated. The models were trained to solve the Reynolds-averaged Navier-Stokes equations incorporating Prandtl's mixing length turbulence model. No training data were deployed to train the models. The superiority of the mixed-variable approach was confirmed for the investigated high Reynolds number flow. The mixed-variable models were sensitive to the network shape. For the two cylinders, differently deep networks showed superior performance. The best fitting models were able to capture important flow phenomena such as stagnation regions, boundary layers, flow separation, and recirculation. We also encountered difficulties when predicting high Reynolds number flows without training data

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.

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.

This research provides a computationally efficient numerical tool for predicting 3-D turbulent flows over hydrofoils. It addresses a gap in existing numerical methods by extending the capabilities of the laminar VIScous Vorticity Equation (VISVE) solver to handle 3-D turbulent flow scenarios. The implementation integrates the k − ω SST model into the 3-D VISVE solver using the Finite Volume Method (FVM), extending its capability to handle turbulent flows. The turbulence model is also parallelized to improve computational efficiency. The application of the newly developed model to 3-D hydrofoil cases expands its scope and serves as preparatory work for future rotational propeller applications. Rigorous testing includes a convergence analysis concerning grid numbers and time step sizes, and validation against a Reynolds-Averaged Navier-Stokes (RANS) solver ensures the reliability and accuracy of the enhanced VISVE solver for turbulent flow simulations. The turbulent VISVE solver offers advantages such as reduced computational domain and costs through a vorticity-based approach. Turbulence concentration within boundary layers and free shear flows does not compromise the advantages of the small computational domain. The simplified meshing process, automatically generated by denoting the number of panels on the hydrofoil, makes the solver accessible for researchers and engineers.

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.

An approach based on the OpenFOAM library has been developed to solve a high-speed, multicomponent mixture of a reacting, compressible flow. This work presents comprehensive validation of the newly developed solver called compressibleCentralReactingFoam with different supersonic flows, including shocks, expansion waves, and turbulence-combustion interaction. The comparisons of the simulation results with experimental and computational data confirm the fidelity of this solver for problems involving multicomponent high-speed reactive flows. The gas dynamics of turbulence-chemistry interaction are modeled using a partially-stirred reactor formulation and provide promising results to better understand the complex physics involved in supersonic combustors. A time-scale analysis based on local Damköhler numbers reveals different regimes of turbulent combustion. In the core of the jet flow, the Damköhler number is relatively high indicating that the reaction time scale is smaller than the turbulent mixing time scale. This means that the combustion is controlled by turbulent mixing. In the shear layer, where the heat release rate and the scalar dissipation rate have the highest value, the flame is stabilized due to finite rate chemistry with small Damköhler numbers and a high fraction of fine structure. This solver allows three-dimensional gas dynamic simulation of high-speed multicomponent reactive flows relevant to practical combustion applications.

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.

The devastating effects of COVID-19 permeates virtually all aspects of life across geographic boundaries. As the lockdown is eased, and reinstated, across many parts of the world, attention is now shifting towards how to navigate the ‘now normal’ rather than the ‘next normal,’ of doing business with a view to achieving sustainable development. This study highlights some of the challenges and opportunities for micro, small, and medium enterprises (MSMEs) in times of crisis such as COVID-19. Conceptual by orientation, the study provides a broad overview of trends in crises times, drawing upon a range of qualitative accounts. The overarching message is the current direction of travel for MSMEs in sustaining their bottom-line, not after the crisis, but at this particular moment in time. The study highlights some implications for ‘pivoting,’ in order to cushion the impact of the pandemic on efforts to capacitate MSMEs in leveraging value co-creation with the consumers. The implications of the here and now are discussed.

The computational model that is able to estimate the temperature distribution inside a solid specimen during the film boiling phase of immersion quenching (water entry) process has been presented in paper, and it is based on the prescribed initial temperatures of the solid specimen and the liquid quenchant. In addition, the turbulence effects have to be considered using the assumed turbulence kinetic energy value, i.e., the “frozen turbulence” approach, that remains constant thorough the simulation. The studied material is nickel alloy, Inconel 600, for which the extensive experimental data is available. The work has been carried out using ANSYS Fluent computational fluid dynamics software and the methods for solution of Stefan problem by Eulerian two fluid VOF model.. A satisfactory agreement between the experimental and the calculated data has been achieved, yielding thereby the computationally obtained data that fits to a great extent to the prescribed error band of ± 10 % during the estimated film boiling phase of the immersion quenching process itself.

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.

Fiber coupling efficiency is one of the key technical indicators of atmospheric laser communication, which directly affects the communication system BER. According to the optical fiber coupling efficiency model, the optical fiber coupling model under atmospheric slant link is established, and the relationship between wavelength, zenith angle, receiving aperture, altitude and coupling efficiency is studied through numerical simulation. Based on the principle of adaptive optics turbulence correction, the relationship between fiber coupling efficiency and adaptive optics order is established. A fast method for estimating the coupling efficiency based on Strehl ratio is proposed. According to the principle of laser communication system, the semi-physical simulation experimental system of atmospheric slant-range channel coupling efficiency is constructed, and the coupling efficiency test and analysis are completed.

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.

RANS simulations have been broadly used to investigate turbulence in the oceans and atmosphere. Within these environments there are a multitude of tracers undergoing reactions (eg. phytoplankton growth, chemical reactions). The distribution of these reactive tracers is strongly influenced by turbulent mixing. With a 50 member ensemble of two dimensional Rayleigh-Taylor induced turbulent mixing, we show that the dynamics of a reactive tracer growing according to Fisher’s equation are poorly captured by the mean. A fluctuation-dependent sink introduced by Reynolds averaging Fisher’s equation transfers concentration from the mean to the fluctuations. We compare the dynamics of the reactive tracer with those of a passive tracer. The reaction causes the concentration of reactive tracer to increase thereby increasing Fickian diffusion and allowing the reactive tracer to diffuse into turbulent structures that the passive tracer cannot reach. A positive feedback between turbulent mixing and fluctuation growth is identified. We show that eddy viscosity and diffusivity parameterizations fail to capture the bulk trends of the system and identify a need for negative eddy diffusivities. One must therefore be cautious when interpreting RANS results for reactive tracers.

We have found an infinite dimensional manifold of exact solutions of the Navier-Stokes loop equation for the Wilson loop in decaying Turbulence in arbitrary dimension $d >2$. This solution family is equivalent to a fractal curve in complex space $\mathbb C^d$ with random steps parametrized by $N$ Ising variables $\sigma_i=\pm 1$, in addition to a rational number $\frac{p}{q}$ and an integer winding number $r$, related by $\sum \sigma_i = q r$. This equivalence provides a \textbf{dual} theory describing a strong turbulent phase of the Navier-Stokes flow in $\mathbb R_d$ space as a random geometry in a different space, like ADS/CFT correspondence in gauge theory. From a mathematical point of view, this theory implements a stochastic solution of the unforced Navier-Stokes equations. For a theoretical physicist, this is a \textbf{quantum} statistical system with integer-valued parameters, satisfying some number theory constraints. Its long-range interaction leads to critical phenomena when its size $N \rightarrow \infty$ or its chemical potential $\mu \rightarrow 0$. The system with fixed $N$ has different asymptotics at odd and even $N\rightarrow \infty$, but the limit $\mu \rightarrow 0$ is well defined. The energy dissipation rate is analytically calculated as a function of $\mu$ using methods of number theory. It grows as $\nu/\mu^2$ in the continuum limit $\mu \rightarrow 0$, leading to anomalous dissipation at $\mu \propto \sqrt{\nu} \to 0$. The same method is used to compute all the local vorticity distribution, which has no continuum limit but is renormalizable in the sense that infinities can be absorbed into the redefinition of the parameters. The small perturbation of the fixed manifold satisfies the linear equation we solved in a general form. This perturbation decays as $t^{-\lambda}$, where the anomalous dimensions $\lambda$ satisfy the spectral equation \eqref{spectralEq}. The spectrum becomes a continuum in the statistical limit $N \to \infty$, presumably leading to multifractal phenomena.

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.

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.

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.

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.

We present a mixing length based algebraic turbulence model calibrated to pipe flow; the main purpose of the model is to capture the increasing turbulence production-to-dissipation ratio observed in connection with the high Reynolds number transition region. The model includes the mixing length description by Gersten and Herwig which takes the observed variation of the von Kármán number with Reynolds number into account. Pipe wall roughness effects are included in the model. Results are presented for area-averaged (integral) quantities which can be used both as a self-contained model and as initial inlet boundary conditions for computational fluid dynamics simulations.

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.

The nonlocality (superdiffusion) of turbulence is expressed in the empiric Richardson t3 scaling law for the mean square of the mutual separation of a pair of particles in a fluid or gaseous medium. The development of the theory of nonlocality of various processes in physics and other sciences based on the concept of Lévy flights resulted in the idea of Shlesinger and colleagues about the possibility of describing the nonlocality of turbulence using a linear integro-differential equation with a slowly falling kernel. The close approach developed by us made it possible to establish the closeness of the superdiffusion parameter of plasma density fluctuations moving across a strong magnetic field in a tokamak to the Richardson law. In this paper, we show the possibility of a universal description of the characteristics of nonlocality of transfer in a stochastic medium (including turbulence of gases and fluids) using the Biberman-Holstein approach to the transfer of excitation of a medium by photons, generalized to take into account the finiteness of the velocity of excitation carriers. This approach enables us to propose a scaling that generalizes Richardson's t3 scaling law to the combined regime of Lévy flights and Lévy walks.

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.

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.

A rigid protruding surface, agitator and vortex generator often dissipate more energy than to produce a satisfied turbulence production. In contrast, a passively-oscillating flexible protruding surface has the ability in generating greater turbulence level. In the current study, a circular finite cylinder (cantilever) was used as the geometry of the rigid and protruding surface. They were varied through the aspect ratio (AR) and materials. Also, a local Reynolds number within the subcritical flow range (102<ReD<105) was considered. The results from the rigid protruding surface (finite cylinder) is served as a validation to the published results and benchmark to the improve-ment of the turbulence generated by the flexible protruding surface. It was mentioned in the pre-vious work that the increment of turbulence production by the flexible cylinder is speculated to be oscillation, this work investigates the validity of the speculation through manipulating the struc-tural stiffness of the flexible cylinder by employing different materials. The results have further demonstrated that the flexible cylinder is capable of generating a greater turbulence through the examination of turbulence intensity, turbulence production term and Reynolds stress. All the flex-ible cylinders that oscillate show increment in turbulence production but at a different percentage. The cylinders studied in this work ranged from the least structural stiffness (EVA), moderate (al-uminium) to the highest structural stiffness (carbon steel). Through studying the normalized am-plitude responses graph of the flexible cylinder, it is identified that the oscillating motion does in-deed contribute to increment. A further investigation into the results, it is found that the factor which governs the increment is the structural velocity instead of just oscillating motion alone.

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.

Non-normal mode stability analysis has long been employed to study linear transient characteristics of wall bounded shear flows. However, the application of this approach to investigate shear flows with stratified thermal conditions (STCF) has received limited attention. This paper examines thermal effects on Taylor-Couette flow. The results in this study indicate the possibility of a subcritical bifurcation in stratified Taylor-Couette flow.

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.

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.

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.

The DNS of film boiling requires strong computational resources that are difficult to obtain for a daily CFD use by the practitioners of the industrial R&amp;D experts. On the other hand, the film boiling experiments are associated with the usage of the expensive and highly sophisticated apparatus, and the research to this end is found to be relatively difficult due to high heat flow rates that are present in the process itself. When combined with a transient heat conduction in a solid, the problem becomes significantly difficult. Therefore, a novel method in computation of conjugate heat transfer during film boiling in a quiescent liquid has been proposed in this paper. The method relies on the solution of mass, momentum and energy conservation equations in a two-fluid framework, supplemented with the appropriate closures. Furthermore, the turbulent flow has been found as an important parameter in obtaining the accurate solution of the temperature field evolution in a solid specimen, via the proper modeling of turbulent kinetic energy (TKE) value, that has been imposed as a constant value, i.e., the frozen turbulence approach. It was found, in addition, that the appropriate TKE value can be obtained by use of Kelvin-Helmholtz instability theory in conjunction with the boundary layer theory. The obtained results show excellent agreement with the experimental data within the first 15 s of the experiment, i.e., the first ca. 10 % of the total duration of the film boiling mode of heat transfer. Furthermore, the heat transfer coefficient has matched the error bands prescribed by the authors of the paper that has presented the correlations, whilst the averaged values are far beyond this band, i.e., are slightly more than 30 % higher. The further inspection has revealed a measure of similarity between the computational result of the volume fraction field distribution and the experiment, thus confirming the capability of the method to obtain realistic interface evolution in time. The method has shown full capability for further pursuing the industrial-scale film boiling problems that involve turbulent flow and the conjugate heat transfer approach.

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.

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.

In this work, an analytical and numerical analysis of the transition to chaos in five nonlinear systems of ordinary and partial differential equations, which are models of the autocatalytic chemical processes and the numbers of interacting populations, is carried out. It is analytically and numerically shown that in all considered systems of equations, further complication of the dynamics of solutions and the transition to chemical and biological turbulence is carried out in full accordance with the universal Feigenbaum-Sharkovsky-Magnitskii bifurcation theory through subharmonic and homoclinic cascades of bifurcations of stable limit cycles. In this case, irregular (chaotic) attractors in all cases are exclusively singular attractors in the sense of the FShM theory.

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.

Physics-informed neural networks are a promising method to yield surrogate models of flow fields. We present a metamodeling technique for variable geometries based on physics-informed neural networks. The method was applied to the DU99W350 airfoil at a Reynolds number of 1×10⁵. The model predicted the Reynolds-averaged velocity and pressure field around the airfoil for arbitrary angles of attack between 10.0° and 17.5°. The model was trained with data from CFD simulations for a limited set of angles of attack. Additionally, satisfaction of the a priori known boundary conditions as well as the Reynolds-averaged Navier-Stokes equations were trained. A sensitivity analysis concerning the Reynolds number, the amount and distribution of training data, and the turbulence model was conducted showing the superiority of the pseudo-Reynolds stress method and the demand of labeled training data in the domain. The trained network was capable of predicting the developing flow separation on the suction surface and exhibited excellent agreement with CFD results even in the proximity to the wall for interpolations as well as extrapolations from the labeled data set.

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.

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.

In a power plant that uses seawater as a coolant, a debris filter is required to remove foreign substances from the seawater, and differential pressure leads to a decrease in the coolant flow rate, leading to a decrease in power generation efficiency. In this study, analysis was performed under different initial flow velocity conditions of 1.5m/s, 2.0m/s, and 2.5m/s for a debris filter used in a power plant to understand flow characteristics in CTCS. The flow and differential pressure characteristics of a CTCS with an installed debris filter were considered in a comparative analysis of velocity, pressure, and turbulent kinetic energy distribution. The results confirmed that a vortex was generated in the pipe with the debris filter, apparently due to the collision of the flow with the bracket of the debris filter. As the flow rate increased, the range of the vortex increased, causing a loss in flow.

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.

The fundamental challenge for the structural design of complicated shapes and lightweight structures such as Tensile Membrane Structures (TMS) is a reliable estimation of lateral loading such as wind load. Wind pressure coefficients can’t be correctly calculated by standards for a large range of intricate forms of structures. Computational Fluid Dynamics (CFD), a strong computational technique for evaluating wind pressure distribution on complicated geometry, can be used to determine wind-related parameters. The average wind related-distributions for the tensile membrane surface are validated with a few contributions from prior studies, which are important to validate with benchmarks in experimental wind tunnel experiments. For the double-curved surface, the results show good agreement between experimental tests and CFD simulations. In this study, CFD simulation was used to calculate the mean surface pressure coefficient (Cp) for generic shapes and tensile membrane forms using the steady Reynolds-Averaged Navier-Stokes (RANS) technique. Wind force coefficients and the effect of wind directions are also explored across a wide range of classical and tensile forms. The purpose of the current research is to look at the effect of wind assault angle on the Cp distribution for tensile membrane structures and general shapes.

A formula for the turbulent Prandtl number is derived in terms of local variables available from two-equations turbulence models. The derivation is a direct consequence of the expected square additivity of the molecular and turbulent parameters defining the effective viscosity and the effective conductivity. The formula does not degenerate and leads to a Kays like formulation if approximated. Its utility and its limitation are discussed.

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.

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 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.

This paper presents a novel theoretical and computational methodology for the generation of an onset turbulent field with prescribed properties in the numerical simulation of an arbitrary viscous flow. The methodology is based on the definition of a suitable distribution of volume force terms in the right–hand side of the Navier–Stokes equations. The distribution is represented by harmonic functions that are randomly variable in time and space. The intensity of the distribution is controlled by a simple PID strategy in order to obtain that the generated turbulent flow matches a prescribed turbulence intensity. A further condition is that a homogeneous isotropic flow is es- tablished downstream of the region where volume force terms are imposed. Although it is general, the proposed methodology is primarily intended for the computational modelling of hydrokinetic turbines in turbulent flows representative of tidal or riverine installations. A first numerical applica- tion is presented by considering the injection of homogeneous and isotropic turbulence with 16% intensity into a uniform unbounded flow. The analysis of statistical properties as auto-correlation, power spectral density, probability density functions, demonstrates that the generated flow tends to achieve satisfactory levels of stationarity and isotropy, whereas the simple control strategy used determines overestimated turbulent intensity levels.

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.

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.

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.

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.

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.

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.

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.

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 Finite size Lyapunov exponent (FSLE) has been used extensively since the late 1990’s to diagnose turbulent regimes from Lagrangian experiments and to detect Lagrangian coherent structures in geophysical flows and two-dimensional turbulence. Historically, the FSLE was defined in terms of its computational method rather than via a mathematical formulation, and the behavior of the FSLE in the turbulent inertial ranges is based primarily on scaling arguments. Here we propose an exact definition of the FSLE based on conditional averaging of the finite amplitude growth rate (FAGR) of the particle pair separation. With this new definition, we show that the FSLE is a close proxy for the inverse structural time, a concept introduced a decade before the FSLE. The (in)dependence of the FSLE on initial conditions is also discussed, as well as the links between the FAGR and other relevant Lagrangian metrics, such as the finite time Lyapunov exponent and the second order velocity structure function.

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.

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 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 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.

Spatiotemporal dynamics of the atmospheric kinetic energy and its components caused by the ordered and turbulent motions of air masses are estimated from minisodar measurements of three velocity vector components and their variances within the lowest 5–200 m layer of the atmosphere, with a particular emphasis on the turbulent kinetic energy. The layered structure of the total atmospheric kinetic energy has been established. From the diurnal hourly dynamics of the altitude profiles of the turbulent kinetic energy (TKE) retrieved from minisodar data, four layers are established by the character of the altitude TKE dependence, namely, the near-ground layer, the surface layer, the layer with a linear TKE increase, and the transitive layer above. In the first layer, the most significant changes of the ТКЕ were observed in the evening hours. In the second layer, no significant changes in the TKE values were observed. A linear increase in the TKE values with altitude was observed in the third layer. In the fourth layer, the TKE slightly increased with altitude and exhibited variations during the entire observation period. The altitudes of the upper boundaries of these layers depended on the time of day. The MKE values were much less than the corresponding TKE values, they did not exceed 50 m²/s². From two to four MKE layers were distinguished based on the character of its altitude dependence. The two-layer structures were observed in the evening and at night (under conditions of the stable atmospheric boundary layer). In the morning and daytime, the four-layer MKE structures with intermediate layers of linear increase and subsequent decrease in the MKE values were observed. Our estimates demonstrated that the ТКЕ contribution to the total atmospheric kinetic energy considerably (by a factor of 2.5–3) exceeded the corresponding МКЕ contribution.

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.

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.

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.

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 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.

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.

In line with Pomeau’s conjecture about the relevance of directed percolation (DP) to turbulence onset/decay in wall-bounded flows, we propose a minimal stochastic model dedicated to the interpretation of the spatially intermittent regimes observed in channel flow before its return to laminar flow. Numerical simulations show that a regime with bands obliquely drifting in two stream-wise symmetrical directions bifurcates into an asymmetrical regime, before ultimately decaying to laminar flow. The model is expressed in terms of a probabilistic cellular automaton evolving von Neumann neighbourhoods with probabilities educed from a close examination of simulation results. It implements band propagation and the two main local processes: longitudinal splitting involving bands with the same orientation, and transversal splitting giving birth to a daughter band with orientation opposite to that of its mother. The ultimate decay stage observed to display one-dimensional DP properties in a two-dimensional geometry is interpreted as resulting from the irrelevance of lateral spreading in the single-orientation regime. The model also reproduces the bifurcation restoring the symmetry upon variation of the probability attached to transversal splitting, which opens the way to a study of the critical properties of that bifurcation, in analogy with thermodynamic phase transitions.

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.

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.

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.

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.

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.

In the present paper, an aerodynamic investigation of a high-speed train is performed. In the first section of this article, a generic high-speed train against a turbulent flow is simulated, numerically. The Reynolds-Averaged Navier-Stokes (RANS) equations combined with the SST turbulence model are applied to solve incompressible turbulent flow around a high-speed train. Flow structure, velocity and pressure contours and streamlines at some typical wind directions are the most important results of this simulation. The maximum and minimum values are specified and discussed. Also, the pressure coefficient for some critical points on the train surface is evaluated. In the following, the wind direction influence the aerodynamic key parameters as drag, lift, and side forces at the mentioned wind directions are analyzed and compared. Moreover, the effects of velocity changes (50, 60, 70, 80 and 90 m/s) are estimated and compared on the above flow and aerodynamic parameters. In the second section of the paper, various data-driven methods including Gene Expression Programming (GEP), Gaussian Process Regression (GPR), and random forest (RF), are applied for predicting output parameters. So, drag, lift and side forces and also minimum and a maximum of pressure coefficients for mentioned wind directions and velocity are predicted and compared using statistical parameters. Obtained results indicated that RF in all coefficients of wind direction and most coefficients of free stream velocity provided the most accurate predictions. As a conclusion, RF may be recommended for the prediction of aerodynamic coefficients.

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.

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.

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.

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.

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.

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.

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.

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.

A practical wind farm controller for production maximisation based on coordinated control is presented. The farm controller emphasises computational efficiency without compromising accuracy. The controller combines Particle Swarm Optimisation (PSO) with a turbulence intensity based Jensen wake model (TI-JM) for exploiting the benefits of either curtailing upstream turbines using coefficient of power ($C_P$) or deflecting wakes by applying yaw-offsets for maximising net farm production. First, TI-JM is evaluated using convention control benchmarking WindPRO and real time SCADA data from three operating wind farms. Then the optimized strategies are evaluated using simulations based on TI-JM and PSO. The innovative control strategies can optimise a medium size wind farm, Lillgrund consisting of 48 wind turbines, requiring less than 50 seconds for a single simulation, increasing farm efficiency up to a maximum of 6% in full wake conditions.

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.

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%.