We propose a strategy for optimizing distribution of flow in a typical benchtop microfluidic chamber for dielectrophoretic application. It is aimed at encouraging uniform flow velocity along the whole analysis chamber in order to ensure DEP force is evenly applied to biological particle. Via the study, we have come up with a constructive strategy in improving the design of microfluidic channel which will greatly facilitate the use of DEP system in laboratory and primarily focus on the relationship between architecture and cell distribution, by resorting to the tubular structure of blood vessels. The design was validated by hydrodynamic flow simulation using COMSOL Multiphysics v4.2a software. Simulations show that the presence of 2-level bifurcation has developed portioning of volumetric flow which produced uniform flow across the channel. However, further bifurcation will reduce the volumetric flow rate, thus causing undesirable deposition of cell suspension around the chamber. Finally, an improvement of microfluidic design with rounded corner is proposed to encourage a uniform cell adhesion within the channel.
Understanding fluid dynamics under extreme confinement, where device and intrinsic fluid length scales become comparable, is essential to successfully develop the coming generations of fluidic devices. Here we report measurements of advancing fluid fronts in such a regime, which we dub superconfinement. We find that the strong coupling between contact-line friction and geometric confinement gives rise to a new stability regime where the maximum speed for a stable moving front exhibits a distinctive response to changes in the bounding geometry. Unstable fronts develop into drop-emitting jets controlled by thermal fluctuations. Numerical simulations reveal that the dynamics in superconfined systems is dominated by interfacial forces. Henceforth, we present a theory that quantifies our experiments in terms of the relevant interfacial length scale, which in our system is the intrinsic contact-line slip length. Our findings show that length-scale overlap can be used as a new fluid-control mechanism in strongly confined systems.
This study offers the numerical solutions for the problem of mixed convection stagnation-point flow along a permeable
vertical flat plate in an Oldroyd-B fluid. The present investigation considers the effects of thermal radiation and heat
generation/absorption in the fluid flow. The similarity transformation simplifies the complex model and the bvp4c function
generates the numerical solutions according to the variations in the governing parameters. A higher degree of shrinking
hastens flow separations. The dual solutions are visible in the range of buoyancy opposing flow. The results from this study
may be useful for the scientist to understand the behaviour of the dilute polymer solutions in the industrial applications,
for example, the drag reduction in pipe flows.
In this study, the effects of suction and injection on the mixed convection flow of a nanofluid, over a moving permeable
vertical plate were discussed. A similarity variable was used to transform the governing equations to the ordinary
differential equations, which were then solved numerically using the bvp4c programme from MATLAB. Dual solutions
(upper and lower branches) were found within a certain range of the mixed convection parameter in assisting and
opposing flow regions. A stability analysis was implemented to confirm that the upper branch solution was stable, while
the lower branch solution was unstable.
The present work was carried out to investigate the blood flow behavior and the severity of blockage caused in the
arterial passage due to the different geometries such as elliptical, trapezium and triangular shapes of stenosis. The study
was conducted with respect to various sizes of stenosis in terms of 70%, 80% and 90% area blockage of the arterial
blood flow. The study was carried out numerically with the help of advance computational fluid dynamic software. It
was found that the shape of the stenosis plays an important role in overall pressure drop across the blockage region
of artery. The highest level of pressure drop was observed for trapezoidal shape of stenosis followed by elliptical and
then by triangular shaped stenosis. The wall shear stress across the stenosis is great for trapezoidal shape followed by
triangular and elliptical stenosis for same blockage area in the artery.
A computational fluid dynamic analysis (CFD) is presented in the study of low Reynolds number fluid flow moving past bluff bodies. The study is focusing on the understanding of the effects of the apex-angles orientation on the flow structure and related occurring force. The apex-angle both facing upstream and downstream were computationally investigated. The simulation results of the cylinder solid are compared with available experimental data to justify the results and the model used. Results obtained in the present work were Strouhal number, drag coefficient, and Fast Fourier Transform (FFT). The study had found that the value of the drag force is increasing directly proportional to the apex angle. In contrast, the value of Strouhal number inversely proportional to the increasing of the apex angle. This was due to the flow over a cylinder creating a vortex shedding in the wake region which influenced the flow separation of fluid. Through the changing on orientation of the apex angle, it was also found that the characteristic linear dimension of the geometry will also be changed, thus affecting the flow pattern.
According to numerous studies, rowing performance is influenced by several factors including rower's biomechanics, rower's physiology, the force generated and stroke style. However, there is a missing gap linking such factors with rowing performance in the available literature. This paper aims to investigate the rowing mechanism in terms of rower anthropometry and physiology, which can impact its biomechanics and performance. The corresponding hydrodynamic force generated by the oar blade to accelerate the boat is also considered in the current study. To test the objectives, systematical online searching was conducted in search of the inclusion literature criteria. All included studies used Preferred Reporting item for Systematic Review and Meta-analysis (PRISMA) guidelines to obtain the final collection of articles for this review. In order to rate the quality of the articles, risk bias assessment was performed. A total of 35 studies were included in the assessment. The studies discussed the aspects of anthropometry and physiological of the rower, the biomechanics of the rower, corresponding hydrodynamic force on the oar blade and the rowing mechanism concerning boat performance. Based on the information obtained, an understanding of the important aspects of the rowing mechanism was achieved to provide an update for comprehensive improvement.
Differential equations are commonly used to model various types of real life applications. The complexity of these models may often hinder the ability to acquire an analytical solution. To overcome this drawback, numerical methods were introduced to approximate the solutions. Initially when developing a numerical algorithm, researchers focused on the key aspect which is accuracy of the method. As numerical methods becomes more and more robust, accuracy alone is not sufficient hence begins the pursuit of efficiency which warrants the need for reducing computational cost. The current research proposes a numerical algorithm for solving initial value higher order ordinary differential equations (ODEs). The proposed algorithm is derived as a three point block multistep method, developed in an Adams type formulae (3PBCS) and will be used to solve various types of ODEs and systems of ODEs. Type of ODEs that are selected varies from linear to nonlinear, artificial and real life problems. Results will illustrate the accuracy and efficiency of the proposed three point block method. Order, stability and convergence of the method are also presented in the study.
The drying of Piper betle Linn (betel) leaf extract using a lab scale spray dryer was simulated using Computational Fluid Dynamics (CFD). Three different turbulent models (standard k-ε, RNG k-ε and realizable k-ε) were used in the present study to determine the most suitable model for predicting the flow profile. Parametric studies were also conducted to evaluate the effect of process variables on the final moisture content. Four different initial droplet sizes (36, 79, 123 and 166 μm) were tested with four sets of combination of hot air temperature (140 and 160°C) and feed rate (4, 9.5 and 15 ml/min). It was found that standard k-ε is the most suitable turbulent model to predict the flow behaviour Moreover, the lowest final moisture content present in samples was obtained at 140°C and a feed rate of 15.0 ml/min.
In this study, the numerical simulation in a mixing vessel agitated by a six bladed Rushton turbine has
been carried out to investigate the effects of effective parameters to the mixing process. The study is intended to screen the potential parameters which affect the optimization process and to provide the detail insights into the process. Three-dimensional and steady-state flow has been performed using the fully predictive Multiple Reference Frame (MRF) technique for the impeller and tank geometry. Process optimization is always used to ensure the optimum conditions are fulfilled to attain industries’ satisfaction or needs (ie; increase profit, low cost, yields, etc). In this study, the range of recommended speed to accelerate optimization is 100, 150 and 200rpm respectively and the range of recommended clearance is 50, 75 and 100mm respectively for dual Rushton impeller. Thus, the computer fluid dynamics (CFD) was introduced in order to screen the suitable parameters efficiently and to accelerate optimization. In this study,
In this paper, the problem of steady laminar boundary layer flow of an incompressible viscous fluid over a moving thin needle is considered. The governing boundary layer equations were first transformed into non-dimensional forms. These non-dimensional equations were then transformed into similarity equations using the similarity variables, which were solved numerically using an implicit finite-difference scheme known as the Keller-box method. The solutions were obtained for a blunt-nosed needle. Numerical computations were carried out for various values of the dimensionless parameters of the problem which included the Prandtl number Pr and the parameter a representing the needle size. It was found that the heat transfer characteristics were significantly
influenced by these parameters. However, the Prandtl number had no effect on the flow characteristics due to the decoupled boundary layer equations.
Chemical interesterification of rubber seed oil has been investigated for four different designed orifice devices in a pilot scale hydrodynamic cavitation (HC) system. Upstream pressure within 1-3.5bar induced cavities to intensify the process. An optimal orifice plate geometry was considered as plate with 1mm dia hole having 21 holes at 3bar inlet pressure. The optimisation results of interesterification were revealed by response surface methodology; methyl acetate to oil molar ratio of 14:1, catalyst amount of 0.75wt.% and reaction time of 20min at 50°C. HC is compared to mechanical stirring (MS) at optimised values. The reaction rate constant and the frequency factor of HC were 3.4-fold shorter and 3.2-fold higher than MS. The interesterified product was characterised by following EN 14214 and ASTM D 6751 international standards.
This paper presents a numerical analysis of a stagnation-point flow towards a nonlinearly stretching/shrinking sheet immersed in a viscous fluid. The stretching/shrinking velocity and the external flow velocity impinges normal to the stretching/shrinking sheet are assumed to be in the form U ~ xm, where m is a constant and x is the distance from the stagnation point. The governing partial differential equations are converted into ordinary ones by a similarity transformation, before being solved numerically. The variations of the skin friction coefficient and the heat transfer rate at the surface with the governing parameters are graphed and tabulated. Different from a stretching sheet, it is found that the solutions for a shrinking sheet are non-unique for m > 1/3.
The paper reconsiders the problem of the mixed convection boundary layer flow near the lower stagnation point of a horizontal circular cylinder with a second order slip velocity model and a constant surface heat flux studied recently by RoKa et al. (2013). The ordinary (similarity) differential equations are solved numerically using the function bvp4c from Matlab for different values of the governing parameters. It is found that the similarity equations have two branches, upper and lower branch solutions, in a certain range of the mixed convection parameters. A stability analysis has been performed to show that the upper branch solutions are stable and physically realizable, while the lower branch solutions are not stable and therefore, not physically possible. This stability analysis is different by that presented by RoKa et al. (2013), who have presented a time-dependent analysis to determine the stability of the solution branches.
The magnetohydrodynamic (MHD) boundary-layer flow and heat transfer due to a shrinking sheet in a porous medium is considered for the first time. The Navier-Stokes equations and the heat equation are reduced to two nonlinear ordinary differential equations via similarity transformations. The transformed equations are solved by a semi-analytic method. The effects of the suction and porosity parameters, the Prandtl and Hartmann numbers on the skin friction, heat transfer rate, velocity and temperature profiles are discussed and presented, respectively.
Low-gradient magnetic separation (LGMS) of magnetic nanoparticles (MNPs) has been proven as one of the techniques with great potential for biomedical and environmental applications. Recently, the underlying principle of particle capture by LGMS, through a process known as magnetophoresis, under the influence of hydrodynamic effect has been widely studied and illustrated. Even though the hydrodynamic effect is very substantial for batch processes, its impact on LGMS operated at continuous flow (CF) condition remained largely unknown. Hence, in this study, the dynamical behaviour of LGMS process operated under CF was being studied. First, the LGMS experiments using poly(sodium 4-styrenesulfonate)-functionalized-MNP as modelled particle system were performed through batchwise (BW) and CF modes at different operating conditions. Here BW operation was used as a comparative study to elucidate the transport mechanism of MNP under the similar environment of CF-LGMS process, and it was found out that the convection induced by magnetophoresis (timescale effective is ∼1200 s) is only significant at far-from-magnet region. Hence, it can be deduced that forced convection is more dominant on influencing the transport behaviour of CF-LGMS (with resident time ≤240 s). Moreover, we found that the separation efficiency of CF-LGMS process can be boosted by the higher number of magnets, the higher MNP concentration and the lower flowrate of MNP solution. To better illustrate the underlying dynamical behaviour of LGMS process, a mathematical model was developed to predict its kinetic profile and separation efficiency (with average error of ∼2.6% compared to the experimental results).
In attached microalgae cultivation systems, cell detachment due to fluid hydrodynamic flow is not a subject matter that is commonly looked into. However, this phenomenon is of great relevance to optimizing the operating parameters of algae cultivation and feasible reactor design. Hence, this current work miniaturizes traditional benchtop assays into a microfluidic platform to study the cell detachment of green microalgae, Chlorella vulgaris, from porous substrates during its early cultivation stage under precisely controlled conditions. As revealed by time lapse microscopy, an increase in bulk flow velocity facilitated nutrient transport but also triggered cell detachment events. At a flow rate of 1000 μL min-1 of growth medium for 120 min, the algal cell coverage was up to 5% lower than those at 5 μL min-1 and 50 μL min-1. In static seeding, the evolution of attached cell resistance toward liquid flows was dependent on hydrodynamic zones. The center zone of the microchannel was shown to be a "comfortable zone" of the attached cells to sequester nutrients effectively at lower medium flow rates but there was a profile transition where outlet zones favored cell attachment the most at higher flow rates (1.13 times higher than the center zone for 1000 μL min-1). Besides, computational fluid dynamics (CFD) simulations illustrated that the focusing band varied between cross-sections and depths, while the streamline was the least concentrated along the side walls and bottom plane of the microfluidic devices. It was intriguing to learn that cell detachment was not primarily happening along the symmetry streamline. Insight gained from this study could be further applied in the optimization of operating conditions of attached cultivation systems whilst preserving laminar flow conditions.
Despite the advancement of cardiac imaging technologies, these have traditionally been limited to global geometrical measurements. Computational fluid dynamics (CFD) has emerged as a reliable tool that provides flow field information and other variables essential for the assessment of the cardiac function. Extensive studies have shown that vortex formation and propagation during the filling phase acts as a promising indicator for the diagnosis of the cardiac health condition. Proper setting of the boundary conditions is crucial in a CFD study as they are important determinants, that affect the simulation results. In this article, the effect of different transmitral velocity profiles (parabolic and uniform profile) on the vortex formation patterns during diastole was studied in a ventricle with dilated cardiomyopathy (DCM). The resulting vortex evolution pattern using the uniform inlet velocity profile agreed with that reported in the literature, which revealed an increase in thrombus risk in a ventricle with DCM. However the application of a parabolic velocity profile at the inlet yields a deviated vortical flow pattern and overestimates the propagation velocity of the vortex ring towards the apex of the ventricle. This study highlighted that uniform inlet velocity profile should be applied in the study of the filling dynamics in a left ventricle because it produces results closer to that observed experimentally.
Influences of river hydrodynamic behaviours on hydrochemistry (salinity, pH, dissolved oxygen saturations and dissolved phosphorus) were evaluated through high spatial and temporal resolution study of a sandbar-regulated coastal river. River hydrodynamic during sandbar-closed event was characterized by minor dependency on tidal fluctuations, very gradual increase of water level and continual low flow velocity. These hydrodynamic behaviours established a hydrochemistry equilibrium, in which water properties generally were characterized by virtual absence of horizontal gradients while vertical stratifications were significant. In addition, the river was in high trophic status as algae blooms were visible. Conversely, river hydrodynamic in sandbar-opened event was tidal-controlled and showed higher flow velocity. Horizontal gradients of water properties became significant while vertically more homogenised and with lower trophic status. In essence, this study reveals that estuarine sandbar directly regulates river hydrodynamic behaviours which in turn influences river hydrochemistry.
Taking into account the effect of constant convective thermal and mass boundary conditions, we present numerical solution of the 2-D laminar g-jitter mixed convective boundary layer flow of water-based nanofluids. The governing transport equations are converted into non-similar equations using suitable transformations, before being solved numerically by an implicit finite difference method with quasi-linearization technique. The skin friction decreases with time, buoyancy ratio, and thermophoresis parameters while it increases with frequency, mixed convection and Brownian motion parameters. Heat transfer rate decreases with time, Brownian motion, thermophoresis and diffusion-convection parameters while it increases with the Reynolds number, frequency, mixed convection, buoyancy ratio and conduction-convection parameters. Mass transfer rate decreases with time, frequency, thermophoresis, conduction-convection parameters while it increases with mixed convection, buoyancy ratio, diffusion-convection and Brownian motion parameters. To the best of our knowledge, this is the first paper on this topic and hence the results are new. We believe that the results will be useful in designing and operating thermal fluids systems for space materials processing. Special cases of the results have been compared with published results and an excellent agreement is found.