Flow of micropolar fluid through a tube with stenosis

mathematical model for the steady flow of non-Newtonian fluid through a stenotic region is presented. The results indicate that the general shape and size of the stenosis together with rheological properties of blood are important in understanding the flow characteristics and the presence of flow separation.

Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

Styrene Oxidation to Benzaldehyde by Polymer Anchored Wilkinson Complex

Wilkinson complex, insolubilized by anchoring to polymeric Amberlite beads, had been used for the liquid-phase catalytic oxidation of styrene to benzaldehyde and formaldehyde in toluene medium. Styrene conversion was followed by measuring the oxygen volume in contact with the reaction mixture in a specially designed closed batch apparatus. Styrene conversion depended upon catalyst loading and distribution inside the porous beads, while temperature had little effect on it. The internal diffusional effects on the conversion process have been taken into consideration by a mathematical model which allowed calculation of effectiveness factors for various catalyst loadings and corresponding catalyst distributions. The influence of external diffusion was separately determined by plotting initial rate versus catalyst loading. The proposed method can be readily extended to immobilized enzymes in porous matrices.

Styrene Oxidation to Benzaldehyde by Polymer Anchored Wilkinson Complex

Wilkinson complex, insolubilized by anchoring to polymeric Amberlite beads, had been used for the liquid-phase catalytic oxidation of styrene to benzaldehyde and formaldehyde in toluene medium. Styrene conversion was followed by measuring the oxygen volume in contact with the reaction mixture in a specially designed closed batch apparatus. Styrene conversion depended upon catalyst loading and distribution inside the porous beads, while temperature had little effect on it. The internal diffusional effects on the conversion process have been taken into onsideration by a mathematical model which allowed calculation of effectiveness factors for various catalyst loadings and corresponding catalyst distributions. The influence of external diffusion was separately determined by plotting initial rate versus catalyst loading. The proposed method can be readily extended to immobilized enzymes in porous matrices.

A study on bullock carts. Part 1. Engineering analysis of the two-wheel bullock cart design

An engineering analysis of the design of two-wheel bullock carts has been carried out with the aid of a mathematical model. Non-dimensional expressions for the pull and the neck load have been developed. In the first instance, the cart is assumed to be cruising at constant velocity on a terrain with the effective coefficient of rolling friction varying over a wide range (0.001 to 0.5) and the gradient varying between +0.2 to −0.2. Subsequently, the effect of inertia force due to an acceleration parallel to the ground is studied. In the light of this analysis, two modifications to the design of the cart have been proposed and the relative merits of the current designs and the proposed designs are discussed.

Stochastic Modeling of an Aircraft Traversing a Runway Using Time Series Analysis

Time series, from a narrow point of view, is a sequence of observations on a stochastic process made at discrete and equally spaced time intervals. Its future behavior can be predicted by identifying, fitting, and confirming a mathematical model. In this paper, time series analysis is applied to problems concerning runwayinduced vibrations of an aircraft. A simple mathematical model based on this technique is fitted to obtain the impulse response coefficients of an aircraft system considered as a whole for a particular type of operation. Using this model, the output which is the aircraft response can be obtained with lesser computation time for any runway profile as the input.

On the axial vibrations of rotating bars

One of the important developments in rotary wing aeroelasticity in the recent past has been the growing awareness and acceptance of the fact that the problem is inherently non-linear and that correct treatment of aeroelastic problems requires the development of a consistent mathematical model [l]. This has led to a number of studies devoted to the derivation of a consistent set of “second order” non-linear equations, for example, those of Hodges and Dowel1 [2], of Rosen and Friedmann [3], and of Kvaternik, White and Kaza [4], each of which differs from the others on the question of the inclusion of certain terms in the equations of motion. The final form of the equations depends first upon the ordering scheme used for characterizing the displacements and upon the consistency with which this is applied in omitting terms of lower order. The ideal way of achieving this would be to derive the equations of motion with all the terms first included regardless of their relative orders of magnitude and then to apply the ordering scheme.

Dynamics of a class of social interaction systems

A mathematical model of social interaction in the form of two coupler! first-order non-linear differential equations, forms the topic of this study. This non-conservative model io representative of such varied social interaction problems as coexisting sub-populations of two different species, arms race between two rival countries and the like. Differential transformation techniques developed elsewhere in the literature are seen to be effective tools of dynamic analysis of this non-linear non-conservative mode! of social interaction process.

Ferrous iron oxidation by foam immobilized Acidithiobacillus ferrooxidans: Experiments and modeling

Ferrous iron bio-oxidation by Acidithiobacillus ferrooxidans immobilized on polyurethane foam was investigated. Cells were immobilized on foams by placing them in a growth environment and fully bacterially activated polyurethane foams (BAPUFs) were prepared by serial subculturing in batches with partially bacterially activated foam (pBAPUFs). The dependence of foam density on cell immobilization process, the effect of pH and BAPUF loading on ferrous oxidation were studied to choose operating parameters for continuous operations. With an objective to have high cell densities both in foam and the liquid phase, pretreated foams of density 50 kg/m3 as cell support and ferrous oxidation at pH 1.5 to moderate the ferric precipitation were preferred. A novel basket-type bioreactor for continuous ferrous iron oxidation, which features a multiple effect of stirred tank in combination with recirculation, was designed and operated. The results were compared with that of a free cell and a sheet-type foam immobilized reactors. A fivefold increase in ferric iron productivity at 33.02 g/h/L of free volume in foam was achieved using basket-type bioreactor when compared to a free cell continuous system. A mathematical model for ferrous iron oxidation by Acidithiobacillus ferrooxidans cells immobilized on polyurethane foam was developed with cell growth in foam accounted by an effectiveness factor. The basic parameters of simulation were estimated using the experimental data on free cell growth as well as from cell attachment to foam under nongrowing conditions. The model predicted the phase of both oxidation of ferrous in shake flasks by pBAPUFs as well as by fully activated BAPUFs for different cell loadings in foam. Model for stirred tank basket bioreactor predicted within 5% both transient and steady state of the experiments closely for the simulated dilution rates. Bio-oxidation at high Fe2+ concentrations were simulated with experiments when substrate and product inhibition coefficients were factored into cell growth kinetics.

Modeling of heat and mass transfer for solid state fermentation process in tray bioreactor

A mathematical model is developed to simulate oxygen consumption, heat generation and cell growth in solid state fermentation (SSF). The fungal growth on the solid substrate particles results in the increase of the cell film thickness around the particles. The model incorporates this increase in the biofilm size which leads to decrease in the porosity of the substrate bed and diffusivity of oxygen in the bed. The model also takes into account the effect of steric hindrance limitations in SSF. The growth of cells around single particle and resulting expansion of biofilm around the particle is analyzed for simplified zero and first order oxygen consumption kinetics. Under conditions of zero order kinetics, the model predicts upper limit on cell density. The model simulations for packed bed of solid particles in tray bioreactor show distinct limitations on growth due to simultaneous heat and mass transport phenomena accompanying solid state fermentation process. The extent of limitation due to heat and/or mass transport phenomena is analyzed during different stages of fermentation. It is expected that the model will lead to better understanding of the transport processes in SSF, and therefore, will assist in optimal design of bioreactors for SSF.

An adaptive drug delivery design using neural networks for effective treatment of infectious diseases: A simulation study

An adaptive drug delivery design is presented in this paper using neural networks for effective treatment of infectious diseases. The generic mathematical model used describes the coupled evolution of concentration of pathogens, plasma cells, antibodies and a numerical value that indicates the relative characteristic of a damaged organ due to the disease under the influence of external drugs. From a system theoretic point of view, the external drugs can be interpreted as control inputs, which can be designed based on control theoretic concepts. In this study, assuming a set of nominal parameters in the mathematical model, first a nonlinear controller (drug administration) is designed based on the principle of dynamic inversion. This nominal drug administration plan was found to be effective in curing "nominal model patients" (patients whose immunological dynamics conform to the mathematical model used for the control design exactly. However, it was found to be ineffective in curing "realistic model patients" (patients whose immunological dynamics may have off-nominal parameter values and possibly unwanted inputs) in general. Hence, to make the drug delivery dosage design more effective for realistic model patients, a model-following adaptive control design is carried out next by taking the help of neural networks, that are trained online. Simulation studies indicate that the adaptive controller proposed in this paper holds promise in killing the invading pathogens and healing the damaged organ even in the presence of parameter uncertainties and continued pathogen attack. Note that the computational requirements for computing the control are very minimal and all associated computations (including the training of neural networks) can be carried out online. However it assumes that the required diagnosis process can be carried out at a sufficient faster rate so that all the states are available for control computation.

An optimal dynamic inversion-based neuro-adaptive approach for treatment of chronic myelogenous leukemia

Combining the advanced techniques of optimal dynamic inversion and model-following neuro-adaptive control design, an innovative technique is presented to design an automatic drug administration strategy for effective treatment of chronic myelogenous leukemia (CML). A recently developed nonlinear mathematical model for cell dynamics is used to design the controller (medication dosage). First, a nominal controller is designed based on the principle of optimal dynamic inversion. This controller can treat the nominal model patients (patients who can be described by the mathematical model used here with the nominal parameter values) effectively. However, since the system parameters for a realistic model patient can be different from that of the nominal model patients, simulation studies for such patients indicate that the nominal controller is either inefficient or, worse, ineffective; i.e. the trajectory of the number of cancer cells either shows non-satisfactory transient behavior or it grows in an unstable manner. Hence, to make the drug dosage history more realistic and patient-specific, a model-following neuro-adaptive controller is augmented to the nominal controller. In this adaptive approach, a neural network trained online facilitates a new adaptive controller. The training process of the neural network is based on Lyapunov stability theory, which guarantees both stability of the cancer cell dynamics as well as boundedness of the network weights. From simulation studies, this adaptive control design approach is found to be very effective to treat the CML disease for realistic patients. Sufficient generality is retained in the mathematical developments so that the technique can be applied to other similar nonlinear control design problems as well.

Testing and Modeling of MR Damper and Its Application to SDOF Systems Using Integral Backstepping Technique

Magnetorheological dampers are intrinsically nonlinear devices, which make the modeling and design of a suitable control algorithm an interesting and challenging task. To evaluate the potential of magnetorheological (MR) dampers in control applications and to take full advantages of its unique features, a mathematical model to accurately reproduce its dynamic behavior has to be developed and then a proper control strategy has to be taken that is implementable and can fully utilize their capabilities as a semi-active control device. The present paper focuses on both the aspects. First, the paper reports the testing of a magnetorheological damper with an universal testing machine, for a set of frequency, amplitude, and current. A modified Bouc-Wen model considering the amplitude and input current dependence of the damper parameters has been proposed. It has been shown that the damper response can be satisfactorily predicted with this model. Second, a backstepping based nonlinear current monitoring of magnetorheological dampers for semi-active control of structures under earthquakes has been developed. It provides a stable nonlinear magnetorheological damper current monitoring directly based on system feedback such that current change in magnetorheological damper is gradual. Unlike other MR damper control techniques available in literature, the main advantage of the proposed technique lies in its current input prediction directly based on system feedback and smooth update of input current. Furthermore, while developing the proposed semi-active algorithm, the dynamics of the supplied and commanded current to the damper has been considered. The efficiency of the proposed technique has been shown taking a base isolated three story building under a set of seismic excitation. Comparison with widely used clipped-optimal strategy has also been shown.

Free convection heat transfer between vertical parallel plates

A new mathematical model for the solution of the problem of free convection heat transfer between vertical parallel flat isothermal plates under isothermal boundary conditions, has been presented. The set of boundary layer equations used in the model are transformed to nonlinear coupled differential equations by similarity type variables as obtained by Ostrach for vertical flat plates in an infinite fluid medium. By utilising a parameter ηw* to represent the outer boundary, the governing differential equations are solved numerically for parametric values of Pr = 0.733. 2 and 3, and ηw* = 0.1, 0.5, 1, 2, 3, 4, ... and 8.0. The velocity and temperature profiles are presented. Results indicate that ηw* can effectively classify the system into (1) thin layers where conduction predominates, (2) intermediate layers and (3) thick layers whose results can be predicted by the solutions for vertical flat plates in infinite fluid medium. Heat transfer correlations are presented for the 3 categories. Several experimental and analytical results available in the literature agree with the present correlations.

Compliance of rough cylinders in compression

Compression of a rough turned cylinder between two hard, smooth, flat plates has been analysed with the aid of a mathematical model based on statistical analysis. It is assumed that the asperity peak heights follow Gaussian or normal and beta distribution functions and that the loaded asperities comply as though they are completely isolated from the neighbouring ones. Equations have been developed for the loadcompliance relation of the real surface using a simplified relation of the form W0 = K1δn for the load-compliance of a single asperity. Parameters K1 and n have considerable influence on the load-compliance curve and they depend on the material, tip angle of the asperity, standard deviation of the asperity peak height distribution and the density of the asperities.