Full-scale self-propulsion simulation with a discretized propeller

A comparison between towing tank testing and full-scale CFD simulations is presented at three different target speeds. For the current self-propulsion simulation, the self-propulsion point was obtained using polynomial interpolation. The studies of boundary layer thickness, a basic grid uncertainty assessment and verification were performed to give some confidence of grid application to current self-propulsion simulation. All simulations are performed using a commercial CFD software STAR-CCM+. It is concluded that with high-fidelity numerical methods, it’s possible to treat hull roughness and directly calculate full-scale flow characteristics, including the effects of the free surface, none-linearity, turbulence and the interaction between propeller, hull and the flow field.

Influence of the Filler Particles’ Surface Morphology on the Polyurethane Matrix’s Structure Formation in the Composite

This article presents the surface morphology effect of silicon carbide (SiC) particles on the polyurethane binder’s structure formation in a dispersed-filled composite. The difference in the morphology and surface relief of filler particles was ensured by the implementation of plasma chemical modification. As a result of this modification, the filler consisted of core-shell particles characterized by a SiC core and a carbon shell (SiC@C), as well as a carbon shell decorated with silicon nanoparticles (SiC@C/SiNP) or nanos (SiC@C/SiNW). The study of the relaxation properties of polyurethane composites has shown that the strongest limiting effect on the molecular mobility of boundary layer’s chain segments is exerted by a highly developed surface with a complex relief of SiC@C/SiNP and SiC@C/SiNW particles. An empirical method was proposed to find the polymer fractions spent on the formation of the boundary, transition and bulk layers of the polymer matrix in the composite. It was shown that the morphology of the filler particles’ surface does not affect the dependence of the boundary layer thickness on the filler’s volume fraction. However, with an increase in the degree of surface development, the boundary layer thickness decreases.

Numerical investigations of flow around subsea covers at high Reynolds numbers

Two-dimensional (2D) numerical simulations of flow over wall-mounted rectangular and trapezoidal ribs subjected to a turbulent boundary layer flow with the normalized boundary layer thickness of δ/D = 0.73,1.96,2.52 (D is the height of the ribs) have been carried out by using the Reynolds-averaged Navier-Stokes (RANS) equations combined with the k – ω SST (Shear Stress Transport) turbulence model. The angles of the two side slopes of trapezoidal rib varies from 0° to 60°. The Reynolds number based on the free-stream velocity U ∞ and D are 1 × 106 and 2 × 106. The results obtained from the present numerical simulations are in good agreement with the published experimental data. Furthermore, the effects of the angle of the two side slopes of the trapezoidal ribs, the Reynolds number and the boundary layer thickness on the hydrodynamic quantities are discussed.

On mesh requirements for Large Eddy Simulation with Wall Functions

The results of using Large Eddy Simulation with Wall Functions (WFLES) in application to basic wall-bounded flows, such as turbulent boundary layer and channel flow cases are presented. In particular, it is shown that WFLES is suitable for predicting wall bounded flows and provides reasonable accuracy when using appropriate grids. The grid for WFLES should have isotropic cells with size smaller than 10% of the boundary layer thickness. Using coarser grids or anisotropic cells leads to significant reduction of accuracy.

Spectral computation of reactive bi-directional hydromagnetic non-Newtonian convection flow from a stretching upper parabolic surface in non-Darcyporous medium

This paper presents a mathematical model for bi-directional convection magnetohydrodynamic (MHD) tangent hyperbolic nanofluid flow from the upper horizontal subsurface of a stretching parabolic surface to a non-Darcian porous medium, as a simulation of nanocoating. Chemical reaction, activation energy and thermo solutal buoyancy effects are included. The Darcy–Brinkman–Forchheimer model is deployed which permits the analysis of inertial (second order) porous drag effects. The Buongiorno nanoscale model is deployed which includes Brownian motion and thermophoresis effects. The dimensionless, transformed, nonlinear, coupled ordinary differential equations are solved by implementing the spectral relaxation method (SRM). Validation with previous studies is included. The numerical influence of key parameters on transport characteristics is evaluated and visualized graphically. Velocity is elevated (and momentum boundary layer thickness is reduced) with increasing wall thickness parameter, permeability parameter, Forchheimer parameter, Weissenberg (rheological) parameter and modified Hartmann (magnetic body force) number. Velocity enhancement is also computed with increment in stretching rate parameter, rheological power-law index, thermal Grashof number, and species (solutal) Grashof number, and momentum boundary layer thickness diminishes. Temperature is suppressed with increasing stretching rate index and Prandtl number whereas it is substantially elevated with increasing Brownian motion and thermophoresis parameters. Velocity and temperature profiles are reduced adjacent to the parabolic surface with larger wall thickness parameter for stretching rate index [Formula: see text]1, whereas the reverse behavior is observed for stretching rate index [Formula: see text]1. Nanoparticle concentration magnitude is depleted with larger numeric of Lewis number and the Brownian motion parameter, whereas it is enhanced with greater values of the stretching index and thermophoresis parameter. The nanoparticle concentration magnitude is reduced with an increase in chemical reaction rate parameter whereas it is boosted with activation energy parameter. Skin friction, Nusselt number and Sherwood number are also computed. The study is relevant to electromagnetic nanomaterials coating processes with complex chemical reactions.

The structure of a highly decelerated axisymmetric turbulent boundary layer

Experiments were performed over a body of revolution at a length-based Reynolds number of 1.9 million. While the lateral curvature parameters are moderate ( $\delta /r_s < 2, r_s^+>500$ , where $\delta$ is the boundary layer thickness and r s is the radius of curvature), the pressure gradient is increasingly adverse ( $\beta _{C} \in [5 \text {--} 18]$ where $\beta_{C}$ is Clauser’s pressure gradient parameter), representative of vehicle-relevant conditions. The mean flow in the outer regions of this fully attached boundary layer displays some properties of a free-shear layer, with the mean-velocity and turbulence intensity profiles attaining self-similarity with the ‘embedded shear layer’ scaling (Schatzman & Thomas, J. Fluid Mech., vol. 815, 2017, pp. 592–642). Spectral analysis of the streamwise turbulence revealed that, as the mean flow decelerates, the large-scale motions energize across the boundary layer, growing proportionally with the boundary layer thickness. When scaled with the shear layer parameters, the distribution of the energy in the low-frequency region is approximately self-similar, emphasizing the role of the embedded shear layer in the large-scale motions. The correlation structure of the boundary layer is discussed at length to supply information towards the development of turbulence and aeroacoustic models. One major finding is that the estimation of integral turbulence length scales from single-point measurements, via Taylor's hypothesis, requires significant corrections to the convection velocity in the inner 50 % of the boundary layer. The apparent convection velocity (estimated from the ratio of integral length scale to the time scale), is approximately 40 % greater than the local mean velocity, suggesting the turbulence is convected much faster than previously thought. Closer to the wall even higher corrections are required.

Numerical Simulation of the Motion of a Large Scale Unmanned Surface Vessel in High Sea State Waves

The motion stability of the Unmanned Surface Vessel (USV) is threatened by the action of waves under a rough sea state. In the present paper, the motion of a large-scale USV is numerically simulated under high sea state of level 5 and 7. The overset grid method and Reynolds Averaged Navier–Stokes (RANS) approach are employed to solve Navier–Stokes (N-S) equations. For the case of wave incident angle 0° and 30°, the heave, pitch and roll motion response of a large scale USV are investigated by using the six Degrees of Freedom (6-DOF) numerical model. The effects of different sea states, as well as different wave directions, on the motion of USV are compared. The comparative results indicate that the response of this USV in waves is the periodic free-motion according to the corresponding amplitude, which does not exceed the stable range, and there are no overturning and other situations that may affect the safety, in the case of level 5 and 7 sea states. The corresponding pressure at the bottom of this USV meets the range of material strength, and no structural damage or injury to the hull occurs, although the pressure varies at different wave periods. For the case of different wave directions, the analysis of the boundary layer thickness shows that the wave direction is of great importance to the boundary layer thickness distribution, both in the level 5 and level 7 sea states.

Mean Current Profile over Rippled-Beds in the Presence of Non-Breaking Waves and Analysis of Its Influencing Factors

Classical eddy viscosity model deviates from the actual mean current profiles, when calculating the mean current profiles over rippled-beds in the presence of non-breaking waves, owing to the neglect of the enhancement of the wave boundary layer thickness by ripples and the wave-induced shear stress (the radiation stress and the wave Reynolds stress). Considering these shortcomings, a semi-empirical one-dimensional vertical (1DV) model is presented in this study. The present model was obtained using the two-dimensional Navier–Stokes equations and eddy viscosity assumptions, which differ from those of previous researchers, while a top-to-bottom sequence was adopted to calculate the mean current profiles. Empirical formulae were derived from the laboratory measurements and used in the present model to accurately predict the wave boundary layer thickness and bed roughness. The present model is in satisfactory agreement with the data from laboratory experiments. The factors influencing the mean current profiles were analyzed also. The wave-induced second-order shear stresses were found to be the principal reason for the deviations of the mean current profiles in the near-surface layer; as the influencing factors of wave-induced shear stress, the intensity of the wave relative to the current, the angle between the wave and current, and the size of ripples can also have a non-negligible effect on the mean current profiles.