Reynolds-averaged Navier-Stokes Equations Research Articles

The role of annular 2-lobe nozzle in a strut injection system on the mechanism of fuel distribution at a scramjet engine

The strut injection system is a conventional technique for the injection of the fuel in combustor in hypersonic vehicles. The present article investigates the usage of the annular 2-lobe injector to improve the fuel jet behind the strut. Computational studies have been employed for the modeling of fuel jets released from 2-lobe injectors at the supersonic free stream of the combustion chamber. Simulation of the compressible flow is done by solving Reynolds-Averaged Navier-Stokes (RANS) equations with the turbulence model of SST. The roles of the proposed injection configuration on the flow and fuel jet penetration are disclosed in this work. The shape of fuel jet plume behind the injector is also investigated to explain the mechanism of fuel jet distribution in this technique. Presented data disclosed that the usage of the annular 2-lobe nozzle instead of a single circular jet improves fuel mixing by about 50 % behind the strut model.

Data-assisted training of a physics-informed neural network to predict the separated Reynolds-averaged turbulent flow field around an airfoil under variable angles of attack

AbstractPhysics-informed neural networks are a promising method to yield surrogate models of flow fields. We present a metamodeling technique for variable geometries based on physics-informed neural networks. The method was applied to the DU99W350 airfoil at a Reynolds number of $$1\times 10^{5}$$ 1 × 10 5 . Using our technique, the angle of attack was introduced as an additional input parameter of the network and the model was trained to predict the Reynolds-averaged velocity and pressure fields around the airfoil for arbitrary angles of attack between 10.0 and 17.5$$^{\circ }$$ ∘ . Furthermore, we present an effective method to generate the training points for the parameterized geometry. The model was trained with data from simulations for a limited set of angles of attack. Additionally, satisfaction of the a priori known boundary conditions as well as the Reynolds-averaged Navier–Stokes equations was attained. A sensitivity analysis concerning the Reynolds number, the amount and distribution of training data, and the turbulence model was conducted showing the superiority of the pseudo-Reynolds stress method and the demand for labeled training data in the domain. The trained network was capable of predicting the flow separation progressing with angle of attack on the suction surface and exhibited excellent agreement with numerically simulated results, even in the proximity of the wall for interpolations as well as extrapolations from the labeled data set. Our study demonstrates that physics-informed neural networks can be used to obtain accurate flow field surrogate models of variable geometries.

Design and Analysis of a Base Bleed Unit for the Drag Reduction of a High-Power Rocket Operating at Transonic Speeds

In the present study, a passive flow device is considered for drag reduction purposes through implementation in a transonic high-power rocket. The high-power rocket serves as a reference platform that, apart from the operating conditions, enforces several constraints in terms of available volume and placement locations. A step-by-step methodology is suggested, where the unit is initially broken down into an inlet and an outlet component. The flow field is investigated by means of computational modeling (CFD), where the Reynolds-averaged Navier–Stokes equations are solved coupled with turbulence models that vary depending on the design phase and the individual component. In the first design phase, the best alternative configuration is selected for each component by comparing mass flow rates and discharge coefficients. In the second design phase, each component is analyzed in greater detail based on the first phase results. Indicatively, the protruding inlet diffuser-type channel is converted into a protruding inlet nozzle-type channel to avoid choked flow phenomena, and a nozzle geometry is selected as the outlet amongst the other considered scenarios. The two components are eventually integrated into a common base bleed unit and a final assessment is made. The computational results are used to predict the performance and trajectory of the rocket through a well-established trajectory software. The overall methodology is validated against full-scale test flight data. The results show that the base bleed unit developed in the framework of this study yields a drag reduction of approximately 15% at transonic speeds without impacting the rocket mass and stability.

Novel adjustable backflow structure in electric coolant pump for balancing hydraulic performance and reducing temperature rise

An electric coolant pump serves as the source for coolant circulation in the thermal management of electric vehicles. To address the thermal issues for the printed circuit board integrated within the electric coolant pump, a novel backflow structure with an adjustable orifice diameter is designed from the impeller outlet to the printed circuit board to promote convective heat dissipation. Eventually, this flow returns to the impeller inlet, thereby altering the hydraulic performance of pump. The balance between hydraulic performance and temperature rise is achieved by adjusting the orifice diameter. Considering the variations in coolant properties with respect to temperature, numerical simulations of the electric coolant pump are performed using the unsteady compressible Reynolds-Averaged Navier-Stokes Equations. There is good agreement between experimental and numerical simulations results in both cases with and without an orifice. On an average, by increasing the orifice diameter by 1.5 mm which results in a 1.11 % decrease in efficiency, a 0.008 decrease in pressure coefficient, and a reduction of temperature by 2.98 K on printed circuit board. As the orifice diameter increases, leakage losses in the rear side chamber increases with the improvement of heat dissipation. The clash between the backflow jet and the inflow from the suction pipe results in significant entropy generation, forming jet vortices and increasing in pressure fluctuation. However, when the orifice diameter exceeds 3 mm, the cooling advantages no longer justify the trade-off in hydraulic performance. Thus, 3 mm orifice diameter is an optimal compromise between hydraulic performance and temperature control.

Broadband Noise Reduction of a Two-Stage Fan with Wavy Trailing-Edge Blades

In this paper, a numerical investigation is performed to study the broadband noise of a fan stage with wavy trailing-edge blades. A study of the wavelength and ratio of amplitude to wavelength (H/L) is conducted to better understand the noise reduction effect of wavy trailing-edge blades. A rotor–stator interaction broadband noise prediction method based on the result of a Reynolds-averaged Navier–Stokes equation is used. The results show that all wavy trailing-edge configurations reduce the sound power level of the fan stage. The noise reduction effect of H20L10 is the best among all the wavy trailing-edge configurations, and the sound power level is reduced by 2.4 dB at 1000 Hz. When the H/L remains unchanged, the noise reduction effect of the wavy trailing-edge configuration increases with the increase in wavelength. When the wavelength remains unchanged, the noise reduction effect of the wavy trailing-edge configuration with an H/L of 2 is the best. The use of wavy trailing-edge configurations reduces the turbulent kinetic energy and turbulent integral length scale upstream of the stator by changing the wake of the rotor, thereby reducing the rotor–stator interaction broadband noise of the fan stage.

Numerical Tools For Wave Overtopping At Rubble Mound Breakwaters With Submerged Berms

In the era of climate change, the adaptation of existing coastal structures as a response to potential increasing loads, has become a trending topic in coastal engineering. The combination of sea level rise, storm surge, tides and waves is going to dramatically increase the wave overtopping. In this context, this research is focused on the evaluation of wave overtopping for conventional rubble mound breakwaters, modified by the introduction of a submerged berm. Recent studies illustrated how the introduction of emerged berms on the seaward slope reduces the wave overtopping. The goal then becomes to turn the spotlight on the submerged berms, whose effectiveness has yet to be attested, although it has proven useful in reducing wave loads on the armour layer and seabed pressure under the structure. Therefore, the reliability of two different numerical models in detecting overtopping phenomenon have been assessed. IHFOAM and SWASH have been used (see Figure 1). The first one can solve both Reynolds Averaged Navier Stokes and the Volume Averaged Reynolds Averaged Navier–Stokes equations. The second one solves the non-linear shallow water equations with a non-hydrostatic pressure term, representing a simplified form of the Navier-Stokes equations, with associated limited computational cost. The numerical models have been validated based on experimental tests carried out in a wave flume at the University of L’Aquila.

Influence Evaluation of Open Propellers with Boss Cap Fins: Case Studies on Types B4-70 and Ka4-70

Numerical analysis of fins effect on propeller performance was conducted, specifically using the B4-70 and Ka4-70 propellers. The study investigated different types of fins, including bare fins and PBCF (Propeller Boss Cap Fins) using computational fluid dynamics (CFD) simulations. The explicit algebraic stress model (EASM) based on Reynolds-Averaged Navier-Stokes (RANS) equations and turbulence modeling was employed to determine the optimal results. The main objective of this research was to enhance energy efficiency in ships by examining various open propeller configurations. The CFD simulation results for open propellers B4-70 and Ka4-70, with the addition of boss cap fins, revealed interesting phenomena. When the open propellers B4-70 and Ka4-70 were equipped with PBCF, they would experience an increase in efficiency (η0). This was because the performance of the fins functioned optimally when the advance ratio (J) is high, as evident from the high velocity values. Thus, with higher velocity and lower pressure in the boss cap region at high J values, there was an elevation in thrust force due to the reduction of hub vortex. In the case of open propeller B4-70 with added PBCF, there was an increase in the efficiency value (η0) ranging from 3% to 5% when J varied from 0 to 0.7. Similarly, for propeller Ka4-70 with the addition of PBCF, there was an increase in the efficiency value (η0) ranging from 1% to 3% when J varied from 0 to 0.7. Notably, the use of an Energy-Saving Device (ESD) in the form of PBCF can increase the efficiency of ship propeller, as reported in this paper. Consequently, these findings affirmed the reliability of the overall calculations using the CFD approach.

Numerical study of wave run-up on sea dikes with vegetated foreshores

Integrating coastal vegetation into sea dikes is a nature-based approach aimed at combining disaster prevention with ecological sustainability in coastal areas. This study investigates the impact of vegetation on the wave run-up on dikes through numerical analysis. The numerical model used in the study solves the Reynolds-averaged Navier–Stokes equations by adding a vegetation resistance force to account for momentum loss. A stabilized k–ω shear stress transport model considering the vegetation effect was adopted for turbulence closure. A series of numerical simulations was carried out on the wave run-up (Ru) on dikes, focusing on the effects of different vegetation heights, densities, zone lengths, and dike slopes under various wave conditions. The results indicate that vegetation can significantly decrease Ru and may cause the wave to change from breaking to nonbreaking on dikes. The Ru behaviors depend on whether waves break and can be well characterized by the Iribarren number and dimensionless wave momentum flux parameter under breaking and nonbreaking conditions, respectively. Finally, the multivariate non-linear regression (MNLR) and artificial neural network (ANN) methods were adopted to explore a prediction model for evaluating Ru. Comparisons showed that the prediction performance of the ANN model is superior to that of the MNLR model. The ANN model has the potential as a promising predictive tool for obtaining wave run-up on dikes with vegetated foreshores under breaking and nonbreaking conditions.

Suppression of the vortex-induced vibration of the cylinder by the baffle plates at different submerged depths

In this paper, the suppression effect of baffle on VIV of cylinder at different depths is studied. According to two-dimensional transient model and RANS equation, VIV responses of a cylinder with two baffle plates at a fixed distance and a fixed Angle at different depths are numerically simulated. The submersion depth radio h* is mainly 1.2, 1.5, 2.5 and 5. The experimental data and predicting results are compared to display the accuracy of numerical method. When h* is 2.5, the baffle plates have the best suppression effect on VIV. Compared to the bare cylinder, the amplitude is decreased by 69%. Meanwhile, the vibration of the cylinder at different depths shows instability, and the motion response has great jitter. The lower baffle plate is subject to a unidirectional force when h* = 1.5 and 1.2.

Integrating Computational Fluid Dynamics for Maneuverability Prediction in Dual Full Rotary Propulsion Ships: A 4-DOF Mathematical Model Approach

To predict the maneuverability of a dual full rotary propulsion ship quickly and accurately, the integrated computational fluid dynamics (CFD) and mathematical model approach is performed to simulate the ship turning and zigzag tests, which are then compared and validated against a full-scale trial carried out under actual sea conditions. Initially, the RANS equations are solved, employing the Volume of Fluid (VOF) method to capture the free water surface, while a numerical simulation of the captive model test is conducted using the rigid body motion module. Secondly, hydrodynamic derivatives for the MMG model are obtained from the CFD simulations and empirical formula. Lastly, a four-degree-of-freedom mathematical model group (MMG) maneuvering model is proposed for the dual full rotary propulsion ship, incorporating full-scale simulations of turning and zigzag tests followed by a full-scale trial for comparative validation. The results indicate that the proposed method has a high accuracy in predicting the maneuverability of dual full-rotary propulsion ships, with an average error of less than 10% from the full-scale trial data (and within 5% for the tactical diameters in particular) in spite of the influence of environmental factors such as wind and waves. It provides experience in predicting the maneuverability of a full-scale ship during the ship design stage.

Numerical modeling of turbulent air-water mixture flows in spilling breaking waves

This paper presents a numerical study of turbulent air-water mixture flows under spilling breaking waves by using a mixture-flow model. The model is extended from the earlier 2D numerical framework NEWFLUME, which solves the Reynolds-averaged Navier-Stokes (RANS) equations for mean flow motions. The turbulence closure is accomplished by a modified k−ε two-equations model, in which the effect of entrained bubbles on turbulence production and dissipation is considered by incorporating additional buoyancy terms. A transient equation based on two-fluid formulation is solved for air bubble transport. The present model is validated by comparing the simulated free surface displacement, void fraction, mean velocity and turbulence kinetic energy against existing experimental data and the previous numerical results. The effects of the additional convective term in transient equation of void fraction and the bubble-induced turbulence term in turbulence transport equation are further explored. The results reveal that the additional convective term in the transient equation of void fraction plays an essential role of suppressing unphysical turbulence growth and bringing correct amount of air entrainment into water during wave breaking processes. The additional bubble-induced turbulence production and dissipation terms become effective after wave breaks when the void fraction has been appropriately modeled, and they can further reduce the simulated turbulence intensities in spilling breaking wave, rendering a better agreement with the experimental measurements.

Effect of aircraft landing angle on hydrodynamics characteristics during wave landing

In order to study the effect of different courses on the landing performance of aircraft when landing in waves, a Computational Fluid Dynamics (CFD) method is developed to investigate the landing performance of the NACA TN 2929 A aircraft model. The solution of fluid flow is based on the Reynolds-averaged Navier-Stokes equations with the shear stress transport (SST) k−ω turbulence model, and the numerical scheme employed is the Finite Volume Method (FVM). The fifth-order Stokes wave is generated by using the velocity inlet wave-making method. In the work, the slamming load, pitch change and motion response of the aircraft landing at different wave positions are analyzed. The landing characteristics of aircraft landing parallel to the wave, perpendicular to the wave and at the yaw angle are studied. It is found that the landing characteristics of the aircraft at the 1/2 π phase are better than those at other ones, with the maximum load at the π phase. When the aircraft lands parallel to the wave, yaw displacement generated at 1/2 π phase and 3/2 π phase is equal in magnitude and opposite in direction in the initial short time. These results are important for reducing slamming loads and improving the airworthiness of the aircraft structure when landing on a wave surface.

Transformation of internal solitary waves at the edge of ice cover

Abstract. Internal wave-driven mixing is an important factor in the balance of heat and salt fluxes in the polar regions of the ocean. Transformation of internal waves at the edge of the ice cover can enhance the mixing and melting of ice in the Arctic Ocean and Antarctica. In the polar oceans, internal solitary waves (ISWs) are generated by various sources, including tidal currents over bottom topography, the interaction of ice keels with tides, time-varying winds, vortices, and lee waves. In this study, a numerical investigation of the transformation of ISWs propagating from open water in the stratified sea under the edge of the ice cover is carried out to compare the depression ISW transformation and loss of energy on smooth ice surfaces, including those on the ice shelf and glacier outlets, with the processes beneath the ridged underside of the ice. They were carried out using a non-hydrostatic model that is based on the Reynolds-averaged Navier–Stokes equations in the Boussinesq approximation for a continuously stratified fluid. The Smagorinsky turbulence model extended for stratified fluid was used to describe the small-scale turbulent mixing explicitly. Two series of numerical experiments were carried out in an idealized 2D setup. The first series aimed to study the processes of the ISWs of depression transformation under an ice cover of constant submerged ice thickness. Energy loss was estimated based on a budget of depth-integrated pseudoenergy before and after the wave transformation. The transformation of ISWs of depressions is controlled by the blocking parameter β, which is the ratio of the minimum thickness of the upper layer under the ice cover to the incident wave amplitude. The energy loss was relatively small for large positive and large negative values of β. The maximal value of energy loss was about 38 %, and it was reached at β≈0 for ISWs. In the second series of experiments, a number of keels were located on the underside of the constant-thickness ice layer. The ISW transformation under ridged ice also depends on the blocking parameter β. For large keels (β<0), more than 40 % of energy is lost on the first keel, while for relatively small keels (β>0.3), the losses on the first keel are less than 6 %. Energy losses due to all keels depend on the distance between them, which is characterized by the parameter μ, i.e. the ratio of keel depth to the distance between keels. In turn, for a finite length of the ice layer, the distance between keels depends on the keel quantity.

Numerical study on stem-generated turbulence due to emergent rigid vegetation in water waves

Based on the RANS equations and applying the waves2Foam solver in OpenFOAM, the study numerically investigated the turbulent kinetic energy (TKE) generated due to emergent rigid vegetation in water waves. The numerical results show that the TKE generated due to emergent vegetation in water waves increases with decreasing water depth and increasing wave height and period, as well as vegetation density and stem diameter. Moreover, compared to wave height, the influence of wave period on the stem-generated TKE in waves is more significant. Furthermore, the ratio of the amplitude of wave motion trajectory to stem spacing has a significant impact on the stem-generated TKE in waves. When the amplitude of wave motion trajectory is less than half of the stem spacing, the stem-generated TKE in waves is similar to the TKE of the bare bed in waves. In the case where the amplitude of wave motion trajectory is greater than half of the stem spacing, the stem-generated TKE in waves in the front of the vegetation will be significantly enhanced. Based on the study, an improved formula is proposed for predicting the stem-generated turbulence in waves.

A robust computational framework for variational data assimilation of mean flows with sparse measurements corrupted by strong outliers

A variational framework for the reconstruction of time-averaged mean flows using a sparse set of observations with large magnitudes of noise (referred to as outliers), is presented. The observations constitute a set of point-wise measurements of the mean flow with outliers at certain measurement locations and are incorporated into a numerical simulation governed by the two-dimensional, incompressible Reynolds-averaged Navier-Stokes (RANS) equations with an unknown momentum forcing. This forcing, which corresponds to the divergence of the Reynolds stress tensor, is calculated from a direct-adjoint optimization procedure to reduce the deviation between the measured and estimated mean velocities. ℓ2, ℓ1, Huber, and hybrid loss functions are used to represent the discrepancy in the mean velocity field between the measurements and the predictions. A variety of algorithms are considered to solve the optimization problem with these loss functions and a performance comparison in terms of the quality and physical features of the recovered mean flow is presented. The Huber loss function performed best as it remained robust to strong outliers in the measurements with its ℓ1 contribution and also ensured the uniqueness of the optimal solution with its ℓ2 contribution. Huber loss functions restrict the effect of outliers at the local measurement locations, thereby not affecting the quality of the high-dimensional reconstructed mean flow field. The hybrid loss function, a modified form of the continuous Huber loss function, also recovered the mean flow with high accuracy. We demonstrate the performance of the data assimilation framework for the case of two-dimensional laminar flow around a circular cylinder at Re=100. We then extend the analysis to the case of two-dimensional laminar flow over a backward-facing step at Re=500, to further assess the efficacy and robustness of the data assimilation framework.

Integrated Modeling of Coastal Processes Driven by an Advanced Mild Slope Wave Model

Numerical modeling of wave transformation, hydrodynamics, and morphodynamics in coastal regions holds paramount significance for combating coastal erosion by evaluating and optimizing various coastal protection structures. This study aims to present an integration of numerical models to accurately simulate the coastal processes with the presence of coastal and harbor structures. Specifically, integrated modeling employs an advanced mild slope model as the main driver, which is capable of describing all the wave transformation phenomena, including wave reflection. This model provides radiation stresses as inputs to a hydrodynamic model based on Reynolds-averaged Navier–Stokes equations to simulate nearshore currents. Ultimately, these models feed an additional model that can simulate longshore sediment transport and bed level changes. The models are validated against experimental measurements, including energy dissipation due to bottom friction and wave breaking; combined refraction, diffraction, and breaking over a submerged shoal; wave transformation and wave-generated currents over submerged breakwaters; and wave, currents, and sediment transport fields over a varying bathymetry. The models exhibit satisfactory performance in simulating all considered cases, establishing them as efficient and reliable integrated tools for engineering applications in real coastal areas. Moreover, leveraging the validated models, a numerical investigation is undertaken to assess the effects of wave reflection on a seawall on coastal processes for two ideal beach configurations—one with a steeper slope of 1:10 and another with a milder slope of 1:50. The numerical investigation reveals that the presence of reflected waves, particularly in milder bed slopes, significantly influences sediment transport, emphasizing the importance of employing a wave model that takes into account wave reflection as the primary driver for integrated modeling of coastal processes.

NUMERICAL SIMULATION OF TURBULENT FLOW AROUND DARRIEUS AND SAVONIUS ROTORS

The study of unsteady processes during the flow around the rotors of vertical-axis wind turbines was carried out based on the Navier-Stokes equations. General definition of the coupled problem on vertical-axis wind turbines rotors aerodynamics and dynamics is formulated. The solution to the system of initial equations is obtained using an implicit finite-volume numerical algorithm based on the method of artificial compressibility and multi-block structured grids. A code for the calculation of vertical-axis wind turbines rotors aerodynamic and power characteristics is developed on the base of Reynolds averaged Navier-Stokes equations for incompressible fluids. Two types of Darrieus and Savonius rotors with different geometric parameters and number of blades were calculated. The influence of Reynolds numbers, tip-speed-ratios and solidity on the power characteristics of the rotors has been established.

КОМП'ЮТЕРНЕ МОДЕЛЮВАННЯ НЕСТАЦІОНАРНОГО ПРОЦЕСУ РОЗВИТКУ НАДЗВУКОВОГО ПЕРЕРОЗШИРЕНОГО СТРУМЕНЮ

Questions of numerical simulation of non-stationary process of the development of a supersonic overexpanded jet are considered. Unsteady Reynolds-averaged Navier-Stokes equations (URANS) of compressible flows, written in an arbitrary coordinate system, along with the Spalart-Allmaras one-equation differential turbulence model are applied. The calculations of non-stationary formation of an air cold supersonic jet during the blowing through a conical Laval nozzle (simulation of engine start-up) were made. The physical features of the shock and acoustic waves generation and propagation are discussed. In the stationary phase an overexpanded free supersonic turbulent jet is formed, which includes a Mach disc in its structure. The obtained results were compared with experimental data and calculations by other authors based on the ANSYS Fluent and FloEFD CFD tools.

Modelling Turbulent Flows near Coastal Structures Using Various Turbulence Models

Introduction. The reduction in beach width due to erosion is a significant issue that can either be mitigated or exacerbated by coastal protection structures. Modelling breaking waves near the coast and around coastal structures can be used to determine their impact on the dynamics of the coastal zone. The objective of this study is to model and analyze the dynamics of turbulent structures around a single breakwater, obtained using two turbulence modelling schemes: RANS and LES.Materials and Methods. Turbulence induced by breaking waves was investigated. The modelling was based on bathymetric measurements conducted along the Azov Sea coast and a three-dimensional wave hydrodynamics model supplemented with various turbulence calculation configurations.Results. Modelling results of wave processes generating turbulent flows in the presence of coastal protection structures using different turbulence models were obtained. Results obtained based on Reynolds-averaged Navier-Stokes (RANS) equations are compared with the results of Large Eddy Simulation (LES) approach with Smagorinsky dynamic subgrid-scale model (DSM).Discussion and Conclusions. The results showed that wave heights simulated by LES were higher than those simulated by RANS in the front and leeward regions of the coastal protection structure and were lower in its upper part. Thus, according to LES, a greater amount of wave energy was preserved after passing over the breakwater. Velocity vectors of the water medium showed the formation of a vortex when LES was used, whereas no evidence of such turbulent vortices was detected in the case of RANS, confirming the better performance of LES for turbulence modelling in the coastal zone. According to the presented results, LES is the preferred tool for generating turbulence under incoming wave conditions in engineering practices.

Study and optimization of a deformable trimaran amphibious transport vehicle's resistance characteristic

Amphibious transport vehicles have increasingly become an important tool for connecting water and land transportation, serving both water and land transport missions. However, due to the unique shape of amphibious vehicles, their sailing resistance is high, which severely restricts their widespread use. In this paper, a deformable trimaran amphibious transport vehicle with retractable side hulls is designed. Based on the Reynolds-averaged Navier–Stokes equation, a resistance prediction method for amphibious vehicles is established and verified through towing experiments. The resistance characteristic of the trimaran vehicle is studied and optimized. First, the effects of the side hull and its position on resistance performance are analyzed, and the sensitivity of position parameters is evaluated. Then, the bow plate angle on resistance is analyzed. Finally, a surrogate model of the side hull position parameters and the bow plate angle with resistance is established using the Kriging method. The Particle Swarm Optimization is used to search for the minimum resistance configuration Opt 2, which achieves a drag reduction effect of 19.84% at the design speed.