Computing and Visualizing Solution Sets of Interval Linear Systems

The paper has been presented at the 12th International Conference on Applications of Computer Algebra, Varna, Bulgaria, June, 2006

A point interpolation method with locally smoothed strain field (PIM-LS2) for mechanics problems using triangular mesh

A point interpolation method with locally smoothed strain field (PIM-LS2) is developed for mechanics problems using a triangular background mesh. In the PIM-LS2, the strain within each sub-cell of a nodal domain is assumed to be the average strain over the adjacent sub-cells of the neighboring element sharing the same field node. We prove theoretically that the energy norm of the smoothed strain field in PIM-LS2 is equivalent to that of the compatible strain field, and then prove that the solution of the PIM- LS2 converges to the exact solution of the original strong form. Furthermore, the softening effects of PIM-LS2 to system and the effects of the number of sub-cells that participated in the smoothing operation on the convergence of PIM-LS2 are investigated. Intensive numerical studies verify the convergence, softening effects and bound properties of the PIM-LS2, and show that the very ‘‘tight’’ lower and upper bound solutions can be obtained using PIM-LS2.

The Galerkin finite element approximation of the fractional cable equation

The cable equation is one of the most fundamental equations for modeling neuronal dynamics. Cable equations with a fractional order temporal derivative have been introduced to model electrotonic properties of spiny neuronal dendrites. In this paper, the fractional cable equation involving two integro-differential operators is considered. The Galerkin finite element approximations of the fractional cable equation are proposed. The main contribution of this work is outlined as follow: • A semi-discrete finite difference approximation in time is proposed. We prove that the scheme is unconditionally stable, and the numerical solution converges to the exact solution with order O(Δt). • A semi-discrete difference scheme for improving the order of convergence for solving the fractional cable equation is proposed, and the numerical solution converges to the exact solution with order O((Δt)2). • Based on the above semi-discrete difference approximations, Galerkin finite element approximations in space for a full discretization are also investigated. • Finally, some numerical results are given to demonstrate the theoretical analysis.

Renormalization group interpretation of the Born and Rytov approximations

In this paper the method of renormalization group (RG) [Phys. Rev. E 54, 376 (1996)] is related to the well-known approximations of Rytov and Born used in wave propagation in deterministic and random media. Certain problems in linear and nonlinear media are examined from the viewpoint of RG and compared with the literature on Born and Rytov approximations. It is found that the Rytov approximation forms a special case of the asymptotic expansion generated by the RG, and as such it gives a superior approximation to the exact solution compared with its Born counterpart. Analogous conclusions are reached for nonlinear equations with an intensity-dependent index of refraction where the RG recovers the exact solution. © 2008 Optical Society of America.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

Static analysis of infinite laminated cylindrical shell panel

This paper presents an approximate three-dimensional elasticity solution for an infinitely long, cross-ply laminated circular cylindrical shell panel with simply supported boundary conditions, subjected to an arbitrary discontinuous transverse loading. The solution is based on the principal assumption that the ratio of the thickness of the lamina to its middle surface radius is negligible compared to unity. The validity of this assumption and the range of application of this approximate solution have been established through a comparison with an exact solution. Results of classical and first-order shear deformation shell theories have been compared with the results of the present solution to bring out the accuracy of these theories. It is also shown that for very shallow shell panels the definition of a thin shell should be based on the ratio of thickness to chord width rather than the ratio of thickness to mean radius.

Flow and heat transfer in a stagnation-point flow over a stretching sheet with a magnetic field

The flow and heat transfer problem in the boundary layer induced by a continuous moving surface is important in many manufacturing processes in industry such as the boundary layer along material handling conveyers, the aerodynamic extrusion of plastic sheet, the cooling of an infinite metalic plate in a cooling bath (which may also be electrolyte). Glass blowing, continuous casting and spinning of fibres also involve the flow due to a stretching surface. Sakiadis [1] was the first to study the flow induced by a semi-infinite moving wall in an ambient fluid. On the other hand, Crane [2] first studied the flow over a linearly stretching sheet in an ambient fluid. Subsequently, Crane [3] also investigated the corresponding heat transfer problem. Since then several authors [4-8] have studied various aspects of this problem such as the effects of mass transfer, variable wall temperature, constant heat flux, magnetic field etc. Recently, Andersson [9] has obtained an exact solution of the Navier-Stokes equations for the MHD flow over a linearly stretching sheet in an ambient fluid. Also Chiam [10] has studied the heat transfer with variable thermal conductivity on a stretching sheet when the velocities of the sheet and the free stream are equal.

MHD Flow of a 2nd Grade Fluid in an Orthogonal Rheometer

The magnetohydrodynamics (MHD) flow of a conducting, homogeneous incompressible Rivlin-Ericksen fluid of second grade contained between two infinite, parallel, insulated disks rotating with the same angular velocity about two noncoincident axes, under the application of a uniform transverse magnetic field, is investigated. This model represents the MHD flow of the fluid in the instrument called an orthogonal rheometer, except for the fact that in the rheometer the rotating plates are necessarily finite. An exact solution of the governing equations of motion is presented. The force components in the x and y directions on the disks are calculated. The effects of magnetic field and the viscoelastic parameter on the forces are discussed in detail.

Survival probability for a diffusive process on a growing domain

We consider the motion of a diffusive population on a growing domain, 0 < x < L(t ), which is motivated by various applications in developmental biology. Individuals in the diffusing population, which could represent molecules or cells in a developmental scenario, undergo two different kinds of motion: (i) undirected movement, characterized by a diffusion coefficient, D, and (ii) directed movement, associated with the underlying domain growth. For a general class of problems with a reflecting boundary at x = 0, and an absorbing boundary at x = L(t ), we provide an exact solution to the partial differential equation describing the evolution of the population density function, C(x,t ). Using this solution, we derive an exact expression for the survival probability, S(t ), and an accurate approximation for the long-time limit, S = limt→∞ S(t ). Unlike traditional analyses on a nongrowing domain, where S ≡ 0, we show that domain growth leads to a very different situation where S can be positive. The theoretical tools developed and validated in this study allow us to distinguish between situations where the diffusive population reaches the moving boundary at x = L(t ) from other situations where the diffusive population never reaches the moving boundary at x = L(t ). Making this distinction is relevant to certain applications in developmental biology, such as the development of the enteric nervous system (ENS). All theoretical predictions are verified by implementing a discrete stochastic model.

Finite-amplitude love-waves on an isotropic layered half-space

A pair of semi-linear hyperbolic partial differential equations governing the slow variations in amplitude and phase of a quasi-monochromatic finite-amplitude Love-wave on an isotropic layered half-space is derived using the method of multiple-scales. The analysis of the exact solution of these equations for a signalling problem reveals that the amplitude of the wave remains constant along its characteristic and that the phase of the wave increases linearly behind the wave-front.

Axisymmetric stress distribution in a transversely isotropic thick plate with a circular hole

An exact solution for the stresses in a transversely isotropic infinite thick plate having a circular hole and subjected to axisymmetric uniformly distributed load on the plane surfaces has been given. The solution is in the form of Fourier-Bessel series and integrals. Numerical results for the stresses are given using the elastic constants for magnesium, and are compared with the isotropic case.

Vibration of Non-Uniform Thin-Walled Beams of Arbitrary Shape

A generalised theory for the natural vibration of non-uniform thin-walled beams of arbitrary cross-sectional geometry is proposed. The governing equations are obtained as four partial, linear integro-differential equations. The corresponding boundary conditions are also obtained in an integro-differential form. The formulation takes into account the effect of longitudinal inertia and shear flexibility. A method of solution is presented. Some numerical illustrations and an exact solution are included.

Stiff-String Basis Functions for Vibration Analysis of High Speed Rotating Beams

A new rotating beam finite element is developed in which the basis functions are obtained by the exact solution of the governing static homogenous differential equation of a stiff string, which results from an approximation in the rotating beam equation. These shape functions depend on rotation speed and element position along the beam and account for the centrifugal stiffening effect. Using this new element and the Hermite cubic finite element, a convergence study of natural frequencies is performed, and it is found that the new element converges much more rapidly than the conventional Hermite cubic element for the first two modes at higher rotation speeds. The new element is also applied for uniform and tapered rotating beams to determine the natural frequencies, and the results compare very well with the published results given in the literature.

Approximate solutions for the structure of shock waves

Approximate solutions of the B-G-K model equation are obtained for the structure of a plane shock, using various moment methods and a least squares technique. Comparison with available exact solution shows that while none of the methods is uniformly satisfactory, some of them can provide accurate values for the density slope shock thickness delta n . A detailed error analysis provides explanations for this result. An asymptotic analysis of delta n for largeMach numbers shows that it scales with theMaxwell mean free path on the hot side of the shock, and that their ratio is relatively insensitive to the viscosity law for the gas.

On the design of quadratic weirs

This paper considers the problem of the design of the quadratic weir notch, which finds application in the proportionate method of flow measurement in a by-pass, such that the discharge through it is proportional to the square root of the head measured above a certain datum. The weir notch consists of a bottom in the form of a rectangular weir of width 2W and depth a over which a designed curve is fitted. A theorem concerning the flow through compound weirs called the “slope discharge continuity theorem” is discussed and proved. Using this, the problem is reduced to the determination of an exact solution to Volterra's integral equation in Abel's form. It is shown that in the case of a quadratic weir notch, the discharge is proportional to the square root of the head measured above a datum Image a above the crest of the weir. Further, it is observed that the function defining the shape of the weir is rapidly convergent and its value almost approximates to zero at distances of 3a and above from the crest of the weir. This interesting and significant behaviour of the function incidentally provides a very good approximate solution to a particular Fredholm integral equation of the first kind, transforming the notch into a device called a “proportional-orifice”. A new concept of a “notch-orifice” capable of passing a discharge proportional to the square root of the head (above a particular datum) while acting both as a notch, and as an orifice, is given. A typical experiment with one such notch-orifice, having A = 4 in., and W = 6 in., shows a remarkable agreement with the theory and is found to have a constant coefficient of discharge of 0.61 in the ranges of both notch and orifice.