Numerical and series solutions of the peristaltic motion of an Oldroyd 8-constant fluid in an endoscope

In the present article, we have presented the peristaltic flow of an Oldroyd 8-constant fluid in an endoscope. The governing equations for the flow problem are simplified using long wavelength and low Reynold's number approximations. The solutions of the simplified problem are calculated using (i) Homotopy analysis method and (ii) Shooting method. The comparison of both the solutions shows a very good agreement between the results. The graphical results for the velocity field and stresses are presented to show the physical behaviour of all the parameters appearing in the problem.     

Model based calculation of substrate/inducer feed-rate profiles in fed-batch processes for recombinant protein production

A model-based feed-rate profile optimization problem is discussed for the fed-batch recombinant protein production. Two optimization procedures, an evolutionary programming technique and a simplified method using the dynamic programming concept, are discussed and compared. Modeling as well as experimental results are presented.     

Effect of ocular shape and vascular geometry on retinal hemodynamics: a computational model

A computational model for retinal hemodynamics accounting for ocular curvature is presented. The model combines (i) a hierarchical Darcy model for the flow through small arterioles, capillaries and small venules in the retinal tissue, where blood vessels of different size are comprised in different hierarchical levels of a porous medium; and (ii) a one-dimensional network model for the blood flow through retinal arterioles and venules of larger size. The non-planar ocular shape is included by (i) defining the hierarchical Darcy flow model on a two-dimensional curved surface embedded in the three-dimensional space; and (ii) mapping the simplified one-dimensional network model onto the curved surface. The model is solved numerically using a finite element method in which spatial domain and hierarchical levels are discretized separately. For the finite element method, we use an exterior calculus-based implementation which permits an easier treatment of non-planar domains. Numerical solutions are verified against suitably constructed analytical solutions. Numerical experiments are performed to investigate how retinal hemodynamics is influenced by the ocular shape (sphere, oblate spheroid, prolate spheroid and barrel are compared) and vascular architecture (four vascular arcs and a branching vascular tree are compared). The model predictions show that changes in ocular shape induce non-uniform alterations of blood pressure and velocity in the retina. In particular, we found that (i) the temporal region is affected the least by changes in ocular shape, and (ii) the barrel shape departs the most from the hemispherical reference geometry in terms of associated pressure and velocity distributions in the retinal microvasculature. These results support the clinical hypothesis that alterations in ocular shape, such as those occurring in myopic eyes, might be associated with pathological alterations in retinal hemodynamics.     

Measurement of blood flow from an assist ventricle by computation of pneumatic driving parameters

The measurement of blood flow from an assist ventricle is important but sometimes difficult in artificial heart experiments. Along with the development of a pneumatic cylinder-piston driver coupled with a ventricular assist device, a simplified method for measuring pump flow was established. From driving parameters such as the piston (or cylinder) displacement and air pressure, the pump flow could be calculated by the use of the equation of state for an ideal gas. The results of this method are broadly in agreement with electromagnetic and Doppler measurements.     

On the pressure and flow-rate distributions in tree-like and arterial-venous networks

A solution algorithm yielding the pressure and flow-rate distributions for steady flow in an arbitrary, tree-like network is provided. Given the tree topology, the conductance of each segment and the pressure distribution at the boundary nodes, the solution is obtained from a simple recursion based on perfect Gauss elimination. An iterative solution method using this algorithm is suggested to solve for the pressure and flow-rate distributions in an arbitrary diverging-converging (arterial-venous) network consisting of two tree-like networks which are connected to each other at the capillary nodes. A number of special solutions for tree-like networks are obtained for which the general algorithm is either simplified or can be replaced by closed form solutions of the pressure and flow-rate distributions. These special solutions can also be obtained for each tree of diverging-converging networks having particular topologies and conductance distributions. Sample numerical results are provided.     

A mathematical model for the peristaltic flow of chyme movement in small intestine

A mathematical model based on viscoelastic fluid (fractional Oldroyd-B model) flow is considered for the peristaltic flow of chyme in small intestine, which is assumed to be in the form of an inclined cylindrical tube. The peristaltic flow of chyme is modeled more realistically by assuming that the peristaltic rush wave is a sinusoidal wave, which propagates along the tube. The governing equations are simplified by making the assumptions of long wavelength and low Reynolds number. Analytical approximate solutions of problem are obtained by using homotopy analysis method and convergence of the obtained series solution is properly checked. For the realistic values of the emerging parameters such as fractional parameters, relaxation time, retardation time, Reynolds number, Froude number and inclination of tube, the numerical results for the pressure difference and the frictional force across one wavelength are computed and discussed the roles played by these parameters during the peristaltic flow. On the basis of this study, it is found that the first fractional parameter, relaxation time and Froude number resist the movement of chyme, while, the second fractional parameter, retardation time, Reynolds number and inclination of tube favour the movement of chyme through the small intestine during pumping. It is further revealed that size of trapped bolus reduces with increasing the amplitude ratio whereas it is unaltered with other parameters.     

Three-Dimensional Flow of an Oldroyd-B Nanofluid towards Stretching Surface with Heat Generation/Absorption

This article addresses the steady three-dimensional flow of an Oldroyd-B nanofluid over a bidirectional stretching surface with heat generation/absorption effects. Suitable similarity transformations are employed to reduce the governing partial differential equations into coupled nonlinear ordinary differential equations. These nonlinear ordinary differential equations are then solved analytically by using the homotpy analysis method (HAM). Graphically results are presented and discussed for various parameters, namely, Deborah numbers and , heat generation/absorption parameter Prandtl parameter , Brownian motion parameters, thermophoresis parameter and Lewis number . We have seen that the increasing values of the Brownian motion parameter and thermophoresis parameter leads to an increase in the temperature field and thermal boundary layer thickness while the opposite behavior is observed for concentration field and concentration boundary layer thickness. To see the validity of the present work, the numerical results are compared with the analytical solutions obtained by Homotopy analysis method and noted an excellent agreement for the limiting cases.     

Numerical Modelling of Electroosmotic Driven Flow in Nanoporous Media by Lattice Boltzmann Method

The lattice Boltzmann method was employed to simulate electroosmotic driven flow and Debye layer screening in conducting electrolyte around a porous structure with average size of 40 nm. The charge screening around the nanopores was investigated by solving the vector-superpositioned potential equilibrium distribution function and adding electro-kinetic force term to the evolution equation. In this intermediate case of moderate Debye length, the electrophoresis problem becomes complicated. The motion of the particles distorts the screening cloud, which becomes asymmetric, resulting in very complex interactions between the electrolyte, the screening cloud and the particle; but the Electroosmotic Flow (EOF) behaviour was still considered based on the Helmoholtz-Smoluchowski model with adaptation to fit nanoporous flow in the porous structure. In the present approach, the flow in the nanopores is directly modelled; the detailed flow information can be obtained by simplifying the repeated macrostructure. Due to the symmetry of the domain, the size of computational domain can be largely reduced by less repeated spherical nanoparticles. Each pore of the medium contains several lattice nodes on the simplified curved edges and potential gradients are produced by adjusting the zeta potential value. The velocity results for pressure-driven and EOF flows agree well with the analytical solutions and recent experimental results. In particular, the interface between solid particles and fluids, the influences of porosity, solid particle diameter, yield stress and electric parameters in EOF were investigated. The anti-adhesion effect of electroosmotic pumping effect was evaluated, and the pulsed DC was applied in order to enhance the performance of the electroosmotic pumping. The results demonstrate that the present lattice Boltzmann model is capable of modelling flow through nanoporous media at certain restrictions while some results deviate from the predictions based on the macroscopic theories.     

Numerical analysis of flow in an elastic artery model.

Oscillatory and pulsatile flows of Newtonian fluids in straight elastic tubes are simulated numerically with the aid of Ling and Atabek's "local flow" assumption for the nonlinear convective acceleration terms. For the first time, a theoretical assessment of the local flow assumption is presented, and the range of validity of the assumption is estimated by comparison with perturbation solutions of the complete flow problem. Subsequent simulations with the local flow model indicate that the flow field and associated wall shear stress are extremely sensitive to the phase angle between oscillatory pressure and flow waves (impedance phase angle). This phase angle, which is a measure of the wave reflection present in the system, is known to be altered by arterial disease (e.g., hypertension) and vasoactive drugs. Thus, the paper elucidates a mechanism by which subtle changes in systemic hemodynamics (i.e., phase angles) can markedly influence local wall shear stress values.     

Chess games: a model for RNA based computation

Here we develop the theory of RNA computing and a method for solving the 'knight problem' as an instance of a satisfiability (SAT) problem. Using only biological molecules and enzymes as tools, we developed an algorithm for solving the knight problem (3 x 3 chess board) using a 10-bit combinatorial pool and sequential RNase H digestions. The results of preliminary experiments presented here reveal that the protocol recovers far more correct solutions than expected at random, but the persistence of errors still presents the greatest challenge.     

Towards accurate numerical method for monodomain models using a realistic heart geometry

The simulation of cardiac electrophysiological waves are known to require extremely fine meshes, limiting the applicability of current numerical models to simplified geometries and ionic models. In this work, an accurate numerical method based on a time-dependent anisotropic remeshing strategy is presented for simulating three-dimensional cardiac electrophysiological waves. The proposed numerical method greatly reduces the number of elements and enhances the accuracy of the prediction of the electrical wave fronts. Illustrations of the performance and the accuracy of the proposed method are presented using a realistic heart geometry. Qualitative and quantitative results show that the proposed methodology is far superior to the uniform mesh methods commonly used in cardiac electrophysiology.     

Radiative Hydromagnetic Flow of Jeffrey Nanofluid by an Exponentially Stretching Sheet

Two-dimensional hydromagnetic flow of an incompressible Jeffrey nanofluid over an exponentially stretching surface is examined in the present article. Heat and mass transfer analysis is performed in the presence of thermal radiation, viscous dissipation, and Brownian motion and thermophoresis effects. Mathematical modelling of considered flow problem is developed under boundary layer and Rosseland’s approximations. The governing nonlinear partial differential equations are converted into ordinary differential equations via transformations. Solution expressions of velocity, temperature and concentration are presented in the series forms. Impacts of physical parameters on the dimensionless temperature and concentration are shown and discussed. Skin-friction coefficients are analyzed numerically. A comparison in a limiting sense is provided to validate the present series solutions.     

Simplified partial digest problem: enumerative and dynamic programming algorithms

We study the Simplified Partial Digest Problem (SPDP), which is a mathematical model for a new simplified partial digest method of genome mapping. This method is easy for laboratory implementation and robust with respect to the experimental errors. SPDP is NP-hard in the strong sense. We present an $O(n2;n)$ time enumerative algorithm and an O(n(2q)) time dynamic programming algorithm for the error-free SPDP, where $n$ is the number of restriction sites and n is the number of distinct intersite distances. We also give examples of the problem, in which there are 2(n+2)/(3)-1 non-congruent solutions. These examples partially answer a question recently posed in the literature about the number of solutions of SPDP. We adapt our enumerative algorithm for handling SPDP with imprecise input data. Finally, we describe and discuss the results of the computer experiments with our algorithms.     

On Comparison of Series and Numerical Solutions for Flow of Eyring-Powell Fluid with Newtonian Heating And Internal Heat Generation/Absorption

In this paper, we have investigated the combined effects of Newtonian heating and internal heat generation/absorption in the two-dimensional flow of Eyring-Powell fluid over a stretching surface. The governing non-linear analysis of partial differential equations is reduced into the ordinary differential equations using similarity transformations. The resulting problems are computed for both series and numerical solutions. Series solution is constructed using homotopy analysis method (HAM) whereas numerical solution is presented by two different techniques namely shooting method and bvp4c. A comparison of homotopy solution with numerical solution is also tabulated. Both solutions are found in an excellent agreement. Dimensionless velocity and temperature profiles are plotted and discussed for various emerging physical parameters.     

Rapid detection of mutagen induced micronucleated erythrocytes by flow cytometry

Summary The micronucleus test is a cytogenetic method for the screening of mutagens and carcinogens which exhibit clastogenic mechanisms of action. After application of clastogenic agents chromosomal fragments or even whole chromsomes can remain as conspicuous structures (micronuclei) in a small fraction of anucleated polychromatic erythrocytes which can be visually scored using a microscope following staining with May-Grünwald-Giemsa solutions. These time-consuming, painstaking, and tedious manual evaluations are often sources of unreliability and uncertainty. Here, a fluorescence technique is presented which applies DNA and protein fluorochromes to discriminate normal anucleated erythrocytes from micronucleated erythrocytes using a fluorescence microscope. This method is particularly tailored to be applied to flow cytometric instrumentation. Data obtained manually and automatically in flow show a strong linear correlation with high significance (r=0.96) as far as the percentage of micronucleated erythrocytes as an indicator for the mutagenicity of a given drug is concerned. These results have been obtained by means of the established clastogens cyclophosphamide and mitomycin C.     

On the possibility for computing the transmembrane potential in the heart with a one shot method: an inverse problem

We analyze the possibility for using body surface potential maps (BSPMs), a priori information about the voltage distribution in the heart and the bidomain equations to compute the transmembrane potential throughout the myocardium. Our approach is defined in terms of an inverse problem for elliptic partial differential equations (PDEs). More precisely, we formulate it in terms of an output least squares framework in which a goal functional is minimized subject to suitable PDE constraints. The problem is highly unstable and, even under optimal recording conditions, it does not have a unique solution. We propose a methodology for stabilizing and enforcing uniqueness for this inverse problem. Moreover, a fully implicit method for solving the involved minimization problem is presented. In other words, we show how one may solve it in terms of a system consisting of three linear elliptic PDEs, i.e. we derive a so-called one shot method (also commonly referred to as an all-at-once method). Finally, our theoretical findings are illuminated by a series of numerical experiments. These examples indicate that, in the presence of regional ischemia, it might be possible to approximately recover the transmembrane potential during the resting and plateau phases of the heart cycle. This is probably due to the fact that rather accurate a priori information is available during these time intervals. The problem of computing the transmembrane potential at an arbitrary time instance during a heart beat is still an open problem.     

Human stick balancing: an intermittent control explanation

There are two issues in balancing a stick pivoting on a finger tip (or mechanically on a moving cart): maintaining the stick angle near to vertical and maintaining the horizontal position within the bounds of reach or cart track. The (linearised) dynamics of the angle are second order (although driven by pivot acceleration), and so, as in human standing, control of the angle is not, by itself very difficult. However, once the angle is under control, the position dynamics are, in general, fourth order. This makes control quite difficult for humans (and even an engineering control system requires careful design). Recently, three of the authors have experimentally demonstrated that humans control the stick angle in a special way: the closed-loop inverted pendulum behaves as a non-inverted pendulum with a virtual pivot somewhere between the stick centre and tip and with increased gravity. Moreover, they suggest that the virtual pivot lies at the radius of gyration (about the mass centre) above the mass centre. This paper gives a continuous-time control-theoretical interpretation of the virtual-pendulum approach. In particular, by using a novel cascade control structure, it is shown that the horizontal control of the virtual pivot becomes a second-order problem which is much easier to solve than the generic fourth-order problem. Hence, the use of the virtual pivot approach allows the control problem to be perceived by the subject as two separate second-order problems rather than a single fourth-order problem, and the control problem is therefore simplified. The theoretical predictions are verified using the data previously presented by three of the authors and analysed using a standard parameter estimation method. The experimental data indicate that although all subjects adopt the virtual pivot approach, the less expert subjects exhibit larger amplitude angular motion and poorly controlled translational motion. It is known that human control systems are delayed and intermittent, and therefore, the continuous-time strategy cannot be correct. However, the model of intermittent control used in this paper is based on the virtual pivot continuous-time control scheme, handles time delays and moreover masquerades as the underlying continuous-time controller. In addition, the event-driven properties of intermittent control can explain experimentally observed variability.     

Robust stabilization and H(∞) controllers design for stochastic genetic regulatory networks with time-varying delays and structured uncertainties

This paper deals with the problem of stabilization design and H(∞) control for a class of genetic regulatory networks (GRNs) with both intrinsic perturbation and extrinsic perturbation. Some delay-dependent mean-square stabilization criteria are put forward for nominal systems and uncertain systems by using an improved free-weighting matrix approach. As a result, the corresponding stabilization controllers and H(∞) controllers of GRNs are constructed with time delays compensated and suboptimal solutions are obtained via exploiting an iterative procedure together with the linear matrix inequality (LMI) method and the cone complementarity liberalization (CCL) algorithm. Finally, three numerical examples are presented to illustrate the effectiveness of the proposed theoretical results.     

Dual Solutions for Nonlinear Flow Using Lie Group Analysis

`The aim of this analysis is to investigate the existence of the dual solutions for magnetohydrodynamic (MHD) flow of an upper-convected Maxwell (UCM) fluid over a porous shrinking wall. We have employed the Lie group analysis for the simplification of the nonlinear differential system and computed the absolute invariants explicitly. An efficient numerical technique namely the shooting method has been employed for the constructions of solutions. Dual solutions are computed for velocity profile of an upper-convected Maxwell (UCM) fluid flow. Plots reflecting the impact of dual solutions for the variations of Deborah number, Hartman number, wall mass transfer are presented and analyzed. Streamlines are also plotted for the wall mass transfer effects when suction and blowing situations are considered.     

Unsteady Boundary Layer Flow and Heat Transfer of a Casson Fluid past an Oscillating Vertical Plate with Newtonian Heating

In this paper, the heat transfer effect on the unsteady boundary layer flow of a Casson fluid past an infinite oscillating vertical plate with Newtonian heating is investigated. The governing equations are transformed to a systems of linear partial differential equations using appropriate non-dimensional variables. The resulting equations are solved analytically by using the Laplace transform method and the expressions for velocity and temperature are obtained. They satisfy all imposed initial and boundary conditions and reduce to some well-known solutions for Newtonian fluids. Numerical results for velocity, temperature, skin friction and Nusselt number are shown in various graphs and discussed for embedded flow parameters. It is found that velocity decreases as Casson parameters increases and thermal boundary layer thickness increases with increasing Newtonian heating parameter.