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 the results of direct numerical simulations of power spectral densities for kinetic energy, convective entropy and heat flux for unsteady Rayleigh-Bénard magnetoconvection in the frequency space. For larger values of frequency, the power spectral densities for all the global quantities vary with frequency f as f⁻². The scaling exponent is independent of Rayleigh number, Chandrasekhar’s number and thermal Prandtl number.

Polymeric drag reducers have been developed over many years due to the great number of practical applications. In all of them, the molecular stability is an essential requirement. Usually, polymers break down under turbulent flows, which causes a decrease in their efficiency as drag reducers. Besides that, some specific applications, in agro and biomedical fields, impose a specific requirement that must be fulfilled, which is the use of non-toxic materials. A suitable stable material that is elected to accomplish this necessity is the mucilage of aloe vera, which is a bio-polymer that can be used as an alternative to the synthetic ones. Here, we investigate the role played by the aging of aloe vera on its capacity to reduce drag. The results obtained by ¹H nuclear magnetic resonance indicate that the compositions of young and mature leaves of aloe vera are different and such a difference plays an important role on their efficiency as drag reducers. Tests were performed to analyse the drag reduction in a rotating apparatus and in a pipeline system and the efficiencies of leaves of different ages were compared to their composition. The main conclusion of these experiments is that the young mucilage samples, which are richer in complex polysaccharides and exhibit lower acid contents, are more efficient drag reducers.

Numerical simulations of turbulent flows in a stirred dead-end membrane bioreactor are performed by the RNG k-ε model based on finite volume method using Fluent codes. Comparisons of numerical and experimental results confirm the reliability and feasibility of the constructed model. Flow structures such as wake flows and circulation loops in stirred flows were well simulated. An increase of stirring speed is proposed to use to minimize the low velocity region. The single vane stirrer is found to be beneficial for biological separations. Results reveal that the increase of vane number can enhance the mixing effect in flow domains. However, a circular disk stirrer goes against the formation of vertical circulations. The six-vane stirrer is found to be able to provide a uniform distribution of high shear stress.

Nowadays, rotor-stator mixers are widely used in the chemical process industry to carry out the dispersion and de-agglomeration processes for manufacturing the homogenous electrode slurries. Development of the process simulation models for the turbulent mixing of these slurries helps to optimize the mixing parameters such as rotational speed and dispersion time. In this work, computational fluid dynamics (CFD) approach through ANSYS Fluent are explored to determine and analyze the flow patterns, velocity fields and viscosity distribution during the dispersion of anode slurries. It utilizes a non-homogeneous Eulerian multiphase approach and the typical k-ε model to create turbulence fields and compute the mixing rate. This work significantly improves our understanding of the relationship between rotor speed, dispersion time, viscosity distribution, and mixer geometry in rotor-stator mixers.

A century-old scientific conundrum is solved in this paper. The Prandtl mixing length modelled plane boundary turbulent flow is described by: $\frac{du^+}{d{y^+} }+\kappa^2{y^+} ^2\left( \frac{du^+}{d{y^+} }\right)^2=1$, together with boundary condition $ {y^+} =0:\, u^+=0$. Only an approximate solution to this nonlinear ordinary differential equation (ODE) has been sought so far, however, the exact solution to this ODE has not been obtained. By introducing a transformation, $2\kappa y^+=\sinh \xi$, I successfully find the exact solution of the ODE as follows: $u^+=\frac{1}{\kappa}\ln(2\kappa {y^+} +\sqrt{1+4\kappa^2{y^+} ^2})-\frac{2{y^+} }{1+\sqrt{1+4\kappa^2{y^+} ^2}}$.

Computational fluid dynamics (CFD) modeling has emerged as a valuable tool for investigating complex processes like microencapsulation. This paper aims to validate the ability of CFD simulations to predict particle size distribution in a polymer microencapsulation process. The CFD modeling approach employed an Eulerian multiphase framework, incorporating a discrete population balance model to track the evolution of the droplet population. A realizable k-ε turbulence model and a multiple reference frame strategy were utilized to capture the system’s flow dynamics. The results reveal that while the CFD simulations align well with experimental data at higher agitation speeds (>10000 rpm), discrepancies arise at lower speeds (<7500 rpm), indicating a challenge in accurately capturing turbulent viscous regimes. Despite these challenges, the CFD model demonstrates robust predictive capabilities for droplet formation and distribution in microencapsulation processes, validated by error margins within acceptable limits. The validated model can be used as a reliable tool to guide experimental efforts and optimize process parameters, contributing to an enhanced understanding and control of microencapsulation processes.

Wind loads have become one of the key influence factors for the running safety of vehicles and comfort of passengers. Investigation on the wind speed spectrum characteristics of a moving vehicle under turbulent crosswinds is of great influence. Expressions of the wind speed spectrum of a moving vehicle was obtained from the von Kármán spectrum based on Taylor’s frozen flow hypothesis. The influence factors, including the ratio of the vehicle speed to the wind speed and the wind yaw angle from 15° to 175°, were analyzed. The maximum value of the wind speed spectrum and the corresponding frequency were studied as well. The results show that the maximum values of the wind speed spectrum of the moving vehicle were larger than those of the static vehicle. The maximum value of the wind speed spectrum corresponding to the moving vehicle first increased and then decreased as the wind yaw angle increased. Some of the frequencies corresponding to the longitudinal wind speed spectrum values of the moving vehicles were smaller than those of the static vehicle. For moving vehicles, the frequency values corresponding to the maximum values of the longitudinal wind speed spectrum first increased and then decreased as the ratio of the vehicle speed to the wind speed and the wind yaw angle increased.

In Shark Bay, a large hypersaline bay in Western Australia, longitudinal density gradients force gravitational circulation that is important for Bay-ocean exchange. Observations of vertical stratification and velocity are presented, confirming a steady near-bed dense water outflow from Shark Bay’s northern Geographe Channel that persisted through all stages of the tide. Outflow velocities were 2–3 times stronger than the outflows recorded previously in Naturaliste Channel (in the west) and were more resistant to breakdown by tidal mixing. Estimates of turbulent kinetic energy production derived from ADCP data using the variance method showed a more complex structure in Geographe Channel, due to shear associated with the stratified conditions, with peak levels of turbulence occurring during reversal of tidal flows. For both channels the main source of turbulence was tidal flow along the seabed, with the bottom current speed cubed, |U_b³|, providing a reasonable proxy for tidal mixing and prediction of dense water outflows from Shark Bay the majority of the time. Orientation and deeper water of Geographe Channel along the main axis of the longitudinal density gradient provided an explanation for the predominant outflow from the Bay’s northern entrance. These density-driven currents could potentially influence recruitment of commercially fished scallops and prawns through dispersal and flushing of larvae.

The current study primarily aimed to simulate detonation initiation via turbulent jet flame acceleration in partial-premixed H2-air mixtures. Different vertical concentration gradients were generated by varying the duration of hydrogen injection (referred to as diffusion time) within an enclosed channel filled with air. H2-air mixtures with average hydrogen concentrations of 22.5% (lean mixture) and 30% (near stoichiometric mixture) were investigated at diffusion times of 3, 5, and 60 seconds. Numerical results show that the vertical concentration gradient has a major influence on the early-stage of flame acceleration (FA). In the stratified lean mixture, detonation began in all the diffusion times, and comparing the flame-speed graphs showed that a decrease in the diffusion time and an increase in the mixture inhomogeneity speeded up the flame propagation and the jet flame to detonation transition occurrence in the channel. In the stratified H2-air mixture with an average hydrogen concentration of 30%, transition from a turbulent jet flame to detonation occurred in all the cases, and the mixture inhomogeneity weakened the FA and delayed the detonation initiation.

The pre-chamber technologies can improve the ignition performance of IC engines by more than two orders in magnitude and thereby substantial economic benefits. Compared with the common pre-chamber, a plug-and-play passive scheme is suitable for quick retrofit, which is getting more attention from the automobile industry. Good scavenging is the precondition for improving turbulent jet ignition performance for a passive pre-chamber. Therefore, detailed evaluations of the scavenging process and turbulent jet ignition deserve investigations for new pre-chamber schemes. In this paper, the effects of design parameters on ignition processes of plug-and-play passive pre-chamber were numerically studied, allowing for the lateral angle, orifice diameter, and vertical angle design. Seven pre-chamber schemes were evaluated, and four optimal ones were selected for bench tests. The characteristics of the scavenging process, turbulent jet ignition, and main-chamber combustion were investigated and analyzed. The results show that allowing for the trade-off between ignition energy and scavenging efficiency, the volume ratio of the pre-chamber to clearance is recommended to be 0.2~0.7%, and the corresponding area-volume ratio is 0.003~0.006 mm-1. Compared with the original natural gas engine, the pre-chamber retrofit can save up to 13.2% fuel consumption, which presents a significant improvement in fuel economy.

In the presence of aquatic plants in streams, a variety of flow structures have been reported. Re-sults of earlier investigations indicated that the hydrodynamic features of flows in rivers, and streams are influenced by a wide range of submerged and non-submerged vegetation. In the pre-sent study, eight experiment runs were conducted in two artificial pools built in a laboratory flume to assess the flow structures across vegetation cover. The slopes for entrance and exit sec-tions of these two pools are of 5 and 10 degrees, respectively. A patch of artificial grass covered the entire pool area over the gravel bed. The effect of stream flow velocity was investigated by considering two aspect ratios of 2 and 2.7. Also, distributions of flow velocity, Reynolds normal and shear stresses, turbulence intensities, turbulent kinetic energy (TKE), and higher-order mo-ments of velocity fluctuations have been investigated. By means of quadrant analysis, it is found that the skewness coefficients of velocity profiles are significant parameters to evaluate. Based on skewness values for the bursting processes, along the pool entrance section with a smaller slope, the contribution of sweep motion occurs in a more restricted zone above the canopy, and this im-plies that the ejection motion occurs in a pool entrance section with a higher slope. Immediately above the vegetation canopy, the dominant process is the sweep motion of bursting events. Whereas, at the upstream boundary of the vegetation patch, an outward motion with slightly pos-itive value occurs; and around the downstream boundary of the vegetation patch, the isotropic turbulence appears at the boundary of the vegetation and gravel bed. Over the vegetation canopy in the pool, less mixing eddies result in lower Reynolds shear stress values, and thus the turbulent flow becomes weak with lower turbulent kinetic energy (TKE). Because of the shallow flow in an artificial channel with an aspect ratio less than five and the availability of vertical non-zero ve-locity components, secondary currents become stronger in the channel. Consequently, there is a divergence from a linear dispersion of Reynolds shear stress.

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.

Predictions of mean values of statistical variables of large scale turbulent flow by the widely used basic k-ϵ model and a new fundamentally based model are verified against published results of Direct Numerical Simulations DNS of the Navier-Stokes equations. The verification concerns turbulent channel flow at shear Reynolds numbers of 950, 2000 and 104. The basic k-ϵ model is largely based on empirical formulations accompanied by calibration constants. This contrasts with the new model where descriptions of leading statistical quantities are based on the general principles of statistical turbulence at large Reynolds number. Predicted values of major output variables such as turbulent viscosity, diffusivity of passive admixture, temperature, and fluid velocities compare well with DNS in case of the new model. Significant differences are seen in case of the basic k-ϵ model.

Widely used in hydraulics, the Colebrook equation for flow friction relates implicitly to the input parameters; the Reynolds number, Re and the relative roughness of inner pipe surface, ε/D with the output unknown parameter; the flow friction factor, λ; λ=f(λ, Re, ε/D). In this paper, a few explicit approximations to the Colebrook equation; λ≈f(Re, ε/D), are generated using the ability of artificial intelligence to make inner patterns to connect input and output parameters in explicit way not knowing their nature or the physical law that connects them, but only knowing raw numbers, {Re, ε/D}→{λ}. The fact that the used genetic programming tool does not know the structure of the Colebrook equation which is based on computationally expensive logarithmic law, is used to obtain better structure of the approximations which is less demanding for calculation but also enough accurate. All generated approximations are with low computational cost because they contain a limited number of logarithmic forms used although for normalization of input parameters or for acceleration, but they are also sufficiently accurate. The relative error regarding the friction factor λ, in best case is up to 0.13% with only two logarithmic forms used. As the second logarithm can be accurately approximated by the Padé approximation, practically the same error is obtained also using only one logarithm.

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.

The turbulent Poiseuille flow between two parallel plates is one of the simplest possible physical situations, and it has been studied intensively. In this paper, we propose a modified Prandtl-van Driest mixing length that satisfies both boundary conditions and wall damping effects. With our new formulations, we numerically solve the problem and, moreover, propose an approximate analytical solution of mean velocity. As applications of our solution, an approximate analytical friction coefficient of turbulent Poiseuille flow is proposed.

The problem of fluid dynamics can be greatly simplified if, for every point in space, the strain-rate tensor is diagonalized. This tensor is introduced into the Navier-Stokes equations via material law and divergence of the stress tensor. This article shows that local SO(3)xU(1) gauge fields can be used to locally diagonalize the diffusion components of the strain-rate tensor. The gauge fields resulting from the connection can be interpreted as convection components of the flow, they show properties of quasiparticles and can be interpreted as elementary vortices. Thus, the proposed approach not only offers new insights for the solution and situative simplification of the Navier-Stokes equations, it also uncovers hidden symmetries within the flow convection, allowing - depending on boundary conditions - further interpretation.

The aim of work is to derive an explicit expression for a function of vertical mixing induced by wind-waves. To this end, in the Navier-Stokes equations, a current is decomposed into four constituents: the mean flow, the wave-orbital motion, the wave-induced turbulent and the background turbulent currents. This decomposition allows separating the wave-induced Reynolds stress, Rw, from the background one, Rb. To make a statistical closure for Rw, the Prandtl approach for the background turbulent fluctuations is used that results in an implicit expression for the wave-induced vertical mixing function, Bv. Expression for Bv is specified based on the author’s results for the eddy viscosity found earlier in the frame of the three-layer concept for a wavy air–sea interface, used for modelling wind-drift currents [1]. Finally, the explicit parameterization for Bv(a, u*, z) is found as a linear function in both the wave amplitude at depth z, a(z), and the friction velocity in the air, u*. The linear dependence of function Bv(a) on the wave amplitude provides the enhanced vertical mixing induced by wind–waves in comparison with function Bv(a) having the cubic dependence found in [2], as far as the wind-wave amplitude a(z) decays exponentially with depth.

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.

In this work, we introduce a scalable and efficient GPU-accelerated methodology for volumetric particle advection and finite-time Lyapunov exponent (FTLE) calculation, focusing on the analysis of Lagrangian Coherent Structures (LCS) in large-scale Direct Numerical Simulation (DNS) datasets across incompressible, supersonic, and hypersonic flow regimes. LCS play a significant role in turbulent boundary layer analysis, and our proposed methodology offers valuable insights into their behavior in various flow conditions. Our novel owning-cell locator method enables efficient, constant-time cell search, and the algorithm draws inspiration from classical search algorithms and modern multi-level approaches in numerical linear algebra. The proposed method is implemented for both multi-core CPUs and Nvidia GPUs, demonstrating strong scaling up to 32,768 CPU cores and up to 62 Nvidia V100 GPUs. By decoupling particle advection from other problems, we achieve modularity and extensibility, resulting in consistent parallel efficiency across different architectures. Our methodology was applied to calculate and visualize the FTLE on four turbulent boundary layers at different Reynolds and Mach numbers, revealing that coherent structures grow more isotropic proportional to the Mach number, and their inclination angle varies along the streamwise direction. We also observed increased anisotropy and FTLE organization at lower Reynolds numbers, with structures retaining coherency along both spanwise and streamwise directions. Additionally, we demonstrated the impact of lower temporal frequency sampling by upscaling with an efficient linear upsampler, preserving general trends with only 10% of the required storage. In summary, we present a particle search scheme for particle advection workloads in the context of visualizing LCS via FTLE that exhibits strong scaling performance and efficiency at scale. Our proposed algorithm is applicable across various domains requiring efficient search algorithms in large structured domains. While this manuscript focuses on the methodology and its application to LCS, an in-depth study of the physics and compressibility effects in LCS candidates will be explored in a future publication.

new subsonic blower wind tunnel design has been studied both numerically and experimentally, it is also referred to as “blower” wind tunnel. This paper, is initially aimed to address each sequential stage of the wind tunnel design process. Rather than applying the standard method of modelling solely the flow in the test section, a large-scale CFD model of the whole wind tunnel was employed. The loss of every constituent element was calculated and then all the losses are added up to determine the power needed for the wind tunnel operation which is used as “intake fan” boundary conditions in the CFD model. Then, flow uniformity and turbulent intensity measurements in an empty test section using a Pitot-static tube and hot wire anemometer (HWA) were introduced to validate the CFD results. The results showed that flow quality was significantly affected by boundary layer controllers (honeycomb and mesh screens) in the settling chamber and wide-angle diffuser. Investigations were also conducted to evaluate the flow deficit in the wake area behind a convex hump model using HWA. This was additional experimental tests carried out to validate the suitability of the wind tunnel flow aerodynamic research.

The exact similarity solutions (also called as special exact solutions) of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, no similarity solutions for two dimensional turbulent boundary layers have been reported in the literature. With the help of dimensional analysis and invariance principle, Prandtl mixing length $\ell=\kappa y$ for one dimensional turbulent boundary layer is extended to $\ell(x,y)=\kappa y (1-\frac{y}{\delta})\sqrt{\frac{\nu}{U\delta}} $ for the two dimensional turbulent boundary layers, furthermore, with a similarity transformation, we successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation $f'''+ ff''+\beta(1-f'^2)+\kappa^2[\eta^2(1-\eta)^2f''|f''|]'=0$. As an application, similarity solutions of the two dimensional turbulent boundary layer on a flat plate at zero incidence have been studied in detail. To solve the ordinary differential equation numerically, a complete Maple code is provided.

The exact similarity solutions of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, for two dimensional turbulent boundary layers, no similarity solutions (special exact solutions) have ever been found. In the light of Blasius' pioneer works, we extend Blasius similarity transformation to the two dimensional turbulent boundary layers, and successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation. By author's Maple code, we numerically solve the ordinary differential equation and produce some useful quantities.

Vegetation has a significant influence on velocity distribution and turbulent energy in a confluence channel. A laboratory measurement with ADV was used to investigate the flow through a Y-shaped confluence channel partially covered with rigid vegetation on its inner bank. In this study, the flow velocities in cases with and without vegetation were measured by the ADV in a Y-shaped confluence channel. The results clearly show that the existence of non-submerged rigid plants has changed the internal flow structure, that the velocity in the non-vegetated area is greater than in the vegetated area, and that there is a large exchange of mass and momentum between the vegetated and non-vegetated areas. The velocity on both sides is significantly reduced when vegetation is present. In the vicinity of tributaries, due to the presence of vegetation, the high-velocity area moved rapidly to the middle of the non-vegetated area, and the secondary flow phenomenon disappeared. In the mainstream, when vegetation was present, circulation disappeared, and the degree of lateral mixing decreased. The presence of vegetation caused a great change in the internal flow structure and made the flow in non-vegetated areas more intense.

In today's volatile and competitive business landscape, achieving synergy between supply chain agility and marketing flexibility is imperative for organizations striving to enhance resilience and maintain competitiveness. This qualitative study explores the integration of these strategic elements and their impact on organizational performance in turbulent markets. Semi-structured interviews were conducted with key stakeholders from diverse industries to capture insights into adaptation strategies, challenges, and outcomes associated with aligning supply chain agility and marketing flexibility. Findings reveal that supply chain agility, characterized by rapid response capabilities to disruptions and fluctuating demands, is essential for optimizing operational efficiency and maintaining continuity in uncertain environments. Meanwhile, marketing flexibility enables organizations to adjust strategies swiftly based on real-time market insights, enhancing customer engagement and market responsiveness. The study identifies organizational silos, resistance to change, and technological limitations as primary barriers to achieving seamless integration. Benefits of synergy include improved customer satisfaction, accelerated time-to-market for new products, and enhanced competitive advantage through personalized marketing strategies and optimized resource allocation. Moreover, organizations that effectively integrate supply chain agility with marketing flexibility reported greater innovation, employee satisfaction, and sustainable growth.

The Colebrook equation ζ is implicitly given in respect to the unknown flow friction factor λ ; λ=ζ(Re,ε*,λ) which cannot be expressed explicitly in exact way without simplifications and use of approximate calculus. Common approach to solve it is through the Newton-Raphson iterative procedure or through the fixed-point iterative procedure. Both requires in some case even eight iterations. On the other hand numerous more powerful iterative methods such as three-or two-point methods, etc. are available. The purpose is to choose optimal iterative method in order to solve the implicit Colebrook equation for flow friction accurately using the least possible number of iterations. The methods are thoroughly tested and those which require the least possible number of iterations to reach the accurate solution are identified. The most powerful three-point methods require in worst case only two iterations to reach final solution. The recommended representatives are Sharma-Guha-Gupta, Sharma-Sharma, Sharma-Arora, Džunić-Petković-Petković; Bi-Ren-Wu, Chun-Neta based on Kung-Traub, Neta, and Jain method based on Steffensen scheme. The recommended iterative methods can reach the final accurate solution with the least possible number of iterations. The approach is hybrid between iterative procedure and one-step explicit approximations and can be used in engineering design for initial rough, but also for final fine calculations.

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.

Transient conjugate heat transfer measurements under varying temperature and velocity inlet boundary conditions at incompressible flow conditions were performed for flat plate and ribbed channel geometries. Therefrom local adiabatic wall temperatures and heat transfer coefficients are determined. Those data are analyzed using typical heat transfer correlations, e.g. Nu=CRe^mPr^n determining the local distribution of C and m. It is shown that they are closely linked. A relationship lnC=A−mB is observed with A and B as modeling parameters. They could be related to parameters in log law or power law representations for turbulent boundary layer flows. The parameter m is shown to have a close link to local pressure gradients and therewith near wall streamlines as well as friction factor distributions. A normalization of the C parameter allows to derive a Reynolds analogy factor and therefrom local wall shear stresses.

In this paper, we introduce a comprehensive solution aimed at enhancing the understanding of three-dimensional atmospheric pollutant dispersion. This innovative solution develops a generalized model that extends previous research and is applicable to all parameterization schemes of the equation, including wind speed profiles and turbulent diffusion coefficients, while incorporating the dry deposition criterion. Our methodology involves subdividing the atmospheric boundary layer into distinct sublayers, which facilitates a detailed examination of pollutant dispersion dynamics. Extensive validation with data from the Hanford experiment has demonstrated the accuracy of this solution in simulating pollutant concentrations. The results exhibit a strong correlation between projected and observed concentrations, underscoring the statistical reliability of our approach. This validation situates the statistical indices of our solution within an acceptable range, confirming its accuracy in predicting atmospheric pollutant dispersion. These findings thus establish our solution as a valid and effective method for studying complex environmental phenomena.

The Colebrook equation is a popular model for estimating friction loss coefficients in water and gas pipes. The model is implicit in the unknown flow friction factor . To date, the captured flow friction factor can be extracted from the logarithmic form analytically only in the term of the Lambert -function. The purpose of this study is to find an accurate and computationally efficient solution based on the shifted Lambert -function also known as the Wright -function. The Wright -function is more suitable because it overcomes the problem with the overflow error by switching the fast growing term of the Lambert -function to the series expansions that further can be easily evaluated in computers without causing overflow run-time errors. Although the Colebrook equation transformed through the Lambert -function is identical to the original expression in term of accuracy, a further evaluation of the Lambert -function can be only approximate. Very accurate explicit approximations of the Colebrook equation that contains only one or two logarithms are shown. The final result is an accurate explicit approximation of the Colebrook equation with the relative error of no more than 0.0096%. The presented approximations are in the form suitable for everyday engineering use, they are both accurate and computationally efficient.

The eighty years old empirical Colebrook function widely used as an informal standard for hydraulic resistance relates implicitly the unknown flow friction factor , with the known Reynolds number and the known relative roughness of a pipe inner surface ; . It is based on logarithmic law in the form that captures the unknown flow friction factor in a way from which it cannot be extracted analytically. As an alternative to the explicit approximations or to the iterative procedures that require at least a few evaluations of computationally expensive logarithmic function or non-integer powers, this paper offers an accurate and computationally cheap iterative algorithm based on Padé polynomials with only one -call in total for the whole procedure (expensive -calls are substituted with Padé polynomials in each iteration with the exception of the first). The proposed modification is computationally less demanding compared with the standard approaches of engineering practice, but does not influence the accuracy or the number of iterations required to reach the final balanced solution.

This paper gives a new unified formula for the Newtonian fluids valid for all pipe flow regimes from laminar to the fully rough turbulent. It includes laminar, unstable sharp jump from laminar to turbulent, and all types of the turbulent regimes: smooth turbulent regime, partial non-fully developed turbulent and fully developed rough turbulent regime. The formula follows the inflectional form of curves as suggested in Nikuradse’s experiment rather than monotonic shape proposed by Colebrook and White. The composition of the proposed unified formula consists of switching functions and of the interchangeable formulas for laminar, smooth turbulent and fully rough turbulent flow. The proposed switching functions provide a smooth and a computationally cheap transition among hydraulic regimes. Thus, the here presented formulation represents a coherent hydraulic model suitable for engineering use. The model is compared to existing literature models, and shows smooth and computationally cheap transitions among hydraulic regimes.

Today, Colebrook’s equation is mostly accepted as an informal standard for modeling of turbulent flow in hydraulically smooth and rough pipes including transient zone in between. The empirical Colebrook’s equation relates the unknown flow friction factor (λ) with the known Reynolds number (R) and the known relative roughness of inner pipe surface (ε/D). It is implicit in unknown friction factor (λ). Implicit Colebrook’s equation cannot be rearranged to derive friction factor (λ) directly and therefore it can be solved only iteratively [λ=f(λ, R, ε/D)] or using its explicit approximations [λ≈f(R, ε/D)]. Of course, approximations carry in certain error compared with the iterative solution where the highest level of accuracy can be reached after enough number of iterations. The explicit approximations give a relatively good prediction of the friction factor (λ) and can reproduce accurately Colebrook’s equation and its Moody’s plot. Usually, more complex models of approximations are more accurate and vice versa. In this paper, numerical values of parameters in various existing approximations are changed (optimized) using genetic algorithms to reduce maximal relative error. After this improvement computational burden stays unchanged while accuracy of approximations increases in some of the cases very significantly.

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

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.

An overview is presented of several diverse branches of work in the area of effectively 2D fluid equilibria which have in common that they are constrained by an infinite number of conservation laws. Broad concepts, and the enormous variety of physical phenomena that can be explored, are highlighted. These span, roughly in order of increasing complexity, Euler flow, nonlinear Rossby waves, 3D axisymmetric flow, shallow water dynamics, and 2D magnetohydrodynamics. The classical field theories describing these systems bear some resemblance to perhaps more familiar fluctuating membrane and continuous spin models, but the fluid physics drives these models into unconventional regimes exhibiting large scale jet and eddy structures. From a dynamical point of view these structures are the end result of various conserved variable forward and inverse cascades. The resulting balance between large scale structure and small scale fluctuations is controlled by the competition between energy and entropy in the system free energy, in turn highly tunable through setting the values of the conserved integrals. Although the statistical mechanical description of such systems is fully self-consistent, with remarkable mathematical structure and diversity of solutions, great care must be taken because the underlying assumptions, especially ergodicity, can be violated or at minimum lead to exceedingly long equilibration times. Generalization of the theory to include weak driving and dissipation (e.g., non-equilibrium statistical mechanics and associated linear response formalism) could provide additional insights, but has yet to be properly explored.

Empirical Colebrook equation from 1939 is still accepted as an informal standard to calculate friction factor during the turbulent flow (4000 < Re < 10⁸) through pipes from smooth with almost negligible relative roughness (ε/D→0) to the very rough (up to ε/D = 0.05) inner surface. The Colebrook equation contains flow friction factor λ in implicit logarithmic form where it is, aside of itself; λ, a function of the Reynolds number Re and the relative roughness of inner pipe surface ε/D; λ = f (λ, Re, ε/D). To evaluate the error introduced by many available explicit approximations to the Colebrook equation, λ ≈ f(Re, ε/D), it is necessary to determinate value of the friction factor λ from the Colebrook equation as accurate as possible. The most accurate way to achieve that is using some kind of iterative methods. Usually classical approach also known as simple fixed point method requires up to 8 iterations to achieve the high level of accuracy, but does not require derivatives of the Colebrook function as here presented accelerated Householder’s approach (3^rd order, 2^nd order: Halley’s and Schröder’s method and 1^st order: Newton-Raphson) which needs only 3 to 7 iteration and three-point iterative methods which needs only 1 to 4 iteration to achieve the same high level of accuracy. Strategies how to find derivatives of the Colebrook function in symbolic form, how to avoid use of the derivatives (Secant method) and how to choose optimal starting point for the iterative procedure are shown. Householder’s approach to the Colebrook’s equations expressed through the Lambert W-function is also analyzed. One approximation to the Colebrook equation based on the analysis from the paper with the error of no more than 0.0617% is shown.

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.

This paper attempts to develop a Dynamic Mode Decomposition (DMD) based Reduced Order Model (ROMs) which can quickly but accurately predict the forces and moments experienced by a road vehicle such that they be used by an on-board controller to determine the vehicle's trajectory. DMD can linearize a large dataset of high-dimensional measurements by decomposing them into low-dimensional coherent structures and associated time-dynamics. This ROM then also can be applied to predict the future state of the fluid flow. Existing literature on DMD is limited to low Reynolds number applications. This paper presents DMD analyses of the flow around an idealized road vehicle, called the Ahmed body, at a Reynolds number of 2.7E6. The high dimensional dataset used in this paper was collected from a computational fluid dynamics (CFD) simulation performed using Menter's shear stress transport (SST) turbulence model within the context of Improved Delayed Detached Eddy Simulations (IDDES). The DMD algorithm, as available in literature, was found to suffer nonphysical dampening of the medium-to-high frequency modes. Enhancements to the existing algorithm were explored, and a modified DMD approach is presented in this paper which includes: (a) a requirement of higher sampling rate to obtain higher resolution of data, and (b) a custom filtration process to remove spurious modes. The modified DMD algorithm thus developed was applied to the high Reynolds number, separation dominated flow past the idealized ground vehicle. The effectiveness of the modified algorithm was tested by comparing future predictions of force and moment coefficients as predicted by the DMD-based ROM to the reference CFD simulation data, and found to offer significant improvement.

Downwind odor characteristics can be very different depending upon the size of the upwind point-source, interim topography, and wind conditions. At one extreme, the downwind odor plume from a relatively large, confined animal feeding operation (CAFO), located on a flat open plain and under stable, near-straight-line wind conditions can be rather broad, sustained and predictable relative to a fixed receptor site downwind. In contrast, the plume from a small point-source (e.g., a roof vent stack) located on irregular topography and under rapidly shifting wind conditions can be intermittent and fleeting. These transient odor events can be surprisingly intense and offensive, in spite of their fleeting occurrence and perception. This work reports on efforts to optimize an environmental odor sampling strategy, which is optimized for the challenges of (1) sampling of such transient odor 'spikes' and (2) the prioritization of odors/odorants from multiple, closely co-located point-sources, under such transient event conditions. Protocol refinement has emerged by way of 2 environmental odor assessment projects which have been undertaken on behalf of the Missouri Department of Natural Resources. The challenge of transient odor events has been mitigated utilizing rapid, odor cued whole-air grab capture sampling into metalized-FEP gas sampling bags, followed by immediate adsorption transfer onto SPME fibers or sorbent tubes for stabilization during the shipment and storage interval between collection and final analysis. Initial results demonstrated approximately 11 fold increases in target odorant yields for 900 mL sorbent tube transfers from 2-3 second 'burst' odor event bag-captures, as compared to equivalent direct collections at the same downwind receptor location but during perceived (stable) odor 'lull' periods. Results-to-date targeting refinement and field trials of this integrated environmental odor assessment strategy are presented. Preliminary application targeting general odor sampling and point-source differentiation utilizing tracer gases is also presented.

The aim of the studies was to test the hydraulic pressure losses along the length, local pressure losses, absolute and relative roughness values of the pipeline walls made of modern polymer material. The article presents the results of the experimental studies of hydraulic resistance of pipelines from modern polymer materials, the method of theoretical determination of resistance coefficients in the VALTEC ISO 21003 metal polymer pipeline, the classification and method of determining the values of local hydraulic resistance coefficients of T-Junctions are given, the dependence for determining the Darcy λ hydraulic friction coefficient associated with the value of the equivalent hydraulic roughness coefficient of the inner surface of pipes and the Reynolds number. The experimental data of the hydraulic coefficient of friction Darcy λ was compared with the known theoretical relationships of A.D. Altschul, Colbrooke for smooth pipes, Colbrooke and White, an alternative to the Colbrooke-White, Offengenden equation for low-shear pipes and pipes from polymer materials with Reynolds numbers Re>10^4. Dependence of hydraulic coefficient of friction on relative equivalent roughness of pipeline walls and Reynolds number, specific energy loss per one linear meter on total flow rate in system of hydraulically short pipeline is obtained. The empirical dependence of values of coefficients of local resistance ζ=f(Q2/Q3) from a ratio of the divided expenses of Q2/Q3 of the equal-pass T-Junction α=90 the pipeline to a branch up, on pass in the direct T-Junction and division of streams at various internal pipeline diameters of VALTEC.