فهرست مطالب

Applied Fluid Mechanics - Volume:15 Issue: 4, Jul-Aug 2022

Journal Of Applied Fluid Mechanics
Volume:15 Issue: 4, Jul-Aug 2022

  • تاریخ انتشار: 1401/01/20
  • تعداد عناوین: 25
|
  • H. Schulz * Pages 959-972
    A second order differential equation for the energy dissipation rate of turbulence is presented. The derivation procedure is explained. The obtained governing equation is a Euler equation, which integration naturally conduces to power laws for the energy dissipation rate as a function of the wavenumber, a result that is extended to the energy spectrum of turbulence. Power laws are obtained for the cases of two equal and two different real roots. For the case of two conjugate complex roots, the solution is a sum of sine and cosine functions of the normal logarithm of the wavenumber. The differential equation accrues from a more basic equation obtained through thermodynamic-type steps that joint part of already consolidated empirical and semi-empirical information on turbulence existing in the literature, and is formally analogue to the Thermodynamics equation of thermal radiation. It is also shown that parameters of turbulence like length and velocity scales may be related to this formulation.
    Keywords: Turbulence spectra, Thermodynamic analogy, Euler differential equation, Empirical laws, Constitutive equations
  • S. H. Delbari, A. Nejat *, A. Hajinezhad Pages 973-984
    This study carries out computational aeroacoustic calculations around the tip region of the CART II wind turbine blade. Additionally, two modified tip designs, namely; tip O and tip R, are further investigated to determine how the geometry of the tip region affects the noise emission characteristics of the wind turbine. The study focuses on the tip vortex noise mechanism using hybrid computational aeroacoustic to tackle the issue of the enormous computational power required for direct noise simulation. Improved Delayed Detached Eddy Simulation (IDDES) technique is used to calculate the instantaneous turbulent flow field near the sound source region, and the noise prediction in far-field is performed using the Ffowcs Williams and Hawking’s (FW-H) acoustic analogy. The method visualizes the flow field near the blades’ tip, assisting researchers to have an accurate understanding of aerodynamically induced noise mechanisms in that highly complex flow region, thus being able to modify tip design in a way that contributes to lower overall noise emission.  The results for the outboard section of the CART II wind turbine’s blade are validated with experimental data. Broadband noise sources such as turbulent-boundary-layer trailing-edge (TBL-TE) noise and the tip vortex noise mechanisms are investigated for the base case as well as tip O and tip R. The results show that the overall sound pressure level (OASPL) and the generated torque of tip R and tip O, are 2.0 %, 5.0 % and 0.8 %, 2.2 % lower than the base case, respectively.
    Keywords: Horizontal axis wind turbine, Aerodynamic noise, Turbulence, Ffowcs Williams, Hawkings acoustic analogy, Improved delayed detached eddy simulation
  • K. Roy *, B. Das, K. K. Pathak, A. Giri Pages 985-998
    The influence of slight inclination, ‘α’ (i.e., 10°, 15°, and 20°) of the channel comprised of shrouded vertical rectangular non-isothermal fin array has been computationally investigated. Simulations are performed to obtain the convective coefficient of heat transfer for the different dimensionless fin spacing (S*= 0.2, 0.3, 0.5), non-dimensional fin tip to shroud clearances (C*= 0.1, 0.2 and 0.3), Grashof numbers (Gr= 1.08×105, 4.42×105 and 11.5×105), fin lengths (L= 0.25m and 0.5m) and fin heights (H= 0.025m, 0.04m, and 0.055m). Hydrodynamic behavior of the fluid indicates that a significant amount of flow reversal occurs near the entrance of the channel at very low inclination which vanishes with the increase in inclination. Further, at higher length, reverse flow is found to occur only in the clearance zone. An increase in the value of ‘α’ from 10° to 20°, results in enhancement of convective coefficient up to a maximum of 161%. An increase in the value of fin spacing from 0.2 to 0.5 results in the enhancement of the convective heat transfer coefficient. At α= 15° and 20°, the heat transfer is enhanced by 84.1% and 101.6%, respectively, while at α=10°, the same is reduced by 33.3%. At lower fin spacing (S*<0.3) increase in C* tends to reduce the efficiency. Further, the efficiency of the finite conductive fin is found to reduce by almost 8% with an increase in Grashof number.
    Keywords: Grashof number, Inclination, Clearance, Convective coefficient, Fin efficiency
  • G. Ou, X. Cao, C. Wang, A. Duan, H. Jin * Pages 999-1015
    The regulating valve applied in coal liquefaction systems becomes seriously worn and its service life decreases because of a significant difference in pressure and solid–liquid flow. This study proposes a wear-resistant multistage pressure relief string-regulating valve, and examines the characteristics of depressurization of flow and the erosion of its throttle element by using computational fluid dynamics (CFD) simulations and the discrete element method (DEM) combined with the E/CRC (Erosion/Corrosion Center) erosion model of the WC-Co coating. The influence of different sizes of the opening of the valve and varying differences in pressure on the characteristics of erosion is analyzed. Moreover, properties of collisions between particles (such as the particle impact velocity and mass per second) are used to represent and analyze the characteristics of erosion. The numerical results show that the cylinder, flat, and bevel of the valve are at high risk of erosion. The characteristics of erosion of the multistage pressure relief string-regulating valve studied in this paper can provide a reference for optimizing methods of erosion prevention.
    Keywords: Multistage depressurization, Regulating valve, String valve, CFD-DEM, Erosion characteristic, Solid–liquid flow
  • T. M. Arispe Angulo *, R. G. Ramirez Gustavo, W. De Oliveira, E. R. Da Silva, G. E. Niño Del Río Pages 1017-1034
    The sound generated by axial flow fans has become increasingly important in several industrial areas. These devices represent considerable noise sources that must be considered from the design stage. The scope of this research work is to present an axial fan design methodology based on the theory of the wing lift and to compare the aerodynamic and aeroacoustics behavior with the free vortex and the non-free vortex conditions. The effect of forward circumferential sweep on the rotor blades is analyzed in both conditions, which consists of a displacement in the tangential direction in the direction of rotor rotation. A cubic polynomial function is proposed for the construction of the blade geometry with circumferential sweep. It defines the center line of the blade and it is taken at 50% of the leading edge of each blade profile. The methodology for analyzing aeroacoustics behavior in fans is based on the integration of Computational Aeroacoustics and Computational Fluid Dynamics in order to analyze aerodynamic behavior, sound power level and sound pressure level at the preliminary design stage. It was verified that although both rotors present a similar aerodynamic behavior (the NFV rotor has a slightly lower performance), according to the analysis of local sound sources, the NFV rotor presents a slight decrease in the sound sources. It was also verified that the NFV rotor presents 2.2 dB less in sound pressure level than the FV rotor.
    Keywords: CFD analysis, Design, Axial flow fan, Circumferential sweep, Aeroacoustics, Computational fluid Dynamic, free vortex, Non-Free Vortex
  • A. Nazarian Shahrbabaki, R. Khoshkhoo *, M. Bazazzadeh Pages 1035-1047
    The current paper aimed to investigate the primary control mechanism of various Nanosecond Dielectric Barrier Discharge (NS-DBD) plasma actuators on the Shock Wave/Boundary Layer Interaction (SWBLI). For this purpose, the effects of the NS-DBD actuator have been investigated on an M=2.8 supersonic flow numerically. The Reynolds Averaged Navier-Stokes (RANS) equations and  SST turbulence model were used as the governing equations to simulate the supersonic flow characteristics. The numerical simulation of the baseline flow (without plasma actuator) was verified using an investigation on wall pressure distribution and the size of SWBLI. Then, NS-DBD phenomenological model based on the energy deposition model in accordance with experimental data was applied to the baseline simulation. Moreover, various stream-wise and span-wise NS-DBD plasma actuator models were used to investigate the actuator effects on the studied flow’s low-density separation zone. Comparing the numerical results of the stream-wise and span-wise actuations revealed that both actuator types cause a momentum transferred to the flow, consequently decreasing the SWBLI region’s size and the boundary layer’s thickness. The results showed that the presence of the NS-DBD actuator increased the local temperature of flow over the insulated electrode. In this regard, a stream-wise NS-DBD actuator with a length of 90 mm upstream of the SWBLI increased the separation flow velocity by 33.7% and decreased the length of the separation region by 5 mm. Also, in this case, after 170 microseconds from the start of actuation, the size of SWBLI decreased by 4.2 mm. Therefore, it can be concluded that the stream-wise type of actuation was more effective in reducing the flow separation and SWBLI size than the span-wise type due to vortex generation into the inlet flow and suppressing the SWBLI region. The proposed NS-DBD actuators were mainly capable of applying the momentum to the boundary layer and reducing the velocity of separated flow in the SWBLI zone. The micro shock wave propagation through the flow associated with the NS-DBD discharge of the actuators can produce more effective high-speed flow control.
    Keywords: Span-wise, Stream-wise actuation, Nanosecond DBD, Plasma Actuator, Flow separation, Shock Wave, Boundary Layer Interaction
  • W. Shang *, Z. Huang, X. Zu, Q. Xiao, X. Jia Pages 1049-1059
    Energy spectrum is an important tool for reflecting time–frequency characteristics of signals. The energy spectra of air shock wave at different explosion distances ware analyzed to investigate the energy evolution mechanism in the propagation and attenuation process. Waveforms of air shock wave were obtained from explosion tests and numerical simulations. Energy propagation and attenuation mechanisms of air shock wave were discussed on the basis of the wave theory of fluid. Instantaneous and marginal energy spectra of air shock wave were calculated using variational modal decomposition (VMD) and Hilbert transform (HT). Energy evolution laws of air shock wave with time, frequency, and explosion distance were analyzed according to the statistical results of the energy spectra. Results showed that the instantaneous energy peak of air shock wave is directly proportional to the square of its pressure peak while inversely proportional to the third power of propagation distance. Nonlinear attenuation of air shock wave will cause frequency dispersion and decelerate the attenuation rate of the total energy of air shock wave. The energy evolution laws of air shock wave with time and frequency reflected by instantaneous and marginal energy spectra were consistent with the theoretical analysis results.
    Keywords: Air shock wave, Energy evolution, Wave propagation, attenuation, Variational modal decomposition (VMD), Hilbert transform (HT), Energy spectrum
  • P. V. Hudson Dharmaraj, T. H. Kim, A. Suryan, H. D. Kim * Pages 1061-1071
    Wet compression has been well known as a promising technology to effectively increase the power output of the gas turbine engine by introducing fine droplets of water into the compressor stages. This lessens the temperature rise throughout the compression process of gas, thus leading to the reduction in the compressor work. A considerable amount of research on the wet compression process has been performed to date, but the detailed compression process of the two-phase gas media and fine droplets is not well understood yet. In the present study, the thermo-fluid analysis has been done on a simple flow model that has a spherical droplet of water inside a cylinder-piston system. The model is validated by using a quasi-steady D2 – Law of evaporation model. The compression rate is varied by employing the piston movement under various flow conditions such as percentage of overspray, water droplet diameter, relative humidity, and temperature. The results obtained show that for a higher percentage of overspray, smaller droplet diameters and a slower compression rate the efficiency obtained is high. Which results in lesser compressor work. The effects of compression rate, droplet diameter, overspray, and the efficiency of the wet compression have been explained and analyzed in detail.
    Keywords: Overspray, Latent heat, Evaporation rate, Wet compression, Compressor power saving, Compression rate
  • M. M. Noroozi, R. Maddahian *, M. H. Ramezani, M. R. Ansari Pages 1073-1085
    The multiscale multiphase flow contains both small-scale (dispersed phase) and large-scale (continuous phase) structures. Standard interface-averaging multiphase models are appropriate for the simulation of flows including small-scale structures. Standard interface-resolving multiphase models are commonly used for the simulation of flow regimes containing large-scale structures. The accurate simulation of different regimes has a crucial role to investigate the physics of multiphase flows. To cover the inability of standard models to simulate multiscale multiphase flows, various generalized hybrid models have been developed. The present research aims to present an LES-like approach to identify the large-scale structures by comparing the equivalent diameter of structures and the averaging length scale. The main difference between the presented model and the models available in the literature is the independency of the model to the thresholds of the local volume fraction to recognize the flow regime. The switching criterion is set based on the cell size and the physical phenomena including the break-up and coalescence mechanisms. To assess the capabilities of the presented multiscale model, four different benchmark cases including the bubble column, the impinging jet, the dam break, and the Rayleigh-Taylor instability are investigated. The physical behavior of the flow is considered as a reference and compared with numerical results. It is demonstrated that the present multifluid model is capable to capture the physical characteristics of both dispersed and segregated flow regimes, and it is a forward step to develop a generalized multiscale hybrid multiphase model.
    Keywords: Multiphase flow, Computational multiphase fluid dynamics, Transient flow regime, Hybrid model, Multiscale model
  • T. Lei, J. Han *, H. Liu Pages 1087-1097
    Flow dynamics accurate prediction is critical for the early-stage design and performance optimization of the piezoelectric (PZT) printhead. To achieve this, the Computational Fluid Dynamics (CFD) method has been widely used. However, for accurate fluid simulation of the PZT printhead, the optimal parameters settings still need to be clarified, which will be discussed in full in this paper. The modelling work is divided into two sub-parts, namely, a three-dimensional (3-D) modelling for the ink chamber and a two-dimensional (2-D) modelling for the nozzle-air domain. Simulations of the 3-D ink chamber were carried out firstly, thereby the transient mass flow rate of the ink outflow from the chamber could be obtained, which will be set as the inlet boundary condition of the 2-D nozzle-air simulations. To ensure accuracy and convergence of 2-D simulations, independence tests of the mesh grid and time step were performed, where the Fine level mesh density and 1e-8 s time step were identified as the optimal choice. And this combination was adopted in the transient simulations of the droplet ejection process. For the model validation purpose, an experimental test rig was developed, and comparisons between the simulations and experimental tests show a good agreement, verifying the accuracy of the developed model. In addition, to validate the feasibility of the developed model, the effect of the ink viscosity on the droplet ejection process was tested, and the results were consistent with those produced by published literature, confirming the feasibility of the CFD model developed in this paper.
    Keywords: 3-D printing, PZT printhead, CFD, Two-phase flow
  • C. Kang *, C. Lu, K. G. Seah, W. Zhang Pages 1099-1109
    The present study aims to describe characteristics of cavitation during the startup process of a condensate pump. The pump is featured by an impeller equipped with five splitter blades. A computational fluid dynamics (CFD) work was conducted to plumb the evolution of cavitation in the pump. Effect of the volumetric flow rate on instantaneous cavitation patterns as the rotational speed of the pump increased was analyzed. The results show that high resistance to cavitation of the pump depends greatly on large area of the impeller eye, which is related to the deployment of the splitter blades. The splitter blades are insignificantly affected by cavitation. During the startup process, both the pump head and the pump efficiency vary drastically, which is insensitive to the flow rate. At a net positive suction head (NPSH) of 2.0 m, high flow rates are responsible for intensified cavitation. High volume fraction of cavitation arises near the inlet of long blades. As the rotational speed increases, the evolution of cavitation is featured by intermittency and diversified cavity patterns. Furthermore, the sum of the volume fraction of cavitation fluctuates with continuously increasing rotational speed.
    Keywords: Cavitation, Condensate pump, Startup, Flow rate, Net positive suction head
  • L. J. Zhai, H. X. Chen *, Z. Ma Pages 1111-1124
    A new turbulent model based on Delayed Detached Eddy Simulation (DDES) with non-linear eddy viscosity model (NLEVM) was developed to predict the complex turbulent flow. The numerical simulation of the triangular cylinder and the centrifugal pump was carried out to investigate the ability and applicability of the DDES model based on NLEVM (DDES_NL). Compared to the turbulent model based on the eddy viscosity model, the computational results of the triangular cylinder showed the advantage of the non-linear eddy viscosity modification in the DDES_NL model which can improve the accuracy of the prediction in the flow phenomenon with a relatively simple turbulence structure. Regrettably, some small-scale turbulent structures among those still cannot be captured accurately. The numerical simulation of the centrifugal pump predicted by the DDES_NL model shows more abundant flow structures and gets close to the realistic statistical characteristics. It also proves the good applicability of the DDES_NL model in complex flow. This study aims to contribute to the growing area of turbulence modeling by exploring it.
    Keywords: Turbulence model, Non-linear eddy viscosity model, Delayed detached eddy simulation model
  • G. Lin, C. Geng, L. Zhang, F. Liu * Pages 1125-1136
    The lattice Boltzmann models, especially the pseudopotential models, have been developed to investigate multicomponent multiphase fluids in presence of phase change process. However, the interparticle force between different components causes compressibility error in the non-phase-change component. This restricts the model capability in quantitative analysis of the physical foaming process, such as expansion rate and decay time. In the present study, a multicomponent multiphase pseudopotential phase change model (the MMPPCM) is improved by introducing an effective mass form of high-pressure-difference multicomponent model in the non-phase-change component. The improved model is compared with the MMPPCM based on simulations of the phase change process of static and moving fluids, as well as the physical foaming process. Density variation of non-phase-change component and its effect on flow field characteristics are analyzed during the phase change process. Simulation results of physical foaming process lead to about 10% ~ 20% reduction of the compressibility error for the improved model as compared with the results of MMPPCM. The improved model also enhances the computational stability of phase change simulation of the static droplets.
    Keywords: Lattice Boltzmann method, Multicomponent, Multiphase, Physical foaming, Compressibility error
  • M. R. Salimi *, R. Askari, M. Hasani Pages 1137-1153
    Dual-throat Nozzle (DTN) is known as one of the most effective approaches of fluidic thrust-vectoring. It is gradually flourishing into a promising technology to implement supersonic and hypersonic thrust-vector control in aircrafts. The main objective of the present study is numerical investigation of the effects of secondary injection geometry on the performance of a fuel-injected planar dual throat thrust-vectoring nozzle. The main contributions of the study is to consider slot and circular geometries as injector cross-sections for injecting four different fuels; moreover, the impact of center-to-center distance of injection holes for circular injector is examined. Three-dimensional compressible reacting simulations have been conducted in order to resolve the flowfield in a dual throat nozzle with pressure ratio of 4.0. Favre-averaged momentum, energy and species equations are solved along with the standard  k- ε model for the turbulence closure, and the eddy dissipation model (EDM) for the combustion modelling. Second-order upwind numerical scheme is employed to discretize and solve governing equations. Different assessment parameters such as discharge coefficient, thrust ratio, thrust-vector angle and thrust-vectoring efficiency are invoked to analyze the nozzle performance. Computationally predicted data are agreed well with experimental measurements of previous studies. Results reveal that a maximum vector angle of 17.1 degrees is achieved via slot injection of methane fuel at a secondary injection rate equal to 9% of primary flow rate. Slot injection is performing better in terms of discharge coefficient, thrust-vector angle and thrust-vectoring efficiency, whereas circular injection provides higher thrust ratio. At 2% secondary injection for methane fuel, vector angle and vectoring efficiency obtained by slot injector is 8% and 34% higher than the circular injector, respectively. Findings suggest that light fuels offer higher thrust ratio, vector angle and vectoring efficiency, while heavy fuels have better discharge coefficient. Increasing center-to-center distance of injector holes improves thrust ratio, while having a negative effect on discharge coefficient, vector angle and vectoring efficiency. A comparison between fuel injectant of current study and inert injectant in the previous studies indicates that fuel reaction could exhibit substantial positive effects on vectoring performance. Secondary-to-primary momentum flux ratio is found to play a crucial role in nozzle performance.
    Keywords: Dual throat nozzle, Fuel injection, Injection geometry, Thrust-vector angle, Thrust-vectoring efficiency
  • X. Y. Zhang *, D. Li, H. J. Yan, Z. Li, Y. H. Li, X. B. Zhang Pages 1155-1164
    Wake generated by wind turbine can greatly influence the performance of downstream turbine. To better understand the wake self-similarity characteristics of vertical axis wind turbine (VAWT), the shear stress transport (SST) turbulence model with the addition of the γ-Reθ transition model is performed to model a two-blade VAWT at different operating conditions. The simulated blade surface pressure and torque are compared with existing experimental results for validation. Results show that, the simulated results after considering the transition model are more consistent with the experimental results. Analysis of the flow field shows that the average streamwise velocity of the wake in the horizontal plane under different tip speed ratios is asymmetry, but symmetric in the vertical plane. Further analysis indicates that, at different downstream positions, the non-dimensional streamwise velocity deficit in the vertical plane remains self-similarity and basically coincides with the Gaussian distribution curve exclude the wake edges. In addition, the larger the tip speed ratio, the easier streamwise velocity deficit reach self-similarity state at downstream of the VAWT. The results of this study will be helpful to establish the wake model of the VAWT.
    Keywords: Vertical axis wind turbine, Wake, Self-similarity, CFD, Tip speed ratio
  • A. Bel Hadj Taher *, H. Kanfoudi, R. Zgolli Pages 1165-1177
    In industry, the phenomenon of cavitation erosion can reduce the lifetime of the components of hydraulic machines. In this article, we present a new numerical approach to predict the mechanical impact resulting from the implosion of a cloud of bubbles, based on an energy approach. The objective of this approach is to determine the main damage mechanisms and to estimate the intensity of the impact pressure near the surface. The large eddy simulation (LES) approach is coupled with a homogeneous cavitation model to assess the risk of erosion around the hydrofoil NACA0009. Indeed, three functions, namely the Pressure Intensity Function (PIF), the steam intensity function and the Erosive Power Function (EPF), are applied to assess the spatial distribution of eroded areas. The calculations show that the functions based on the pressure term are in good agreement with the experiments, namely: the PIF and EPF functions. On the other hand, we assume that the implosion of the cloud of bubbles produces a pressure wave, which in turn causes the implosion of small bubbles near the wall. Then the erosion will be the result of these secondary implosions and not of the cloud of bubbles. Therefore, we vary the degree of proximity of these micro-bubbles near the wall to choose either the shock wave or the micro-jet to extrapolate the pressure field. We can compare these estimates with the existing erosion measurements and we can conclude that the calculations respond more to the probability of the presence of a micro-jet than to the presence of a shock wave.
    Keywords: Cavitation erosion, Shock wave, Micro-jet, 3D unsteady simulation, Ansys-CFX
  • Z. Li, Y. Liu *, J. G. Yang, H. Q. Wang Pages 1179-1192
    The study of intelligent design methods is becoming a hot topic for the design of turbine cascades. This paper proposes a data-based policy model to achieve intelligent design. To gain a high-quality policy model, empirical equations and "space extending + elitism" are adopted to dynamically optimize database. This guarantees the quality of the model. Compared to traditional optimization design approaches, the proposed method relies on less human experience to design a turbine cascade. Ten different turbine cascades are used to verify this method. Results show that, aerodynamic performance of the cascades redesigned is either the same as or better than that of the traditional cascades. The computing time is reduced by more than one order of magnitude compared to a "CFD + optimization algorithm" or "surrogate model + optimization algorithm" method. With the advantages in computing time and intelligence, the proposed novel method shows the possibility of replacing traditional design methods.
    Keywords: Cascade design, Surrogate model, Artificial intelligence, Optimization design, Neural network
  • M. H. Zhang *, H. J. Shi Pages 1193-1206
    To solve the problem of the strong noise generated in the galvanizing process on the surface of the guardrail board, optimal design of the outlet structure of the blowing device is carried out according to the sound absorption and noise reduction theory of microperforated plates. The aerodynamic characteristics and aerodynamic noise analysis of the blowing device are investigated by large eddy simulation with dynamic grid technology. The oblique surface of the outlet is processed with blind holes, and then the influence of blind holes on the aerodynamic noise of the blowing device is explored, including different shapes, porosities and depths. The spectral study reveals that when the guardrail board just enters the blowing device, there is greater noise compared to other working conditions. The place with the highest noise sound pressure level (SPL) is at the outlet of the blowing device at the monitoring point of R=1 m and the direction of 90°. The SPLs of the monitoring points at 0° and 180° are smaller than those in other directions, while the SPL distribution of the monitoring points in other directions is relatively even. Compared with the original blowing device, the best noise reduction performance is achieved when the blowing device has cylindrical holes, with a porosity of 10% and a hole depth of 3 mm. The noise reduction value reaches up to 28.4 dB. In addition, an aerodynamic noise test was carried out on the blowing device in the corrugated board galvanizing workshop to demonstrate the correctness of the results of the numerical simulation.
    Keywords: Blowing device, Impinging jet, Aerodynamic noise, Porous parameters, Brand noise reduction
  • N. S. Wahid *, N. Md Arifin, N. S. Khashi’Ie, I. Pop, N. Bachok, M. E. Hafidz Hafidzuddin Pages 1207-1220
    Magnetic nanofluids (MNFs) have been the focus of extensive research nowadays owing to their potential usefulness as a transfer medium. This study is concerned with the boundary layer flow and heat transfer of MNF past a rotating vertical cone with the embedment of the porosity regime and mixed convection. The buoyancy opposing flow on the combined free and forced convection is being emphasized in this study to evaluate the behavior of the fluid within this region and predict the point of the boundary layer transition. The initial formulation of the model is simplified appropriately by employing the suitable similarity transformation. The package of bvp4c MATLAB is employed to execute the numerical solutions. Analysis of stability is also reported. Due to the mixed convection parameter, the opposing flow contributes towards two different alternative solutions, but the second solution is not stable. A higher local Nusselt number are achieved by increasing the concentration of magnetic nanofluid up to 2% and enlarging the mixed convection parameter under the influence of the porosity regime in the vertical rotating cone. It has been established in this study that the addition of cobalt ferrite as the magnetic nanoparticles (MNPs) is proven to have the ability in enhancing the thermal performance of the fluid.
    Keywords: Mixed convection, Nanofluid, Porous medium, Cone, Stability analysis
  • Y. Q. Zhang, J. P. Hu * Pages 1221-1229
    High-order nonlinearities may be an important cause of freak wave generation; however, it is still unclear how to stimulate the generation of freak waves in deep-water random waves. This study employs the modified fourth-order nonlinear Schrӧdinger equation (mNLSE) to simulate the occurrence of freak waves and analyses the influence of high-order nonlinearities on the evolution of random wave trains described initially by the JONSWAP spectrum. In the evolution of freak wave generation, variations in the linear and nonlinear terms of the mNLSE are displayed with the nonlinear growth of surface elevations. For comparison, the corresponding results from the cubic nonlinear Schrӧdinger equation (CSE) and the linear Schrӧdinger equation (LSE) are also obtained. Power spectra and spectral peakedness curves in the evolution of the wave train are also given to analyze the potential mechanism of freak wave formation. Additionally, the probabilities of freak wave appearances are estimated for different initial parameters and different governing equations. The results show that the fourth-order nonlinearity plays an important role in the generation of freak waves, but this single factor is not enough to generate freak waves, and freak wave occurrence is the contribution of multiple factors to the unstable evolution of the wave train. The higher-order nonlinearity, concentrated initial random phases, larger wave steepness, narrower initial spectral width, and smaller sideband instability parameter can increase the probability of freak wave generation.
    Keywords: Freak waves, High-order nonlinearity, mNLSE, CSE, LSE, Random waves
  • J. Ahn, K. Ahn *, H. Y. Lim Pages 1231-1241
    The application of three-dimensional (3D) printing technology in rocket engine development has numerous benefits considering cost and time. Among the various techniques of 3D printing, the selective laser melting method has the advantage of being able to manufacture complex structures and process multiple materials. In this study, five types of coaxial injectors with different internal configurations were manufactured using metal 3D printing. To confirm the atomization and mixing performance according to the structure of each injector's oxidizer and fuel post, a cold flow test using water and air was performed under a wide range of experimental conditions. As a result of analyzing the injection pressure drop, discharge coefficient, spray pattern, breakup length, and spray angle, the shape of the oxidizer post had a significant influence on the performance of the injector. In comparison, the effect of the fuel post structure was relatively small; however, there was a meaningful difference in the breakup length and spray angle depending on the direction of rotation.
    Keywords: 3D printing, Cold flow test, Shear coaxial injector, swirl coaxial injector, spray characteristics
  • R. Harbi Monfared, M. Taeibi Rahni *, M. Zareh, G. Ahmadi, S. Etemadi Haghighi Pages 1243-1253
    With the advent of various advanced materials, the idea of flying like birds has attracted considerable attention in recent years. In addition, aeroacoustics has become an important issue and is being widely studied. In this work, based on the shape of Owls’ wings, an attempt was made to improve the aeroacoustic and aerodynamic performances of conventional aircraft wings. For this purpose, wings with different elements, namely, square, triangular, and semicircular, on their top surface were examined. In addition, three different spatial distributions of the elements according to the Owl’s wings shape were considered. For incompressible airflow, aerodynamic and aeroacoustic parameters of wing with structured surfaces were investigated. Also, a wing with serrations was examined. The results indicate that wings with elements distributed starting from maximum section thickness and continuing up to the trailing edge are the most suitable case for both aerodynamic and aeroacoustic improvements. On the other hand, a two-sided serrated wing and a serrated wing in the trailing edge reduce the sound level significantly. In addition, the use of both elements and serrations delays wing stall and thus markedly increases the maximum lift coefficient.
    Keywords: Aeroacoustics, Aerodynamics, UAVs, Owl, Wing with structured surfaces
  • S. Ma *, F. Li, L. Gui, D. Chen, T. Duan Pages 1255-1268
    To improve the compressor performances at variable conditions, in this research, we used numerical simulation software to analyze and discuss the effects of three kinks of a wavy blade-end on an axial compressor cascade. One of the wavy blade-ends was a traditional structure with a wavy leading edge and a trailing edge (WBE_P). The other two wavy blades were new structures with a straight trailing edge and a wavy leading edge (WBE_R, WBE_S). The difference between WBE_R and WBE_S was that the blade profile used to form the wavy blade was different. For WBE_R, we obtained the wavy leading edge by rotating the element at different heights, and the trailing edge was the center of a circle. For WBE_S, we obtained the wavy leading edge by using a new element with a different camber angle. Because of the new geometry involving some discussions of a variety of parameters, we further discussed the geometric parameters and structure to guide the design and optimization of the blade. We drew some conclusions at the stable condition (6° incidence). For the first case (WBE_P), when the wavelength of the wavy blade-end was increased to 0.3 times the blade height, the total pressure loss was enhanced. When the wavelength was 0.2 times the blade height, the amplitude was 0.06 times the pitch, the minimum area of blocked range (Ab) and the total pressure loss coefficient ζ could be obtained, and the ζ and Ab were decreased by 5.49%, and 12.72% compared with the baseline. Furthermore, to weaken the influence of changed blade profile on the exit airflow angle, we proposed WBE_P and WBE_R with new structures. We further improved the flow-field characteristics of WBE_P by WBE_R, and the ζ and Ab decreased by 5.84% and 13.79% compared with the baseline. Additionally, the stall point was delayed from 7.5° to 7.7° incidence. Therefore, the wavy blade-end applied in the leading edge had a greater contribution to improving the flow characteristic, and the straight trailing edge was beneficial to maintaining the exit flow angle. WBE_S with a smoother blade surface did not show an advantage compared with WBE_R. WBE_R with a wavy blade-end leading edge and a straight trailing edge showed an advantage in improving the cascade performance, which was suitable for use in the stator from the last stage of a high-load axial compressor.
    Keywords: Wavy blade-end, Corner separation, Axial compressor cascade, Secondary flow, Flow control
  • H. Sevil Ergur * Pages 1269-1280
    This study aims the use of an intelligent system to analyze and enhance the performance of a single stage centrifugal pump (SSCP) with respect to the blade number of impeller, which is crucial for the centrifugal pump design. In general, maximizing the efficiency (h) is the most common performance cost. Moreover, maximizing the pump head (H) and minimizing the power (P) are the other important criteria that should be considered for the optimal blade number. These goals are not simultaneously considered usually in the studies on this topic. The motivation of this study is to expand the system management perspective by evaluating the performance of the system with these multiple criteria. The centrifugal pump design is a typical Multi-Objective Optimization (MOO) problem. The MOO approach consists of five decision variables and three objective functions with different blade numbers for impellers. To determine the optimal solution of this MOO problem, the experimental data is initially expanded using a learning by examples methods based on an intelligent network model training. Then, extended data is used for calculating performance measures at each input-output pairs. In order to evaluate the results of the proposed approach on a real-time system, B50-200/100 type model pump with the speed of 2950 rpm and an outlet angle of 22.150 was tested at State Sugar Machine Plant, Eskisehir (Turkey). Since the main purpose is to determine the optimal blade number, different blade numbers has been used while other geometric parameters were kept constant. To determine the optimal solution, experimental data has used to train a selected soft computing model known as FWNN model. This model has advanced to express the relationship between the inlet and outlet values of the centrifugal pump. FWNN model achieves the general characteristics of the performance measure. The analysis with the use of FWNN model shows that, the optimum number of blades by considering the specified performance parameters for the centrifugal pump design is seven. Comparing with pump impeller with number of blades 5, 6, 8 and 9 increases in efficiency rates are 0.4%, 2.57%, 6.4% and 7.2%, respectively. The FWNN model over the performance analysis algorithm completes the missing data if exists and indicates the best performance solution as given.
    Keywords: Centrifugal pump, Impeller blade, Blade number, Optimization, FWNN model
  • X. Gu, C. Liu *, Y. Zhang, H. Jiang Pages 1281-1289
    To achieve an automatic technology for over-the-shoulder (OTS) launching of air-to-air missiles, this study numerically simulated the overturning process of a slender body by using the dynamic mesh method in the ANSYS Fluent 2021 software. Motion trends and force conditions during the self-turning process were obtained for different center of gravity positions. This investigation showed that a proper center of gravity position was essential for achieving the self-turning of a slender body at high and extra-wide angles of attack. The pressure center of the slender body jumped (discontinuously changed) during the overturning process. The change in the relative position between the pressure center and the center of gravity caused the angular velocity of the slender body to first increase, then decrease and gradually stabilize. These results can be used as a reference for designing the structures of self-turning slender bodies and to realize a new technology for the OTS launching of air-to-air missiles.
    Keywords: Self-turning, Dynamic mesh, Pressure center, Missile, Center of gravity