Traditional numerical techniques such as sliding mesh, dynamic grid mesh, and others, have many limitations in dealing with flow simulation with large-scale movement of solid boundaries; which is the case for complex-shaped hydraulic turbomachinery such as propellers, pumps and turbines. The immersed boundary (IB) method provides a new approach to solve the above limitations. Therefore, this study proposes an sharp-interface IB method based on the level-set function that is suitable for simulating the flow through turbomachinery with complex geometries. This method is applied to actual three-dimensional numerical simulations of high-Reynolds number propellers using an in-house computational fluid dynamics solver. The results show that the proposed method can provide comparatively accurate predictions of unsteady load coefficients within the propeller flow passage, and capture the correct propeller wake characteristics as well as the interaction between the propeller wake and free surface. This study aims at providing a theoretical basis and engineering reference for the application of the IB method in engineering numerical simulations.

Wave is a common environmental load that often causes serious damages to offshore structures. In addition, the stability for the submarine artificial slope is also affected by the wave loading. Although the landslide of submarine slopes induced by the waves has received wide attention, the research on the influence of solitary-wave is rare. In this study, a 2-D integrated numerical model is developed to investigate the stability of the foundation trench under the solitary wave loading. The Reynolds-Averaged-Stokes (RANS) equations are used to simulate the propagation of a solitary wave, while the current is realized by setting boundary inlet/outlet velocity. The pore pressure induced by the solitary wave is calculated by Darcy’s law and the seabed is characterized by Mohr-Coulomb constitutive model. Firstly, the wave model is validated through the comparison between analytical solution and experimental data. The initial consolidation state of slope under hydrostatic pressure is achieved as the initial state. Then, the factor of stability (FOS) for the slope corresponding to different distance between wave crest and slope top is calculated with the strength reduction method. The minimum of FOS is defined as the stability index for the slope with specific slope ratio during the process of dynamic wave loading. The parametric study is conducted to examine the effects of soil strength parameters, slope ratio and current direction. At last, the influence of upper slope ratio in a two-stage slope is also discussed.

Metachrony is defined as coordinated asynchronous movement throughout multiple appendages, such as the cilia of cells and swimmerets of crustaceans. We investigate the process of metachronal propulsion, notably used by species of crustaceans and microscopic cells to move through fluid. A rigid crustacean model with paddles moving in symmetric strokes was created to simulate metachronal movement. Coupled with the surrounding fluid domain, we employ the immersed boundary method to analyze fluid-structure interactions. To explore the effect of nonlinear morphology on the efficiency of metachronal propulsion, we generate and simulate a range of crustacean body shapes, from upward curves to downwards curves. We found that the highest propulsion velocity is achieved when the crustacean model morphology follows a downwards curve, specifically a parabola of leading coefficient k=−0.4. This curved morphology results in 4.5% higher velocity compared to the linear model. As k deviates from −0.4, the propulsion velocity decreases with increasing magnitude, forming a concave downward trend. The impact of body shape on propulsion velocity is shown by how the optimal velocity with k=−0.4 is 71.5% larger than the velocity at k=1. Overall, this study suggests that morphology has a significant impact on metachronal propulsion.

Magnetorheological fluids (MRFs) are smart materials consisting of micro-scale magnetizable particles suspended in a carrier fluid. The rheological properties of a MRF can be changed from a fluid-state to a solid-state upon the application of an external magnetic field. This study reports the development of a particle-level simulation code for magnetic solid spheres moving through an incompressible Newtonian carrier fluid. The numerical algorithm is implemented within an open-source finite-volume solver coupled with an immersed boundary method (FVM-IBM) to perform fully-resolved simulations. The particulate phase of the MRF is modeled using the discrete element method (DEM). The resultant force acting on the particles due to the external magnetic field (i.e., magnetostatic polarization force) is computed based on the Clausius-Mossotti relationship. The fixed and mutual dipole magnetic models are then used to account for the magnetic (MAG) interactions between particles. Several benchmark flows were simulated using the newly-developed FVM-IBM-DEM-MAG algorithm to assess the accuracy and robustness of the calculations. First, the sedimentation of two spheres in a rectangular duct containing a Newtonian fluid is computed without the presence of an external magnetic field, mimicking the so-called drafting-kissing-tumbling (DKT) phenomenon. The numerical results obtained for the DKT case study are verified against published data from the scientific literature. Second, we activate both the magnetostatic polarization and the dipole-dipole forces and resultant torques between the spheres for the DKT case study. Next, we study the robustness of the FVM-IBM-DEM-MAG solver by computing multi-particle chaining (i.e., particle assembly) in a two-dimensional (2D) domain for area volume fractions of 20% (260 particles) and 30% (390 particles) under vertical and horizontal magnetic fields. Finally, the fourth computational experiment describes the multi-particle chaining in a three-dimensional (3D) domain allowing to study fully-resolved MRF simulations of 580 magnetic particles under vertical and horizontal magnetic fields.

Oil immersed paper insulation condition is a crucial aspect of power transformer’s life condition diagnostic. The measurement testing database collected over the years made it possible for researchers to implement classification analysis to in-service power transformer. This article presents classification analysis of transformer oil-immersed paper insulation condition. The measurements data (dielectric characteristics, dissolved gas analysis, and furanic compounds) of 149 transformers with primary voltage of 150 kV had been gathered and analyzed. The algorithm used for developing classification model is Support Vector Machine (SVM). The model has been trained and tested using different datasets. Different models have been created and the best chosen, resulting in 90.63% accuracy in predicting the oil-immersed paper insulation condition. Further implementation was executed to classify oil-paper condition of 19 Transformers which Furan data is not available. The classification results combined, reviewed, and compared to conventional assessment methods and standards, confirming that the model developed has the ability to do classification of current oil-paper condition based on Dissolved Gasses and Dielectric Characteristics.

A numerical study for the effect of crest width, breaking parameter and trunk permeability on hydrodynamics and flow behavior in the vicinity of rubble-mound, permeable, zero-freeboard breakwaters (ZFBs) is presented. The modified two dimensional Navier-Stokes equations for two-phase flows in porous media with a Smagorinsky model for the subgrid scale stresses were solved numerically. An immersed-boundary/level-set method was used. The numerical model was validated for the cases of wave propagation over a submerged impermeable trapezoidal bar and over a low-crested permeable breakwater. Five cases of breakwaters were examined, and the main results are: (a) The size of the crest width, B, does not notably affect the wave reflection, vorticity and currents in the seaward region of ZFBs, while wave transmission, currents in the leeward side, and mean overtopping discharge, all decrease with increasing B. A non-monotonic behavior of the wave setup is also observed. (b) As the breaking parameter decreases, wave reflection, transmission, currents, mean overtopping discharge, and wave setup decrease. This observation is also verified by relevant empirical formulas. (c) As the ZFB trunk permeability decreases, an increase of the wave reflection, currents, wave setup, and a decrease of wave transmission and mean overtopping discharge is observed.