A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Generalized integral transform solution for hydrodynamically developing non-Newtonian flows in circular tubes

ABSTRACT: The Generalized Integral Transform Technique (GITT) is applied to the solution of the momentum equations in a hydrodynamically developing laminar flow of a non-Newtonian power-law fluid inside a circular duct. A primitive variables formulation is adopted in order to avoid the singularity of the auxiliary eigenvalue problem in terms of Bessel functions at the centerline of the duct when the GITT approach is applied. Results for the velocity field and friction factor-Reynolds number product are computed for different power-law indices, which are tabulated and graphically presented as functions of the dimensionless coordinates. Critical comparisons with previous results in the literature are also performed, in order to validate the numerical codes developed in the present work and to demonstrate the consistency of the final results.

Integral transform solution for the forced convection of herschel-bulkley fluids in circular tubes and parallel-plates ducts

ABSTRACT: The thermal entry region in laminar forced convection of Herschel-Bulkley fluids is solved analytically through the integral transform technique, for both circular and parallel-plates ducts, which are maintained at a prescribed wall temperature or at a prescribed wall heat flux. The local Nusselt numbers are obtained with high accuracy in both developing and fully-developed thermal regions, and critical comparisons with previously reported numerical results are performed.

Integral transform method for laminar heat transfer convection of herschel-bulkley fluids within concentric annular ducts

ABSTRACT: Related momentum and energy equations describing the heat and fluid flow of Herschel-Bulkley fluids within concentric annular ducts are analytically solved using the classical integral transform technique, which permits accurate determination of parameters of practical interest in engineering such as friction factors and Nusselt numbers for the duct length. In analyzing the problem, thermally developing flow is assumed and the duct walls are subjected to boundary conditions of first kind. Results are computed for the velocity and temperature fields as well as for the parameters cited above with different power-law indices, yield numbers and aspect ratios. Comparisons are also made with previous work available in the literature, providing direct validation of the results and showing that they are consistent.

Approximate analytical solutions for porous catalysts with nonuniform activity

Using a novel finite integral transform technique, the problem of diffusion and chemical reaction in a porous catalyst with general activity profile is investigated theoretically. Analytical expressions for the effectiveness factor are obtained for pth order and Michaelis-Menten kinetics. Perturbation methods are employed to provide useful asymptotic solutions for large or small values of Thiele modulus and Biot number.

Finite element homogenization technique for the characterization of d(15) shear piezoelectric macro-fibre composites

A finite element homogenization method for a shear actuated d(15) macro-fibre composite (MFC) made of seven layers (Kapton, acrylic, electrode, piezoceramic fibre and epoxy composite, electrode, acrylic, Kapton) is proposed and used for the characterization of its effective material properties. The methodology is first validated for the MFC active layer only, made of piezoceramic fibre and epoxy, through comparison with previously published analytical results. Then, the methodology is applied to the seven-layer MFC. It is shown that the packaging reduces significantly the shear stiffness of the piezoceramic material and, thus, leads to significantly smaller effective electromechanical coupling coefficient k(15) and piezoelectric stress constant e(15) when compared to the piezoceramic fibre properties. However, it is found that the piezoelectric charge constant d(15) is less affected by the softer layers required by the MFC packaging.

Multicomponent reversible reactions inside a porous catalyst pellet

This work deals with a solution method to handle multicomponents reversible reactions occurring inside a porous catalyst pellet. The complexity of this problem arises from the fact that the effective diffusivities and Biot number, which characterizes the external mass transfer, are different for each chemical species. In mathematical terms, this means that each chemical species has its own subspace and, therefore, when the technique of finite integral transform is applied to solve this multicomponent problem, each chemical species is associated with its own integral transform kernel. The analytical solutions obtained for this problem are compact and simple for any further manipulation. Application of this result to the catalytic reforming of C7 hydrocarbon system is shown in this paper.

Covalidation of Integral Transforms and Method of Lines in Nonlinear Convection-Diffusion with Mathematica

The Mathematica system (version 4.0) is employed in the solution of nonlinear difusion and convection-difusion problems, formulated as transient one-dimensional partial diferential equations with potential dependent equation coefficients. The Generalized Integral Transform Technique (GITT) is first implemented for the hybrid numerical-analytical solution of such classes of problems, through the symbolic integral transformation and elimination of the space variable, followed by the utilization of the built-in Mathematica function NDSolve for handling the resulting transformed ODE system. This approach ofers an error-controlled final numerical solution, through the simultaneous control of local errors in this reliable ODE's solver and of the proposed eigenfunction expansion truncation order. For covalidation purposes, the same built-in function NDSolve is employed in the direct solution of these partial diferential equations, as made possible by the algorithms implemented in Mathematica (versions 3.0 and up), based on application of the method of lines. Various numerical experiments are performed and relative merits of each approach are critically pointed out.

Generating finite integral relation algebras

Relation algebras and categories of relations in particular have proven to be extremely useful as a fundamental tool in mathematics and computer science. Since relation algebras are Boolean algebras with some well-behaved operations, every such algebra provides an atom structure, i.e., a relational structure on its set of atoms. In the case of complete and atomic structure (e.g. finite algebras), the original algebra can be recovered from its atom structure by using the complex algebra construction. This gives a representation of relation algebras as the complex algebra of a certain relational structure. This property is of particular interest because storing the atom structure requires less space than the entire algebra. In this thesis I want to introduce and implement three structures representing atom structures of integral heterogeneous relation algebras, i.e., categorical versions of relation algebras. The first structure will simply embed a homogeneous atom structure of a relation algebra into the heterogeneous context. The second structure is obtained by splitting all symmetric idempotent relations. This new algebra is in almost all cases an heterogeneous structure having more objects than the original one. Finally, I will define two different union operations to combine two algebras into a single one.

Análise da magnetohidrodinâmica com transferência de calor em canais de placas paralelas via transformação integral

The main goal of the present work is related to the dynamics of the steady state, incompressible, laminar flow with heat transfer, of an electrically conducting and Newtonian fluid inside a flat parallel-plate channel under the action of an external and uniform magnetic field. For solution of the governing equations, written in the parabolic boundary layer and stream-function formulation, it was employed the hybrid, numericalanalytical, approach known as Generalized Integral Transform Technique (GITT). The flow is sustained by a pressure gradient and the magnetic field is applied in the direction normal to the flow and is assumed that normal magnetic field is kept uniform, remaining larger than any other fields generated in other directions. In order to evaluate the influence of the applied magnetic field on both entrance regions, thermal and hydrodynamic, for this forced convection problem, as well as for validating purposes of the adopted solution methodology, two kinds of channel entry conditions for the velocity field were used: an uniform and an non-MHD parabolic profile. On the other hand, for the thermal problem only an uniform temperature profile at the channel inlet was employed as boundary condition. Along the channel wall, plates are maintained at constant temperature, either equal to or different from each other. Results for the velocity and temperature fields as well as for the main related potentials are produced and compared, for validation purposes, to results reported on literature as function of the main dimensionless governing parameters as Reynolds and Hartman numbers, for typical situations. Finally, in order to illustrate the consistency of the integral transform method, convergence analyses are also effectuated and presented

Heat transfer in laminar flow of non-Newtonian fluids in ducts of elliptical section

Laminar forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumption used in this work is a laminar flow of a power flow inside elliptical tube, under a boundary condition of first kind with constant physical properties and negligible axial heat diffusion (high Peclet number). To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number and the average Nusselt number for various cross-section aspect ratios. (C) 2006 Elsevier. SAS. All rights reserved.

Thermally developing forced convection of non-Newtonian fluids inside elliptical ducts

Laminar-forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumptions used in this work are a non-Newtonian fluid, laminar flow, constant physical properties, and negligible axial heat diffusion (high Peclet number). Most of the previous research in elliptical ducts deal mainly with aspects of fully developed laminar flow forced convection, such as velocity profile, maximum velocity, pressure drop, and heat transfer quantities. In this work, we examine heat transfer in a hydrodynamically developed, thermally developing laminar forced convection flow of fluid inside an elliptical tube under a second kind of a boundary condition. To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform, where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number, and the average Nusselt number for various cross-section aspect ratios.

Aplicação da transformada integral generalizada no escoamento potencial em contrações

Pós-graduação em Engenharia Mecânica - FEIS

Aplicação da transformada integral e da transformação conforme na solução de uma classe de problemas difusivo-convectivos em domínios de geometrias não-convencionais

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Levy stable distributions via associated integral transform

We present a method of generation of exact and explicit forms of one-sided, heavy-tailed Levy stable probability distributions g(alpha)(x), 0 <= x < infinity, 0 < alpha < 1. We demonstrate that the knowledge of one such a distribution g a ( x) suffices to obtain exactly g(alpha)p ( x), p = 2, 3, .... Similarly, from known g(alpha)(x) and g(beta)(x), 0 < alpha, beta < 1, we obtain g(alpha beta)( x). The method is based on the construction of the integral operator, called Levy transform, which implements the above operations. For a rational, alpha = l/k with l < k, we reproduce in this manner many of the recently obtained exact results for g(l/k)(x). This approach can be also recast as an application of the Efros theorem for generalized Laplace convolutions. It relies solely on efficient definite integration. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4709443]