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

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.

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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

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

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.

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

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

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

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

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

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.

In this paper, for the Reynolds-averaged Navier-Stokes equations, a self-closed turbulence model without any adjustable parameter is formulated. 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.

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.

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

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

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.

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.

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.

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.

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

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

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

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

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

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

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

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

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

Wave steepness is presented as an extension and a valuable add-on to the conventional set of sea state parameters retrieved from satellite altimetry data. Following physical model based on recent advances of weak turbulence theory wave steepness is estimated from directly measured spatial gradient of wave height. In this way the method works with altimetry trajectories rather than with point-wise data. Moreover, in contrast to widely used parametric models this approach provides us with instantaneous values of wave steepness and period. Relevance of single-track estimates of wave steepness (period) is shown for wave climate studies and confirmed by a simple probabilistic model. The approach is verified via comparison against buoy and satellite data including crossover points for standard 1 second data of Ku-band altimeters. High quality of the physical model and robustness of the parametric ones are examined in terms of global wave statistics. Prospects and relevance of both approaches in the ocean wave climate studies are discussed.

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.

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

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

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

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

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

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

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

To gain insight into the differences between onshore and offshore atmospheric turbulence, 2 pressure fluctuations were measured for offshore wind under different environmental conditions. 3 A durable piezo-electric sensor was used to sample turbulent pressure data at 50 kHz. Offshore 4 measurements were performed at 100 m height on Germany’s FINO3 offshore platform in the 5 German Bight together with additional meteorological data provided by Deutscher Wetterdienst 6 (DWD). The statistical evaluation revealed that the stability state in the atmospheric boundary has a 7 large impact on turbulent fluctuations. Therefore, we used higher statistical properties (described 8 by so-called shape factors) to the stability state. Data was classified to be either within the unstable, 9 neutral or stable stratification. We found that in case of stable stratification, the shape factor is 10 mostly close to zero, indicating that a thermally stable environment produces closer-to Gaussian 11 distributions. Non-Gaussian distributions were found in unstable and neutral boundary layer states 12 and an occurrence probability was estimated. Possible impact on laminar-turbulent transition on the 13 blade is discussed with application of so-called laminar aerofoils on wind turbine blades. Use of a 14 cut-off frequency to separate load and aerodynamic turbulence is proposed.

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

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

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

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

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

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

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.

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.

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.

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.

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.

Based on the insight gained from single jet analysis performed earlier [1], CFD analysis on turbulent jets impinging on one another at an angle was performed. Multiple impingement angles were considered for this study to gain better understanding of the parameters affecting resultant jet growth and velocity distribution. From the study of single jet, it was concluded that the SST k-ω turbulence model was the ideal turbulence model capable of accurately predicting the flow physics of the jets exiting a fully developed pipe at low Reynolds number. Hence, for the study of impinging jets, SST k-ω turbulence model was used to study the velocity and jet growth characteristics. It became evident that the mesh alignment with the velocity vector at exit of the pipe domain plays a crucial role in the accuracy of the results. The parameter used to evaluate this condition was identified as False Diffusion and was observed to affect the TKE parameter significantly. Methods to reduce false diffusion are also discussed in this article. Based on the mesh obtained from the grid sensitivity study, jets impinging at 30, 45 and 60 degrees at Reynolds number of 7500 were numerically analyzed. It was observed that the profile of the resultant jet closely matched with the prediction of elliptical profile predicted by past researchers, [2] and [3]. Also, it was seen that higher jet growth was predicted in case of jets impinging at a higher impingement angle.

Tidal stream energy is a low carbon energy source. Tidal stream turbines operate in a turbulent environment, and the effect of the structure between the turbine and seabed on this environment is not fully understood. An experimental study using 1:72 scale models based on a commercial turbine design was carried out to study the support structure influence on turbulence intensity around turbine blades. The study was conducted using the wave-current tank at LABIMA, University of Florence. A realistic flow environment (ambient turbulent intensity = 11%) was established. Turbulence intensity was measured upstream and downstream of a turbine mounted on two different support structures (one resembling a commercial design, the other the same with an additional vertical element), in order to quantify any variation in turbulence and performance between the support structures. Turbine drive power was used to calculate power generation. Acoustic Doppler Velocimetry was used to record and calculate upstream and downstream turbulence intensity. In otherwise identical conditions, performance variation of only 4% was observed between two support structures. Turbulent intensity at 1, 3 and 5 blade diameters, both upstream and downstream, showed variation up to 21% between the two cases. The additional turbulent structures generated by the additional element of the second support structure appears to cause this effect, and the upstream propagation of turbulent intensity is believed to be permitted by surface waves. This result is significant for the prediction of turbine array performance.

The three-dimensional single fin configuration finds application in an intake geometry where the cowl-shock wave interacts with the side-wall boundary-layer. Accurate numerical simulation of such three-dimensional shock/turbulent boundary-layer interaction flows, which are characterized by the appearance of strong crossflow separation, is a challenging task. Reynolds-averaged Navier-Stokes computations using the shock-unsteadiness modified Spalart-Allmaras model is carried out at Mach of 5 at large fin angle of 23^◦. The computed results using the modified model are compared to the standard Spalart-Allmaras model and validated against the experimental data. The focus of work is to implement the modified model and to study the flow physics in detail in the complex region of swept-shock-wave turbulent boundary-layer interaction in terms of the shock structure, expansion fan, shear layer and the surface streamlines. The flow structure is correlated to the wall pressure and skin friction in detail. It is observed that the standard model predicts an initial pressure location downstream of the experiments. The modified model reduces the eddy viscosity at the shock and predicts close to the experiments. Overall, the surface pressure using modified model is predicted accurately at all the locations. The skin friction is under predicted by both the models in the reattachment region and is attributed to the poor performance of turbulence models due to flow laminarization.

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.

In this paper, spatial pulse position modulation (SPPM) is used as a case study to investigate the performance of the optical spatial modulation (SM) technique in outdoor atmospheric turbulence (AT). A closed-form expression for the upper bound on the asymptotic symbol error rate (SER) of SPPM in AT is derived and validated by closely-matching simulation results. The error performance is evaluated in weak to strong AT conditions. As the AT strength increases from the weak to strong, the channel fading coefficients become more dispersed and differentiable. Thus, a better error performance is observed under moderate-to-strong AT compared to weak AT. The performance in weak AT can be improved by applying unequal power allocation to make FSO links more distinguishable at the receiver. Receive diversity is considered to mitigate irradiance fluctuation and improve the robustness of the system to turbulence-induced channel fading. The diversity order is computed as half of the number of detectors. Performance comparisons, in terms of energy and spectral efficiencies, are drawn between the SPPM scheme and conventional MIMO schemes such as repetition coding and spatial multiplexing.

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.