Resolution Of Boundary Layers Research Articles

An algebraic global linelet preconditioner for incompressible flow solvers

Given the wide range of scales present in fluid dynamics, highly anisotropic meshes are often employed to resolve boundary layer flow at high Reynolds numbers. However, such meshes result in distorted elements, which can deteriorate the convergence of the preconditioned conjugate gradient solver (PCG) used to solve the discrete pressure equations. To address this issue, the Linelet Preconditioner is commonly used to accelerate the convergence rate of the PCG by building line segments (linelets) along the direction of the strongest couplings and applying a specific operation to each linelet. Under the current state of the art, linelets are constrained to operate independently within each domain partition. In these scenarios, PCG convergence was observed to deteriorate when the number of domain partitions is increased. This work proposes a communication step to be integrated into the preconditioning step, enabling different linelets to exchange information. The proposed method is observed to improve the convergence rate of linear systems with high stiffness due to highly stretched elements. Also, an algorithm is developed to generate the preconditioning matrix by purely algebraic considerations, in contrast to the usual approach of building linelets based on geometrical considerations. Furthermore, since the current approach is agnostic to the domain partition, it does not impose any constraint on the domain decomposition. Hence, the proposed method is a robust preconditioner with a numerical performance almost independent of the quality of the domain decomposition, providing a way for a more balanced load distribution.

Flow around a pair of 2D cylinders using a hybrid Eulerian-Lagrangian solver

The field of external aerodynamics encompasses various engineering disciplines with a significant impact on wind energy technology. Aerodynamic investigations provide insights not only into the characteristics of individual blades or standalone wind turbines but also into entire wind farms. As advancements in wind turbine design continue, understanding the interactions between turbines in close proximity becomes crucial, presenting a multi-body problem. Researchers require efficient and accurate tools to comprehensively study such dynamics. This paper presents a hybrid Eulerian-Lagrangian solver designed to leverage the strengths of Eulerian solvers in resolving boundary layers and Lagrangian solvers in convecting wakes downstream without introducing significant numerical diffusion. The solver adeptly handles multi-body simulations, allowing the construction of independent Eulerian meshes that communicate seamlessly through Lagrangian particles. In this way, the computational study of multibody problems does not require very large and dense meshes. Validation in single-body cases has already been conducted, with this paper demonstrating the solver’s application to a pair of cylinders in different configurations. A comparative performance analysis is carried out against pure Eulerian solvers. The results highlight that the hybrid solver efficiently reproduces the accuracy of the Eulerian solver, demonstrating its effectiveness in handling complex aerodynamic simulations.

Improving Stratocumulus Cloud Amounts in a 200‐m Resolution Multi‐Scale Modeling Framework Through Tuning of Its Interior Physics

AbstractHigh‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation.

Measurements and Prediction of Ash Deposition in a Cyclone-Fired Boiler Operating under Variable Load Conditions

Measurements of ash deposition rates were made between the secondary superheater and reheater sections of a 450 MW cyclone-fired lignite boiler as the operational load varied from 33 to 100%. Significant reductions in deposition rates with a decrease in operational load were observed. To uncover the causative mechanisms behind these observations, operational data from the power plant were used to carry out computational fluid dynamic (CFD) simulations of the boiler. After ascertaining that the gas temperatures and velocities at various sections within the boiler were being represented adequately, decoupled simulations of the ash deposition process on the deposit probe were carried out using a finely resolved boundary layer mesh. Fly ash particle size distribution (PSD) and its concentration for the decoupled calculations were determined from stand-alone cyclone barrel simulations. The ash partitioning (mass %) between the fly ash and slag was found to be ~50:50, which was in line with previous field observations, and it did not vary significantly across different cyclone loads. The predicted PSD of the deposit ash was concentrated in the size range 10–30 microns, which was in agreement with cross-sectional images of the deposit obtained from the measurements. At lower loads, sharp variations in the deposition rates were predicted in the gas temperature range 950–1150 K. The particle kinetic energy—particle viscosity-based capture methodology utilized in this study in conjunction with appropriate ash compositions, ash viscosity models and gas temperature estimates can help estimate slagging propensities at different loads reasonably well in these systems.

An improved slip‐wall model for large eddy simulation and its implementation in the local domain‐free discretization method

AbstractIn this paper, a slip‐wall model for large eddy simulation (LES) is improved and implemented in an immersed boundary method named the local domain‐free discretization (DFD) method. Considering that the matching location may be in the viscous sublayer, the physics‐based interpretation of a Robin‐type wall closure is complemented. Then, the slip‐wall wall model is improved, in which the slip length is redefined and the Robin boundary condition is imposed at the solid wall. The improved slip‐wall model is implemented in the local DFD method to evaluated the tangential velocity at an exterior dependent node, and then the requirement on high resolution of boundary layers can be alleviated. The non‐equilibrium effects are accounted for by adding an explicit correction to the wall shear stress. In order to validate the present wall‐modeled LES/DFD method, a series of turbulent channel flows at various Reynolds numbers, the flow over periodic hills and the flows over a NACA 4412 airfoil at a high Reynolds number are simulated. The predicted results agree well with the referenced experimental data and numerical results. Especially, the results of the separated flow over the airfoil at a near‐stall condition demonstrate the performance of the present wall‐modeled LES/DFD method for complex flows.

Surface Roughness in RANS Applied to Aircraft Ice Accretion Simulation: A Review

Experimental and numerical fluid dynamics studies highlight a change of flow structure in the presence of surface roughness. The changes involve both wall heat transfer and skin friction, and are mainly restricted to the inner region of the boundary layer. Aircraft in-flight icing is a typical application where rough surfaces play an important role in the airflow structure and the subsequent ice growth. The objective of this work is to investigate how surface roughness is tackled in RANS with wall resolved boundary layers for aeronautics applications, with a focus on ice-induced roughness. The literature review shows that semi-empirical correlations were calibrated on experimental data to model flow changes in the presence of roughness. The correlations for RANS do not explicitly resolve the individual roughness. They principally involve turbulence model modifications to account for changes in the velocity and temperature profiles in the near-wall region. The equivalent sand grain roughness (ESGR) approach emerges as a popular metric to characterize roughness and is employed as a length scale for the RANS model. For in-flight icing, correlations were developed, accounting for both surface geometry and atmospheric conditions. Despite these research efforts, uncertainties are present in some specific conditions, where space and time roughness variations make the simulations difficult to calibrate. Research that addresses this gap could help improve ice accretion predictions.

Finite element and isogeometric stabilized methods for the advection-diffusion-reaction equation

We develop two new stabilized methods for the steady advection-diffusion-reaction equation, referred to as the Streamline GSC Method and the Directional GSC Method. Both are globally conservative and perform well in numerical studies utilizing linear, quadratic, cubic, and quartic Lagrange finite elements and maximally smooth B-spline elements. For the streamline GSC method we can prove coercivity, convergence, and optimal-order error estimates in a strong norm that are robust in the advective and reactive limits. The directional GSC method is designed to accurately resolve boundary layers for flows that impinge upon the boundary at an angle, a long-standing problem. The directional GSC method performs better than the streamline GSC method in the numerical studies, but it is not coercive. We conjecture it is inf-sup stable but we are unable to prove it at this time. However, calculations of the inf-sup constant support the conjecture. In the numerical studies, B-spline finite elements consistently perform better than Lagrange finite elements of the same order and number of unknowns.

Isoparametric Finite Element Method to Generate Structured Grid for Numerical Flow Simulation

Eight-nodded quadrilateral finite element is used to generate structured grid for a numerical simulation. The computational domain is sub-divided into number of quadrilateral regions that display the geometry of each sub-domain. The inner and the outer boundaries of the computational domain are described in terms of the surface coordinates. An algebraic homotopy procedure is used to generate grid clustering in order to resolve boundary layer. The structured grid is linked with the flow solver based on finite volume of space discretiztion scheme with multi-stage Runge-Kutta time stepping technique. Examples are illustrated to demonstrate the grid generation procedure and data processing for a forward facing aero-disc spike attached to a hemispherical blunt body at Mach 6 and over a heat shield of a satellite launch vehicle at Mach 0.8 at an angle of attack 5 deg. The present grid generation method is convenient for checking the grid independence test by varying the stretching factor. The quadrilateral grid generated by finite element, vector plot of velocity and contour plots of the computed flow field data are easily drawn with the help of MATLAB.

Accurate RBF-FD meshless numerical simulation of thermo-fluid problems for generic 3D geometries

In contrast to traditional mesh-based methods for the numerical solution of boundary value problems, e.g., Finite Element (FEM) and Finite Volume (FVM), in the recent period many meshfree approaches have been proposed in order to avoid those typical issues due to the mesh. For example, the quality of the mesh greatly affects the reliability of the final solution in the case of CFD problems and the human intervention of a professional is often still needed when dealing with complex-shaped domains. This in turn increases both cost and time required for the reliable simulation of problems of engineering relevance. Meshless methods, on the other side, usually rely on a simpler distribution of nodes and do not require the storage of connectivity information. Among others, one of the most promising meshless methods in terms of accuracy and flexibility is the one based on the Radial Basis Function – Finite Difference (RBF-FD) scheme. RBF-FD methods, however, are usually affected by severe ill conditioning issues when Neumann boundary conditions are employed. This fact is the main responsible for the appearance of large discretization errors near the boundary and for the lack of stability of traditional time integration schemes. In order to address this issue, some new algorithms for the robust treatment of boundary conditions have been developed and successfully employed to solve fluid flow problems with heat transfer. Furthermore, it is well acknowledged that the efficient resolution of boundary layers arising in this class of problems requires an adequate spatial discretization in the neighbourhood of the boundary, i.e., increased node/mesh density along the direction of large gradients only. This result is achieved by employing anisotropic node distributions, which is a novelty in the context of the RBF-FD method to the best of the authors’ knowledge. The method described above is successfully employed for the accurate solution of a representative 3D heat transfer problem with incompressible fluid flow.

Multiscale Analysis of Surface Heterogeneity–Induced Convection on Isentropic Coordinates

Abstract This study uses semi-idealized simulations to investigate multiscale processes induced by the heterogeneity of soil moisture observed during the 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign. The semi-idealized simulations have realistic land heterogeneity, but large-scale winds are removed. Analysis on isentropic coordinates enables the tracking of circulation that transports energy vertically and facilitates the identification of the primary convective processes induced by realistic land heterogeneity. The isentropes associated with upward motion are found to connect the ground characterized by high latent heat flux to cloud bases directly over the ground with high sensible heat flux, while isentropes associated with downward motion connect precipitation to the ground characterized by high sensible heat fluxes. The mixing of dry, warm parcels ascending from the ground with high sensible heat fluxes and moist parcels from high latent heat regions leads to cloud formation. This new mechanism explains how soil moisture heterogeneity provides the key ingredients such as buoyancy and moisture for shallow cloud formation. We also found that the submesoscale dominates upward energy transport in the boundary layer, while mesoscale circulations contribute to vertical energy transport above the boundary layer. Our novel method better illustrates and elucidates the nature of land atmospheric interactions under irregular and realistic soil moisture patterns. Significance Statement Models that resolve boundary layer turbulence and clouds have been used extensively to understand processes controlling land–atmosphere interactions, but many of their configurations and computational expense limit the use of variable land properties. This study aims to understand how heterogeneous land properties over multiple spatial scales affect energy redistribution by moist convection. Using a more realistic land representation and isentropic analyses, we found that high sensible heat flux regions are associated with relatively higher vertical velocity near the surface, and the high latent heat flux regions are associated with relatively higher moist energy. The mixing of parcels rising from these two regions results in the formation of shallow clouds.

RANS-Based Modelling of Turbulent Flow in Submarine Pipe Bends: Effect of Computational Mesh and Turbulence Modelling

Pipe bend is a critical integral component, widely used in slurry pipeline systems involving various engineering applications, including natural gas hydrate production. The aim of this study is to assess the capability of RANS-based CFD models to capture the main features of the turbulent single-phase flow in pipe bends, in view of the future investigation of the hydrate slurry flow in the same geometry. This is different from the available literature in which only a few accounted for the effects of a combination of computational mesh, turbulence model, and near-wall treatment approach. In this study, three types of mesh configuration were adopted to carry out the computations, namely unstructured mesh and two structured meshes with a uniform and nonuniform inflation layer, respectively. To explore the influence of the turbulence model, standard k-ε, low-Reynolds k-ε, and nonlinear eddy viscosity turbulence model were selected to close RANS equations. Pressure coefficient, mean axial velocity, turbulence intensity, secondary flow velocity, and magnitude of secondary flow were regarded as the critical variables to make a comprehensive sensitivity analysis. Predicted results suggest that turbulent kinetic energy is the most sensitive variable to the computational mesh while others tend to stabilize. The largest difference of turbulence kinetic energy was around 26% between unstructured mesh and structured mesh with a nonuniform inflation layer. Additionally, a fully resolved boundary layer can reduce the sensitivity of mesh on turbulent kinetic energy, especially for a nonlinear turbulence model. However, the large gradient and peak value of turbulence intensity near the inner wall of the bend was not captured by the case with a fully resolved boundary layer, compared with that of the wall function used. Furthermore, it has been confirmed that the same rule was detected also for different curvature ratios, Reynolds numbers, and dimensionless wall distance y+.

A posteriori boundary layer detection and resolution in hpq-adaptive finite element methods for 3D-based hierarchical plate and shell models

A posteriori boundary layer detection and resolution in hpq-adaptive finite element methods for 3D-based hierarchical plate and shell models

Assessment of CUPID code used for condensation heat transfer analysis under steam-air mixture conditions

In this study, three condensation models of the CUPID code, i.e., the resolved boundary layer approach (RBLA), heat and mass transfer analogy (HMTA) model, and an empirical correlation, were tested and validated against the COPAIN and CAU tests. An improvement on HMTA model was also made to use well-known heat transfer correlations and to take geometrical effect into consideration. The RBLA was a best option for simulating the COPAIN test, having mean relative error (MRE) about 0.072, followed by the modified HMTA model (MRE about 0.18). On the other hand, benchmark against CAU test (under natural convection and occurred on a slender tube) indicated that the modified HMTA model had better accuracy (MRE about 0.149) than the RBLA (MRE about 0.314). The HMTA model with wall function and the empirical correlation underestimated significantly, having MRE about 0.787 and 0.55 respectively. When using the HMTA model, consideration of geometrical effect such as tube curvature was essential; ignoring such effect leads to significant underestimation. The HMTA and the empirical correlation required significantly less computational resources than the RBLA model. Considering that the HMTA model was reasonable accurate, it may be preferable for large-scale simulations of containment.

High-Order Accurate Numerical Simulation of Supersonic Flow Using RANS and LES Guided by Turbulence Anisotropy

This paper discusses accuracy improvements to Reynolds-Averaged Navier–Stokes (RANS) modeling of supersonic flow by assessing a wide range of factors for physics capture. Numerical simulations reveal complex flow behavior resulting from shock and expansion waves and so, a supersonic jet emanating from rectangular nozzle is considered. PIV based experimental data for the jet is available from literature and is used for validation purposes. Effect of various boundary conditions and turbulence modeling approaches is assessed qualitatively and quantitatively. Of particular interest are the inlet conditions considering the turbulence intensity and the effect of upstream air supply duct, the effect of nozzle wall surface roughness on nozzle internal flow and downstream, wall y+ sensitivity for boundary layer resolution and laminar to turbulent transition modeling. In addition to mesh sensitivity, domain dependency is conducted to evaluate the appropriate domain size to capture the kinetic energy dissipation downstream of the nozzle. To further improve the flow characteristics, accounting for the anisotropy of Reynolds stresses is also one of the focuses. Therefore, non-linear eddy viscosity-based two-equation model and Reynolds stress transport model are also investigated. Additionally, the results of baseline linear (Boussinesq) RANS are compared. Corresponding comparisons with high-fidelity LES are presented. Jet self-similar behavior resulting from all simulation fidelities is assessed and it appears that turbulent flow in LES becomes self-similar, but not in RANS. Finally, various factors such as the nozzle geometry and numerical modeling choices influencing the anisotropy in jet turbulence are discussed.

The interaction of longitudinal vortex pairs with a turbulent boundary layer

In this study, we demonstrate an efficient approach to investigating the interaction of vortex pairs on a turbulent boundary layer. Our aim is to assess how vortex characteristics impact the downstream flow. Wall-modelled large eddy simulations are used together with an inlet defined by a Batchelor vortex model superimposed on a turbulent boundary layer profile, generated using synthetic turbulence and a precursor Reynolds-averaged Navier–Stokes calculation. This set-up allows for the efficient testing of multiple configurations whilst providing adequate resolution of outer boundary layer and the vortex–vortex interactions. After validating the methodology, we report a set of simulations for both co- and counter-rotating vortex pairs at different separations and asymmetric strengths. The separation distance between the vortices was found to have a significant effect upon the merging of the vortices. Asymmetric strength vortex pairs, analogous to vortex generators in yaw, demonstrate performance independent of flow direction for small angles. Detailed analysis of this flow provides insight into turbulence generation mechanisms and Reynolds stress anisotropy – valuable reference data for the development of lower-order models. Skin-friction enhancement is shown to be more effective for counter-rotating vortex pairs than co-rotating pairs of the same strength and spacing. Additionally, a wider spacing between the initial vortex positions results in a faster rise in the skin-friction coefficient.

Computational framework for resolving boundary layers in electrochemical systems using weak imposition of Dirichlet boundary conditions

We present a finite element based computational framework to model electrochemical systems. The electrochemical system is represented by the coupled Poisson–Nernst–Planck (PNP) and Navier–Stokes (NS) equations. The key quantity of interest in such simulations is the current (flux) at the system boundaries. Accurately computing the current flux is challenging due to the small critical dimension of the boundary layers (small Debye layer) that require fine mesh resolution at the boundaries. We present a numerical framework which resolves this challenge by utilizing a weak imposition of Dirichlet boundary conditions for Poisson–Nernst–Planck equations. In this numerical framework we utilize a block iterative strategy to solve NS and PNP equations. This allows us to efficiently and easily implement the weak imposition of Dirichlet boundary conditions. The results from our numerical framework shows excellent agreement when compared to strong imposition of boundary conditions (strong imposition requires a much finer mesh). Furthermore, we show that the weak imposition of the boundary conditions allows us to resolve the fluxes in the boundary layers with much coarser meshes compared to strong imposition. We also show that the method converges as we refine the mesh near the boundaries at a much faster rate compared to strong imposition of the boundary layer. We present multiple test cases with varying boundary layer thickness to illustrate the utility of the numerical framework. We illustrate the approach on canonical 3D problems that otherwise would have been computationally intractable to solve accurately. Lastly, we simulate electrokinetic instabilities near a perm-selective membrane with weakly imposed boundary conditions on the membrane. This approach substantially reduces the computational cost of modeling thin boundary layers in electrochemical systems.

Resolving Boundary Layers with Harmonic Extension Finite Elements

In recent years, the standard numerical methods for partial differential equations have been extended with variants that address the issue of domain discretisation in complicated domains. Sometimes similar requirements are induced by local parameter-dependent features of the solutions, for instance, boundary or internal layers. The adaptive reference elements are one way with which harmonic extension elements, an extension of the p-version of the finite element method, can be implemented. In combination with simple replacement rule-based mesh generation, the performance of the method is shown to be equivalent to that of the standard p-version in problems where the boundary layers dominate the solution. The performance over a parameter range is demonstrated in an application of computational asymptotic analysis, where known estimates are recovered via computational means only.

Impedance boundary conditions for acoustic time‐harmonic wave propagation in viscous gases in two dimensions

We present impedance boundary conditions for the viscoacoustic equations for approximative models that are in terms of the acoustic pressure or in terms of the macroscropic acoustic velocity. The approximative models are derived by the method of multiple scales up to order 2 in the boundary layer thickness. The boundary conditions are stable and asymptotically exact, which is justified by a complete mathematical analysis. The models can be discretized by finite element methods without resolving boundary layers. In difference to an approximation by asymptotic expansion for which for each order 1 PDE system has to be solved, the proposed approximative are solutions to one PDE system only. The impedance boundary conditions for the pressure of first and second orders are of Wentzell type and include a second tangential derivative of the pressure proportional to the square root of the viscosity and take thereby absorption inside the viscosity boundary layer of the underlying velocity into account. The conditions of second order incorporate with curvature the geometrical properties of the wall. The velocity approximations are described by Helmholtz‐like equations for the velocity, where the Laplace operator is replaced by , and the local boundary conditions relate the normal velocity component to its divergence. The velocity approximations are for the so‐called far field and do not exhibit a boundary layer. Including a boundary corrector, the so‐called near field, the velocity approximation is accurate even up to the domain boundary. The results of numerical experiments illustrate the theoretical foundations.

3D unsteady and steady modeling of heat and mass transfer during Cz Si crystal growth with a horizontal magnetic field

A combined LES/RANS approach is applied to modeling Czochralski (Cz) Si growth with a strong horizontal magnetic field stabilizing the melt flow. A model furnace geometry with a cylindrical crucible is considered, which can be suggested for benchmarking Cz Si melt convection of industrial scale. Grid resolution in boundary layers is carefully adjusted and the 4th approximation order for diffusion terms is shown to provide sufficient computation accuracy even with comparatively rough grids. The computed flow structure, temperature distribution on the melt free surface, and RMS temperature fluctuations of the melt are in good agreement with published experimental data and modeling results for industrial furnaces. The effect of the crystal rotation rate on the interface shape and oxygen incorporation into the crystal is analyzed. Due to high thermal conductivity (or low Pr number) of the melt, the unsteady oscillations do not exert a significant effect on the temperature and heat flux distribution in the melt, which allows quite accurate predictions of the interface shape, even neglecting the flow unsteadiness. In contrast, due to a low diffusivity (or high Sc number) of atomic oxygen in the melt, the unsteady oscillations quite significantly affect its transport through the melt and incorporation into the crystal, which results in remarkable overestimation of the oxygen concentration without consideration of unsteady melt flow evolution. The formulation of the LES/RANS approach used in the benchmark problem was verified on experimental data for an industrial furnace for Magnetic Cz (MCz) growth of 300 mm silicon crystals.

An alternative Vorticity based Adaptive Mesh Refinement (V-AMR) technique for tip vortex cavitation modelling of propellers using CFD methods

ABSTRACT This study focuses on the investigation of cavitating flow around the benchmark INSEAN E779A model propeller with the main aim of further improving the computational efficiency of the tip vortex cavitation (TVC) modelling by using a commercial CFD solver. Also, the effects of various key computational parameters including, numerical modelling, grid size, timestep, water quality and boundary layer resolution, on the TVC formation and its extension in the propeller slipstream are investigated systematically. The numerical simulations are conducted in uniform and open water conditions using RANS, DES and LES solvers implemented in the commercial CFD code, Start CCM+. In order to achieve the aim of the study, an alternative and new Vorticity-based Adaptive Mesh Refinement (V-AMR) technique is introduced for enhanced modelling of the TVC on the blades and downstream. For the CFD modelling of cavitation, the Schneer Sauer cavitation model based on the reduced Rayleigh Plesset equation is used for the sheet, tip and hub vortex cavitation. The hydrodynamic results and cavity patterns are validated with the experimental data. The results show that the application of the V-AMR technique further improves the representation of the TVC with minimal increase in computational cost. However, the eddy viscosity at the propeller blade tips increases with applying the V-AMR technique using the RANS solver due to its inherent modelling errors for the solution of the flow inside the tip vortex. This consequently results in an insufficient extension of TVC in the propeller slipstream compared to the predictions by the DES and LES based numerical solvers. Also, the evolution of the TVC is found to be sensitive to the boundary layer resolution when the standard RANS solver is used. The study will help to widen further applications of the CFD methods involving TVC, particularly for propeller induced underwater noise prediction and analysis.