Periodic solution of the generalized Rayleigh equation

The periodic solutions of a strongly cubic nonlinear oscillator whose motion is described with the generalized Rayleigh equation are studied. Approximate analytic solving methods are introduced. A new method based on homotopy and averaging is developed to determine the limit cycle motion. The obtained analytical solutions are compared with those calculated by the elliptic harmonic balance method with generalized Fourier series and Jacobian elliptic functions. Three types of cubic nonlinearity are considered: the coefficients of the linear and cubic terms are positive, the coefficient of the linear term is positive and that of the cubic term is negative and the opposite case. Comparisons of the analytical solution and numerical solution, obtained by using the Runge-Kutta method, are illustrated with examples.     

激光衍射微粒测量中的独立模式反演算法

在用激光衍射法测量微粒直径及其分布时得到的矩阵方程组是一组病态方程 ,无法直接求解。目前大多数求解方法是采用已知模式的拟合算法。但由于在实际中的微粒直径的分布是多种多样的 ,因此难以用一种或几种分布函数去近似描述。试图从微粒衍射的特征着手提出一种未知分布的独立模式的求解方法     

Numerical study of a stochastic particle algorithm solving a multidimensional population balance model for high shear granulation

This paper is concerned with computational aspects of a multidimensional population balance model of a wet granulation process. Wet granulation is a manufacturing method to form composite particles, granules, from small particles and binders. A detailed numerical study of a stochastic particle algorithm for the solution of a five-dimensional population balance model for wet granulation is presented. Each particle consists of two types of solids (containing pores) and of external and internal liquid (located in the pores). Several transformations of particles are considered, including coalescence, compaction and breakage. A convergence study is performed with respect to the parameter that determines the number of numerical particles. Averaged properties of the system are computed. In addition, the ensemble is subdivided into practically relevant size classes and analysed with respect to the amount of mass and the particle porosity in each class. These results illustrate the importance of the multidimensional approach. Finally, the kinetic equation corresponding to the stochastic model is discussed.     

A Method for Solving Initial-value Problem and Finding Exact Solutions of Nonlinear Partial Differential Equations

Homogeneous balance method for solving nonlinear partial differential equation(s) is extended to solving initial-value problem and getting new solution(s) from a known solution of the equation(s) under consideration. The approximate equations for long water waves are chosen to illustrate the method, infinitely many simple-solitary-wave solutions and infinitely many rational function solutions, especially the closed form of the solution for initial-value problem, are obtained by using the extended homogeneous balance method given here.     

A fluctuation problem in particle transport theory II

A balance equation is formulated for the probability that a particle injected into an infinite, amorphous medium will have suffered N collisions and have given rise to n new particles in a given energy range at time t. The method of regeneration points has been employed and this leads, in the case of two particle production, to a non-linear, integro-differential equation for the probability generating function. This equation is solved for the case of foreign particles slowing down, in which case it becomes linear and results are obtained which include the effects of electronic stopping and absorption, thus generalizing the work in part I. In the cascade problem, a single particle gives rise to two new particles in every collision and it is shown, for a simple hard-sphere model with 1/v scattering and absorption, how the non-linear equation may be solved. The probability for the number of particles and the number of collisions suffered to absorption is obtained in the case of zero absorption, the probability law is shown to obey a Furry distribution. The limitations of the method described in part I for dealing with cascades are highlighted.     

A reliable treatment for nonlinear Schrödinger equations

Exp-function method is used to find a unified solution of nonlinear wave equation. Nonlinear Schrödinger equations with cubic and power law nonlinearity are selected to illustrate the effectiveness and simplicity of the method. It is shown that the Exp-function method, with the help of symbolic computation, provides a powerful mathematical tool for solving nonlinear equation.     

The shallow water equations in Lagrangian coordinates

Recent advances in the collection of Lagrangian data from the ocean and results about the well-posedness of the primitive equations have led to a renewed interest in solving flow equations in Lagrangian coordinates. We do not take the view that solving in Lagrangian coordinates equates to solving on a moving grid that can become twisted or distorted. Rather, the grid in Lagrangian coordinates represents the initial position of particles, and it does not change with time. We apply numerical methods traditionally used to solve differential equations in Eulerian coordinates, to solve the shallow water equations in Lagrangian coordinates. The difficulty with solving in Lagrangian coordinates is that the transformation from Eulerian coordinates results in solving a highly nonlinear partial differential equation. The non-linearity is mainly due to the Jacobian of the coordinate transformation, which is a precise record of how the particles are rotated and stretched. The inverse Jacobian must be calculated, thus Lagrangian coordinates cannot be used in instances where the Jacobian vanishes. For linear (spatial) flows we give an explicit formula for the Jacobian and describe the two situations where the Lagrangian shallow water equations cannot be used because either the Jacobian vanishes or the shallow water assumption is violated. We also prove that linear (in space) steady state solutions of the Lagrangian shallow water equations have Jacobian equal to one. In the situations where the shallow water equations can be solved in Lagrangian coordinates, accurate numerical solutions are found with finite differences, the Chebyshev pseudospectral method, and the fourth order Runge–Kutta method. The numerical results shown here emphasize the need for high order temporal approximations for long time integrations.     

A fluctuation problem in particle transport theory I

The fluctuation in the number of collision suffered by particles as they slow down in a moderating medium is studied via a probability balance equation. The equation describes the collision history of foreign particles slowing down in a host medium and also accounts for the recoil particles produced in the collision. The equations are solved by the introduction of a generating function from which the space and time dependent probability distributions are obtained. That is, the probability that a particle will suffer just N collisions to reach energy E at a time t after injection. In the space dependent case it is the probability that a particle suffers just N collisions to travel a given path length before coming to rest.Explicit expressions for the means and variances are obtained by solving a difference equation. From this solution it has been possible to obtain exact expressions for hard spheres and for a variety of models based on the inverse power law approximation. A number of new results are presented and some old ones rederived in a more efficient and general manner. The results are of value in the understanding of radiation damage cascades and in neutron slowing down in moderating materials.     

蒙特卡罗方法在解微分方程边值问题中的应用

介绍了蒙特卡罗方法的基本原理以及随机数的产生方法。基于蒙特卡罗方法的思想,结合有限差分方法,建立了求解微分方程边值问题的随机概率模型,并以第一类边界条件的拉普拉斯方程和一个给定初值及边界条件的非稳态热传导方程为数值算例,研究了蒙特卡罗方法在求解微分方程边值问题中的应用。结果表明:利用蒙特卡罗方法,不仅可以有效解决给定边界条件的微分方程,对于给定初值条件的微分方程,也可以从时域有限差分方程出发,采用蒙特卡罗方法进行求解。数值模拟和对误差的理论分析均表明,增加蒙特卡罗试验中的模拟粒子点数,可以提高计算结果的精度。     

The radial Marchenko-Blažek integral equation for complete spectra

The same simple and standard way, by means of which an inverse problem for scattering of spinless particles by central potentials is solved in the Gelfand-Levitan method, is applied to the Marchenko method for general angular momenta including the bound states. We first derive an integral equation for the kernel with a triangularity property, which relates it to a potential and then the other one, which connects the kernel with spectral data. A solution corresponding to the general Yukawa potential is found and some formulae are checked by solving the problem of the restrained phase equivalent potentials.     

The general solution of the binding mean spherical approximation for pairing ions

The mean spherical approximation (MSA) for an arbitrary mixture of charged hard spheres with saturating bonds is solved in the Wertheim formalism. Only pairs are allowed. It is shown that the general solution is given in terms of two scaling parameters and the set of binding fractions. One of the scales is the screening parameter, and the other is a coupling parameter analogous to that of the simple MSA, but that now is found solving a cubic equation. Therefore the full solution requires solvingm+2 nonlinear algebraic equations for a system withm components. A brief discussion of the thermodynamics is given.     

Modified Kudryashov method via new exact solutions for some conformable fractional differential equations arising in mathematical biology

In this study, the modified Kudryashov method is used to construct new exact solutions for some conformable fractional differential equations. By implementing the conformable fractional derivative and compatible fractional complex transforms, the fractional generalized reaction duffing (RD) model equation, the fractional biological population model and the fractional diffusion reaction (DR) equation with quadratic and cubic nonlinearity are discussed. As an outcome, some new exact solutions are formally established. All solutions have been verified back into its corresponding equation with the aid of maple package program. We assure that the employed method is simple and robust for the estimation of the new exact solutions, and practically capable for reducing the size of computational work for solving a various class of fractional differential equations arising in applied mathematics, mathematical physics and biology.     

Tau method approximation for radiative transfer problems in a slab medium

An approximate method for solving the radiative transfer equation in a slab medium with an isotropic scattering is proposed. The method is based upon constructing the double Legendre series to approximate the required solution using Legendre tau method. The differential and integral expressions which arise in the radiative transfer equation are converted into a system of linear algebraic equations which can be solved for the unknown coefficients. Numerical examples are included to demonstrate the validity and applicability of the method and a comparison is made with existing results.     

Solution procedure of residue harmonic balance method and its applications

This paper presents a simple and rigorous solution procedure of residue harmonic balance for predicting the accurate approximation of certain autonomous ordinary differential systems. In this solution procedure, no small parameter is assumed. The harmonic residue of balance equation is separated in two parts at each step. The first part has the same number of Fourier terms as the present order of approximation and the remaining part is used in the subsequent improvement. The corrections are governed by linear ordinary differential equation so that they can be solved easily by means of harmonic balance method again. Three kinds of different differential equations involving general, fractional and delay ordinary differential systems are given as numerical examples respectively. Highly accurate limited cycle frequency and amplitude are captured. The results match well with the exact solutions or numerical solutions for a wide range of control parameters. Comparison with those available shows that the residue harmonic balance solution procedure is very effective for these autonomous differential systems. Moreover, the present method works not only in predicting the amplitude but also the frequency of bifurcated period solution for delay ordinary differential equation.     

Lattice Boltzmann Study of Mixed Convection in aCubic Cavity

The problem of the mixed convection in a cubic cavity is studied with lattice Boltzmann method. A multiple-relaxation-time lattice Boltzmann model for incompressible flow in the cubic cavity and another thermal lattice Boltzmann model for solving energy/temperature equation are proposed. The present models are first validated through a comparison with some available results, and then, we present a detailed parameter study on the mixed
convection in the cubic cavity. The numerical results show that the flow and temperature patterns change greatly with variations of the Reynolds and Richardson numbers.     

Electromagnetic wave scattering by small bodies

A reduction of the Maxwell's system to a Fredholm second-kind integral equation with weakly singular kernel is given for electromagnetic (EM) wave scattering by one and many small bodies. This equation is solved asymptotically as the characteristic size of the bodies tends to zero. The technique developed is used for solving the many-body EM wave scattering problem by rigorously reducing it to solving linear algebraic systems, completely bypassing the usage of integral equations. An equation is derived for the effective field in the medium, in which many small particles are embedded. A method for creating a desired refraction coefficient is outlined.     

Helmholtz and parabolic equation solutions to a benchmark problem in ocean acoustics

The Helmholtz equation (HE) describes wave propagation in applications such as acoustics and electromagnetics. For realistic problems, solving the HE is often too expensive. Instead, approximations like the parabolic wave equation (PE) are used. For low-frequency shallow-water environments, one persistent problem is to assess the accuracy of the PE model. In this work, a recently developed HE solver that can handle a smoothly varying bathymetry, variable material properties, and layered materials, is used for an investigation of the errors in PE solutions. In the HE solver, a preconditioned Krylov subspace method is applied to the discretized equations. The preconditioner combines domain decomposition and fast transform techniques. A benchmark problem with upslope-downslope propagation over a penetrable lossy seamount is solved. The numerical experiments show that, for the same bathymetry, a soft and slow bottom gives very similar HE and PE solutions, whereas the PE model is far from accurate for a hard and fast bottom. A first attempt to estimate the error is made by computing the relative deviation from the energy balance for the PE solution. This measure gives an indication of the magnitude of the error, but cannot be used as a strict error bound.     

Particle synthesis in flames

From the view point of combustion science, the fundamentals of particle formation in a flame environment are discussed. The principles of converting a gas phase precursor dopant into particles and also their growth are addressed. Various experimental methods and examples for the synthesis of particles of both single and mixed oxides are reviewed. First attempts to tune the stoichiometry of oxide particles by varying the combustion parameters of premixed flames are demonstrated. Some aspects of modelling and the problems which still need to be solved are illustrated by means of the particles’ population balance equation.     

Thermal analysis for transient radiative cooling of a conducting semitransparent layer of ceramic in high-temperature applications

A method is developed for obtaining transient temperature distribution in a cooling semitransparent layer of ceramic. The layer is emitting, absorbing, isotropically scattering and heat conducting with a refractive index ranging from 1 to 2. The solution involves solving simultaneously the energy equation and the integral equation for the radiative flux gradient. The energy equation is solved using an implicit finite volume scheme and the integral equation of radiative heat transfer is solved using the singularity technique and Gaussian integration. The effects of scattering are investigated. It is shown that scattering has a significant effect on the transient temperature distribution and the transient mean temperature of the layer.     

Time shift and superposition method for solving transient radiative transfer equation

Due to the long-time transient response of pulse irradiation, the computational time required for solving transient radiative transfer (TRT) is often very long, especially for the case in which the boundary is subjected to continuous pulse train and the geometry is complicated. In addition, sometimes, before actual experiments are carried out, a suitable pulse shape or type often needs to be selected by numerical simulation and comparison. Because the numerical solution of TRT needs to be repeated many times, the selection processes is very time-consuming. In this paper, by considering that the TRT equation and its initial and boundary conditions are linear, a time shift and superposition method is developed for solving TRT equation in non-emitting media, in which only the transient response of a short square pulse needs to be solved, and the solution of TRT under any pulse shape can be constructed by time shift and then superposition using the basis solution of the short square pulse. Three numerical examples are studied to illustrate the peformance of the superposition method in solving TRT problems. The results show that the superposition is effective, accurate and very suitable for solving TRT in the medium subjected to a series of pulse train.