0
Research Papers: Fundamental Issues and Canonical Flows

# Effect of Microstructure Geometric Form on Surface Shear StressOPEN ACCESS

[+] Author and Article Information
Kaushik K. Rangharajan

Department of Mechanical and
Aerospace Engineering,
The Ohio State University,
Columbus, OH 43210
e-mail: rangharajan.1@osu.edu

Matthew J. Gerber

Department of Mechanical and
Aerospace Engineering,
The Ohio State University,
Columbus, OH 43210
e-mail: gerber211@g.ucla.edu

Shaurya Prakash

Department of Mechanical and
Aerospace Engineering,
The Ohio State University,
Columbus, OH 43210
e-mail: prakash.31@osu.edu

1K. K. Rangharajan and M. J. Gerber contributed equally to this work.

2Present address: Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, Los Angeles, CA 90095.

3Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received October 8, 2015; final manuscript received July 21, 2016; published online September 20, 2016. Assoc. Editor: Daniel Maynes.

J. Fluids Eng 139(1), 011201 (Sep 20, 2016) (8 pages) Paper No: FE-15-1728; doi: 10.1115/1.4034363 History: Received October 08, 2015; Revised July 21, 2016

## Abstract

Low Reynolds number flow of liquids over micron-sized structures and the control of subsequently induced shear stress are critical for the performance and functionality of many different microfluidic platforms that are extensively used in present day lab-on-a-chip (LOC) domains. However, the role of geometric form in systematically altering surface shear on these microstructures remains poorly understood. In this study, 36 microstructures of diverse geometry were chosen, and the resultant overall and facet shear stresses were systematically characterized as a function of Reynolds number to provide a theoretical basis to design microstructures for a wide array of applications. Through a set of detailed numerical calculations over a broad parametric space, it was found that the top facet (with respect to incident flow) of the noncylindrical microstructures experiences the largest surface shear stress. By systematically studying the variation of the physical dimensions of the microstructures and the angle of incident flow, a comprehensive regime map was developed for low to high surface shear structures and compared against the widely studied right circular cylinder in cross flow.

<>

## Introduction

Low Reynolds number flow of viscous fluids over solid surfaces in either a confined or an open environment is a ubiquitous phenomenon that can be commonly observed. In this flow regime with Reynolds number, Re < 100, and in many LOC practical applications with Re < 1, the fluid motion at all the physical solid–fluid interfaces is usually considered to follow the no-slip boundary condition [1]. Therefore, the resulting transverse, nonzero velocity gradients contribute toward fluid shear which has been successfully exploited in numerous practical applications over the years, including biosensors, advanced health care diagnostics, materials processes, and energy-related technology advancements [27].

The operational efficiency for these diverse applications is critically impacted by the ability to maintain specific regimes of shear stress on the functional surfaces within the underlying devices. For example, lower surface shear was beneficial in improving the functionality of microfluidic chips with enhanced capture yield of circulating tumor cells [3], while in the case of artificial heart valves, reducing high-shear-induced coagulation of blood on the valve surface minimized the risk of thrombogenesis and thromboembolism [4]. For energy-related applications, lower surface shear increased microbial fuel cell sensitivity toward detecting Cu(II) toxicity [2].

On the other hand, flows with higher shear rate were beneficial in inducing microbial enrichment of microbial fuel cell anodes, resulting in a threefold increase in power output [5]. Furthermore, customized microfluidic cell culture chips were shown to exhibit superior gene expression under increased fluid shear [8]. However, in most low Re number (Re < 100) flow applications, shear stress is manipulated solely by modulating the flow rate due to the ease in controlling external pressure or potential gradients that act as common driving forces for these flows, especially in LOC environments. In contrast, an open question remains: What is the extent to which shear stress can be controlled by engineering passive structures that can systematically manipulate the shear without the need to actively adjust flow rate?

Studies on the effect of surface geometries on surface shear stress have mainly focused on considerations of surface roughness. This includes studies on developing a correlation between the effect of actual roughness to that of closely packed sand grains [9], an evaluation of the skin friction coefficient and equivalent sand roughness data on various rough surfaces [10], and an analysis on randomly placed, nonuniform, three-dimensional roughness with irregular geometry and arrangement [11]. Canonical flows over solid surfaces have predominantly focused on cylinders, which have been extensively studied for shear, both numerically and experimentally, for many years, and comprehensive reviews exist [1215]. However, studies that investigate low Reynolds number flows (Re < 100) for surface shear on noncylindrical, three-dimensional structures are scarce—many studies that consider noncylindrical forms are in two-dimensional flows or focus on the fluid flow and do not mention surface shear effects. Indeed, the study of surface shear is limited to specific ad hoc applications or targeted operating conditions [16,17]. Some simple radial geometries such as cylindrical disks and spheroids have been investigated and have shown that high-pressure stagnation zones correspond to areas of low surface shear stress, but no conclusions were formed regarding the relationship between geometric form and surface shear stress [18]. Similarly, in a computational fluid dynamics study that used images to align a competitive swimmer with the flow field, it was found that larger surface shear stresses were observed on areas of the body that presented a complex surface geometry, such as the head, shoulders, and heels [19].

Due to the scarcity of studies that investigate surface shear on distinguishable three-dimensional structural forms in low Reynolds number flows, no generalizations exist on how the surface structure actually manipulates or modulates shear stress under these conditions. Therefore, no generalizations exist on how the surface structure actually manipulates or modulates shear stress. However, in many applications, the geometry and layout of the physical structure in the flow path are critical for the device or system operation [20,21]. For example, many bioprocessing devices and bioreactors rely on the integrity of a biofilm adhered to their functional surfaces, and past works [22,23] have shown the dependence of biofilm morphology on surface shear stress. Surface structure geometry in several practical applications is typically chosen based on considerations such as ease of fabrication [24] or on what materials are commercially available [25] without much design focus on the flow–structure interaction, which actually governs the operation of these devices. Therefore, the resulting systems are inherently unoptimized for intended, shear-dependent functions because specific effects of geometry on surface shear stress were not explicitly considered.

Therefore, the purpose of this work is to provide a systematic analysis for low Reynolds number flow of water (a common working fluid in many LOC flows) over various microstructures on a flat, solid surface and subsequently characterize the induced fluid shear stress on the microstructure surface. Specifically, a regime map that directly correlates various microstructure shapes to regimes of surface shear stress has been developed. The regime map comprises a wide array of microstructures, such as cones, pyramids, rectilinear prisms, and other forms commonly found in engineering applications to provide a direct correlation between structure geometry, inlet flow conditions, and surface shear stress. Such regime maps can potentially provide starting point data to engineer higher operational efficiencies for applications due to the ability to now explicitly incorporate the geometric dependence of shear stress. In addition, the contribution of various microstructure facets (front, rear, side, and top of a structure) with respect to incoming flow toward overall fluid shear was also quantified as the structure was exposed to fully developed low Reynolds number flows. It is important to recognize that the focus of this work is on the study of external flows over microstructures where the wall constraints can be considered negligible.

## Methodology

In this paper, 36 geometric forms, each with the same height of 100 μm, were heuristically chosen based on common geometries found in a variety of engineering applications (as discussed in Sec. 1) and compared for relative shear stress, based on their detailed facet geometries with respect to incident or incoming flow (Fig. 1). The total wetted surface area for all the geometric shapes is plotted in Fig. S1, which is available under the “Supplemental Materials” tab on the ASME Digital Collection. Velocity gradients on individual facets are discussed for a cube, a common geometric form as a representative case. Using a cube of edge length of 100 μm as a base structure, the surface areas, edge lengths, and facet angles were systematically altered in incremental steps to study the effect of gradual changes in geometrical form over the base structure. In order to observe trends in surface shear stress based on the physical orientation of a structure with respect to the incident flow, a right pyramid and a cube were altered in systematic steps through one full rotation as observed with respect to the incident flow.

All the structures were modeled in a three-dimensional computational domain to simulate external flows over the structures' surfaces (see Fig. S2, which is available under the “Supplemental Materials” tab on the ASME Digital Collection). Initial coarse meshes were successively refined until mesh-independent solutions were achieved [26]. To test against the presence of numerical artifacts in the computed solutions [27], each simulation was performed with two iterative numerical methods (GMRES [28] and Bi-CGSTAB [29]) along with two methods of mesh generation Delaunay [30] and advancing front [31]. The results from the methods were compared to confirm agreement of their solutions to within 1% [27]. To eliminate the effect of singularities in the models [32], geometrically sharp features (e.g., edges and corners) were modified to fillets (radius of 5 μm) [32], and the surface shear rate was evaluated at least 3–5 grid elements [32], corresponding to a length of 5 μm away from such features. It is worth noting that the critical dimension is 100 μm for all the features of interest (FOI) (as discussed later), and therefore, the presence of a 5 μm fillet at the edges and corners of structures did not yield any significant changes in the reported results. All the governing equations were solved under steady-state, isothermal, and laminar flow conditions with Reynolds numbers of 0.001, 0.1, and 100 and with water modeled as an incompressible, Newtonian working fluid [33]. Calculation of Re was based on a microstructure critical dimension of 100 μm. The steady-state continuity and Navier–Stokes equations for these conditions reduce to Display Formula

(1)$∇ ⋅ V=0$
Display Formula
(2)$ρ V⋅∇ V=−∇ p+μ ∇ 2V$
Display Formula
(3)$Re=ρ uin hμ$

where $∇$ denotes the gradient operator, $V$ is the fluid velocity vector, $uin$ is the inlet velocity, $ρ$ is the fluid density, $p$ is the fluid pressure, h is the microstructure height (100 μm for all the structures), and $μ$ is the fluid viscosity of water at room temperature. All the numerical calculations were performed with comsol multiphysics v4.4 [34]. Average surface shear stress, $τ¯s$, was obtained from Newton's law of viscosity [33] by determining local shear rate magnitude, $γ˙$ from each grid point and averaging over the domain. The shear rate magnitude is calculated as Display Formula

(4)$γ˙=12Γ.:Γ. Γ˙=∇V+(∇V)T$

where: is the tensor contraction operator, and the superscript $T$ denotes the matrix transpose. Throughout this paper, “front facet” refers to the structural component facing the incident flow. Additionally, the overall and facet shear stress of all the microstructures discussed in this paper were nondimensionalized with respect to a reference case explicitly listed in Sec. 3. All the models were solved using the supercomputing cluster at the Ohio Supercomputer Center that employs 8328 HP Intel Xeon x5650 central processing units with 12 cores and 48 GB of memory per HP SL390 G7 node. The numerical model was validated with the past work [3], and good agreement was found on the spatial variation of velocity around the cylinders and the x-component of shear stress on the periphery of cylinder, as shown in Fig. S3, which is available under the “Supplemental Materials” tab on the ASME Digital Collection.

## Results and Discussion

###### Regime Map.

The 36 distinct geometric shapes considered in this study are divided into four categories based on traditional geometric definitions, namely: (1) rectilinear prisms (eight shapes with noncircular top), (2) radial prisms (11 shapes with top facet having a radial profile), (3) nonvertical prisms (ten shapes with top facet having varying heights), and (4) apex structures (seven shapes with no-top facet) as summarized in Fig. 1.

For a given Re, the overall shear stress of a microstructure, $(τs)all$, was divided by the overall shear stress of a cylinder and reported as $(τ¯s)all$, thus nondimensionalizing the stress and providing a direct comparison to the cylinder, which has been one of the most common structures used in a variety of flow configurations as discussed in Sec. 1. Moreover, comparing the shear stress of the 36 cases chosen here with respect to the cylinder also facilitates ready comparison with the published data [35] and provides easy visualization to existing engineering applications that employ cylinders. Due to the nondimensionalization, $(τ¯s)all$ =1 for a cylinder as shown in the regime map (Fig. 2). $(τs)all$ for a 100 μm diameter (and height) cylinder was calculated to be 7.5 × 10−5 N/m2 (Re = 0.001), 7.5 × 10−3 N/m2 (Re = 0.1), and 16.1 N/m2 (Re = 100). As the nondimensionalized shear stress for all the microstructures reported in this study was similar for Re = 0.001 and 0.1, results for Re = 0.1 will be discussed as a representative condition. The nondimensionalized overall shear stress, $(τ¯s)all$, for each structure (Fig. 2) as well as differences in front, top, side, and rear facets were examined (Figs. 3 and 4). The three-dimensional orientation of all the microstructures with respect to a global Cartesian coordinate frame and also the direction of incident flow are explicitly shown in Fig. 2.

In Fig. 2, a shear stress regime map for the overall shear stress, $(τ¯s)all$, as a function of geometry is presented for Re = 0.1 (with Re = 0.001 being similar to Re = 0.1) and 100. Rectilinear prisms have $(τ¯s)all$ varying in the range of 0.73–0.88 for Re = 0.1, suggesting that the overall shear on a cubic prism is only 73% of the shear stress seen by a cylinder under similar flow conditions and critical dimensions. In rectilinear prisms, when the fluid contacts the front facet, the velocity gradient is larger but limited only to the edges and small everywhere else, in contrast to curvilinear prisms, where the velocity drop is extended out over nearly the entire surface area of the structure, resulting in higher surface shear stress. This is the first of many results that suggest a fluid flow over sharply angled features (such as rectilinear prisms) will, in general, result in lower surface shear stress when compared to curvilinear surfaces. At Re = 100, rectilinear prisms have $(τ¯s)all$ varying over a broader range of 0.62–0.97. Nonvertical prisms exhibited an even wider spread in the magnitude of $(τ¯s)all$ (0.38–0.80 compared to a cylinder) for Re = 0.1, but was confined between 0.37 and 0.72 at Re = 100. Radial prisms exhibited larger regimes of $(τ¯s)all$ in magnitude (0.82–1.1 at Re = 0.1 and 0.74–1.2 at Re = 100); however, apex forms exhibited lower regimes of $(τ¯s)all$ with respect to the rectilinear and radial prisms with the exception of the cone at Re = 0.1, where $(τ¯s)all$ was estimated to be 0.92 or nearly the same as the cylinder, despite not having a top facet. Apex structures due to minimal cross-sectional area (normal to the flow field) would allow for increased momentum of the fluid flow over the entire structure reducing the velocity gradients over their surface. If true, then structures 19, 23, 25, and 26 (all forms with cross-sectional area less than 1 μm2) should also experience less surface shear stress in the Re = 100 case, which is demonstrated by the results.

It is important to note that $(τ¯s)all$ is an area-averaged value of shear stress over all the exposed facets as shown in Eq. (5). Therefore, though individual facets exhibit varied magnitudes of shear stress, given by $(τ¯s)i$ which is independent of the facet area $Ai$, the overall shear stress, $(τ¯s)all$, is influenced by both $(τ¯s)i$ and $Ai$Display Formula

(5)$(τ¯s)all=∑i(τ¯s)iAi∑iAi$

The regime map study thus clearly indicates that although the overall shear stress experienced by a microstructure is similar at Re = 0.001 and 0.1, further extrapolation of shear stress at higher Re (100) is not obvious and therefore was explicitly evaluated (Fig. 2). In short, the regime map provides data trends for low Reynolds number flow that can be used for selection of geometric categories for a variety of applications as discussed in Sec. 1.

###### Velocity Gradients and Impact on Shear Stress.

Nondimensionalized shear stress experienced by individual facets (front, rear, side, and top facets) for all the 36 microstructures is plotted at Re = 0.1, 100 in Figs. 3, and 4, respectively. Shear stress experienced at each facet was nondimensionalized with $(τs)all$ of a cylinder at the respective Re, as discussed in Sec. 3.1.

In Fig. 3, critical facet stress $(τ¯s)i$ of 0.5 (after nondimensionalization, i.e., 50% of the overall shear on a cylinder) was benchmarked for front and rear facets, as the transition between high $((τ¯s)i>0.5)$ and low $((τ¯s)i≤0.5)$ shear stress structures. Similarly, as the top facet was found to experience higher shear stress in comparison to all other facets, a critical facet stress $(τ¯s)i$ of 1.5 was chosen to differentiate between high shear and low shear structures. The $(τ¯s)i$ of side facets, though greater than 0.5, varied over a narrow range for all the microstructures (average $(τ¯s)i≈$ 0.79 ± 0.13). Similarly, critical shear stresses of 1 for the front facet, 0.3 for the rear facet, 1 for the side facet, and 1.5 for the top facet were chosen to differentiate high shear and low shear structures at Re = 100, as shown in Fig. 4. This classification between “high shear” and “low shear” structures in Figs. 3 and 4 was implemented to aid in the selection of microstructures for shear-dependent applications. Nondimensionalized facet shear stress as a function of Re is plotted for all the microstructures in Fig. S4, which is available under the “Supplemental Materials” tab on the ASME Digital Collection.

Next, the role of velocity gradients in contributing to the overall shear was evaluated by considering a cube of side 100 μm as a representative structure. The velocity gradients were nondimensionalized with respect to the overall shear rate of a cylinder at Re = 0.1 (7.89 s−1). The nondimensionalized velocity gradients that contribute to shear in the X, Y, and Z directions across the top facet were calculated and are shown in Fig. 5 for a cube of side 100 μm at Re = 0.1 as a representative case. As discussed in Sec. 2, data were plotted beginning at 5 μm from each edge in order to eliminate the effects of singularities. Therefore, each cross section of the cube reported in Fig. 5 has an area of 90 μm × 90 μm instead of 100 μm × 100 μm. The effect of this strategy is discussed and shown to be valid in the supporting information section (Fig. S5, which is available under the “Supplemental Materials” tab on the ASME Digital Collection).

The velocity vector is given by $V=u x̂+vŷ+w ẑ$, where $x̂$, $ŷ$, and $ẑ$ are the unit vectors along the coordinate axes (shown in Fig. 5). Since the top facet is located in the X–Y plane, and due to the symmetric structure of the cube, the nondimensionalized velocity gradient is expected to dominate along Z (velocity gradient∼O(10 deg)) and was found to be negligible across the X,Y direction (velocity gradient ∼ O(10−15)) as shown in Fig. 5. The gradient of the velocity along Z is attributed to the no-slip boundary condition at the top facet, which leads to the development of a transverse boundary layer. This transverse boundary layer and the resulting velocity gradient are the primary causes of the surface shear stress calculated on the facet. Since the incident flow is in the X direction, $|u|$ (magnitude of $u$) is greater than $|v|$, and therefore, the magnitude of $∂v/∂z$ is smaller by an order of magnitude compared to $∂u/∂z$, as also shown in Fig. 5. In addition, contribution of $∂v/∂z$ toward $(τ¯s)top$ can be ignored when averaged, given the symmetry involved in the structure. Similarly, the average of $∂w/∂z$ across the top facet was estimated to be zero (not shown) and therefore does not contribute toward the shear stress experienced by the top facet. From Fig. 5, $∂u/∂z$ exhibits the maximum magnitude compared to all other velocity gradients that contribute toward shear stress. The negligible magnitude of all other velocity gradients, $(∂v/∂x), (∂w/∂x), (∂u/∂y),(∂w/∂y)$ ∼ O(10−15) as observed in Fig. 5 suggests that the variation of velocity $V$ in a direction normal to a particular facet predominantly influences the magnitude of shear stress at that facet. Therefore, the surface shear stress is largest along the edges of the facet (∼1.5 times larger than the average over the facet) and peaks near the corners of the facet (∼2.5 times larger than the average over the facet). Thus, one would hypothesize that structures having multiple edges and corners should experience higher overall shear stress, but as can be seen from Figs. 3 and 4, this is not the case. This discrepancy arises as the actual surface area of edges and corners is negligible compared to the overall facet area and therefore contributes minimally (<2%) to the overall average shear stress.

Figure 6 shows the nondimensionalized velocity gradients that primarily influence $(τ¯s)i$ plotted across the front ($−x̂$ being the normal unit vector), side ($−ŷ$ being the normal unit vector), and rear ($x̂$ being the normal unit vector) facets for a cube with edge length 100 μm at Re = 0.1. Orientation of the cube with respect to the three-dimensional axis and incoming flow is as shown in Fig. 5. Since $|u|$ > $| v,w |$, the gradient of $v, w$ (experienced only at the front and rear facets) is expected to be less than the gradient of $u$ (experienced only at the side facet). Therefore, as shown in Fig. 6, the shear rate experienced by the side facet is about 2.5 times greater than the front and rear facet shear stresses $((τ¯s)side >(τ¯s)front,(τ¯s)rear)$ for a cube. This trend is in agreement with the results reported for Re = 0.1 for a cube and suggests that, in any structure with facets of similar orientation and shape to those of a cube, the side facets will strongly contribute to the overall shear stress for the structure (and referring to Figs. 2 and 3, this appears to be the case: structures with overall form similar to that of the cube experience comparable magnitudes of surface shear stress). However, for noncubic geometries, it is clear from the regime map (Fig. 2) that the overall geometric form affects the shear stress on individual facets differently (Figs. 3 and 4) and therefore justifies the need for such a regime map in designing next-generation LOC devices with embedded microstructures [3,3639].

###### Parametric Manipulation of Physical Dimensions.

As shown in Fig. 7, four physical dimensions of a cube were incremented to progressively alter its shape, to achieve the final structures shown by the third image in each sequence at the top of each figure panel. The shear stress on each individual facet in this discussion is now nondimensionalized with respect to $(τs)front$ or the shear stress on the front facet of a cube, i.e., the facet of the geometry facing the incoming flow experienced by an unaltered cube at Re = 100, and reported as $(τ¯s)i$ in Figs. 7(a)7(d), to facilitate easy comparison to the unaltered case. For Re = 100, $(τ¯s)front$ of the unaltered cube was calculated to be 12.5 N/m2.

In Fig. 7(a), the results of incrementing 0 deg ≤  $θp$  ≤ 26.6 deg are shown ($θp$ is the angle of tilt for each side of the cube), which transforms a cube into a right pyramid as it increases. It was found that the shear stress at the top facet, $(τ¯s)top$, increased with a quadratic dependence on $θp$ ($R2$  = 0.98, where $R2$ is the coefficient of determination) for the systematic translation of a cube to a pyramid. This scaling is expected because as the flat faces of the cube were systematically altered to approach triangular cross sections, the surface area of the top facet continued to decrease.

Moreover, the front facet exhibits a strong linear correlation ($r$ =−0.91, where $r$ is the linear correlation coefficient) between decreasing $(τ¯s)front$ and increasing $θp$. Though the front and rear facets have the same area, the shear stress experienced by the rear facet is an order of magnitude lower than that of the front facet. Additionally, $(τ¯s)rear$ was found to be invariant with a change in surface area, unlike $(τ¯s)front$, which implies that the shear stress experienced by each facet is strongly influenced by the overall geometric form and Re. At $θp$  = 0, all the facets of a cube have equal surface area and from Eq. (4), the overall stress is dominated by $(τ¯s)top$, which has the maximum magnitude, which is in agreement with the discussion in Sec. 3.2. However, at $θp$  = 26.6 deg, the top surface vanishes (or collapses to a point as the pyramid tip) and therefore does not contribute to overall stress. Unlike the side and rear facets whose shear stress is invariant to $θp$, $(τ¯s)front$ decreases linearly and therefore from Eq. (4), $(τ¯s)all$ decreases linearly as observed in Fig. 7(a), suggesting that the front facet is instrumental in dictating $(τ¯s)all$ with an increase in $θp$.

In Fig. 7(b), results are shown for incrementing 0 ≤  $Rv$  ≤ 50 μm, i.e., increasing cube edge curvature ($Rv$) to round-out the cube and eventually reach a cylinder. The results for the shear stress on the front curves suggest that the “sharper” a geometric feature is (in this case, smaller values of $Rv$), the more surface shear stress it will exhibit, as expected in Ref. [40] and from previous discussions, above. Also, with increase in $Rv$, the surface area of the top facet decreases by 21.5% from a square ($Rv$  = 0 μm) to a circle ($Rv$  = 50 μm), compared to a 100% decrease in the case discussed in Fig. 7(a). In addition, the top facet experiences the highest transverse velocity gradients when $Rv$  > 25 μm, as indicated by a higher value of $(τ¯s)top$ in comparison to other facets in Fig. 7(b). Therefore, coupled with a significant contributed area and shear stress toward the estimation of $(τ¯s)all$ ($(τ¯s)top$ contributes 36.0–41.5% to overall shear with increase in $Rv$), the top facet dictates the variation of $(τ¯s)all$ with $Rv$.

In Fig. 7(c), the results of incrementing the angle of tilt, $θw$, of the front face of a cube are shown. Therefore, the area of the front and side facets varies with $θw$, and the area of the rear facet remains fixed as shown in Fig. 7(c). Again, $(τ¯s)top$ shows a quadratic increase (R2 = 0.97) for the evaluated values of $θw$. Despite an overall increase in $θw$ from 0 deg to 45 deg, thereby increasing the front facet area by 41%, $(τ¯s)front$ only increases by 0.86%. In comparison, the percent contribution of $(τ¯s)sides$ to $(τ¯s)all$ is 16.9 ± 0.25% for all the values of $θw$, suggesting $(τ¯s)sides$ is independent of $θw$. Together, these results suggest that the top and front facets dictate $(τ¯s)all$ for all the values of $θw$.

In Fig. 7(d), results are shown for incrementing 0 ≤  $Rt$  ≤ 50 μm, i.e., increasing the curvature ($Rt$) of the cube's top facet. As $Rt$ is increased, the height (and therefore area) of the front facet was systematically reduced to maintain the total height of the structure at 100 μm. The results for the front curve and the top facets confirmed previous results that sharper geometric features result in larger values of $(τ¯s)i$ along expected trends [40]. $(τ¯s)all$ exhibits only a 3.4% overall decrease, while $(τ¯s)front$ shows a strong linear correlation (r=−0.99) with $Rt$, decreasing 41.5% overall, and remaining within 21.5% of $(τ¯s)all$. It is important to note that the area of the newly formed front curve increases as $Rt$ increases.

A common trend observed in all the cases was that the magnitude of shear stress on the rear facet, $(τ¯s)rear$, remains small (∼3%) and is mostly unaltered by the changes to structure morphology, in agreement with the regime map (Re = 100, Fig. 2(b)). Similarly, minimal changes in the value of $(τ¯s)all$ suggest that in most low Reynolds number based flows, the surfaces exposed to the incoming fluid (typically the front and top facets) dictates the performance characteristics. Also, a decrease in the area of the front facet when transformed from a cube (Figs. 7(a), 7(b), and 7(d)) results in a decrease in $(τ¯s)front$. However, decreasing the surface area of the top facet (Figs. 7(a)7(d)) leads to a quadratic increase in the magnitude of $(τ¯s)top$ most likely due to the increase in the amount of area exposed to a large velocity gradient near the edges of the facet. Overall, the results suggest that to produce a structure with larger overall surface shear stress, smaller top facets, rounder top facets, and flatter front facets normal to the incident flow are preferred, but of most importance is the form of the top facet. Thus, the trends previously observed and discussed in relation to Fig. 2 are further illuminated. In conclusion, the overall results in Fig. 7 suggest that the shear stress experienced by various facets is strongly influenced by the overall geometric form of microstructures.

###### Effect of Altering Angle of Rotation.

The effect of changing the angle of orientation with respect to the incident flow for a pyramid, $ϕp$, and cube, $ϕc$, was investigated (Fig. 8). The axis of rotation for both the pyramid and cube is shown in Fig. 8. To facilitate comparison between the rotated structures and their nonrotated forms, the reported values of shear stress were nondimensionalized with $(τs)front$ at $ϕp$  = 0 deg and $ϕc$  = 0 deg.

In Fig. 8(a), the pyramid's feature of interest (FOI = front facet at $ϕp$  = 0 deg) exhibits a global minimum at $ϕp$  = 180 deg. The maximum values of $(τ¯s)FOI$ occur at $ϕp$  = 60 deg and 300 deg. Therefore, the magnitude of shear stress at the front is maximum when angled at ± 60 deg from the incident flow. From Fig. 2, we see that, in general, structures 5 and 8 (which have angled, square front facets) experience larger surface shear stress than the cube, whose front facet is normal to the flow. This result indicates that the surface shear stress on a structural form can be manipulated and achieved purely by changing its orientation with respect to the incoming flow. The shear stress exhibits a smooth, sinusoidal distribution of both $(τ¯s)all$ and $(τ¯s)FOI$ as $ϕp$ is varied, i.e., the overall and facet shear stress is symmetric as the orientation of the structure with respect to incoming flow is varied. Figure 8(b) shows the effect of changing $ϕc$ for a cube. The maximum values of $(τ¯s)FOI$ occur at $(τ¯s)FOI$  = 50 deg and 310 deg (i.e., ±50 deg with respect to the incident flow, which agrees closely with the results for the pyramid). For the top facet, the maximum values of $(τ¯s)top$ occur at $ϕc$  = 45 deg, 135 deg, 225 deg, and 315 deg, corresponding to orientations that exhibit a maximum gradient in velocity for the same top facet surface area.

## Summary and Conclusions

A systematic numerical analysis for the 36 diverse geometries was analyzed to estimate overall and individual facet shear stress at Reynolds numbers (Re) of 0.001, 0.1, and 100, which spans five orders of magnitude for relatively low Re flows typically seen in many viscous flow and LOC applications. The structures were divided into four categories based on the geometry of the top facet, namely: rectilinear prisms, radial prisms, nonvertical prisms, and apex structures.

Overall shear stress of a microstructure was nondimensionalized with the overall shear stress of a cylinder at a given Re and reported in a regime map. The results indicate that the nondimensionalized facet and overall shear stress for all the 36 structures were found to be the same at Re = 0.001 and 0.1. However, as Re was increased to 100, facet and overall shear stress was found to vary in comparison to Re = 0.1. Since the low Reynolds number flow is incident in one-direction, the magnitude of the transverse velocity components is lower compared to the axial component (which is along the direction of flow). Therefore, facets that include the gradient of axial velocity toward the estimation of shear rate experience a higher magnitude of shear stress compared to facets that take into account the transverse component of velocity. In addition, by changing the angle of rotation, critical angles were found for a cube and pyramid for which the surface shear stress was maximum in magnitude. The results of this work are expected to provide a broad basis for choosing microstructures that are essential for shear-dependent applications.

## Acknowledgements

The computational facilities at the Ohio Supercomputing Center (OSC) are acknowledged for the support. The U.S. Department of Energy (DOE) is acknowledged for the partial funding support for the personnel through ARPA-E under Grant No. DE-AR0000282. The discussions with Tong Lin are also acknowledged.

## Nomenclature

• $Ai$ =

area of facet of interest

• $p$ =

fluid pressure (components vary spatially)

• $r$ =

linear correlation coefficient

• $Rt$ =

radius of curvature of top facet

• $Rv$ =

• $R2$ =

coefficient of determination

• Re =

Reynolds number

• $u$ =

X component of $V$

• $uin$ =

inlet velocity

• $v$ =

Y component of $V$

• $V$ =

Eulerian fluid velocity (components vary spatially)

• $w$ =

Z component of $V$

• X =

global Cartesian coordinate in the direction of inlet flow (see Fig. 2)

• $x̂$ =

unit vector along X

• Y =

global Cartesian coordinate perpendicular to the flow

• $ŷ$ =

unit vector along Y

• Z =

global Cartesian coordinate perpendicular to the flow

• $ẑ$ =

unit vector along Z

• $γ.$ =

average shear rate

• $θp$ =

angle of inclination of each facet of cube

• $θw$ =

angle of incidence of front facet normal to the incident flow

• $μ$ =

fluid viscosity (water in this study)

• $ρ$ =

fluid density (water in this study)

• $(τs)i$ =

average shear stress of facet i (N/m2)

• $(τs)i$ =

shear stress of facet i (N/m2)

• $(τ¯s)i$ =

nondimensionalized average shear stress of facet i

• $(τ¯s)all$ =

nondimensionalized overall shear stress of microstructure

• $(τ¯s)front$ =

nondimensionalized average shear stress of front facet

• $(τ¯s)rear$ =

nondimensionalized average shear stress of rear facet

• $(τ¯s)side$ =

nondimensionalized average shear stress of side facet

• $(τ¯s)top$ =

nondimensionalized average shear stress of top facet

• $ϕc$ =

angle of rotation of cube

• $ϕp$ =

angle of rotation of pyramid

## References

Prakash, S. , and Yeom, J. , 2014, Nanofluidics and Microfluidics: Systems and Applications, Elsevier, New York.
Shen, Y. , Wang, M. , Chang, I. S. , and Ng, H. Y. , 2013, “ Effect of Shear Rate on the Response of Microbial Fuel Cell Toxicity Sensor to Cu(II),” Bioresour. Technol., 136, pp. 707–710. [PubMed]
Nagrath, S. , Sequist, L. V. , Maheswaran, S. , Bell, D. W. , Irimia, D. , Ulkus, L. , Smith, M. R. , Kwak, E. L. , Digumarthy, S. , Muzikansky, A. , Ryan, P. , Balis, U. J. , Tompkins, R. G. , Haber, D. A. , and Toner, M. , 2007, “ Isolation of Rare Circulating Tumour Cells in Cancer Patients by Microchip Technology,” Nature, 450(7173), pp. 1235–1239. [PubMed]
Dasi, L. P. , Murphy, D. W. , Glezer, A. , and Yoganathan, A. P. , 2008, “ Passive Flow Control of Bileaflet Mechanical Heart Valve Leakage Flow,” J. Biomech., 41(6), pp. 1166–1173. [PubMed]
Pham, H. T. , Boon, N. , Aelterman, P. , Clauwaert, P. , De Schamphelaire, L. , Van Oostveldt, P. , Verbeken, K. , Rabaey, K. , and Verstraete, W. , 2008, “ High Shear Enrichment Improves the Performance of the Anodophilic Microbial Consortium in a Microbial Fuel Cell,” Microbiol. Biotechnol., 1(6), pp. 487–496.
Prakash, S. , Pinti, M. , and Bhushan, B. , 2012, “ Theory, Fabrication and Applications of Microfluidic and Nanofluidic Biosensors,” Philos. Trans. R. Soc. London, Ser. A, 370(1967), pp. 2269–2303.
Kellie, B. M. , Silleck, A. C. , Bellman, K. , Snodgrass, R. , and Prakash, S. , 2013, “ Deposition of Few-Layered Graphene in a Microcombustor on Copper and Nickel Substrates,” RSC Adv., 3(19), pp. 7100–7105.
Hattori, K. , Munehira, Y. , Kobayashi, H. , Satoh, T. , Sugiura, S. , and Kanamori, T. , 2014, “ Microfluidic Perfusion Culture Chip Providing Different Strengths of Shear Stress for Analysis of Vascular Endothelial Function,” J. Biosci. Bioeng., 118(3), pp. 327–332. [PubMed]
Sigal, A. , and Danberg, J. E. , 1990, “ New Correlation of Roughness Density Effect on the Turbulent Boundary Layer,” AIAA J., 28(3), pp. 554–556.
Coleman, H. W. , Hodge, B. K. , and Taylor, R. P. , 1984, “ A Re-Evaluation of Schlichting's Surface Roughness Experiment,” ASME J. Fluids Eng., 106(1), pp. 60–65.
van Rij, J. A. , Belnap, B. J. , and Ligrani, P. M. , 2002, “ Analysis and Experiments on Three-Dimensional, Irregular Surface Roughness,” ASME J. Fluids Eng., 124(3), pp. 671–677.
Zdravkovich, M. M. , 1997, Flow Around Circular Cylinders, Vol. 1: Fundamentals, Oxford University Press, Oxford, UK.
Beaudan, B. , and Moin, P. , 1994, “ Numerical Experiments on the Flow Past a Circular Cylinder at Sub-Critical Reynolds-Number,” Report No. TF-62, Center for Turbulence Research, Stanford University, Stanford, CA.
Williamson, C. H. K. , 1996, “ Vortex Dynamics in the Cylinder Wake,” Annu. Rev. Fluid Mech., 28(1), pp. 477–539.
Rosetti, G. F. , Vaz, G. , and Fujarra, A. L. C. , 2012, “ URANS Calculations for Smooth Circular Cylinder Flow in a Wide Range of Reynolds Numbers: Solution Verification and Validation,” ASME J. Fluids Eng., 134(12), p. 121103.
Agimelen, O. S. , and Olmsted, P. D. , 2013, “ Apparent Fracture in Polymeric Fluids Under Step Shear,” Phys. Rev. Lett., 110(20), p. 204503. [PubMed]
Gao, T. , Hu, H. H. , and Castañeda, P. P. , 2012, “ Shape Dynamics and Rheology of Soft Elastic Particles in a Shear Flow,” Phys. Rev. Lett., 108(5), p. 058302. [PubMed]
Gutierrez, R. A. , and Crumpler, E. T. , 2007, “ Potential Effect of Geometry on Wall Shear Stress Distribution Across Scaffold Surfaces,” Ann. Biomed. Eng., 36(1), pp. 77–85. [PubMed]
Popa, C. V. , Zaidi, H. , Arfaoui, A. , Polidori, G. , Taiar, R. , and Fohanno, S. , 2011, “ Analysis of Wall Shear Stress Around a Competitive Swimmer Using 3D Navier–Stokes Equations in CFD,” Acta Bioeng. Biomech., 13(1), pp. 3–11. [PubMed]
Choi, S. , Karp, J. M. , and Karnik, R. , 2012, “ Cell Sorting by Deterministic Cell Rolling,” Lab Chip, 12(8), pp. 1427–1430. [PubMed]
Loutherback, K. , Chou, K. , Newman, J. , Puchalla, J. , Austin, R. , and Sturm, J. , 2010, “ Improved Performance of Deterministic Lateral Displacement Arrays With Triangular Posts,” Microfluid. Nanofluid., 9(6), pp. 1143–1149.
Huang, Z. , McLamore, E. S. , Chuang, H. S. , Zhang, W. , Wereley, S. , Leon, J. L. C. , and Banks, M. K. , 2013, “ Shear-Induced Detachment of Biofilms From Hollow Fiber Silicone Membranes,” Biotechnol. Bioeng., 110(2), pp. 525–534. [PubMed]
Kim, J. , Kim, H.-S. , Han, S. , Lee, J.-Y. , Oh, J.-E. , Chung, S. , and Park, H.-D. , 2013, “ Hydrodynamic Effects on Bacterial Biofilm Development in a Microfluidic Environment,” Lab Chip, 13(10), pp. 1846–1849. [PubMed]
Langhamer, O. , Wilhelmsson, D. , and Engström, J. , 2009, “ Artificial Reef Effect and Fouling Impacts on Offshore Wave Power Foundations and Buoys—A Pilot Study,” Estuarine Coastal Shelf Sci., 82(3), pp. 426–432.
Logan, B. , Cheng, S. , Watson, V. , and Estadt, G. , 2007, “ Graphite Fiber Brush Anodes for Increased Power Production in Air-Cathode Microbial Fuel Cells,” Environ. Sci. Technol., 41(9), pp. 3341–3346. [PubMed]
Patanakar, S. V. , 1980, Numerical Heat Transfer and Fluid Flow, McGraw-Hill, New York.
Sprittles, J. E. , and Shikhmurzaev, Y. D. , 2011, “ Viscous Flow in Domains With Corners: Numerical Artifacts, Their Origin and Removal,” Comput. Methods Appl. Mech. Eng., 200(9–12), pp. 1087–1099.
Saad, Y. , and Schultz, M. , 1986, “ GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems,” SIAM J. Sci. Stat. Comput., 7(3), pp. 856–869.
Van der Vorst, H. A. , 1992, “ Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems,” SIAM J. Sci. Stat. Comput., 13(2), pp. 631–644.
Meng, Q. , Thanh-Tung, C. , and Tiow-Seng, T. , 2013, “ Computing 2D Constrained Delaunay Triangulation Using the GPU,” IEEE Trans. Visual. Comput. Graphics, 19(5), pp. 736–748.
Löhner, R. , 2014, “ Recent Advances in Parallel Advancing Front Grid Generation,” Arch. Comput. Methods Eng., 21(2), pp. 127–140.
Sprittles, J. E. , and Shikhmurzaev, Y. D. , 2011, “ Viscous Flows in Corner Regions: Singularities and Hidden Eigensolutions,” Int. J. Numer. Methods Fluids, 65(4), pp. 372–382.
Kundu, P. K. , Cohen, I. M. , and Dowling, D. R. , 2012, Fluid Mechanics, 5th ed., Academic Press, Boston, MA, p. iii.
COMSOL, 2013, “ comsol Multiphysics Release Notes,” COMSOL Multiphysics, Burlington, MA, accessed May 2013,
Tamayol, A. , Yeom, J. , Akbari, M. , and Bahrami, M. , 2013, “ Low Reynolds Number Flows Across Ordered Arrays of Micro-Cylinders Embedded in a Rectangular Micro/Minichannel,” Int. J. Heat Mass Transfer, 58(1–2), pp. 420–426.
Chien-Yu, F. , Sheng-Yang, T. , Shih-Mo, Y. , Long, H. , Cheng-Hsien, L. , and Hwan-You, C. , 2014, “ A Microfluidic Chip With a U-Shaped Microstructure Array for Multicellular Spheroid Formation, Culturing and Analysis,” Biofabrication, 6(1), p. 015009. [PubMed]
Stone, H. A. , Stroock, A. D. , and Ajdari, A. , 2004, “ Engineering Flows in Small Devices,” Annu. Rev. Fluid Mech., 36(1), pp. 381–411.
Gogolides, E. , Ellinas, K. , and Tserepi, A. , 2015, “ Hierarchical Micro and Nano Structured, Hydrophilic, Superhydrophobic and Superoleophobic Surfaces Incorporated in Microfluidics, Microarrays and Lab on Chip Microsystems,” Microelectron. Eng., 132, pp. 135–155.
Sugioka, K. , Masuda, M. , Hongo, T. , Cheng, Y. , Shihoyama, K. , and Midorikawa, K. , 2004, “ Three-Dimensional Microfluidic Structure Embedded in Photostructurable Glass by Femtosecond Laser for Lab-on-Chip Applications,” Appl. Phys. A, 79(4–6), pp. 815–817.
Park, H. , Park, H. , and Kim, J. , 2013, “ A Numerical Study of the Effects of Superhydrophobic Surface on Skin-Friction Drag in Turbulent Channel Flow,” Phys. Fluids, 25(11), p. 110815.
View article in PDF format.

## References

Prakash, S. , and Yeom, J. , 2014, Nanofluidics and Microfluidics: Systems and Applications, Elsevier, New York.
Shen, Y. , Wang, M. , Chang, I. S. , and Ng, H. Y. , 2013, “ Effect of Shear Rate on the Response of Microbial Fuel Cell Toxicity Sensor to Cu(II),” Bioresour. Technol., 136, pp. 707–710. [PubMed]
Nagrath, S. , Sequist, L. V. , Maheswaran, S. , Bell, D. W. , Irimia, D. , Ulkus, L. , Smith, M. R. , Kwak, E. L. , Digumarthy, S. , Muzikansky, A. , Ryan, P. , Balis, U. J. , Tompkins, R. G. , Haber, D. A. , and Toner, M. , 2007, “ Isolation of Rare Circulating Tumour Cells in Cancer Patients by Microchip Technology,” Nature, 450(7173), pp. 1235–1239. [PubMed]
Dasi, L. P. , Murphy, D. W. , Glezer, A. , and Yoganathan, A. P. , 2008, “ Passive Flow Control of Bileaflet Mechanical Heart Valve Leakage Flow,” J. Biomech., 41(6), pp. 1166–1173. [PubMed]
Pham, H. T. , Boon, N. , Aelterman, P. , Clauwaert, P. , De Schamphelaire, L. , Van Oostveldt, P. , Verbeken, K. , Rabaey, K. , and Verstraete, W. , 2008, “ High Shear Enrichment Improves the Performance of the Anodophilic Microbial Consortium in a Microbial Fuel Cell,” Microbiol. Biotechnol., 1(6), pp. 487–496.
Prakash, S. , Pinti, M. , and Bhushan, B. , 2012, “ Theory, Fabrication and Applications of Microfluidic and Nanofluidic Biosensors,” Philos. Trans. R. Soc. London, Ser. A, 370(1967), pp. 2269–2303.
Kellie, B. M. , Silleck, A. C. , Bellman, K. , Snodgrass, R. , and Prakash, S. , 2013, “ Deposition of Few-Layered Graphene in a Microcombustor on Copper and Nickel Substrates,” RSC Adv., 3(19), pp. 7100–7105.
Hattori, K. , Munehira, Y. , Kobayashi, H. , Satoh, T. , Sugiura, S. , and Kanamori, T. , 2014, “ Microfluidic Perfusion Culture Chip Providing Different Strengths of Shear Stress for Analysis of Vascular Endothelial Function,” J. Biosci. Bioeng., 118(3), pp. 327–332. [PubMed]
Sigal, A. , and Danberg, J. E. , 1990, “ New Correlation of Roughness Density Effect on the Turbulent Boundary Layer,” AIAA J., 28(3), pp. 554–556.
Coleman, H. W. , Hodge, B. K. , and Taylor, R. P. , 1984, “ A Re-Evaluation of Schlichting's Surface Roughness Experiment,” ASME J. Fluids Eng., 106(1), pp. 60–65.
van Rij, J. A. , Belnap, B. J. , and Ligrani, P. M. , 2002, “ Analysis and Experiments on Three-Dimensional, Irregular Surface Roughness,” ASME J. Fluids Eng., 124(3), pp. 671–677.
Zdravkovich, M. M. , 1997, Flow Around Circular Cylinders, Vol. 1: Fundamentals, Oxford University Press, Oxford, UK.
Beaudan, B. , and Moin, P. , 1994, “ Numerical Experiments on the Flow Past a Circular Cylinder at Sub-Critical Reynolds-Number,” Report No. TF-62, Center for Turbulence Research, Stanford University, Stanford, CA.
Williamson, C. H. K. , 1996, “ Vortex Dynamics in the Cylinder Wake,” Annu. Rev. Fluid Mech., 28(1), pp. 477–539.
Rosetti, G. F. , Vaz, G. , and Fujarra, A. L. C. , 2012, “ URANS Calculations for Smooth Circular Cylinder Flow in a Wide Range of Reynolds Numbers: Solution Verification and Validation,” ASME J. Fluids Eng., 134(12), p. 121103.
Agimelen, O. S. , and Olmsted, P. D. , 2013, “ Apparent Fracture in Polymeric Fluids Under Step Shear,” Phys. Rev. Lett., 110(20), p. 204503. [PubMed]
Gao, T. , Hu, H. H. , and Castañeda, P. P. , 2012, “ Shape Dynamics and Rheology of Soft Elastic Particles in a Shear Flow,” Phys. Rev. Lett., 108(5), p. 058302. [PubMed]
Gutierrez, R. A. , and Crumpler, E. T. , 2007, “ Potential Effect of Geometry on Wall Shear Stress Distribution Across Scaffold Surfaces,” Ann. Biomed. Eng., 36(1), pp. 77–85. [PubMed]
Popa, C. V. , Zaidi, H. , Arfaoui, A. , Polidori, G. , Taiar, R. , and Fohanno, S. , 2011, “ Analysis of Wall Shear Stress Around a Competitive Swimmer Using 3D Navier–Stokes Equations in CFD,” Acta Bioeng. Biomech., 13(1), pp. 3–11. [PubMed]
Choi, S. , Karp, J. M. , and Karnik, R. , 2012, “ Cell Sorting by Deterministic Cell Rolling,” Lab Chip, 12(8), pp. 1427–1430. [PubMed]
Loutherback, K. , Chou, K. , Newman, J. , Puchalla, J. , Austin, R. , and Sturm, J. , 2010, “ Improved Performance of Deterministic Lateral Displacement Arrays With Triangular Posts,” Microfluid. Nanofluid., 9(6), pp. 1143–1149.
Huang, Z. , McLamore, E. S. , Chuang, H. S. , Zhang, W. , Wereley, S. , Leon, J. L. C. , and Banks, M. K. , 2013, “ Shear-Induced Detachment of Biofilms From Hollow Fiber Silicone Membranes,” Biotechnol. Bioeng., 110(2), pp. 525–534. [PubMed]
Kim, J. , Kim, H.-S. , Han, S. , Lee, J.-Y. , Oh, J.-E. , Chung, S. , and Park, H.-D. , 2013, “ Hydrodynamic Effects on Bacterial Biofilm Development in a Microfluidic Environment,” Lab Chip, 13(10), pp. 1846–1849. [PubMed]
Langhamer, O. , Wilhelmsson, D. , and Engström, J. , 2009, “ Artificial Reef Effect and Fouling Impacts on Offshore Wave Power Foundations and Buoys—A Pilot Study,” Estuarine Coastal Shelf Sci., 82(3), pp. 426–432.
Logan, B. , Cheng, S. , Watson, V. , and Estadt, G. , 2007, “ Graphite Fiber Brush Anodes for Increased Power Production in Air-Cathode Microbial Fuel Cells,” Environ. Sci. Technol., 41(9), pp. 3341–3346. [PubMed]
Patanakar, S. V. , 1980, Numerical Heat Transfer and Fluid Flow, McGraw-Hill, New York.
Sprittles, J. E. , and Shikhmurzaev, Y. D. , 2011, “ Viscous Flow in Domains With Corners: Numerical Artifacts, Their Origin and Removal,” Comput. Methods Appl. Mech. Eng., 200(9–12), pp. 1087–1099.
Saad, Y. , and Schultz, M. , 1986, “ GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems,” SIAM J. Sci. Stat. Comput., 7(3), pp. 856–869.
Van der Vorst, H. A. , 1992, “ Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems,” SIAM J. Sci. Stat. Comput., 13(2), pp. 631–644.
Meng, Q. , Thanh-Tung, C. , and Tiow-Seng, T. , 2013, “ Computing 2D Constrained Delaunay Triangulation Using the GPU,” IEEE Trans. Visual. Comput. Graphics, 19(5), pp. 736–748.
Löhner, R. , 2014, “ Recent Advances in Parallel Advancing Front Grid Generation,” Arch. Comput. Methods Eng., 21(2), pp. 127–140.
Sprittles, J. E. , and Shikhmurzaev, Y. D. , 2011, “ Viscous Flows in Corner Regions: Singularities and Hidden Eigensolutions,” Int. J. Numer. Methods Fluids, 65(4), pp. 372–382.
Kundu, P. K. , Cohen, I. M. , and Dowling, D. R. , 2012, Fluid Mechanics, 5th ed., Academic Press, Boston, MA, p. iii.
COMSOL, 2013, “ comsol Multiphysics Release Notes,” COMSOL Multiphysics, Burlington, MA, accessed May 2013,
Tamayol, A. , Yeom, J. , Akbari, M. , and Bahrami, M. , 2013, “ Low Reynolds Number Flows Across Ordered Arrays of Micro-Cylinders Embedded in a Rectangular Micro/Minichannel,” Int. J. Heat Mass Transfer, 58(1–2), pp. 420–426.
Chien-Yu, F. , Sheng-Yang, T. , Shih-Mo, Y. , Long, H. , Cheng-Hsien, L. , and Hwan-You, C. , 2014, “ A Microfluidic Chip With a U-Shaped Microstructure Array for Multicellular Spheroid Formation, Culturing and Analysis,” Biofabrication, 6(1), p. 015009. [PubMed]
Stone, H. A. , Stroock, A. D. , and Ajdari, A. , 2004, “ Engineering Flows in Small Devices,” Annu. Rev. Fluid Mech., 36(1), pp. 381–411.
Gogolides, E. , Ellinas, K. , and Tserepi, A. , 2015, “ Hierarchical Micro and Nano Structured, Hydrophilic, Superhydrophobic and Superoleophobic Surfaces Incorporated in Microfluidics, Microarrays and Lab on Chip Microsystems,” Microelectron. Eng., 132, pp. 135–155.
Sugioka, K. , Masuda, M. , Hongo, T. , Cheng, Y. , Shihoyama, K. , and Midorikawa, K. , 2004, “ Three-Dimensional Microfluidic Structure Embedded in Photostructurable Glass by Femtosecond Laser for Lab-on-Chip Applications,” Appl. Phys. A, 79(4–6), pp. 815–817.
Park, H. , Park, H. , and Kim, J. , 2013, “ A Numerical Study of the Effects of Superhydrophobic Surface on Skin-Friction Drag in Turbulent Channel Flow,” Phys. Fluids, 25(11), p. 110815.

## Figures

Fig. 1

In total, 36 structures were considered for systematic determination of surface shear stress and are sorted according to the geometry of the top facet: rectilinear prisms (eight shapes with nonradial, uniform top facet), radial prisms (11 shapes with radial, uniform top facet), nonvertical prisms (ten shapes with top facet of varying height), and apex structures (seven shapes with no explicit top facet)

Fig. 2

A regime map for the 36 structures separated into four categories studied for the overall surface shear stress on each structure for (a) Re = 0.1 and (b) Re = 100. Overall shear stress is nondimensionalized with respect to (τs)all of a cylinder (100 μm in diameter and height). For Re = 0.1, 100 (τs)all of a cylinder was estimated to be 7.5 mN/m2 and 16.1 N/m2, respectively. The direction of fluid flow is indicated by the gray arrow with the coordinate system shown near the bottom left of the regime map.

Fig. 3

Plot showing (τ¯s)i, average facet shear stress of (a) front, (b) rear, (c) side, and (d) top facets of all the 36 structures shown in the regime map at Re = 0.1. Each microstructure was categorized into high and low shear based on critical shear stress determined for each facet. The classification threshold for defining “high” and “low” shear stress for identifying structures was based on the relative comparison to a cylinder in cross-flow as discussed in the main text. Side facets exhibit minimal changes to shear stress as a function of microstructure morphology. The average side facet shear stress of all the 36 microstructures was estimated to be 0.79 ± 0.13 (minimal standard deviation). In general, the magnitude of shear stress is maximum at the top facet, followed by the side facet.

Fig. 4

Plot showing (τ¯s)i average facet shear stress of (a) front, (b) rear, (c) side, and (d) top facets of all the 36 structures shown in the regime map at Re = 100. The classification threshold for defining “high” and “low” shear stress for identifying structures was based on the relative comparison to a cylinder in cross-flow.

Fig. 5

Variation of velocity gradients at the top facet, nondimensionalized with the overall shear rate of a cylinder (7.89 s−1) at Re = 0.1 with equivalent critical dimensions. The gradient of velocity components in a direction normal to the top facet along ẑ contributes significantly toward (τ¯s)top.

Fig. 6

The gradient of velocity in the normal direction to the front, side, and rear facets of a cube at Re = 0.1. Since the flow is along x̂,  |u| > |v,w  |, and therefore, the side facet experiences a greater shear stress compared to the front and rear facets.

Fig. 7

The results of incrementing four physical dimensions of a cube through a range of values are presented. Distinct points represent discrete model solutions; and lines have been added for clarity. Dimensionless values of (τ¯s)i represent average surface shear stress nondimensionalized with respect to (τ¯s)front of an unaltered cube (12.5 N/m2 and Re = 100) with side length of 100 μm. In (a), the flat facets of the cube were altered to approach triangular cross sections (0 ≤  θp  ≤ 26.6 deg) eventually reaching a right pyramid. In (b), the curvature of the front and side facets was increased (0 ≤  Rv  ≤ 50 μm) to eventually be a cylinder. In (c), the angle of the front facet to the incident flow was increased (0 ≤  θw  ≤ 45.0 deg) until the final structure was a wedge. In (d), the curvature of the top facet was increased (0 ≤  Rt  ≤ 50 μm) until the form shown was reached. The direction of flow is indicated by the gray arrow.

Fig. 8

The effects on (τ¯s)i for different angles of orientation for (a) a right pyramid and (b) a cube at Re = 100. At ϕp  = 0 deg and ϕC  = 0 deg, the FOI is in the front facet position. Dimensionless values of (τ¯s)i represent average surface shear stress nondimensionalized with respect to the pyramid and cube at ϕp  = 0 deg and ϕC  = 0 deg, respectively. The direction of fluid flow is always fixed and indicated by the gray arrow; FOI is indicated by a darkened face.

## Discussions

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections