Transition delay using biomimetic fish scale arrays

Aquatic animals have developed effective strategies to reduce their body drag over a long period of time. In this work, the influence of the scales of fish on the laminar-to-turbulent transition in the boundary layer is investigated. Arrays of biomimetic fish scales in typical overlapping arrangements are placed on a flat plate in a low-turbulence laminar water channel. Transition to turbulence is triggered by controlled excitation of a Tollmien–Schlichting (TS) wave. It was found that the TS wave can be attenuated with scales on the plate which generate streamwise streaks. As a consequence, the transition location was substantially delayed in the downstream direction by 55% with respect to the uncontrolled reference case. This corresponds to a theoretical drag reduction of about 27%. We thus hypothesize that fish scales can stabilize the laminar boundary layer and prevent it from early transition, reducing friction drag. This technique can possibly be used for bio-inspired surfaces as a laminar flow control means.

Fish are one of the oldest evolutionary species contributing to more than half of the living vertebrates distributed almost evenly across seawater and freshwater regions of the world 1,2 . It comprises more than 33,100 species and is larger than the sum of all other vertebrates, with the size range varying from a few millimetres to more than 10 m. Of the 33,100 species, more than 26,000 belong to bony fish. Only 1,000 species belong to cartilaginous fish, such as sharks, and just only 100 species belong to jawless fish (e.g. lamprey and hagfish) 3 . Fish are highly dynamic creatures that persistently travel in water to reproduce, feed on their prey, and evade from their predator 4 . As a result, much of the energy expended is largely for locomotion (against drag from skin friction and pressure drag) within the aquatic environment. This can be still water or turbulent water, which is often the case in river flows 5,6 .
Fish have prodigious features for flow control over their bodies, adapted to the environment and living circumstances. The skin and appendages are essentially the main parts where the flow is most likely to be tweaked to meet the need. For example, sharks swim at speeds ranging from 0.3 to 0.9 body-length/s, reaching Reynolds numbers 7 more than a million, which makes the boundary layer turbulent over most of their body 8 . The placoid scales (similar to a riblet shape) have been proven to reduce the turbulent skin friction drag 9 . Most studies over a decade focused on this specific scale shape 10,11 . Similarly, pectoral flippers on humpback whales have leadingedge tubercles that prevent stalling and give them high manoeuvrability 12 . Dolphins (aquatic mammals) reduce their drag by delaying the transition to turbulent flow on their body due to their anisotropic and compliant skin structure which dampens the flow instabilities 13 .
This study will focus on the hydrodynamics of the skin of bony fish. More than 95% of the existing bony fish belong to teleosts whose skin is covered by leptoid scales that are further classified into scales of ctenoid and cycloid type 3 . Scale classification and morphology from various fishes were determined and inferred for possible hydrodynamic functions 14,15 . Most teleost fish operate at the Reynolds number range where transitional boundary layer flow prevails on the fish surface 5,8 . In addition, their elongated body with an elliptical cross-section resembles that of an hydrofoil. Therefore, flow over these bodies can be closely related to the flow over a flat plate. Very few research on the flow dynamics over typical skins of teleost have been performed so far. For grass carp, (Ctenopharyngodon idellus) some geometric parameters of the scales were scanned and a bionic surface was created with individual, non-overlapping elements resembling those scales. When tested on a flat plate in a towing tank, the results showed a drag reduction of approximately 3% 16 . Recently, some of the present authors have investigated the scale structure of European sea bass (Dicentrarchus labrax) and designed a biomimetic surface, which mimics the realistic features of overlapping scales and their characteristic surface pattern. Computation Fluid Dynamics (CFD) was used to study the flow pattern over the surface and revealed a hitherto unknown effect of the scales as a mechanism to generate a regular pattern of parallel streamwise velocity streaks www.nature.com/scientificreports/ in the boundary layer 17 . To prove their existence also on the real fish skin, oil flow visualisation was done on sea bass and common carp, which indeed confirmed their presence in a regular manner along their real body, with the same arrangement relative to the scale array as observed along the biomimetic surface . These results let the  authors hypothesize about a possible mechanism for transition delay, inspired by various previous fundamental  transition studies, where streaky structures generated by cylindrical roughness elements or vortex generator  arrays have shown a delay of transition [18][19][20] .
The main purpose of this study is to prove this hypothesis with experiments in a well-established low-turbulence water channel facility, where the transition process can be studied under well-defined boundary conditions with controlled excitation of the fundamental instability of the Tollmien-Schlichting waves. In addition, the facility allows optical access to the boundary layer flow and detailed measurements of the instability using hot-wire. The facility is established as a unique low-turbulence laminar water channel and there has been a long history of transition experiments [21][22][23] . The biomimetic scale array used in this study is based on the scale structure of European sea bass (Dicentrarchus labrax), and its details have been published in the authors' preceding paper 17 . Although the geometry was derived from sea bass, the scale pattern is common to most species of seawater and freshwater fishes. Figure 1F displays the scale pattern used in this study where the area with pink colour is the central region and the area with a light blue colour represents the region of overlap between adjacent scales. In a direct comparison with real fish scales, the biomimetic scale pattern in Fig. 1F looks very similar to the natural scales of the species in Fig. 1A-E.

Experimental conditions
The experiments were conducted in the Laminar Water Channel of the Institute of Aerodynamics and Gas Dynamics in Stuttgart, Germany. The test section dimensions are 10 m in length (streamwise direction: X-axis), 1.2 m in width (spanwise direction: Z-axis) and 0.2 m in height (wall-normal direction: Y-axis) with a turbulence intensity value less than 0.05% 20 . Hot-film anemometry was used to measure the velocity of the flow inside the water channel. A flat plate was used to generate a laminar boundary-layer flow over it at a flow velocity ( U ∞ ) of 0.086 m/s corresponding to a Reynolds number of Re L ∼ 5.2 · 10 5 based on the flat plate length ( x = L ) of 6 m, where the Reynolds number is defined by ( Re x = xU ∞ /ν ) and ν is the kinematic viscosity of water. The experimental details are shown in Fig. 2A. The leading edge of fish scale arrays was located at ( X ke ) 0.3 m from the leading edge and two models were used to develop different streak amplitudes. The first model, Setting-1 was 3D printed with three rows of scales with 5 mm height followed by three rows of scales with 3 mm height. The second model, Setting-2 was 3D printed with eight rows of scales with 3 mm height. The trailing edge of the model has a smooth ramp to avoid separation behind the model. The roughness height of the scales ( h s ) on sea bass was around 0.15 mm for a fish length of 350 mm. At a swimming speed of around 1.5 m/s the boundary layer thickness ( δ ) will be about 0.5 mm at the location where the scale array begins, which leads to a δ/h s ratio of 3. Hence, in this study the δ/h s ratio is maintained to match the dynamic similarity after the upscaling of the biomimetic scale array 24 . This is a common practice in fluid mechanics to upscale the model based on geometric and kinematic similarities, for example, see the tests on upscaled shark skin scales here 25 . The CAD drawings of both models are shown in Fig. 1G,H in detail. The roughness Reynolds number defined by R kk = ku k /ν is 334 for Setting-1 and 143 for Setting-2, where u k is the undisturbed velocity without roughness at the height (k) of the maximum roughness 26 . Both models were 3D printed in spanwise segments of 0.2 m so that the length of the total setting in Z-axis was 0.6 m. A standard procedure in transition research along a flat plate is used herein to investigate the response of the boundary layer to the modified surface. The method uses a vibrating wire, which is spanned in the spanwise direction and used to excite a 2D Tollmien-Schlitching wave 27 . Location of the wire within the boundary layer and the vibration frequency are adapted to the theoretical instability diagram of the Blasius solution for a laminar 2D boundary layer 27 . A wire of 0.1 mm diameter is located at X TS = 1.2 m from the leading edge of the flat plate 20 at a wall-normal distance of 5 mm from its surface and vibrating at a physical frequency of f = 0.2 Hz (the corresponding normalised frequency is F = 2πf ν/U 2 ∞ = 166 · 10 −6 ) with an amplitude of 0.25 mm.

Results
The surface, mimicking the array of overlapping scales along the flat plate in the low-turbulence facility, again show the generation of streamwise velocity streaks, similar as detected in our previous study on real fish bodies 17 . Contours of constant streamwise velocity in a cross-section of the boundary layer in Fig. 3A illustrate the velocity variation at 1.2 m from CFD simulation and experiment for one wavelength of the streak ( z = 50 mm ) for Setting-1. Note, that the result is periodic in spanwise direction with each row of scales. The Blasius contour at the right depicts the 2-D mean flow over the same wavelength 28 . A detailed comparison of velocity profiles in the high-and low-velocity regions behind the scales are shown in Fig. 3B,C, respectively. In the region of high velocity and low velocity regions the deviation from the reference Blasius profile depicts velocity deficit or increase behind the scale array. Spanwise averaged streamwise velocity in the streaky flow is compared with Blasius profile in Fig. 3D. The streaky flow produces a fuller velocity profile when compared with the reference flow which results in a shape factor of 2.47 instead of 2.59 which is comparable with the Large Eddy Simulation results reported in literature for the streaky base flow 29 . CFD and experimental results are comparable with only minor variations that may be attributed to the boundary conditions and the experimental uncertainties. Measurements at different locations further downstream prove that the streak persists in streamwise direction (not shown here).
This modulation of the velocity is fundamentally different from the streaky structure generated by the lift-up effect caused by a vortex generator or cylinder array 30  Afterwards it drops as a result of the decelerating trailing ramp. Once again the flow reorganises up to some extent to increase the streak amplitude and then the viscosity causes it to decay continuously downstream. Both models did not induce a bypass transition (instantly tripping laminar flow into turbulent), nor did they induce secondary streak instability as seen from hot-film signals and with flow visualisations.
To investigate the response of the boundary layer to the scaled surface with regard to the laminar-to-turbulent transition process, a controlled transition experiment with a representative Tollmien-Schlichting (TS) wave at a given frequency were performed, following the method invented in 27 . In Fig. 4A the neutral stability curve is shown as a black line for the present free-stream velocity of U ∞ = 0.086 m/s . The neutral stability curve is given along the X-axis and also for comparison with non-dimensional parameters in similar studies based on the Reynolds number ( Re x ). The area within the stability curve is the region in which, according to linear stability theory, infinitesimal disturbances will grow exponentially 7 . The velocity signals were measured for 60 s at a data acquisition frequency of 100 Hz at Y = 10 mm from the wall from X = 1.98 m to 5.92 m for the reference flat plate and fish scale array (Setting-1) cases. In the reference case, the induced small disturbances from the vibrating wire grow in the streamwise direction inside the instability region which can be inferred from the black u rms curve from 2 to 2.6 m in Fig. 4B. Initially, the so-called primary instability mechanism increases the velocity fluctuations until secondary instability mechanism set in, afterwards the fluctuations increase rapidly until they reach a peak around 3.8 m in Fig. 4B. From there on the flow is turbulent as observed from the constant u rms plateau after 4.5 m 31 . However, for the flow with fish scale array (Setting-1), as seen from the red line in Fig. 4B the fluctuation level u rms remains almost constant until 4 m and it increases monotonically but with a lower rate when compared with the reference case. The local flow state can be defined generally by the intermittency  www.nature.com/scientificreports/ Next, we explore the flow by means of temporal velocity signals for the two cases previously discussed. Figure 5A displays the velocity signals subtracted from their mean values at 2, 2.5, 3.0 and 3.5 m for a period of 20 s for the reference case without fish scales. The amplitude of the velocity signals in Fig. 5A increases with streamwise coordinate X. The respective frequency spectra of these velocity signals is given in Fig. 5C, where the abscissa is normalised with respect to the vibration frequency ( F 0 = 0.2 Hz ) of the vibrating ribbon. At X = 2 m the spectrum displays a peak at F/F 0 = 1 which indicates that the fluctuation energy is only from the wire's fundamental vibrating frequency. Higher harmonics and subharmonics of the fundamental frequency F 0 appear further downstream (see the peaks at 3 and 3.5 m) and the disturbance energy also increases compared to the spectrum at 2 m. At 3.5 m the energy is increased over all frequencies given in the plot indicating that the flow is becoming increasingly disturbed resulting in turbulence. On the contrary, the velocity magnitude for the flow with fish scale array remains within 2% at all locations, as shown in Fig. 5B. At the same time the fluctuation energy is very small compared to the flat plate case as shown in Fig. 5D. Additionally, the higher harmonic components in the flow are completely absent in the case of the fish scale array. This reflects the very low level of velocity fluctuations u rms depicted by a red line in in Fig. 4B. The increase of u rms beyond 4 m is due to uncontrolled background oscillations of the water tunnel and not necessarily due to a re-amplification of the TS wave.
Complementing information about the flow states in different cases is obtained from flow visualisation using the method of surface streakline generation. Individual potassium permanganate crystals were placed on the flat plate and, while dissolving as dye in the water, they visualize the flow close to the surface. These streak lines will be visible as compact dye lines if the flow is laminar, while, in contrast, they develop kinks and diffuse very quickly when the flow is turbulent. The locations where the visualisations has been done is shown in Fig. 6A  www.nature.com/scientificreports/ for brevity. Figure 6B(Reference) shows the visualisation picture for the reference flat plate case from X = 3.15 to 3.85 m. Certainly, the streaklines are visible for more than 70% of the picture indicating laminar flow with some instabilities. In Figure 6B Drag estimation. Flow over any body will experience drag that has two components, skin friction and pressure drag. For a flat plate, the drag is only from skin friction with the friction coefficient for laminar and turbulent flow given by Eqs. (2) and (3), respectively 28 . These equations were compared with Direct Numerical Simulation results and found to be comparable with similar kind of TS waves 33 .
Hence, the total drag along a flat plate with laminar and turbulent flow regimes can be approximated by the summation of the drag components by the Eq. (4), where, x L is the location from the leading edge of the flat plate where the flow is assumed to change from laminar to turbulent state.
For the sake of comparing the drag for both cases, the location x L is assumed where the intermittency factor reaches a value of 50%. For the reference flat plate case the location x L is estimated from the hot-film measurements to be at 3.3 m and for the case with the fish scale array (Setting-1) the location x L is placed at 5.3 m. For Setting-2 the location x L is chosen from the flow visualisation of about 4.3 m since we do not observe any turbulence  Fig. 2C). The components of drag for the two cases are given in Table 1. The net drag (D Net ) is reduced by 0.02 N with the fish scale array (Setting-1) and 0.008 N for Setting-2. This results in a net drag reduction of about 27% for Setting-1 and 10.7% for Setting-2 when compared with the reference flat plate case. This result can be understood by using the skin friction plots as depicted in Fig. 7A. The laminar skin friction curve is shown as a dashed pink line and the turbulent skin friction curve as a dashed green line. Typical transition curves appear in all the cases considered here 28 . Generally, if a flow becomes turbulent the skin friction coefficient rises to almost twice its value for laminar flow at a particular location. The total drag of the surface is the area under the skin friction curve, therefore, the area under the curve reduces for fish scale array when compared with the reference flat plate case. Furthermore, the reduction in integral amplitude is proportional to the streak amplitude in the experiments considered here.
We have demonstrated that the fish scale array could delay transition to reduce the net drag. The underlying mechanism is the attenuation of the modulated TS waves due to the streamwise velocity streaks in the base flow. The latter produce a spanwise averaged flow with a steeper velocity gradient than the Blasius solution (reference flat plate case). This leads to a smaller shape factor which is known to stabilize the boundary layer 29,34 . This is Table 1. Comparison of estimated drag for reference flat plate and fish scale array (Setting-1).

Configuration
Laminar www.nature.com/scientificreports/ also seen by the streamwise decay of the observed lambda-vortices for the scales. In comparison, in the classical transition scenario (regular Blasius flow) the two dimensional TS waves grow inside the boundary layer within the linear instability region and grow further until the amplitude of fluctuation increases above a critical amplitude, which is when three-dimensional undulations lead to the formation of strong -vortices 35 ( non-linear flow regime based on H-type or K-type transition 29 ) and ultimately to turbulence. The spanwise wavelength of these -vortices (spacing between the legs of the -vortex) is generally larger than half of the TS wavelength 36 .
Herein, for the streaky base flow, the TS waves reorganise already early in the linear phase into weak -vortices as depicted in Fig. 7B due to the streamwise modulation of the flow. The spanwise wavelength of these weak -vortices is equal to the wavelength of the fish scale array, different from the wavelength on natural transition. In addition, because of the stabilizing effect of the smaller shape-factor of the boundary layer 34 the weak -vortices decay in the downstream direction. This proposed mechanism follows similar arguments given in 29 for the simulation of transition delay due to finite amplitude streaks.

conclusions
In this paper, we have investigated and shown the stabilizing influence of a periodic streaky base flow generated by a biomimetic scale array along a flat plate boundary layer on a Tollmien-Schlitching wave. In the following we shall discuss these results with respect to the existing literature on transition delay after listing out the salient outcomes from this study. The following observations were made: • The biomimetic fish scale array produces stable velocity streaks in the otherwise laminar boundary layer flow. The spanwise wavelength of the streaks can be controlled with the spacing between adjacent rows of the scale array. • The streak amplitude can be increased by increasing the number of rows in streamwise direction or by increasing the thickness of the scale array, as long as the roughness height is not exceeding a critical level to trigger bypass transition. • The velocity profile develops towards a smaller shape factor (on average over the span) which is known to stabilize the boundary layer 34 . As a consequence, early 3D instabilities of the TS waves in form of weak -vortices decay again while convecting downstream, which results in a delay of transition, similar as observed in 29 • The observed transition delay has direct consequences on the skin friction drag, as the laminar part of the boundary layer is extended relative to the reference flat plate case. This reduces the energy consumption to overcome the wall friction loss.
As explained in the introduction, most fish species operate in the transitional laminar-to-turbulent Reynolds number regime which was tested in this study (i.e., Re L ∼ 10 5 ). Experiments revealed that the scale array attenuates the TS wave and hence is able to delay laminar-turbulent transition which results in a maximum net drag reduction of about 27% for the given configuration. The present mechanism to generate the streaks differs from those in previous studies where cylindrical roughness elements and vortex generators were used to delay transition. Cylinder arrays can produce maximum streak amplitudes of about 12%, beyond that, the cylinders will trigger bypass transition because of the absolute instability in the wake 31 . In addition, the roughness elements generate parasitic drag due to the pressure drop around the protruding body. Vortex generators can delay transition with higher streak amplitudes without secondary instabilities, however, they will act as bluff bodies when the flow is not perfectly aligned with the orientation of the vortex generator. Additionally, in the case of purely laminar flow where transition does not occur (i.e. Re L < 10 4 ), these two types of bodies will always generate parasitic drag.
To be effective at all operating Reynolds numbers, it requires a multi-role flow control mechanism. Our previous studies showed that the biomimetic fish scale array already achieved laminar drag reduction in a clean laminar flow (where the ratio of boundary layer thickness to the maximum scale height was greater than 10), thus reducing the drag in the low-velocity regime, too 17 . The mechanism by which the overlapping scales generate the streaks is via producing a spanwise flow near the wall, which is sustained by the repeated overlapping along the rows of the scales. When the swimming speed increases and the flow is likely to be transitional, the streaks from the scales tend to prolong the laminar flow by delaying the transition without adding any parasitic drag, thus minimizing skin friction. Note, that the results further indicate no tendency of the scales to generate a bypass transition. Therefore this mechanism is assumed to be robust against moderate variations in details, like scale shape and height. The results let us speculate that the overlapping scale arrays on most bony fishes are an evolutionary result to minimize friction drag by producing streaky flow which produces a fuller laminar velocity profile on the surface (smaller shape-factor).
Despite the promising results from this study, salient limitations should be mentioned. Primarily, the tests were done on a flat plate ignoring the pressure gradient which is inevitable on the body of the fish. However, recent experiments on an aerofoil with streaky base flow generated by miniature-vortex generators delays transition for a particular configuration which motivates the use of these fish scale array on a surface with imposed pressure gradient 37 . Secondly, the flexibility of the scales and also the undulatory motion of the fish is not considered here, which largely changes the transitional boundary layer 38 due to unsteady effects which were not part of the present study. Therefore, the primary focus was on the performance of the biomimetic fish scale array on a very controlled transition scenario for comparison with previous studies which were successful in delaying flat-plate boundary-layer transition. From the observation that streaks exist, already proven with the help of surface flow visualisations on a real fish, and the observation that scale arrays delay transition without by-pass transition, it is hypothesized that the fish scales are efficient in delaying laminar-turbulent transition on a real www.nature.com/scientificreports/ fish body, as well. Additionally, the performance of the scale array (rigid or flexible) in turbulent or separated boundary layers remains to be studied in future research.

Methods
Biomimetic fish scale models. Biomimetic fish scale models were 3D printed at City University of London and University of Stuttgart using ABS plastic material of density about 1080 kg/m 3 . The array was made in many pieces because of the size restriction in the printer. The models were modeled in CATIA and were 200 mm wide and printed as separate tiles with extra 1 mm thickness at the base for stable print. The leading part with smooth ramp and three scale rows was printed as one piece and the other three scale rows with trailing edge ramp were printed as a separate piece. When placed one behind the other it will be acting as a single array of scale rows. The tile with scale array was glued on the flat plate of the Laminar Water Channel to keep it in place without any movements. Two models have been used in this study as shown in Fig. 1G,H.
Experiments using water channel. The experiments were conducted in the open channel closed loop Laminar Water Channel (Laminarwasserkanal) at University of Stuttgart, Germany. The flow is induced by two axial propellers connected to a frequency control drive to vary the RPM. The flow from the pump is passed through a honeycomb chamber followed by a long diffuser with a series of screens to a settling chamber with a contraction ratio of 7.7:1. Before the contraction three additional sets of screens were used to reduce the turbulence. The turbulence intensity lies below 0.05% within a frequency range of 0.1-10 Hz at 0.145 m/s 21 . The dimensions of the test section are 10 m in length, 1.2 m in width and 0.2 m in height. Inside the test section a very long but segmented glass plate is used to create a two dimensional Blasius boundary layer. The leading edge of the first plate is elliptical to reduce the leading edge separation and getting zero pressure gradient on the flat plate quickly. In the spanwise direction the plates are little less than the widths of the water channel which provides natural suction to prevent the corner flows. Constant temperature hot-film anemometry was used to measure the velocity of the water using DANTEC 55R15 with a 16-bit A/D converter. The overheat ratio for the hot-film probe was set at 8% as given in the manufacturer manual. All the measurements in this study have been acquired at 100 Hz for 60 s. Before starting the measurement the hot-film was calibrated in still water by traversing the probe with a controlled series of constant velocities to acquire the corresponding voltage from the data acquisition system. Finally a correlation graph was used to find the coefficients in King's law to relate the voltage (E) with velocity (U) given by The data was post-processed in MATLAB for filtering the signals and also to find the spectrum of the signals. For details on the experimental facility and measurement equipment, the reader is referred to 39 .

Flow visualisation.
To perform flow visualisations the water channel was emptied below the flat plate to dry the surface. A yellow colour sheet of width 600 mm and length of about 2500 mm was placed on the flat plate for better contrast. Then it was left to glue on the surface for a day. The yellow colour sheet was placed from 3 to 5.5 m on the flat plate where the flow disturbances grow from laminar to turbulent flow. Potassium permanganate crystals (less than 2 mm in average) were placed at the start of the yellow sheet. As the crystals dissolve with water they will colour the water without changing its physical properties. The coloured water will be clearly identified in the regions where the flow is laminar and quickly diffuse and disappear in the turbulent flow regions. The dye flow visualisation was recorded from above the tunnel using digital cameras from three regions. All the video recordings were done with a shutter time of 1/30 s. The first frame was from 3.15 to 3.85 m, the second frame from 3.55 to 4.25 m and the third frame from 4.35 to 5.05 m. The setup of the flow visualisation is shown in Fig. 2B.
CFD methodology. The computational study was done using ANSYS-Fluent 19.0. The CFD domain with fish scale array was modeled in CATIA with a spanwise length equal to two wavelengths of the array 2 Z = 100 mm as shown in Fig. 2C. The leading edge of the fish scale array was placed at X = 300 mm from the leading edge of the flat plate as in the experiments. The length of the domain in the streamwise direction was set to 1,200 mm and in the wall normal direction the domain was 200 mm in length. The inlet was specified with a velocity of U ∞ = 0.086 m/s and the outlet was specified as a pressure outlet. Periodic boundary conditions were used in the spanwise direction and the top domain was specified as free shear boundary with zero normal velocity. The domain was dicretised with 2 mm elements in the streamwise direction, 1.43 mm in the spanwise direction. The first cell height in the wall normal direction was set at 0.035 mm with inflation of 20 cells within 2 mm to capture the near wall gradient and a total of 110 grid points were used to mesh the wall normal direction. The major part of the domain was discretised with a Cartesian, structured mesh except the volume with fish scale array which was meshed with both prism and tetrahedral elements.The total number of elements for Setting-1 is 7.8 million and Setting-2 is 8.5 million. Second order pressure and second order upwinding schemes were used for discretisation with a steady state solver which was used to compute the laminar flow through the domain.