Boundary value problems for third-order nonlinear ordinary differential equations

In this paper, we describe how to analyze boundary value problems for third-order nonlinear ordinary differential equations over an infinite interval. Several physical problems of interest are governed by such systems. The seminumerical schemes described here offer some advantages over solutions obtained by using traditional methods such as finite differences, shooting method, etc. These techniques also reveal the analytic structure of the solution function. For illustrative purposes, several physical problems, mainly drawn from fluid mechanics, are considered; they clearly demonstrate the efficiency of the techniques presented here.

Boundary value problems for third-order nonlinear ordinary differential equations

In this paper, we describe how to analyze boundary value problems for third-order nonlinear ordinary differential equations over an infinite interval. Several physical problems of interest are governed by such systems. The seminumerical schemes described here offer some advantages over solutions obtained by using traditional methods such as finite differences, shooting method, etc. These techniques also reveal the analytic structure of the solution function. For illustrative purposes, several physical problems, mainly drawn from fluid mechanics, are considered; they clearly demonstrate the efficiency of the techniques presented here.

Numerical procedure for second order non-linear ordinary differential equations and application to heat transfer problem

In this paper, we have first given a numerical procedure for the solution of second order non-linear ordinary differential equations of the type y″ = f (x;y, y′) with given initial conditions. The method is based on geometrical interpretation of the equation, which suggests a simple geometrical construction of the integral curve. We then translate this geometrical method to the numerical procedure adaptable to desk calculators and digital computers. We have studied the efficacy of this method with the help of an illustrative example with known exact solution. We have also compared it with Runge-Kutta method. We have then applied this method to a physical problem, namely, the study of the temperature distribution in a semi-infinite solid homogeneous medium for temperature-dependent conductivity coefficient.

A parallel multistep predictor-corrector algorithm for solving ordinary differential equations

In this paper we give a generalized predictor-corrector algorithm for solving ordinary differential equations with specified initial values. The method uses multiple correction steps which can be carried out in parallel with a prediction step. The proposed method gives a larger stability interval compared to the existing parallel predictor-corrector methods. A method has been suggested to implement the algorithm in multiple processor systems with efficient utilization of all the processors.

Spectral Invariance of SG Pseudo-Differential Operators on L-P(R-n)

We prove the spectral invariance of SG pseudo-differential operators on L-P(R-n), 1 < p < infinity, by using the equivalence of ellipticity and Fredholmness of SG pseudo-differential operators on L-p(R-n), 1 < p < infinity. A key ingredient in the proof is the spectral invariance of SC pseudo-differential operators on L-2(R-n).

Adaptive mesh generation for singularly perturbed fourth-order ordinary differential equations

In this paper, we consider a singularly perturbed boundary-value problem for fourth-order ordinary differential equation (ODE) whose highest-order derivative is multiplied by a small perturbation parameter. To solve this ODE, we transform the differential equation into a coupled system of two singularly perturbed ODEs. The classical central difference scheme is used to discretize the system of ODEs on a nonuniform mesh which is generated by equidistribution of a positive monitor function. We have shown that the proposed technique provides first-order accuracy independent of the perturbation parameter. Numerical experiments are provided to validate the theoretical results.

Jacobi forms and differential operators

We affirmatively answer a question due to S. Bocherer concerning the feasibility of removing one differential operator from the standard collection of m + 1 of them used to embed the space of Jacobi forms of weight 2 and index m into several pieces of elliptic modular forms. (C) 2014 Elsevier Inc. All rights reserved.

Unsteady-Flow Of A Viscous-Fluid Between 2 Parallel Disks With A Time-Varying Gap Width And A Magnetic-Field

The unsteady incompressible viscous fluid flow between two parallel infinite disks which are located at a distance h(t*) at time t* has been studied. The upper disk moves towards the lower disk with velocity h'(t*). The lower disk is porous and rotates with angular velocity Omega(t*). A magnetic field B(t*) is applied perpendicular to the two disks. It has been found that the governing Navier-Stokes equations reduce to a set of ordinary differential equations if h(t*), a(t*) and B(t*) vary with time t* in a particular manner, i.e. h(t*) = H(1 - alpha t*)(1/2), Omega(t*) = Omega(0)(1 - alpha t*)(-1), B(t*) = B-0(1 - alpha t*)(-1/2). These ordinary differential equations have been solved numerically using a shooting method. For small Reynolds numbers, analytical solutions have been obtained using a regular perturbation technique. The effects of squeeze Reynolds numbers, Hartmann number and rotation of the disk on the flow pattern, normal force or load and torque have been studied in detail

Unsteady flow over a stretching surface with a magnetic field in a rotating fluid

The effect of the magnetic field on the unsteady flow over a stretching surface in a rotating fluid has been studied. The unsteadiness in the flow field is due to the time-dependent variation of the velocity of the stretching surface and the angular velocity of the rotating fluid. The Navier-Stokes equations and the energy equation governing the flow and the heat transfer admit a self-similar solution if the velocity of the stretching surface and the angular velocity of the rotating fluid vary inversely as a linear function of time. The resulting system of ordinary differential equations is solved numerically using a shooting method. The rotation parameter causes flow reversal in the component of the velocity parallel to the strerching surface and the magnetic field tends to prevent or delay the flow reversal. The surface shear stresses dong the stretching surface and in the rotating direction increase with the rotation parameter, but the surface heat transfer decreases. On the other hand, the magnetic field increases the surface shear stress along the stretching surface, but reduces the surface shear stress in the rotating direction and the surface heat transfer. The effect of the unsteady parameter is more pronounced on the velocity profiles in the rotating direction and temperature profiles.

A pseudo-dynamical systems approach to a class of inverse problems in engineering

A pseudo-dynamical approach for a class of inverse problems involving static measurements is proposed and explored. Following linearization of the minimizing functional associated with the underlying optimization problem, the new strategy results in a system of linearized ordinary differential equations (ODEs) whose steady-state solutions yield the desired reconstruction. We consider some explicit and implicit schemes for integrating the ODEs and thus establish a deterministic reconstruction strategy without an explicit use of regularization. A stochastic reconstruction strategy is then developed making use of an ensemble Kalman filter wherein these ODEs serve as the measurement model. Finally, we assess the numerical efficacy of the developed tools against a few linear and nonlinear inverse problems of engineering interest.

Stability Of Large Damped Flexible Spacecraft With Stored Angular Momentum

Vibrational stability of large flexible structurally damped spacecraft carrying internal angular momentum and undergoing large rigid body rotations is analysed modeling the systems as elastic continua. Initially, analytical solutions to the motion of rigid gyrostats under torque-free conditions are developed. The solutions to the gyrostats modeled as axisymmetric and triaxial spacecraft carrying three and two constant speed momentum wheels, respectively, with spin axes aligned with body principal axes are shown to be complicated. These represent extensions of solutions for simpler cases existing in the literature. Using these solutions and modal analysis, the vibrational equations are reduced to linear ordinary differential equations. Equations with periodically varying coefficients are analysed applying Floquet theory. Study of a few typical beam- and plate-like spacecraft configurations indicate that the introduction of a single reaction wheel into an axisymmetric satellite does not alter the stability criterion. However, introduction of constant speed rotors deteriorates vibrational stability. Effects of structural damping and vehicle inertia ratio are also studied.

Self-similar solution of unsteady compressible three-dimensional stagnation-point boundary layers

The self-similar solution of the unsteady laminar compressible boundary-layer flow with variable properties at a three-dimensional stagnation point with mass transfer has been obtained when the free-stream velocity varies inversely as a linear function of time. The resulting ordinary differential equations have been solved numerically using an implicit finite-difference scheme. The results are found to be strongly dependent on the parameter characterizing the unsteadiness in the free-stream velocity. The velocity profiles show some features not encountered in steady flows.

Oscillatory flow in an elastic tube of variable cross-section

An oscillatory flow of a viscous incompressible fluid in an elastic tube of variable cross section has been investigated at low Reynolds number. The equations governing, the flow are derived under the assumption that the variation of the cross-section is slow in the axial direction for a tethered tube. The problem is then reduced to that of solving for the excess pressure from a second order ordinary differential equation with complex valued Bessel functions as the coefficients. This equation has been solved numerically for geometries of physiological interest and a comparison is made with some of the known theoretical and experimental results.

Mass transfer effects on the generalised vortex flow over a stationary surface with or without magnetic field

The effect of injection and suction on the generalised vortex flow of a steady laminar incompressible fluid over a stationary infinite disc with or without magnetic field under boundary-layer approximations has been studied. The coupled nonlinear ordinary differential equations governing the self-similar flow have been numerically solved using the finite-difference scheme. The results indicate that the injection produces a deeper inflow layer and de-stabilises the motion while suction or magnetic field suppresses the inflow layer and produces stability. The effect of decreasingn, the parameter characterising the nature of vortex flow, is similar to that of increasing the injection rate.