The majority of the eye models developed in the late 90s and early 00s considers only heat conduction inside the eye. This assumption is not entirely correct, since the anterior and posterior chambers are filled aqueous humor (AH) that is constantly in motion due to thermally-induced buoyancy. In this paper, a three-dimensional model of the human eye is developed to investigate the effects AH hydrodynamics have on the human eye temperature under exposure to external heat sources. If the effects of AH flow are negligible, then future models can be developed without taking them into account, thus simplifying the modeling process. Two types of external thermal loads are considered; volumetric and surface irradiation. Results showed that heat convection due to AH flow contributes to nearly 95% of the total heat flow inside the anterior chamber. Moreover, the circulation inside the anterior chamber can cause an upward shift of the location of hotspot. This can have significant consequences to our understanding of heat-induced cataractogenesis.
The heat and mass transfer on time dependent hydrodynamic squeeze flow of Jeffrey nanofluid across two plates over permeable medium in the slip condition with heat generation/absorption, thermal radiation and chemical reaction are investigated. The impacts of Brownian motion and thermophoresis is examined in the Buongiorno's nanofluid model. Conversion of the governing partial differential equations to the ordinary differential equations is conducted via similarity transformation. The dimensionless equations are solved by imposing numerical method of Keller-box. The outputs are compared with previous reported works in the journals for the validation of the present outputs and found in proper agreement. The behavior of velocity, temperature, and nanoparticles concentration profiles by varying the pertinent parameters are examined. Findings portray that the acceleration of the velocity profile and the wall shear stress is due to the squeezing of plates. Furthermore, the velocity, temperature and concentration profile decline with boost in Hartmann number and ratio of relaxation to retardation times. It is discovered that the rate of heat transfer and temperature profile increase when viscous dissipation, thermophoresis and heat source/sink rises. In contrast, the increment of thermal radiation reduces the temperature and enhances the heat transfer rate. Besides, the mass transfer rate decelerates for increasing Brownian motion in nanofluid, while it elevates when chemical reaction and thermophoresis increases.
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.