Unsteady mixed convection flow in stagnation region adjacent to a vertical surface

The unsteady two-dimensional laminar mixed convection flow in the stagnation region of a vertical surface has been studied where the buoyancy forces are due to both the temperature and concentration gradients. The unsteadiness in the flow and temperature fields is caused by the time-dependent free stream velocity. Both arbitrary wall temperature and concentration, and arbitrary surface heat and mass flux variations have been considered. The Navier-Stokes equations, the energy equation and the concentration equation, which are coupled nonlinear partial differential equations with three independent variables, have been reduced to a set of nonlinear ordinary differential equations. The analysis has also been done using boundary layer approximations and the difference between the solutions has been discussed. The governing ordinary differential equations for buoyancy assisting and buoyancy opposing regions have been solved numerically using a shooting method. The skin friction, heat transfer and mass transfer coefficients increase with the buoyancy parameter. However, the skin friction coefficient increases with the parameter lambda, which represents the unsteadiness in the free stream velocity, but the heat and mass transfer coefficients decrease. In the case of buoyancy opposed flow, the solution does not exist beyond a certain critical value of the buoyancy parameter. Also, for a certain range of the buoyancy parameter dual solutions exist.

Nonaxisymmetric unsteady motion over a rotating disk in the presence of free convection and magnetic field

The nonaxisymmetric unsteady motion produced by a buoyancy-induced cross-flow of an electrically conducting fluid over an infinite rotating disk in a vertical plane and in the presence of an applied magnetic field normal to the disk has been studied. Both constant wall and constant heat flux conditions have been considered. It has been found that if the angular velocity of the disk and the applied magnetic field squared vary inversely as a linear function of time (i.e. as (1??t*)?1, the governing Navier-Stokes equation and the energy equation admit a locally self-similar solution. The resulting set of ordinary differential equations has been solved using a shooting method with a generalized Newton's correction procedure for guessed boundary conditions. It is observed that in a certain region near the disk the buoyancy induced cross-flow dominates the primary von Karman flow. The shear stresses induced by the cross-flow are found to be more than these of the primary flow and they increase with magnetic parameter or the parameter ? characterizing the unsteadiness. The velocity profiles in the x- and y-directions for the primary flow at any two values of the unsteady parameter ? cross each other towards the edge of the boundary layer. The heat transfer increases with the Prandtl number but reduces with the magnetic parameter.

Computation of restoration of ligand response in the random kinetics of a prostate cancer cell signaling pathway

Mutation and/or dysfunction of signaling proteins in the mitogen activated protein kinase (MAPK) signal transduction pathway are frequently observed in various kinds of human cancer. Consistent with this fact, in the present study, we experimentally observe that the epidermal growth factor (EGF) induced activation profile of MAP kinase signaling is not straightforward dose-dependent in the PC3 prostate cancer cells. To find out what parameters and reactions in the pathway are involved in this departure from the normal dose-dependency, a model-based pathway analysis is performed. The pathway is mathematically modeled with 28 rate equations yielding those many ordinary differential equations (ODE) with kinetic rate constants that have been reported to take random values in the existing literature. This has led to us treating the ODE model of the pathways kinetics as a random differential equations (RDE) system in which the parameters are random variables. We show that our RDE model captures the uncertainty in the kinetic rate constants as seen in the behavior of the experimental data and more importantly, upon simulation, exhibits the abnormal EGF dose-dependency of the activation profile of MAP kinase signaling in PC3 prostate cancer cells. The most likely set of values of the kinetic rate constants obtained from fitting the RDE model into the experimental data is then used in a direct transcription based dynamic optimization method for computing the changes needed in these kinetic rate constant values for the restoration of the normal EGF dose response. The last computation identifies the parameters, i.e., the kinetic rate constants in the RDE model, that are the most sensitive to the change in the EGF dose response behavior in the PC3 prostate cancer cells. The reactions in which these most sensitive parameters participate emerge as candidate drug targets on the signaling pathway. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Unsteady free convection flow under the influence of a magnetic field

The unsteady free convection boundary layer at the stagnation point of a two-dimensional body and an axisymmetric body with prescribed surface heat flux or temperature has been studied. The magnetic field is applied parallel to the surface and the effect of induced magnetic field has been considered. It is found that for certain powerlaw distribution of surface heat flux or temperature and magnetic field with time, the governing boundary layer equations admit a self-similar solution locally. The resulting nonlinear ordinary differential equations have been solved using a finite element method and a shooting method with Newton's corrections for missing initial conditions. The results show that the skin friction and heat transfer coefficients, and x-component of the induced magnetic field on the surface increase with the applied magnetic field. In general, the skin friction, heat transfer and x-component of the induced magnetic field for axisymmetric case are more than those of the two-dimensional case. Also they change more when the surface heat flux or temperature decreases with time than when it increases with time. The skin friction, heat transfer and x-component of the induced magnetic field are significantly affected by the magnetic Prandtl number and they increase as the magnetic Prandtl number decreases. The skin friction and x-component of the magnetic field increase with the dissipation parameter, but heat transfer decreases.

Variable density approach for rotating shallow shell of variable thickness

A new approach based on variable density in conjunction with shallow shell theory is proposed to analyse rotating shallow shell of variable thickness. Coupled non-linear ordinary differential equations governing shallows shells of variable thickness are first derived before applying the variable density approach. Results obtained from the new approach compare well with FEM calculation for a wide range of profiles considered in this paper.

Singularity structure and chaotic behavior of the homopolar disk dynamo

We study in great detail a system of three first-order ordinary differential equations describing a homopolar disk dynamo (HDD). This system displays a large variety of behaviors, both regular and chaotic. Existence of periodic solutions is proved for certain ranges of parameters. Stability criteria for periodic solutions are given. The nonintegrability aspects of the HDD system are studied by investigating analytically the singularity structure of the system in the complex domain. Coexisting attractors (including period-doubling sequence) and coexisting strange attractors appear in some parametric regimes. The gluing of strange attractors and the ungluing of a strange attractor are also shown to occur. A period of bifurcation leading to chaos, not observed for other chaotic systems, is shown to characterize the chaotic behavior in some parametric ranges. The limiting case of the Lorenz system is also studied and is related to HDD.

Symplectic integration of nonlinear Hamiltonian systems

There exist several standard numerical methods for integrating ordinary differential equations. However, if one is interested in integration of Hamiltonian systems, these methods can lead to wrong results. This is due to the fact that these methods do not explicitly preserve the so-called 'symplectic condition' (that needs to be satisfied for Hamiltonian systems) at every integration step. In this paper, we look at various methods for integration that preserve the symplectic condition.

Nonlinear dynamics and chaotic motions in feedback-controlled two- and three-degree-of-freedom robots

The dynamics of a feedback-controlled rigid robot is most commonly described by a set of nonlinear ordinary differential equations. In this paper we analyze these equations, representing the feedback-controlled motion of two- and three-degrees-of-freedom rigid robots with revolute (R) and prismatic (P) joints in the absence of compliance, friction, and potential energy, for the possibility of chaotic motions. We first study the unforced or inertial motions of the robots, and show that when the Gaussian or Riemannian curvature of the configuration space of a robot is negative, the robot equations can exhibit chaos. If the curvature is zero or positive, then the robot equations cannot exhibit chaos. We show that among the two-degrees-of-freedom robots, the PP and the PR robot have zero Gaussian curvature while the RP and RR robots have negative Gaussian curvatures. For the three-degrees-of-freedom robots, we analyze the two well-known RRP and RRR configurations of the Stanford arm and the PUMA manipulator respectively, and derive the conditions for negative curvature and possible chaotic motions. The criteria of negative curvature cannot be used for the forced or feedback-controlled motions. For the forced motion, we resort to the well-known numerical techniques and compute chaos maps, Poincare maps, and bifurcation diagrams. Numerical results are presented for the two-degrees-of-freedom RP and RR robots, and we show that these robot equations can exhibit chaos for low controller gains and for large underestimated models. From the bifurcation diagrams, the route to chaos appears to be through period doubling.

Right-Definite Multiparameter Sturm–Liouville Problems With Eigenparameter-Dependent Boundary Conditions

We study a system of ordinary differential equations linked by parameters and subject to boundary conditions depending on parameters. We assume certain definiteness conditions on the coefficient functions and on the boundary conditions that yield, in the corresponding abstract setting, a right-definite case. We give results on location of the eigenvalues and oscillation of the eigenfunctions.

Tracking analysis of an LMS decision feedback equalizer for a wireless channel

We consider a time varying wireless fading channel, equalized by an LMS Decision Feedback equalizer (DFE). We study how well this equalizer tracks the optimal MMSEDFE (Wiener) equalizer. We model the channel by an Autoregressive (AR) process. Then the LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs (ordinary differential equations). Using these ODEs, we show via some examples that the LMS equalizer moves close to the instantaneous Wiener filter after initial transience. We also compare the LMS equalizer with the instantaneous optimal DFE (the commonly used Wiener filter) designed assuming perfect previous decisions and computed using perfect channel estimate (we will call it as IDFE). We show that the LMS equalizer outperforms the IDFE almost all the time after initial transience.

Tracking performance of an LMS-Linear Equalizer for fading channels

We consider a time varying wireless fading channel, equalized by an LMS linear equalizer. We study how well this equalizer tracks the optimal Wiener equalizer. We model the channel by an Auto-regressive (AR) process. Then the LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs (ordinary differential equations). Using these ODEs, the error between the LMS equalizer and the instantaneous Wiener filter is shown to decay exponentially/polynomially to zero unless the channel is marginally stable in which case the convergence may not hold.Using the same ODEs, we also show that the corresponding Mean Square Error (MSE) converges towards minimum MSE(MMSE) at the same rate for a stable channel. We further show that the difference between the MSE and the MMSE does not explode with time even when the channel is unstable. Finally we obtain an optimum step size for the linear equalizer in terms of the AR parameters, whenever the error decay is exponential.

Noise reduction with perforated three-duct muffler components

This paper describes the authors’ distributed parameter approach for derivation of closed-form expressions for the four-pole parameters of the perforated three-duct muffler components. In this method, three simultaneous second-order partial differential equations are first reduced to a set of six first-order ordinary differential equations. These equations are then uncoupled by means of a modal matrix. The resulting 6 × 6 matrix is reduced to the 2 × 2 transfer matrix using the relevant boundary conditions. This is combined with transfer matrices of other elements (upstream and downstream of this perforated element) to predict muffler performance like noise reduction, which is also measured. The correlation between experimental and theoretical values of noise reduction is shown to be satisfactory.

Co-evolution of content popularity and delivery in mobile P2P networks

Mobile P2P technology provides a scalable approach for content delivery to a large number of users on their mobile devices. In this work, we study the dissemination of a single item of content (e. g., an item of news, a song or a video clip) among a population of mobile nodes. Each node in the population is either a destination (interested in the content) or a potential relay (not yet interested in the content). There is an interest evolution process by which nodes not yet interested in the content (i.e., relays) can become interested (i.e., become destinations) on learning about the popularity of the content (i.e., the number of already interested nodes). In our work, the interest in the content evolves under the linear threshold model. The content is copied between nodes when they make random contact. For this we employ a controlled epidemic spread model. We model the joint evolution of the copying process and the interest evolution process, and derive joint fluid limit ordinary differential equations. We then study the selection of parameters under the content provider's control, for the optimization of various objective functions that aim at maximizing content popularity and efficient content delivery.

Steady mixed convection flow of Maxwell fluid over an exponentially stretching vertical surface with magnetic field and viscous dissipation

The steady mixed convection flow and heat transfer from an exponentially stretching vertical surface in a quiescent Maxwell fluid in the presence of magnetic field, viscous dissipation and Joule heating have been studied. The stretching velocity, surface temperature and magnetic field are assumed to have specific exponential function forms for the existence of the local similarity solution. The coupled nonlinear ordinary differential equations governing the local similarity flow and heat transfer have been solved numerically by Chebyshev finite difference method. The influence of the buoyancy parameter, viscous dissipation, relaxation parameter of Maxwell fluid, magnetic field and Prandtl number on the flow and heat transfer has been considered in detail. The Nusselt number increases significantly with the Prandtl number, but the skin friction coefficient decreases. The Nusselt number slightly decreases with increasing viscous dissipation parameter, but the skin friction coefficient slightly increases. Maxwell fluid reduces both skin friction coefficient and Nusselt number, whereas buoyancy force enhances them.

Phase space path integral approach to harmonic oscillator with a time-dependent force constant

The quantum statistical mechanical propagator for a harmonic oscillator with a time-dependent force constant, m omega(2)(t), has been investigated in the past and was found to have only a formal solution in terms of the solutions of certain ordinary differential equations. Such path integrals are frequently encountered in semiclassical path integral evaluations and having exact analytical expressions for such path integrals is of great interest. In a previous work, we had obtained the exact propagator for motion in an arbitrary time-dependent harmonic potential in the overdamped limit of friction using phase space path integrals in the context of Levy flights - a result that can be easily extended to Brownian motion. In this paper, we make a connection between the overdamped Brownian motion and the imaginary time propagator of quantum mechanics and thereby get yet another way to evaluate the latter exactly. We find that explicit analytic solution for the quantum statistical mechanical propagator can be written when the time-dependent force constant has the form omega(2)(t) = lambda(2)(t) - d lambda(t)/dt where lambda(t) is any arbitrary function of t and use it to evaluate path integrals which have not been evaluated previously. We also employ this method to arrive at a formal solution of the propagator for both Levy flights and Brownian subjected to a time-dependent harmonic potential in the underdamped limit of friction. (C) 2015 Elsevier B.V. All rights reserved.