Search results

The purpose of this paper is to study the propagation of Rayleigh waves in thermoelastic medium with mass diffusion.

The field equations for the linear theory of homogeneous isotropic thermoelastic diffusion medium are taken into consideration by using dual-phase-lag heat transfer (DPLT) and dual-phase-lag diffusion (DPLD) models. Using the potential functions and harmonic wave solution, three coupled dilatational waves and a shear wave is obtained. After developing mathematical formulation, the dispersion equation is obtained, which results to be complex and irrational. This equation is converted into a polynomial form of higher degree.

From the polynomial equation, Rayleigh wave root is found. The secular equation is resolved into a polynomial form to find the roots and therefore to find the existence and propagation of Rayleigh wave. The existence of Rayleigh wave in the assumed model depends on the values of various parameters involved in the secular equation. These roots are resolved for phase velocity and attenuation of the inhomogeneous propagation of Rayleigh wave. Behavior of particle motion of these waves inside and at the surface of the thermoelastic medium with mass diffusion is studied. Particular cases of the interest are also deduced from the present investigation.

Governing equations corresponding to DPLT and DPLD models of thermoelastic diffusion are formulated to study the wave propagation and their dependence on various material parameters. In this paper effects of thermal and diffusion phase lags on the phase velocity, attenuation and on particle paths are observed and depicted graphically.

The purpose of this paper is to develop a new finite-volume method of lines for one-dimensional reaction-diffusion equations that provides piece-wise analytical solutions in space and is conservative, compare it with other finite-difference discretizations and assess the effects of the nonlinear diffusion coefficient on wave propagation.

A conservative, finite-volume method of lines based on piecewise integration of the diffusion operator that provides a globally continuous approximate solution and is second-order accurate is presented. Numerical experiments that assess the accuracy of the method and the time required to achieve steady state, and the effects of the nonlinear diffusion coefficients on wave propagation and boundary values are reported.

The finite-volume method of lines presented here involves the nodal values and their first-order time derivatives at three adjacent grid points, is linearly stable for a first-order accurate Euler’s backward discretization of the time derivative and has a smaller amplification factor than a second-order accurate three-point centered discretization of the second-order spatial derivative. For a system of two nonlinearly-coupled, one-dimensional reaction-diffusion equations, the amplitude, speed and separation of wave fronts are found to be strong functions of the dependence of the nonlinear diffusion coefficients on the concentration and temperature.

A new finite-volume method of lines for one-dimensional reaction-diffusion equations based on piecewise analytical integration of the diffusion operator and the continuity of the dependent variables and their fluxes at the cell boundaries is presented. The method may be used to study heat and mass transfer in layered media.

The purpose of this paper is to depict the effect of thermal and diffusion phase-lags on plane waves propagating in thermoelastic diffusion medium with different material symmetry. A generalized form of mass diffusion equation is introduced instead of classical Fick's diffusion theory by using two diffusion phase-lags, one phase-lag of diffusing mass flux vector, represents the delayed time required for the diffusion of the mass flux and the other phase-lag of chemical potential, represents the delayed time required for the establishment of the potential gradient. The basic equations for the anisotropic thermoelastic diffusion medium in the context of dual-phase-lag heat transfer (DPLT) and dual-phase-lag diffusion (DPLD) models are presented. The governing equations for transversely isotropic and isotropic case are also reduced. The different characteristics of waves like phase velocity, attenuation coefficient, specific loss and penetration depth are computed numerically. Numerically computed results are depicted graphically for anisotropic, transversely isotropic and isotropic medium. The effect of diffusion and thermal phase-lags are shown on the different characteristic of waves. Some particular cases of result are also deduced from the present investigation.

The governing equations of thermoelastic diffusion are presented using DPLT model and a new model of DPLD. Effect of phase-lags of thermal and diffusion is presented on different characteristic of waves.

The effect of diffusion and thermal phase-lags on the different characteristic of waves is appreciable. Also the use of diffusion phase-lags in the equation of mass diffusion gives a more realistic model of thermoelastic diffusion media as it allows a delayed response between the relative mass flux vector and the potential gradient.

Introduction of a new model of DPLD in the equation of mass diffusion.

The purpose of this paper is to develop an alternative numerical approach for describing fractional diffusion in Cartesian and non‐Cartesian domains using a Monte Carlo random walk scheme. The resulting domain shifting scheme provides a numerical solution for multi‐dimensional steady state, source free diffusion problems with fluxes expressed in terms of Caputo fractional derivatives. This class of problems takes account of non‐locality in transport, expressed through parameters representing both the extent and direction of the non‐locality.

The method described here follows a similar approach to random walk methods previously developed for normal (local) diffusion. The key differences from standard methods are: first, the random shifting of the domain about the point of interest with, second, shift steps selected from non‐symmetric, power‐law tailed, Lévy probability distribution functions.

The domain shifting scheme is verified by comparing predictive solutions to known one‐dimensional and two‐dimensional analytical solutions for fractional diffusion problems. The scheme is also applied to a problem of fractional diffusion in a non‐Cartesian annulus domain. In contrast to the axisymmetric, steady state solution for normal diffusion, a non‐axisymmetric solution results.

This is the first random walk scheme to utilize the concept of allowing the domain to undergo the random walk about a point of interest. Domain shifting scheme solutions of fractional diffusion in non‐Cartesian domains provide an invaluable tool to direct the development of more sophisticated grid based finite element inspired fractional diffusion schemes.

The purpose of this paper is to determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of microfluidic droplet ensembles, water flow in unsaturated flows, infiltration, etc., as functions of the powers of the convection and diffusion fluxes and upstream boundary condition; to study numerically the evolution of the wave for two different initial conditions; and to assess the accuracy of several finite difference methods for the solution of the degenerate, nonlinear, advection--diffusion equation that governs the model.

The theory of ordinary differential equations and several explicit, finite difference methods that use first- and second-order, accurate upwind, central and compact discretizations for the convection terms are used to determine the analytical solution for steadily propagating waves and the evolution of the wave fronts from hyperbolic tangent and piecewise linear initial conditions to steadily propagating waves, respectively. The amplitude and phase errors of the semi-discrete schemes are determined analytically and the accuracy of the discrete methods is assessed.

For non-zero upstream boundary conditions, it has been found both analytically and numerically that the shock wave is smooth and its steepness increases as the power of the diffusion term is increased and as the upstream boundary value is decreased. For zero upstream boundary conditions, smooth, cusped and sharp shock waves may be encountered depending on the powers of the convection and diffusion terms. For a linear diffusion flux, the shock wave is smooth, whereas, for a quadratic diffusion flux, the wave exhibits a cusped front whose left spatial derivative decreases as the power of the convection term is increased. For higher nonlinear diffusion fluxes, a sharp shock wave is observed. The wave speed decreases as the powers of both the convection and the diffusion terms are increased. The evolution of the solution from hyperbolic tangent and piecewise linear initial conditions shows that the wave back adapts rapidly to its final steady value, whereas the wave front takes much longer, especially for piecewise linear initial conditions, but the steady wave profile and speed are independent of the initial conditions. It is also shown that discretization of the nonlinear diffusion flux plays a more important role in the accuracy of first- and second-order upwind discretizations of the convection term than either a conservative or a non-conservative discretization of the latter. Second-order upwind and compact discretizations of the convection terms are shown to exhibit oscillations at the foot of the wave’s front where the solution is nil but its left spatial derivative is largest. The results obtained with a conservative, centered second--order accurate finite difference method are found to be in good agreement with those of the second-order accurate, central-upwind Kurganov--Tadmor method which is a non-oscillatory high-resolution shock-capturing procedure, but differ greatly from those obtained with a non-conservative, centered, second-order accurate scheme, where the gradients are largest.

A new, one-dimensional model for microfluidic droplet transport, water flow in unsaturated flows, infiltration, etc., that includes high-order convection fluxes and degenerate diffusion, is proposed and studied both analytically and numerically. Its smooth, cusped and sharp shock wave solutions have been determined analytically as functions of the powers of the nonlinear convection and diffusion fluxes and the boundary conditions. These solutions are used to assess the accuracy of several finite difference methods that use different orders of accuracy in space, and different discretizations of the convection and diffusion fluxes, and can be used to assess the accuracy of other numerical procedures for one-dimensional, degenerate, convection--diffusion equations.

The paper gives the description of boundary element method (BEM) with subdomains for the solution of convection—diffusion equations with variable coefficients and Burgers’ equations. At first, the whole domain is discretized into K subdomains, in which linearization of equations by representing convective velocity by the sum of constant and variable parts is carried out. Then using fundamental solutions for convection—diffusion linear equations for each subdomain the boundary integral equation (in which the part of the convective term with the constant convective velocity is not included into the pseudo‐body force) is formulated. Only part of the convective term with the variable velocity, which is, as a rule, more than one order less than convective velocity constant part contribution, is left as the pseudo‐source. On the one hand, this does not disturb the numerical BEM—algorithm stability and, on the other hand, this leads to significant improvement in the accuracy of solution. The global matrix, similar to the case of finite element method, has block band structure whereas its width depends only on the numeration order of nodes and subdomains. It is noted, that in comparison with the direct boundary element method the number of global matrix non‐zero elements is not proportional to the square of the number of nodes, but only to the total number of nodal points. This allows us to use the BEM for the solution of problems with very fine space discretization. The proposed BEM with subdomains technique has been used for the numerical solution of one‐dimensional linear steady‐state convective—diffusion problem with variable coefficients and one‐dimensional non‐linear Burgers’ equation for which exact analytical solutions are available. It made it possible to find out the BEM correctness according to both time and space. High precision of the numerical method is noted. The good point of the BEM is the high iteration convergence, which is disturbed neither by high Reynolds numbers nor by the presence of negative velocity zones.

The purpose of this paper is to depict the effect of time and thermal and diffusion phase-lags due to axisymmetric heat supply in a ring. The problem is discussed within the context of dual-phase-lag heat transfer and dual-phase-lag diffusion models. The upper and lower surfaces of the ring are traction free and subjected to an axisymmetric heat supply.

The solution is found by using Laplace and Hankel transform technique and a direct approach without the use of potential functions. The analytical expressions of displacements, stresses and chemical potential, temperature and mass concentration are computed in transformed domain. Numerical inversion technique has been applied to obtain the results in the physical domain. Numerically simulated results are depicted graphically. The effect of time and diffusion and thermal phase-lags are shown on the various components. Some particular cases of result are also deduced from the present investigation.

It is observed that change in time changes the behaviour of deformations of the various components of stresses, displacements, chemical potential function, temperature change and mass concentration. The authors find that for t=0.2, trends are oscillatory in all the cases whereas for t=0.1, trends are quite different. A sound impact of diffusion and thermal phase-lags on the various quantities is observed. A lot of difference in the trends of single phase lag and dual phase lag is observed. The use of diffusion phase-lags in the equation of mass diffusion gives a more realistic model of thermoelastic diffusion media as it allows a delayed response between the relative mass flux vector and the potential gradient.

This problem is totally new because dual phase lag is applied in heat conduction and diffusion equation while considering the problem of plate in axisymmetric heat supply.

The paper's aim is to investigate a two‐dimensional deformation of homogeneous, anisotropic generalized thermoelastic diffusion as a result of an inclined load by applying Laplace and Fourier transforms. The inclined load is assumed to be a linear combination of a normal load and a tangential load.

As an application, concentrated and distributed loads have been taken to illustrate the utility of the approach. The transformed solutions are inverted numerically, using a numerical inversion technique.

The variations of normal displacement, temperature distribution and chemical potential distribution due to different sources for different angle of inclinations with distance have been shown graphically to depict the effect of diffusion and anisotropy. A special case is also deduced from the present investigation.

It can contribute to the theoretical consideration of the seismic and volcanic sources since it can account for the deformation fields in the entire volume surrounding the sources region.

The purpose of this study is to develop a new method of lines for one-dimensional (1D) advection-reaction-diffusion (ADR) equations that is conservative and provides piecewise analytical solutions in space, compare it with other finite-difference discretizations and assess the effects of advection and reaction on both 1D and two-dimensional (2D) problems.

A conservative method of lines based on the piecewise analytical integration of the two-point boundary value problems that result from the local solution of the advection-diffusion operator subject to the continuity of the dependent variables and their fluxes at the control volume boundaries is presented. The method results in nonlinear first-order, ordinary differential equations in time for the nodal values of the dependent variables at three adjacent grid points and triangular mass and source matrices, reduces to the well-known exponentially fitted techniques for constant coefficients and equally spaced grids and provides continuous solutions in space.

The conservative method of lines presented here results in three-point finite difference equations for the nodal values, implicitly treats the advection and diffusion terms and is unconditionally stable if the reaction terms are implicitly treated. The method is shown to be more accurate than other three-point, exponentially fitted methods for nonlinear problems with interior and/or boundary layers and/or source/reaction terms. The effects of linear advection in 1D reacting flow problems indicates that the wave front steepens as it approaches the downstream boundary, whereas its back corresponds to a translation of the initial conditions; for nonlinear advection, the wave front exhibits steepening but the wave back shows a linear dependence on space. For a system of two nonlinearly coupled, 2D ADR equations, it is shown that a counter-clockwise rotating vortical field stretches the spiral whose tip drifts about the center of the domain, whereas a clock-wise rotating one compresses the wave and thickens its arms.

A new, conservative method of lines that implicitly treats the advection and diffusion terms and provides piecewise-exponential solutions in space is presented and applied to some 1D and 2D advection reactions.

Taiwan is one of several leading countries in the mobile music context. Accordingly, Taiwan's experiences in promoting mobile music service diffusion are of importance and interest to international practitioners and researchers. Applying Rogers' innovation diffusion theory, this study aims to employ econometric models to investigate whether the diffusion of mobile music service adoption is affected by external influences (e.g. mass media advertising, salespeople, and service providers), internal influences (e.g. interactions and imitations among acquaintances), or a combination of such influences.

To determine which influence best explained the diffusion of mobile music adoptions, the external, internal, and Von Bertalanffy mixed influence diffusion models were tested in this study. GNUS, a strongly functional language and environment to statistically explore data sets, was used to estimate the parameters of each model. The performance of each diffusion model was then examined using the Akaike AIC and Schwarz BIC statistics.

Findings indicated that the Von Bertalanffy mixed influences model best describes the diffusion pattern of mobile music service adoption and that acquaintances' influence in terms of interactions is the dominant factor influencing mobile music service adoption decision in Taiwan.

How managers of a mobile music service provider can use the internal and external influences interchangeably to effectively accelerate the mobile music diffusion at the different stage of product lifecycle is presented in this study. Indeed, the mobile music service is one of the most important industries worldwide not only because its penetration rate in many countries is over 50 percent, but also because of its killer applications. In light of this, the study contributes highly to theoretical and empirical examinations because the diffusion of the mobile music services within a society is the essence of the development/usage of the m‐commerce or music industries.