An alternative general solution of the steady-state heat diffusion equation

A unified analytical solution of the steady-state atmospheric diffusion equation

We consider the preconditioned Krylov subspace method for linear systems arising from the finite volume discretization method of steady-state variable-coefficient conservative space-fractional diffusion equations. We propose to use a scaled-circulant preconditioner to deal with such Toeplitz-like discretization matrices. We show that the difference between the scaled-circulant preconditioner and the coefficient matrix is equal to the sum of a small-norm matrix and a low-rank matrix. Numerical tests are conducted to show the effectiveness of the proposed method for one- and two-dimensional steady-state space-fractional diffusion equations and demonstrate that the preconditioned Krylov subspace method converges very quickly.

In this paper we describe, analyse and implement a parallel iterative method for the solution of the steady-state drift diffusion equations governing the behaviour of a semiconductor device in two space dimensions. The unknowns in our model are the electrostatic potential and the electron and hole quasi-Fermi potentials. Our discretisation consists of a finite element method with mass lumping for the electrostatic potential equation and a hybrid finite element with local current conservation properties for the continuity equations. A version of Gummel's decoupling algorithm which only requires the solution of positive definite symmetric linear systems is used to solve the resulting nonlinear equations. We show that this method has an overall rate of convergence which only degrades logarithmically as the mesh is refined. Indeed the (inner) nonlinear solves of the electrostatic potential equation converge quadratically, with a mesh independent asymptotic constant. We also describe an implementation on a MasPar MP-1 data parallel machine, where the required linear systems are solved by the preconditioned conjugate gradient method. Domain decomposition methods are used to parallelise the required matrix-vector multiplications and to build preconditioners for these very poorly-conditioned systems. Our preconditioned linear solves also have a rate of convergence which degrades logarithmically as the grid is refined relative to subdomain size, and their performance is resilient to the severe layers which arise in the coefficients of the underlying elliptic operators. Parallel experiments are given.

A potential and stream function analysis of two-dimensional steady-state convective diffusion equations involving Laplace fields

This study provides an approximate analytical solution to the fractional vibration problem of thin plate governing anomalous motion of plate subjected to in-plane loads. The method of variable separable is employed to transform the fractional partial differential equations under consideration into a fractional ordinary differential equation in temporal variable and a bi-harmonic plate equation in spatial variable. The technique of conformable fractional derivative is utilized to solve the resulting fractional differential equation and the approach of finite sine integral transform method is used to solve the accompanying bi-harmonic plate equation. The deflection field which measures the transverse displacement of the plate is expressed in terms of product of Bessel and trigonometric functions via the temporal and spatial variables respectively. The obtained solution reduces to the solution of the free vibration problem of thin plate in literature. This work shows that conformable fractional derivative is an efficient mathematical tool for tracking analytical solution of fractional partial differential equation governing anomalous vibration of thin plates.

Problems involving the dynamics of heating or cooling of physical systems submitted to certain boundary conditions are of great importance in a wide range of situations. Since simple situations involving cooling of food where the temperature of each part of the system needs to be controlled and, in general, does not involve the presence of heat sources, even more complex situations like nuclear reactors, where there is a heat source which is dependent on position and the dependence of the reactor reactivity on the temperature requires a good knowledge of the dynamics of this quantity, both for safety and efficiency. Problems of the steady-state heat diffusion equation for a series of physical systems are easy to solve and can be found in a vast literature on the subject. To find the time evolution of temperature T (x, t) of a system submitted to a heat source, however, is not a simple task and this type of problem is usually treated in an approximate way. In this work, for one-dimensional and homogeneous system, equilibrium T(x) and transient T (x, t) solutions are presented for cases of the existence of sources permeating throughout the medium. Each system considered is supposed to have constant temperature in the edges. Calculations are presented for two cases: one is the case of a source non-dependent on position and time and the other is the case of a source with sinusoidal spatial dependence, a situation close to what would occur in a one-dimensional nuclear reactor. The extension to a three-dimensional case, homogeneous sphere with constant temperature on its surface will also be presented. The case involving the contact of two subsystems of different materials will also be analyzed.

Accurate thermoplasmonic simulation of metallic nanoparticles

The brief presentation of the method of fast expansions is given to solve nonlinear differential equations. Application rules of the operator of fast expansions are specified for solving differential equations. According to the method of fast expansions, an unknown function can be represented as the sum of the boundary function and Fourier series sines and cosines for one variable. The special construction of the boundary functions leads to reasonably fast convergence of the Fourier series, so that for engineering calculations, it is sufficient to consider only the first three members. The method is applicable both to linear and nonlinear integro-differential systems. By means of applying the method of fast expansions to nonlinear Navier-Stokes equations the problem is reduced to a closed system of ordinary differential equations, which solution doesn't represent special difficulties. We can reapply the method of fast expansions to the resulting system of differential equations and reduce the original problem to a system of algebraic equations. If the problem is n -dimensional, then after n -fold application of the method of fast expansions the problem will be reduced to a closed algebraic system. Finally, we obtain an analytic-form solution of complicated boundary value problem in partial derivatives. The flow of an incompressible viscous fluid of Navier–Stokes is considered in a curvilinear pipe. The problem is reduced to solving a closed system of ordinary differential equations with boundary conditions by the method of fast expansions. The article considers peculiarities of finding the coefficients of boundary functions and Fourier coefficients for the zero-order and first-order operators of fast expansions. Obtaining the analytic-form solution is of great interest, because it allows to analyze and to investigate the influence of various factors on the properties of the viscous fluid in specific cases.

PARCS code multi-group neutron diffusion constants generation using Monte Carlo method

Measurements of effective sizes and diffusivities of nano-colloids and micro-particles

Proton radiography is an important diagnostic method for laser plasma experiments and is particularly important in the analysis of magnetized plasmas. The theory of radiographic image analysis has heretofore only permitted somewhat limited analysis of the radiographs of such plasmas. We furnish here a theory that remedies this deficiency. We show that to linear order in magnetic field gradients, proton radiographs are projection images of the MHD current along the proton trajectories. We demonstrate that in the linear regime (i.e., the small image contrast regime), the full structure of the projected perpendicular magnetic field can be reconstructed by solving a steady-state inhomogeneous 2-dimensional diffusion equation sourced by the radiograph fluence contrast data. We explore the validity of the scheme with increasing image contrast, as well as limitations of the inversion method due to the Poisson noise, discretization errors, radiograph edge effects, and obstruction by laser target structures. We also provide a separate analysis that is suited to the inference of isotropic-homogeneous magnetic turbulence spectra. Finally, we discuss extension of these results to the nonlinear regime (i.e., the order unity image contrast regime).

The mathematical model proposed by Chapman and Antano (Electrochimica Acta, 56 (2010), 128–132) for the catalytic electrochemical–chemical (EC’) processes in an irreversible second-order homogeneous reaction in a microelectrode is discussed. The mass-transfer boundary layer neighbouring an electrode can contribute to the electrode’s measured AC impedance. This model can be used to analyse membrane-transport studies and other instances of ionic transport in semiconductors and other materials. Two efficient and easily accessible analytical techniques, AGM and DTM, were used to solve the steady-state non-linear diffusion equation’s infinite layers. Herein, we present the generalized approximate analytical solution for the solute, product, and reactant concentrations and current for the small experimental values of kinetic and diffusion parameters. Using the Matlab/Scilab program, we also derive the numerical solution to this problem. The comparison of the analytical and numerical/computational results reveals a satisfactory level of agreement.

Often in reactor dynamics, higher eigenfunctions of the multigroup diffusion equation must be determined. An algorithm to calculate higher eigenfunctions (modes) of the [lambda]-eigenvalue problem corresponding to the steady-state two-group neutron diffusion equation is presented. The method is based on a special type of subspace iteration for large sparse nonsymmetric eigenvalue problems. Having been tested using an International Atomic Energy Agency benchmark problem and also applied to a VVER-1000 pressurized water reactor assembly, the algorithm was found to work very effectively and reliably. In its application, the algorithm presented is not restricted to the [lambda]-eigenvalue problem only but is also generally applicable to large sparse nonsymmetric eigenvalue problems even with multiple and complex eigenvalues.

The material in this thesis is the result of a year's experience in the solution of problems on the Caltech Electric Analog Computer. Although much work has been done elsewhere, the solution of partial differential equations is a relatively new field for the Caltech Computer. It is natural that such an undertaking should initiate points of view and techniques that differ from those of other investigators., This thesis contains the development of certain ideas that have been useful to the author in the solution of problems on the computer. In Part I, finite difference methods are treated with reference to problems with one space variable. Techniques are developed for the representation of differential operators by means of electrical networks and the question of unequal lumping is discussed. The solutions of the fourth order differential equations of a beam are treated in Part II. Such practical considerations as the effects of parasitic impedances and cell size are investigated. Solutions are presented for the normal modes of a cantilever beam, the transient vibration of a cantilever beam, and the coupled modes of vibration of an airplane wing. Problems involving the scalar Laplacian operator are treated in Part III. A general asymmetric network is developed that is useful for problems with irregular boundaries and for problems where it is desired to have variable cell size. These techniques are illustrated with respect to a cavity resonator problem and an electromagnetic radiation problem. Elastic plate problems are treated in Part IV. The analogy for the elastic plate is an extension of the beam analogy to two dimensions. Here the difficult problems are those relating to the representation of boundary conditions, particularly of boundary conditions along an irregular edge. Some conclusions regarding the construction of a network analyzer designed specifically for the solution of partial differential equations are given in Part V. The chief conclusion is that such a computer must contain a much larger number of electrical elements than are at present available in the Caltech Computer.

A theoretical model of reaction/diffusion within conducting polymer microelectrodes is discussed. The model is based on the steady-state diffusion equation containing a nonlinear term related to the Michaelis-Menten kinetic of the enzymatic reaction. An analytical expression pertaining to the concentration of substrate and current is obtained using homotopy perturbation method for all values of diffusion and the saturation parameter. The substrate concentration profile and current response can be used in a large range of concentrations including the non-linear contributions. These approximate analytical results were found to be in good agreement with the previously reported limiting case results.