The effectiveness of honeycomb, in reducing the swirl and turbulence level, is studied by simulating the flow field using standard k-ε turbulence model in commercial computational fluid dynamics (CFD). CFD is the most precise and economical approach to estimate the effectiveness of a honeycomb.
A computational domain of honeycomb is created as shown in Fig. 1
We know computationally, it is very difficult to provide the realistic non-uniform flow at the entry of honeycomb as experienced in the experiments. Such random inlet conditions would essentially simulate the realistic case in which air can enter the honeycomb from any direction and at any level of turbulence. Therefore, special domain is designed for introducing practical inlet condition
The solid model of honeycomb is meshed in GAMBIT 2.3.16. As shown in Fig. 2. A structured rectangular mesh is used for the simulation with square honeycomb configuration. Governing equations for mass and momentum conservations for subsonic flow along with the equations for turbulence and porous flow are solved for the honeycomb using commercial CFD. RANS type RNG k-ε model is used for the turbulence modeling.
The separate domain created upstream of the honeycomb is provided with various inlet conditions to arrive at the disorderly motion at the exit, which should be given as an inlet to the honeycomb cells. This essentially simulates the more realistic case that the flow can enter into the honeycomb from any direction. Specifications of this inlet along with other necessary boundary conditions are mentioned here. Flow at the inlet of the honeycomb shall necessarily have turbulent and swirling motions. Therefore, in order to incorporate these requirements, a separate fluid domain is constructed.
Top and bottom circular faces are considered as inlet to this domain to get a flow field with higher magnitude of lateral velocity. This domain is provided with vertical and horizontal cylinders as an obstruction to the inlet to produce sufficient swirling at the exit of this section. A tetrahedral mesh as shown in Fig. 3 with tetrahedral elements is generated for this geometry. The number of nodes are 1,47,666. Three faces of this configuration are specified as inlets with velocity boundary conditions. Fluid velocity at these inlet faces has been so taken that averaged mean velocity at the outlet is 1 m/s, which is in the operational wind tunnel.
A pressure outlet boundary condition is used at exit of the settling chamber where pressure at outlet is set to zero for gauge pressure.
It is always possible to predict the entire flow field by meshing whole fluid domain; however simulation for the prediction of entire flow field using symmetry boundary condition. This approach reduces the mesh requirement and computational efforts. Therefore, symmetry boundary is used at the periphery of the computational domain.
All the solid boundaries in the computational domain are specified as viscous walls with no-slip wall boundary condition.
Turbulence intensity profile at the exit of turbulence model is shown in Fig. 4. This figure shows the turbulence intensity and which is maximum at the center (30%) and at the walls is around 16-18%, now this profile is incorporated inside the honeycomb as shown in Fig. 2, the profile of turbulence intensity comes out from the honeycomb is shown in Fig. 5. In this profile we can see that the turbulence intensity is reduced from 30% to 1.2% at the center and 16% to 3.5%, it means the honeycomb effectiveness is very high which is around 96%.