Analytical solutions for the multi-term time-fractional diffusion-wave/diffusion equations in a finite domain

Multi-term time-fractional differential equations have been used for describing important physical phenomena. However, studies of the multi-term time-fractional partial differential equations with three kinds of nonhomogeneous boundary conditions are still limited. In this paper, a method of separating variables is used to solve the multi-term time-fractional diffusion-wave equation and the multi-term time-fractional diffusion equation in a finite domain. In the two equations, the time-fractional derivative is defined in the Caputo sense. We discuss and derive the analytical solutions of the two equations with three kinds of nonhomogeneous boundary conditions, namely, Dirichlet, Neumann and Robin conditions, respectively.

An Analysis of a Boundary Layer Diffusion Flame Over a Porous Flat Plate in a Confined Flow

A numerical analysis of the gas dynamic structure of a two-dimensional laminar boundary layer diffusion flame over a porous flat plate in a confined flow is made on the basis of the familiar boundary layer and flame sheet approximations neglecting buoyancy effects. The governing equations of aerothermochemistry with the appropriate boundary conditions are solved using the Patankar-Spalding method. The analysis predicts the flame shape, profiles of temperature, concentrations of variousspecies, and the density of the mixture across the boundary layer. In addition, it also predicts the pressure gradient in the flow direction arising from the confinement ofthe flow and the consequent velocity overshoot near the flame surface. The results of thecomputation performed for an n-pentane-air system are compared with experimental data andthe agreement is found to be satisfactory.

Heat transfer to power-law fluids in thermal entrance region with viscous dissipation for constant heat flux conditions

An analytical solution of the heat transfer problem with viscous dissipation for non-Newtonian fluids with power-law model in the thermal entrance region of a circular pipe and two parallel plates under constant heat flux conditions is obtained using eigenvalue approach by suitably replacing one of the boundary conditions by total energy balance equation. Analytical expressions for the wall and the bulk temperatures and the local Nusselt number are presented. The results are in close agreement with those obtained by implicit finite-difference scheme. It is found that the role of viscous dissipation on heat transfer is completely different for heating and cooling conditions at the wall. The results for the case of cooling at the wall are of interest in the design of the oil pipe line.

A semi-analytical solution for multilayer diffusion in a composite medium consisting of a large number of layers

Diffusion in a composite slab consisting of a large number of layers provides an ideal prototype problem for developing and analysing two-scale modelling approaches for heterogeneous media. Numerous analytical techniques have been proposed for solving the transient diffusion equation in a one-dimensional composite slab consisting of an arbitrary number of layers. Most of these approaches, however, require the solution of a complex transcendental equation arising from a matrix determinant for the eigenvalues that is difficult to solve numerically for a large number of layers. To overcome this issue, in this paper, we present a semi-analytical method based on the Laplace transform and an orthogonal eigenfunction expansion. The proposed approach uses eigenvalues local to each layer that can be obtained either explicitly, or by solving simple transcendental equations. The semi-analytical solution is applicable to both perfect and imperfect contact at the interfaces between adjacent layers and either Dirichlet, Neumann or Robin boundary conditions at the ends of the slab. The solution approach is verified for several test cases and is shown to work well for a large number of layers. The work is concluded with an application to macroscopic modelling where the solution of a fine-scale multilayered medium consisting of two hundred layers is compared against an “up-scaled” variant of the same problem involving only ten layers.

Effect of boundary condition description on convergence of solution in a boundary-value problem

It is well known that the numerical accuracy of a series solution to a boundary-value problem by the direct method depends on the technique of approximate satisfaction of the boundary conditions and on the stage of truncation of the series. On the other hand, it does not appear to be generally recognized that, when the boundary conditions can be described in alternative equivalent forms, the convergence of the solution is significantly affected by the actual form in which they are stated. The importance of the last aspect is studied for three different techniques of computing the deflections of simply supported regular polygonal plates under uniform pressure. It is also shown that it is sometimes possible to modify the technique of analysis to make the accuracy independent of the description of the boundary conditions.

Outward spherical solidification of a superheated melt with time dependent boundary flux

Short-time analytical solutions of solid and liquid temperatures and freezing front have been obtained for the outward radially symmetric spherical solidification of a superheated melt. Although results are presented here only for time dependent boundary flux, the method of solution can be used for other kinds of boundary conditions also. Later, the analytical solution has been compared with the numerical solution obtained with the help of a finite difference numerical scheme in which the grid points change with the freezing front position. An efficient method of execution of the numerical scheme has been discussed in details. Graphs have been drawn for the total solidification times and temperature distributions in the solid.

Non-Newtonian boundary layers on a moving cylinder

The analysis of steady laminar forced convection boundary layer of power-law non-Newtonian fluids on a continuously moving cylinder with the surface maintained at a uniform temperature or uniform heat flux is presented. Of interest were the effects of power-law viscosity index, transverse curvature, generalized Prandtl number and streamwise coordinate on the local Nusselt number as well as on the velocity and temperature profiles. The two thermal boundary conditions yield quite similar results. Comparison of the calculated results with available series expansion solutions for a Newtonian fluid shows a very good performance of the present numerical procedure.

Interannual variations of Indian summer monsoon in a GCM: External conditions versus internal feedbacks

The potential predictability of the Indian summer monsoon due to slowly varying sea surface temperature (SST) forcing is examined. Factors responsible for limiting the predictability are also investigated. Three multiyear simulations with the R30 version of the Geophysical Fluid Dynamics Laboratory's climate model are carried out for this purpose, The mean monsoon simulated by this model is realistic including the mean summer precipitation over the Indian continent. The interannual variability of the large-scale component of the monsoon such as the "monsoon shear index" and its teleconnection with Pacific SST is well simulated by the model in a 15-yr integration with observed SST as boundary condition. On regional scales, the skill in simulating the interannual variability of precipitation over the Indian continent by the model is rather modest and its simultaneous correlation with eastern Pacific SST is negative but poor as observed. The poor predictability of precipitation over the Indian region in the model is related to the fact that contribution to the interannual variability over this region due to slow SST variations [El Nino-Southern Oscillation (ENSO) related] is comparable to those due to regional-scale fluctuations unrelated to ENSO SST. The physical mechanism through which ENSO SST tend to produce reduction in precipitation over the Indian continent is also elucidated. A measure of internal variability of the model summer monsoon is obtained from a 20-yr integration of the same model with fixed annual cycle SST as boundary conditions but with predicted soil moisture and snow cover. A comparison of summer monsoon indexes between this run and the observed SST run shows that the internal oscillations can account for a large fraction of the simulated monsoon variability. The regional-scale oscillations in the observed SST run seems to arise from these internal oscillations. It is discovered that most of the interannual internal variability is due to an internal quasi-biennial oscillation (QBO) of the model atmosphere. Such a QBO is also found in the author's third 18-yr simulation in which fixed annual cycle of SST as well as soil moisture and snow cover are prescribed. This shows that the model QBO is not due to land-surface-atmosphere interaction. It is proposed that the model QBO arises due to an interaction between nonlinear intraseasonal oscillations and the annual cycle. Spatial structure of the QBO and its role in limiting the predictability of the Indian summer monsoon is discussed.

The Wiener-Hopf solution of a class of mixed boundary value problems arising in surface water wave phenomena

Two mixed boundary value problems associated with two-dimensional Laplace equation, arising in the study of scattering of surface waves in deep water (or interface waves in two superposed fluids) in the linearised set up, by discontinuities in the surface (or interface) boundary conditions, are handled for solution by the aid of the Weiner-Hopf technique applied to a slightly more general differential equation to be solved under general boundary conditions and passing on to the limit in a manner so as to finally give rise to the solutions of the original problems. The first problem involves one discontinuity while the second problem involves two discontinuities. The reflection coefficient is obtained in closed form for the first problem and approximately for the second. The behaviour of the reflection coefficient for both the problems involving deep water against the incident wave number is depicted in a number of figures. It is observed that while the reflection coefficient for the first problem steadily increases with the wave number, that for the second problem exhibits oscillatory behaviour and vanishes at some discrete values of the wave number. Thus, there exist incident wave numbers for which total transmission takes place for the second problem. (C) 1999 Elsevier Science B.V. All rights reserved.

Thermal instability in a three-dimensional rigid container with prescribed heat flux at lower boundary

The Bénard–Marangoni convection is studied in a three-dimensional container with thermally insulated lateral walls and prescribed heat flux at lower boundary. The upper surface of the incompressible, viscous fluid is assumed to be flat with temperature dependent surface tension. A Galerkin–Tau method with odd and even trial functions satisfying all the essential boundary conditions except the natural boundary conditions at the free surface has been used to solve the problem. The critical Marangoni and Rayleigh numbers are determined for the onset of steady convection as a function of aspect ratios x0 and y0 for the cases of Bénard–Marangoni, pure Marangoni and pure Bénard convections. It is observed that critical parameters are decreasing with an increase in aspect ratios. The flow structures corresponding to the values of the critical parameters are presented in all the cases. It is observed that the critical parameters are higher for case with heat flux prescribed than those corresponding to the case with prescribed temperature. The critical Marangoni number for pure Marangoni convection is higher than critical Rayleigh number corresponding to pure Bénard convection for a given aspect ratio whereas the reverse was observed for two-dimensional infinite layer.

Rotating beams and non-rotating beams with shared eigenpair for pinned-free boundary condition

In this paper we look for a rotating beam, with pinned-free boundary conditions, whose eigenpair (frequency and mode-shape) is same as that of a uniform non-rotating beam for a particular mode. It is seen that for any given mode, there exists a flexural stiffness function (FSF) for which the ith mode eigenpair of a rotating beam with uniform mass distribution, is identical to that of a corresponding non-rotating beam with same length and mass distribution. Inserting these derived FSF's in a finite element code for a rotating pinned-free beam, the frequencies and mode shapes of a non-rotating pinned-free beam are obtained. For the first mode, a physically realistic equivalent rotating beam is possible, but for higher modes, the FSF has internal singularities. Strategies for addressing these singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test functions for rotating beam codes and also for targeted destiffening of rotating beams.

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90degrees and 180degrees is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress sigma(infinity) and electric field E-infinity or (sigma(infinity), E-infinity) is replaced by strain e and electric displacement D-infinity or (epsilon(infinity), D-infinity). Mixed conditions (sigma(infinity),D-infinity) and (epsilon(infinity),E-infinity) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E-y(infinity),sigma(y)(infinity)) for negative E-y(infinity) and (D-y(infinity), epsilon(y)(infinity)) for positive D-y(infinity) while the mechanical conditions sigma(y)(infinity) or epsilon(y)infinity are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth. (C) 2002 Published by Elsevier Science Ltd.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

A singularly perturbed linear two-point boundary-value problem

We consider the following singularly perturbed linear two-point boundary-value problem:

Ly(x) ≡ Ω(ε)D_xy(x) - A(x,ε)y(x) = f(x,ε) 0≤x≤1 (1a)

By ≡ L(ε)y(0) + R(ε)y(1) = g(ε) ε → 0^+ (1b)

Here Ω(ε) is a diagonal matrix whose first m diagonal elements are 1 and last m elements are ε. Aside from reasonable continuity conditions placed on A, L, R, f, g, we assume the lower right mxm principle submatrix of A has no eigenvalues whose real part is zero. Under these assumptions a constructive technique is used to derive sufficient conditions for the existence of a unique solution of (1). These sufficient conditions are used to define when (1) is a regular problem. It is then shown that as ε → 0^+ the solution of a regular problem exists and converges on every closed subinterval of (0,1) to a solution of the reduced problem. The reduced problem consists of the differential equation obtained by formally setting ε equal to zero in (1a) and initial conditions obtained from the boundary conditions (1b). Several examples of regular problems are also considered.

A similar technique is used to derive the properties of the solution of a particular difference scheme used to approximate (1). Under restrictions on the boundary conditions (1b) it is shown that for the stepsize much larger than ε the solution of the difference scheme, when applied to a regular problem, accurately represents the solution of the reduced problem.

Furthermore, the existence of a similarity transformation which block diagonalizes a matrix is presented as well as exponential bounds on certain fundamental solution matrices associated with the problem (1).