Saturday, 29 September 2018

Improvement of the Volume Flow Rate Through a Blower Fan

     In this post, an improvement in the volume flow rate through the blower fan assembly made are presented. The only thing changed in the blower fan was the cross section of the fan blades. In the previous version, the blade cross section resembled a flat plate with fillets at the leading edge. The trailing edge in the previous design was blunt. In the modified design, there aero-foils were selected, namely NACA 9410, NACA 9420 and the NACA 9430. All the other parameters were kept the same to the previous case. The CAD models of the modified fan blades are shown in Fig. 1.

Fig. 1, Fan blade geometries.

     The velocity contours are shown in Fig. 2 while the pressure contours are shown in Fig. 3, super imposed with velocity vectors and the computational mesh. The volume flow rate was the most for the fan with blades having cross-section of NACA 9410 aero-foil, followed by the fan with blades having cross-section of NACA 9420 and the NACA 9430 cross sections, respectively.

Fig. 2, Pressure contours. Row 1, L-R; fan with the NACA 9430 and NACA 9420 cross sections. Row 2, fan with NACA 9410 cross sections.

Fig. 3, Velocity contours. Row 1, L-R; fan with the NACA 9430 and NACA 9420 cross sections. Row 2, fan with NACA 9410 cross sections.

     Thank you for reading. If you would like to collaborate, both scientifically and financially, on research projects, please reach out.

Computational Fluid Dynamics Analysis of a Blower/Centrifugal Fan: Update 01

     In this post the results from a CFD analysis of a blower fan are presented. The fan had a diameter of 66 mm and a height of 12.57 mm. The fan's rotational velocity was at 10,000 rpm. The CAD model is shown in Fig. 1.


Fig. 1, CAD Assembly of the Blower Fan.

     The simulations were completed in SolidWorks Flow Simulation Premium code. The code employs immersed boundary method to create a Cartesian mesh. The sliding mesh feature was employed to simulate the rotation of the fan at atmospheric conditions. The code employs κ-ε model with Two-Scales Wall Functions approach as the turbulence model. The numerical algorithm implemented is the SIMPLE-R, modified. The second-order upwind discretization scheme is used to approximate the convective fluxes while the diffusive terms are approximated using the central differencing scheme. The time derivatives are approximated with an implicit first-order Euler scheme.

     The numerical model for the fan had 816,994 cells of which 209,421 cells were at the solid-fluid interface. Two mesh controls were employed to refine the mesh near the blades of the fan and at the boundary of the stationery and the rotating domains. The results were indeed, mesh independent. Due to the fact that this was an internal flow problem, domain independence test was not applicable. The mesh and the computational domain is shown in Fig. 2. The curved teal arrow represents the direction of rotation of the fan. The blue arrows represent the pressure boundary conditions at the inlet and at the outlet of the fan assembly. The straight teal arrow represents the force of gravity (the arrow is inverted).


Fig. 2, The mesh and the computational domain.

     The pressure and velocity plots are shown in Fig. 3-4.

Fig. 3, Pressure contours.

Fig. 4, Velocity contour

     Thank you very much for reading. If you would like to collaborate on research projects or want a tutorial for the setup of the numerical simulations such as this one, please reach out.

Update 01

     CAD files are available here.

Monday, 10 September 2018

Computational Fluid Dynamics Analysis of a Symmetrical Wing, Update 01

     This post is about the computational fluid dynamics analysis of a wing. The wing analyzed employed the NACA 0021 section throughout. The wing had a span of 4 m and a chord length of 1 m. The Reynolds number was kept at 3,000,000. The software employed was SolidWorks Flow Simulation Premium.

     The mesh had a total of 385,064 cells of which 84,826 cells were in contact with the wing surface, as shown in Fig. 1. The results are, indeed, mesh independent. Mesh controls were employed to refine the mesh near the wing surface. The computational domain employed was of cylindrical shape.

 
Fig. 1, The computational mesh around the wing.
 
     The velocity variation at various angles of attack around the wing cross-section is shown in Fig. 3 while the pressure variation on the wing surface is shown in Fig. 4. The results were validated against experiments conducted by [1].

 
Fig. 2, Velocity variation around the wing at 0-25 degree AOA, 5 degree increments.

 
Fig. 3, Pressure variation at the wing surface at 0-25 degree AOA, 5 degree increments.

     The purpose of this blog is maintain my online portfolio. I did this analysis because I realized I haven't written anything of this nature before. All of my previous simulations and/or blog entries were from the propulsion, renewable energy and turbo-machinery areas.
 

     Update 01

     CAD files are available here.
 
    
     Thank you for reading. If you would like to collaborate on research projects, please feel free to contact.

     [1] Fernando A. Rocha, Adson A. de Paula, Marcos d. Sousa, André V. Cavalieri, and Vitor G. Kleine, "Lift enhancement by wavy leading edges at Reynolds numbers between 700,000 and 3,000,000," Proceedings of the 2018 Applied Aerodynamics Conference, AIAA AVIATION Forum, Atlanta, GA, 2018.