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

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.

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

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.

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

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

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

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.

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

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

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.

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.

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.

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

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.

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

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

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

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

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.