The hybrid grid implemented DSMC method used in 2D triangular micro cavity flows

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.In this study a new hybrid grid is implemented in a 2D DSMC solver to be used in 2D triangular micro cavity flows. Currently DSMC is the prominent method to analyze micro scale gas flows which are rarefied. Because of the computational cost, DSMC solvers are generally used in rarefied gas conditions in which continuum based solvers are useless. If the efficiency of DSMC solvers is improved, the application range of these solvers can be increased further where the continuum based solvers dominate. Indexing the particles according to their cells is one of the main steps in the DSMC method. Either the particles are traced cell-by-cell along their trajectories or coordinate transformation techniques are used in this step. The first option requires complex trigonometric operations and search algorithms which are computationally expensive. But it can be used in both structured and unstructured grids. Although the second option is computationally more efficient, it demands specially tailored structured grids which are more geometry dependent compared to the unstructured grids. Here it is shown that a novel hybrid grid structure can be used successfully in 2D DSMC solver to analyze triangular shaped lid-driven micro cavity flows. Hybrid grids used in this study are much less dependent of the geometry like unstructured grids. Additionally, hybrid grids like structured grids facilitate coordinate transformation techniques in order to increase the efficiency of the particle indexing step in the DSMC method

The effect of chemical treatment, fibre length, fibre content and injection moulding parameters to uv irradiation resistance of oil palm fibre reinforced composites

There are many types of polymer used in engineering materials expose to UV irradiation, such as automotive parts (car body, bumper, dashboard etc.) which can cause material degradation. Some polymers are used in pure polymer and some of them in composite material. This study has investigated the composite material degradation. In this study, polypropylene was used as a matrix of the composite material samples, while oil palm fibre as reinforcement. The effect of the fibre length, fibre content, fibre treatment, coupling agent and injection moulding parameter to ultraviolet (UV) light resistance of this composite and also the optimum setting of them were investigated. The UV resistance was examined via the change of mechanical properties after UV exposured in the UV accelerated weathering tester. The Linear regression models were generated for tensile strength, strain at maximum stress, break stress, break strain and Charpy Impact strength based on six different UV exposure time i.e.: 0, 96, 336, 504, 1008, and 1512 hours. The significance of the regression models were tested by Analysis Of Variance (ANOVA) and verified by two expose time i.e: 168 and 672 hours. This study found that all of the mechanical properties decrease after UV irradiated. The largest property decrease was break strain of the samples, which was decrease in the range 44.54% to 79.21% after 1512 hours UV irradiated. The lowest decrease was break stress in the range 27.38 % to 63.82%. It was also found in this study that fibre content, and UV irradiation time, significantly affect all properties. Coupling agent and alkali treatment significantly affect all properties except strain at maximum stress. Whereas fibre length and injection moulding parameter only significantly affect the Impact strength of the specimens. All of the regression models are significant which are signed by the p value of each of regression models were lower than 0.05. The equation for predicting the lifetime of UV exposed of oil palm fibre reinforced composite have been generated. The lifetime UV irradiated specimen can be predicted using this equation

Study on SPH Viscosity Term Formulations

For viscosity-dominated flows, the viscous effect plays a much more important role. Since the viscosity term in SPH-governing (Smoothed Particle Hydrodynamics) equations involves the discretization of a second-order derivative, its treatment could be much more challenging than that of a first-order derivative, such as the pressure gradient. The present paper summarizes a series of improved methods for modeling the second-order viscosity force term. By using a benchmark patch test, the numerical accuracy and efficiency of different approaches are evaluated under both uniform and non-uniform particle configurations. Then these viscosity force models are used to compute a documented lid-driven cavity flow and its interaction with a cylinder, from which the most recommended viscosity term formulation has been identified

A stability condition for turbulence model: From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale (EMMS) model for gas-solid systems, is formulated to close the dynamic constraint equations of turbulence, allowing the heterogeneous structural parameters of turbulence to be optimized. We call this model the `EMMS-based turbulence model', and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room. The numerical results show that the EMMS-based turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.Comment: 26 pages, 13 figures, 2 table

Noise radiated by an open cavity at low Mach number: Effect of the cavity oscillation mode

The present work focuses on the study of noise generation and radiation of an infinite open three-dimensional cavity at low Mach number with laminar upstream conditions that is of interest to understand noise generation mechanisms in wall-bounded separated flows. A particular feature of this configuration is the oscillatory mode: shear layer mode or wake mode. For the parameters considered in the present study it is seen that while in shear layer mode the flow shows a two-dimensional behavior, in the wake mode the flow is three-dimensional, resulting in significantly different sound sources. The influence of the acoustic feedback mechanism in the shear layer mode has also been investigated comparing the results between different momentum thickness values at the cavity inlet. This paper presents results of sound radiated by a three-dimensional infinite open cavity with aspect ratio L/D = 4 at Reynolds number based on the cavity depth of ReD = 1500 and Mach number of M = 0.15, both for shear layer (L/θ = 67) and wake (L/θ = 84) oscillation modes. To do so, Curle integral evaluated as a post-process of an incompressible solution will be used. The results are compared with the resulting Curle post-process of a two-dimensional incompressible simulationPeer ReviewedPostprint (author's final draft

Simplex space-time meshes in thermally coupled two-phase flow simulations of mold filling

The quality of plastic parts produced through injection molding depends on many factors. Especially during the filling stage, defects such as weld lines, burrs, or insufficient filling can occur. Numerical methods need to be employed to improve product quality by means of predicting and simulating the injection molding process. In the current work, a highly viscous incompressible non-isothermal two-phase flow is simulated, which takes place during the cavity filling. The injected melt exhibits a shear-thinning behavior, which is described by the Carreau-WLF model. Besides that, a novel discretization method is used in the context of 4D simplex space-time grids [2]. This method allows for local temporal refinement in the vicinity of, e.g., the evolving front of the melt [10]. Utilizing such an adaptive refinement can lead to locally improved numerical accuracy while maintaining the highest possible computational efficiency in the remaining of the domain. For demonstration purposes, a set of 2D and 3D benchmark cases, that involve the filling of various cavities with a distributor, are presented.Comment: 14 pages, 11 Figures, 4 Table