Simulation of laminar microchannel flows with realistic 3D surface roughness

Akbaş, Batuhan
Effects of flow development and surface roughness on the pressure drop characteristics of laminar liquid flowsinside microchannels are investigated numerically using OpenFOAM. Channels with square cross section of 500 μm×500 μmand length of 80 mm are studied. Top surface of the channels are artificially roughened using thespatial frequency methodto create 8 different roughness profiles. Scaling the relative roughness ({u1D700}) values of each profile to three different values (1.0, 2.5and5.0 %),a total of 24 cases are simulated in the Reynolds number range of 100-1500. Results are presented mainly as the variation of the apparent Poiseuille number (Poapp) with Reynolds number (Re). Flow fields are examined by considering velocity magnitude contours at different cross sections to understand the flow development effects, as well as streamlines. For 1% relative roughness, it is seen that {u1D700}is a representative parameter to understand the pressure drop behavior and other roughness characteristics do not affect the flow field considerably. All {u1D700}=1% cases show Poapp–Re variations that have the same behavior as that of the smooth channel theory, with maximum deviations being 3% at Re = 100 and 7% at Re = 1500. However, for {u1D700}=2.5% and 5% cases, {u1D700}is no longer enough to explain the effect of surface roughness on pressure drop. Although certain roughness profiles provide Poapp–Re variations that are similar to that of the smooth channel theory, other profiles resulted in completely different behaviors with deviations increasing with Re. It is seen that for {u1D700}=2.5% Re = 1500, depending on the surface profile, Poappcan deviate from the smooth theory as low as 15% and as high as 55%. These numbers increased to 37% and 151% for {u1D700}=5% cases.It is concluded that the misunderstandings and the conflicting results seen in the literature about pressure drop in microchannels can be related to the use of relative surface roughness as the only defining parameter of roughness or even completely omitting roughness effects, which is clearly not the case for relative surface roughness values higher than 1%. Other parameters such as number of peaks and waviness also need to be considered and reported.


Computational Modeling of the Effects of Viscous Dissipation on Polymer Melt Flow Behavior During Injection Molding Process in Plane Channels
Tutar, M.; Karakuş, Ali (2013-02-01)
The present finite volume method based fluid flow solutions investigate the boundary-layer flow and heat transfer characteristics of polymer melt flow in a rectangular plane channel in the presence of the effect of viscous dissipation and heat transfer by considering the viscosity and density variations in the flow. For different inflow velocity boundary conditions and the injection polymer melt temperatures, the viscous dissipation effects on the velocity and temperature distributions are studied extensive...
Simulation of crystallization and glass formation processes for binary Pd-Ag metal alloys
Kart, HH; Uludogan, M; Cagin, T; Tomak, Mehmet (2003-09-12)
Glass formation and crystallization process of Pd-Ag metallic alloys are investigated by means of molecular dynamics simulation. This simulation uses the quantum Sutton-Chen (Q-SC) potential to study structural and transport properties of Pd-Ag alloys. Cooling rates and concentration effects on the glass formation and crystallization of binary alloys considered in this work are investigated. Pd-Ag alloys show the glass structure at fast cooling rates while it crystallizes at slow cooling rates. Increment of...
Numerical Simulation of Rarefied Laminar Flow past a Circular Cylinder
Çelenligil, Mehmet Cevdet (2014-07-18)
Numerical simulations have been obtained for two-dimensional laminar flows past a circular cylinder in the transitional regime. Computations are performed using the direct simulation Monte Carlo method for Knudsen numbers of 0.02 and 0.2 and Mach numbers of 0.102 and 0.4. For these conditions, Reynolds number ranges from 0.626 to 24.63 and the flows are steady. Results show that separation occurs in the wake region for the flow with Mach number of 0.4 and Knudsen number of 0.02, but for the other eases flow...
Numerical simulation of transient turbulent flow in a heated pipe
Uygur, Ahmet Bilge; Selçuk, Nevin; Oymak, Olcay; Department of Chemical Engineering (2002)
A computational fluid dynamics (CFD) code based on direct numerical simulation (DNS) and the method of lines MOL approach developed previously for the solution of transient two-dimensional Navier-Stokes equations for turbulent, incompressible, internal, non-isothermal flows with constant wall temperature was applied to prediction of turbulent flow and temperature fields in flows dominated by forced convection in circular tubes with strong heating. Predictive ability of the code was tested by comparing its r...
Modeling of particle-resin suspension impregnation in compression resin transfer molding of particle-filled, continuous fiber reinforced composites
Sas, Hatice Sinem; Erdal Erdoğmuş, Merve (Springer Science and Business Media LLC, 2014-03-01)
A particle-resin suspension impregnation model is used for analyzing the mold filling process in compression resin transfer molding (CRTM) of particle-filled, continuous fiber composites. The model is based on Darcy flow coupled with particle filtration and is applicable to two-dimensional impregnation through isotropic/anisotropic fiber preforms. Comparisons with simple analytical solutions and experimental results from the literature were made to validate the numerical solution. Simulations showed that CR...
Citation Formats
B. Akbaş, “Simulation of laminar microchannel flows with realistic 3D surface roughness,” Thesis (M.S.) -- Graduate School of Natural and Applied Sciences. Mechanical Engineering., Middle East Technical University, 2019.